18 REAL(KIND=8),
DIMENSION(3,3) :: A,B,sqB,sqA
22 a(1,:)=(/9.1784812096335050e-009,5.7810098811039660e-009,3.8077584903965249e-009/)
23 a(2,:)=(/5.7810098811039676e-009,3.6411334819038273e-009,2.3982932421037010e-009/)
24 a(3,:)=(/3.8077584903965258e-009,2.3982932421037010e-009,1.5796758080159289e-009/)
30 print*,
'Schur decomposition' 32 CALL sqrtm(a,sqa,info,info2,2)
45 print*,
'Svd Decomposition'
subroutine, public sqrtm(A, sqA, info, info_triu, bs)
Routine to compute a real square-root of a matrix.
subroutine, public printmat(A)
subroutine, public init_sqrt
program test_sqrtm
Small program to test the matrix square-root module.
The model parameters module.
subroutine init_params
Parameters initialisation routine.
Utility module with various routine to compute matrix square root.
subroutine, public sqrtm_svd(A, sqA, info, info_triu, bs)
Routine to compute a real square-root of a matrix via a SVD decomposition.