Can Sgesv divided in Sgetrf and Sgetrs?

General CULA Dense (LAPACK & BLAS) support and troubleshooting. Use this forum if you are having a general problem or have encountered a bug.

Can Sgesv divided in Sgetrf and Sgetrs?

Postby jpeinado » Mon Jan 04, 2010 4:46 am

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
jpeinado
 
Posts: 37
Joined: Mon Sep 14, 2009 10:48 am

Re:Can Sgesv divided in Sgetrf and Sgetrs?

Postby john » Tue Jan 05, 2010 11:04 am

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?

Postby john » Tue Jan 05, 2010 3:50 pm

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?

Postby jpeinado » Tue Jan 05, 2010 3:58 pm

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?

Postby john » Wed Jan 06, 2010 2:22 pm

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?

Postby jpeinado » Wed Jan 06, 2010 4:37 pm

Thank you very much John

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?

Postby jpeinado » Fri Jan 08, 2010 4:44 am

Update:

Installed 1.1b and worked OK.


Thank you very much

jpeinado
jpeinado
 
Posts: 37
Joined: Mon Sep 14, 2009 10:48 am


Return to CULA Dense Support

Who is online

Users browsing this forum: No registered users and 2 guests

cron