### Sparse matrix eigenvalues in fortran

Posted:

**Tue Aug 22, 2017 11:44 am**Hey all

I've been trying to leverage some new GPUs for a micromagnetics project and I'm currently looking into the 'magma_zlobpcg' routine which looks promising. The project's main computation is repeated calculation of eigenvalues for a large finite element Hamiltonian alternated with modifying said Hamiltonian. The wrinkle is that I'm working in Fortran so Magma usage is somewhat indirect. I've compiled the example_f.F90 and some other stuff with 'gfortran -Dmagma_devptr_t="integer(kind=8)"...' but I'm worried whether I can access other variable types, specifically magma_z_matrix. Also of concern is the usage of CUBLAS in the example rather than native Magma routines which suggests to me they might not be readily available. Though this might very well be my untrained eye and not really a problem.

If someone could suggest a bit of pseudocode that can point me in the right direction it would be much appreciated.

I've been trying to leverage some new GPUs for a micromagnetics project and I'm currently looking into the 'magma_zlobpcg' routine which looks promising. The project's main computation is repeated calculation of eigenvalues for a large finite element Hamiltonian alternated with modifying said Hamiltonian. The wrinkle is that I'm working in Fortran so Magma usage is somewhat indirect. I've compiled the example_f.F90 and some other stuff with 'gfortran -Dmagma_devptr_t="integer(kind=8)"...' but I'm worried whether I can access other variable types, specifically magma_z_matrix. Also of concern is the usage of CUBLAS in the example rather than native Magma routines which suggests to me they might not be readily available. Though this might very well be my untrained eye and not really a problem.

If someone could suggest a bit of pseudocode that can point me in the right direction it would be much appreciated.