From a2360e2fa76492529135e8ec0af87982f629c3ad Mon Sep 17 00:00:00 2001 From: William Daniau Date: Mon, 16 Sep 2019 16:36:56 +0200 Subject: [PATCH] =?UTF-8?q?Correction=20d'un=20bug=20apparu=20avec=20l'upd?= =?UTF-8?q?ate=20d'openblas=20le=20tableau=20rwork=20doit=20=C3=AAtre=20al?= =?UTF-8?q?lou=C3=A9=20lors=20de=20l'appel=20query?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- fvn_linear/fvn_linear.f90 | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/fvn_linear/fvn_linear.f90 b/fvn_linear/fvn_linear.f90 index 52a8acf..3ec48fd 100644 --- a/fvn_linear/fvn_linear.f90 +++ b/fvn_linear/fvn_linear.f90 @@ -1281,12 +1281,12 @@ subroutine fvn_c_matev(d,a,evala,eveca,status,sortval) !call ccopy(d*d,a,1,wc_a,1) wc_a(:,:)=a(:,:) - + ! rwork must be allocated before query + allocate(rwork(2*d)) ! query optimal work size call cgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,twork,-1,rwork,info) lwork=int(twork(1)) - allocate(work(lwork)) - allocate(rwork(2*d)) + allocate(work(lwork)) call cgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,work,lwork,rwork,info) if (info /= 0) then @@ -1335,11 +1335,12 @@ subroutine fvn_z_matev(d,a,evala,eveca,status,sortval) !call zcopy(d*d,a,1,wc_a,1) wc_a(:,:)=a(:,:) + ! rwork must be allocated before query + allocate(rwork(2*d)) ! query optimal work size call zgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,twork,-1,rwork,info) lwork=int(twork(1)) - allocate(work(lwork)) - allocate(rwork(2*d)) + allocate(work(lwork)) call zgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,work,lwork,rwork,info) if (info /= 0) then -- 2.16.4