From c879a461b6e19aa95686381cc8047ba72c984266 Mon Sep 17 00:00:00 2001 From: Loris Burth Date: Wed, 2 Apr 2025 15:42:39 +0200 Subject: [PATCH] debugging --- src/utils/complex_orthogonalization.f90 | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/src/utils/complex_orthogonalization.f90 b/src/utils/complex_orthogonalization.f90 index bb1fea4..215d482 100644 --- a/src/utils/complex_orthogonalization.f90 +++ b/src/utils/complex_orthogonalization.f90 @@ -38,7 +38,12 @@ subroutine complex_orthonormalize(N,vectors,A) ! Copy the input matrix to a temporary matrix allocate(L(N,N),Linv(N,N)) - + write(*,*) "Size A" + write(*,*) size(A,1) + write(*,*) size(A,2) + write(*,*) "Size VECTORS" + write(*,*) size(vectors,1) + write(*,*) size(vectors,2) L = matmul(matmul(transpose(vectors),A),vectors) call complex_cholesky_decomp(N,L) call complex_inverse_matrix(N,L,Linv)