Can Sgesv divided in Sgetrf and Sgetrs?

[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

[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.

[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.

[jpeinado]

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

[john]

This is now fixed in release 1.1b. My apologies for the inconvenience.

[jpeinado]

Thank you very much John

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

jpeinado

[jpeinado]

Update:

Installed 1.1b and worked OK.


Thank you very much

jpeinado