Skip to content

Commit

Permalink
bug fix
Browse files Browse the repository at this point in the history
  • Loading branch information
wxj6000 committed Apr 26, 2024
1 parent ad57ac2 commit c2081e2
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions gpu4pyscf/solvent/hessian/pcm.py
Original file line number Diff line number Diff line change
Expand Up @@ -139,10 +139,10 @@ def pcm_grad_scanner(mol):
for ix in range(3):
dv = numpy.zeros_like(coords)
dv[ia,ix] = eps
mol.set_geom_(coords + dv, unit='Bohr'
mol.set_geom_(coords + dv, unit='Bohr')
g0 = pcm_grad_scanner(mol)

mol.set_geom_(coords - dv, unit='Bohr'
mol.set_geom_(coords - dv, unit='Bohr')
g1 = pcm_grad_scanner(mol)
de[ia,:,ix] = (g0 - g1)/2.0/eps
t1 = log.timer_debug1('solvent energy', *t1)
Expand Down Expand Up @@ -176,10 +176,10 @@ def pcm_vmat_scanner(mol):
for ix in range(3):
dv = numpy.zeros_like(coords)
dv[ia,ix] = eps
mol.set_geom_(coords + dv, unit='Bohr'
mol.set_geom_(coords + dv, unit='Bohr')
vmat0 = pcm_vmat_scanner(mol)

mol.set_geom_(coords - dv, unit='Bohr'
mol.set_geom_(coords - dv, unit='Bohr')
vmat1 = pcm_vmat_scanner(mol)

grad_vmat = (vmat0 - vmat1)/2.0/eps
Expand Down

0 comments on commit c2081e2

Please sign in to comment.