PTEROSOR/QuAcK
10
1
Fork 0
mirror of https://github.com/pfloos/quack synced 2025-05-06 23:24:58 +02:00
Code Releases Activity

debugging

Browse Source
This commit is contained in:
Loris Burth 2025-04-02 15:42:39 +02:00
parent 3d3807dc34
commit c879a461b6
1 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/utils/complex_orthogonalization.f90
View File

@ -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)