5 ! Contributed by Andres Legarra.
10 INTEGER, PARAMETER :: dp
= KIND(1.0D0
)
11 END MODULE LA_PRECISION
15 SUBROUTINE DGESV_F90( A
, B
, IPIV
, INFO
)
16 USE la_precision
, ONLY
: wp
=> dp
18 INTEGER, INTENT(OUT
), OPTIONAL
:: INFO
19 INTEGER, INTENT(OUT
), OPTIONAL
, TARGET
:: IPIV(:)
20 REAL(WP
), INTENT(IN OUT
) :: A(:,:), B(:,:)
21 END SUBROUTINE DGESV_F90
25 SUBROUTINE DGESV_F90( A
, B
, IPIV
, INFO
)
26 USE la_precision
, ONLY
: wp
=> dp
28 INTEGER, INTENT(OUT
), OPTIONAL
:: INFO
29 INTEGER, INTENT(OUT
), OPTIONAL
, TARGET
:: IPIV(:)
30 REAL(WP
), INTENT(IN OUT
) :: A(:,:), B(:,:)
31 END SUBROUTINE DGESV_F90
36 integer, parameter :: r8
= SELECTED_REAL_KIND( 15, 307 )
37 real(r8
)::denseop_tol
=1.d
-50
43 real(r8
),allocatable
::x_o(:,:)
44 allocate(x_o(size(x
,1),size(x
,1)))