make a_copy and b_copy allocatable.

gha3mi · Dec 6, 2023 · 3b31a4f · 3b31a4f
1 parent 8ddf79b
commit 3b31a4f
Showing 1 changed file with 4 additions and 4 deletions.
diff --git 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