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

avoid matmul nesting with large matrices

Browse Source
This commit is contained in:
Loris Burth 2025-04-02 15:50:31 +02:00
parent c879a461b6
commit 8a3c738135
1 changed files with 2 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/utils/complex_orthogonalization.f90
View File

@ -38,13 +38,8 @@ 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)
L = matmul(transpose(vectors),A)
L = matmul(L,vectors)
call complex_cholesky_decomp(N,L)
call complex_inverse_matrix(N,L,Linv)
vectors = matmul(vectors,transpose(Linv))