Can Sgesv divided in Sgetrf and Sgetrs?
7 posts
• Page 1 of 1
Can Sgesv divided in Sgetrf and Sgetrs?
Hi:
As part of debugging my CulaDeviceSgesv problems, see:
http://www.culatools.com/forums/15-gene ... art=12#337
I tried to divide CulaDeviceSgesv in:
- CulaDeviceSgetrf and
- CulaDeviceSgetrs
When executing CulaDeviceSgetrs, I obtain the following runtime error:
"A required feature is not yet implemented"
Am I doing anything bad?
With many thanks in advance
jpeinado
As part of debugging my CulaDeviceSgesv problems, see:
http://www.culatools.com/forums/15-gene ... art=12#337
I tried to divide CulaDeviceSgesv in:
- CulaDeviceSgetrf and
- CulaDeviceSgetrs
When executing CulaDeviceSgetrs, I obtain the following runtime error:
"A required feature is not yet implemented"
Am I doing anything bad?
With many thanks in advance
jpeinado
- jpeinado
- Posts: 37
- Joined: Mon Sep 14, 2009 10:48 am
Re:Can Sgesv divided in Sgetrf and Sgetrs?
Getrs only accepts the 'N' job code presently and will respond with the message you received for any other job code. But you only need the 'N' case to implement GESV. Don't forget that you also need a matrix permute to implement GESV.
- john
- Administrator
- Posts: 587
- Joined: Thu Jul 23, 2009 2:31 pm
Re:Can Sgesv divided in Sgetrf and Sgetrs?
Update: we have discovered an error regarding the detection of the 'N' job code for getrs. We will likely post an update soon to correct this; I apologize for the inconvenience. In the meantime, the gesv routine will function correctly.
- john
- Administrator
- Posts: 587
- Joined: Thu Jul 23, 2009 2:31 pm
Re:Can Sgesv divided in Sgetrf and Sgetrs?
john wrote:Update: we have discovered an error regarding the detection of the 'N' job code for getrs. We will likely post an update soon to correct this; I apologize for the inconvenience. In the meantime, the gesv routine will function correctly.
OK. It was very strange because I have just connected to my machine work in the University and the code was this
- Code: Select all
// We can use only a call to solve the complete system
// status = culaDeviceSgesv(m,n,cu_C22,m,cu_ipiv,cu_C21,m);
// checkStatus(status);
// Or We can use two calls
status = culaDeviceSgetrf(m,n,cu_C22,m,cu_ipiv);
checkStatus(status);
status = culaDeviceSgetrs('N',m,n,cu_C22,m,cu_ipiv,cu_C21,m);
checkStatus(status);
The problem with this is that I need 'sgetrs'. I explain this to you: several friends of mine, working in CUDA/CUBLAS, developed several routines to solve linear systems. But they only solve the LU (or Cholesky) decomposition. They do not solve the triangular systems.
Do you know when an update of this will be available?
Thank you very much
jpeinado
- jpeinado
- Posts: 37
- Joined: Mon Sep 14, 2009 10:48 am
Re:Can Sgesv divided in Sgetrf and Sgetrs?
This is now fixed in release 1.1b. My apologies for the inconvenience.
- john
- Administrator
- Posts: 587
- Joined: Thu Jul 23, 2009 2:31 pm
Re:Can Sgesv divided in Sgetrf and Sgetrs?
Thank you very much John
I will talk with the system administrator of my machine for him to install this version
jpeinado
I will talk with the system administrator of my machine for him to install this version
jpeinado
- jpeinado
- Posts: 37
- Joined: Mon Sep 14, 2009 10:48 am
Re:Can Sgesv divided in Sgetrf and Sgetrs?
Update:
Installed 1.1b and worked OK.
Thank you very much
jpeinado
Installed 1.1b and worked OK.
Thank you very much
jpeinado
- jpeinado
- Posts: 37
- Joined: Mon Sep 14, 2009 10:48 am
7 posts
• Page 1 of 1
Who is online
Users browsing this forum: No registered users and 0 guests