Commit a2360e2fa76492529135e8ec0af87982f629c3ad
1 parent
b5f099f3cf
Exists in
master
and in
1 other branch
Correction d'un bug apparu avec l'update d'openblas
le tableau rwork doit être alloué lors de l'appel query
Showing 1 changed file with 6 additions and 5 deletions Side-by-side Diff
fvn_linear/fvn_linear.f90
... | ... | @@ -1281,12 +1281,12 @@ |
1281 | 1281 | !call ccopy(d*d,a,1,wc_a,1) |
1282 | 1282 | wc_a(:,:)=a(:,:) |
1283 | 1283 | |
1284 | - | |
1284 | + ! rwork must be allocated before query | |
1285 | + allocate(rwork(2*d)) | |
1285 | 1286 | ! query optimal work size |
1286 | 1287 | call cgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,twork,-1,rwork,info) |
1287 | 1288 | lwork=int(twork(1)) |
1288 | - allocate(work(lwork)) | |
1289 | - allocate(rwork(2*d)) | |
1289 | + allocate(work(lwork)) | |
1290 | 1290 | call cgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,work,lwork,rwork,info) |
1291 | 1291 | |
1292 | 1292 | if (info /= 0) then |
1293 | 1293 | |
... | ... | @@ -1335,11 +1335,12 @@ |
1335 | 1335 | !call zcopy(d*d,a,1,wc_a,1) |
1336 | 1336 | wc_a(:,:)=a(:,:) |
1337 | 1337 | |
1338 | + ! rwork must be allocated before query | |
1339 | + allocate(rwork(2*d)) | |
1338 | 1340 | ! query optimal work size |
1339 | 1341 | call zgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,twork,-1,rwork,info) |
1340 | 1342 | lwork=int(twork(1)) |
1341 | - allocate(work(lwork)) | |
1342 | - allocate(rwork(2*d)) | |
1343 | + allocate(work(lwork)) | |
1343 | 1344 | call zgeev('N','V',d,wc_a,d,evala,vl,1,eveca,d,work,lwork,rwork,info) |
1344 | 1345 | |
1345 | 1346 | if (info /= 0) then |