eigenvalues and eigenvectors for Symmetry/Hermitian Matrix
gemm
subroutine
gemm
gemv
subroutine
gemv
eye
subroutine
eye
print
subroutine
print matrix
geut
subroutine
$U^{\dagger}AU$ or $UAU^{\dagger}$
equal
function
equal
program check
use linalg_fortran
real(8)::a(3,3)
real(4)::b(3,3)
callrandom_number(a)
b=a
if(linalg%inv(a))thencall linalg%print(a,"a")
end ifif(linalg%inv(b))thencall linalg%print(b,"b")
end ifend program check
Given the name of the module iso_linalg_mod I wonder if the Fortran standards committee is thinking of adding such a module to the language. Is there a paper about this?
It might be nice to use the preprocessor to support single precision. Users can enable the preprocessor in the top-level fpm project according to their own needs, and choose to use iso_linalg_mod with single or double precision.