MAGMA Forum

Mark Thank you, thank you for the NB. That is indeed a bug in my code and I am astonished it was working before. I will do some more exploration as to the other aspect of my problem. I cannot easily print out the matrix data as it is about 7000 square! This is part of a calculation to locate a small...

This is a question about the assumptions which can be made about the lifetime of memory on a GPU set by calls to cublas_set_matrix. I have been converting some Fortran code from CPU (Lapack) to GPU operation (MAGMA). In the code I call magmaf_zgetrs_gpu twice using the same factorization to solve tw...

I have now implemented a version I have called magma_zcgesv_opt_gpu which gives the option to do either the factorization or the solution or both.

I have it working from Fortran using a modification of the interface code from the other thread.

Thank you for your help.

John

I have built the attached Fortran program testing_zcgesv_gpu_f by extending the Fortran one for zgetrf_gpu with material suggested by the C++ one for zcgesv_gpu. This routine does a similar calculation to the one in the Fortran one for zgetrf_gpu with magma_zcgesv_gpu as well and prints both sets of...

Mark

Thank you for this. I will give these a try out.

John

Mark

Thank you for the quick response. I have had a look at these files. Presumably I just compile and link these files into my code.

John

P.S. I am currently using Magma 1.6.2

I am wishing to run some routines such as magma_zcgesv from FORTRAN. There are wrappers for quite a lot of routines but not for this one. I can see that the control subdirectory contains the files magma_zfortran.F90 and magma_zf77.cpp and it looks as though I can use some of the code there as a mode...

When I wrote the above I had not yet looked at the Magma version of ZCGESV. I have now done so, and find that it does allow for transpose operations in its interface, unlike the LAPACK one. It has another difference from the LAPACK version as follows. It does the SP factorization using magma_cgetrf_...

I am working with various problems which require the repeated solution of N sets of equations Ax = b where there are N copies of A which don't change from one iteration to the next while the b blocks do change. One way of doing this is to run ZGETRF on each A and save the results, so that ZGETRS can...

Hartwig Thank you, I have thought of that. I do need to do that eventually. At the moment I have two computers, chosen to have the same graphics card (GTX 460) when they were new. Both are on Ubuntu 12.04. One dual boots to Windows 7 as well. At the moment the dual boot one will only do low res grap...