From f573241e8ee2f27d51676f2909660d4445280cae Mon Sep 17 00:00:00 2001
From: Antony Lewis <antony@cosmologist.info>
Date: Wed, 18 Sep 2019 17:20:55 +0100
Subject: [PATCH] fix type check compat in matrixUtils

---
 MatrixUtils.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/MatrixUtils.f90 b/MatrixUtils.f90
index 8ca3f5e..6aec732 100644
--- 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)