LCPQ/qp2
9
1
Fork 0
mirror of https://github.com/QuantumPackage/qp2.git synced 2025-01-03 09:05:39 +01:00
Code Releases Activity

more error checking

Browse Source
This commit is contained in:
Kevin Gasperich 2020-05-07 17:27:32 -05:00
parent e2802ea5b9
commit c120ccf523
1 changed files with 24 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
src/utils/linear_algebra.irp.f
View File

@ -224,18 +224,34 @@ subroutine ortho_qr_complex(A,LDA,m,n)
call zgeqrf( m, n, A, LDA, tau, work, lwork, info ) call zgeqrf( m, n, A, LDA, tau, work, lwork, info )
lwork=int(work(1)) lwork=int(work(1))
deallocate(work) deallocate(work)
if (info.ne.0) then
print*,irp_here,' The ',-info,' argument to zgeqrf had an illegal value'
stop 1
endif
allocate(work(lwork)) allocate(work(lwork))
call zgeqrf(m, n, A, LDA, tau, work, lwork, info ) call zgeqrf(m, n, A, LDA, tau, work, lwork, info )
deallocate(work) deallocate(work)
if (info.ne.0) then
print*,irp_here,' The ',-info,' argument to zgeqrf had an illegal value'
stop 2
endif
lwork=-1 lwork=-1
allocate(work(1)) allocate(work(1))
call zungqr(m, n, n, A, LDA, tau, work, lwork, info) call zungqr(m, n, n, A, LDA, tau, work, lwork, info)
lwork=int(work(1)) lwork=int(work(1))
deallocate(work) deallocate(work)
if (info.ne.0) then
print*,irp_here,' The ',-info,' argument to zgeqrf had an illegal value'
stop 3
endif
allocate(work(lwork)) allocate(work(lwork))
call zungqr(m, n, n, A, LDA, tau, work, lwork, info) call zungqr(m, n, n, A, LDA, tau, work, lwork, info)
deallocate(work,tau) deallocate(work,tau)
if (info.ne.0) then
print*,irp_here,' The ',-info,' argument to zgeqrf had an illegal value'
stop 4
endif
end end
subroutine ortho_qr_unblocked_complex(A,LDA,m,n) subroutine ortho_qr_unblocked_complex(A,LDA,m,n)
@ -260,7 +276,15 @@ subroutine ortho_qr_unblocked_complex(A,LDA,m,n)
allocate(tau(n),work(n)) allocate(tau(n),work(n))
call zgeqr2(m,n,A,LDA,tau,work,info) call zgeqr2(m,n,A,LDA,tau,work,info)
if (info.ne.0) then
print*,irp_here,' The ',-info,' argument to zgeqr2 had an illegal value'
stop 1
endif
call zung2r(m,n,n,A,LDA,tau,work,info) call zung2r(m,n,n,A,LDA,tau,work,info)
if (info.ne.0) then
print*,irp_here,' The ',-info,' argument to zung2r had an illegal value'
stop 2
endif
deallocate(work,tau) deallocate(work,tau)
end end