fix type check compat in matrixUtils

cmbant · Sep 18, 2019 · f573241 · f573241
1 parent 7ccba81
commit f573241
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/MatrixUtils.f90 b/MatrixUtils.f90
@@ -1710,7 +1710,7 @@ subroutine Matrix_SVD_VT(Mat,m, n, D, U)
     integer WorkSize, ierr
     integer,allocatable, dimension (:) :: IWork
     real(dm), allocatable, dimension (:) :: rv1
-    real(dm) OptWk
+    real(dm) OptWk(1)
 
     if (n<=m) call MpiStop('Matrix_SVD_VT assumed n>m. ')
 
@@ -1725,7 +1725,7 @@ subroutine Matrix_SVD_VT(Mat,m, n, D, U)
 #else
         call DGESDD('O',m,n, Mat, m ,D,U,m,Mat,n,OptWk,WorkSize,IWork,ierr)
 #endif
-        WorkSize = nint(OptWk)
+        WorkSize = nint(OptWk(1))
         allocate(rv1(WorkSize))
 #ifdef MATRIX_SINGLE
         call SGESDD('O',m,n, Mat, m ,D,U,m,Mat,n,rv1,WorkSize,IWork,ierr)