Page 1 of 1

Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Mon Jan 04, 2010 4:46 am
by jpeinado
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

Re:Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Tue Jan 05, 2010 11:04 am
by john
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.

Re:Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Tue Jan 05, 2010 3:50 pm
by john
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.

Re:Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Tue Jan 05, 2010 3:58 pm
by jpeinado
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

Re:Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Wed Jan 06, 2010 2:22 pm
by john
This is now fixed in release 1.1b. My apologies for the inconvenience.

Re:Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Wed Jan 06, 2010 4:37 pm
by jpeinado
Thank you very much John

I will talk with the system administrator of my machine for him to install this version

jpeinado

Re:Can Sgesv divided in Sgetrf and Sgetrs?

PostPosted: Fri Jan 08, 2010 4:44 am
by jpeinado
Update:

Installed 1.1b and worked OK.


Thank you very much

jpeinado