From 3b31a4f0e06e1209cf842ba439c0af74c4a4450e Mon Sep 17 00:00:00 2001 From: Seyed Ali Ghasemi Date: Wed, 6 Dec 2023 20:29:43 +0100 Subject: [PATCH] make a_copy and b_copy allocatable. --- src/forsolver.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/forsolver.f90 b/src/forsolver.f90 index 1527bd8..5f15cb2 100644 --- a/src/forsolver.f90 +++ b/src/forsolver.f90 @@ -92,8 +92,8 @@ pure function dgesv_rel(A, b) result(x) integer :: info ! result info integer :: n, lda, ldb integer, dimension(size(A, 2)) :: ipiv - real(rk), dimension(size(A,1),size(A,2)) :: a_copy - real(rk), dimension(size(b)) :: b_copy + real(rk), dimension(:,:), allocatable :: a_copy + real(rk), dimension(:), allocatable :: b_copy ! interface for dgels subroutine interface @@ -140,8 +140,8 @@ pure function dgels_rel(A, b) result(x) integer :: m, n, lda, ldb, lwork real(rk), allocatable :: work(:) real(rk) :: work1(1) - real(rk), dimension(size(A,1),size(A,2)) :: a_copy - real(rk), dimension(size(b)) :: b_copy + real(rk), dimension(:,:), allocatable :: a_copy + real(rk), dimension(:), allocatable :: b_copy ! interface for dgels subroutine interface