comp-chem/NWChem
1
0
Fork 0
mirror of https://github.com/nwchemgit/nwchem.git synced 2026-07-28 06:05:44 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 1

Use DSYRK instead of DGEMM

Browse source
This commit is contained in:
Daniel Mejia-Rodriguez 2023-09-28 21:40:17 -07:00 • committed by edoapra
parent edb236eba2
commit a94cd1ecdb
No known key found for this signature in database
GPG key ID: 9E6A0B70826967BA
1 changed files with 1 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/geom/geom_hnd.F
View file

@ -308,16 +308,7 @@ c
write(LuOut,*) ' B matrix '
call output(x, 1, ncart, 1, nzvar, ncart, nzvar, 1)
endif
call dgemm('t', 'n', nzvar, nzvar, ncart, 1d0,
$ x, ncart, x, ncart, 0d0, y, nzvar) ! BT.B -> y ... should symmetrize
do i = 1, nzvar
do j = 1, nzvar
ij = i + (j-1)*nzvar
ji = j + (i-1)*nzvar
y(ij) = 0.5d0*(y(ij) + y(ji))
y(ji) = y(ij)
end do
end do
call dsyrk('u', 't', nzvar, ncart, 1d0, x, ncart, 0d0, y, nzvar)
call dsyev('v', 'u', nzvar, y, nzvar, e, z, ndim**2, info)
if (info .ne. 0) call errquit('geom_binvr: dsyev failed', info,
& MEM_ERR)