Skip to content

Commit

Permalink
bugfix
Browse files Browse the repository at this point in the history
  • Loading branch information
wxj6000 committed Apr 27, 2024
1 parent 1ae630f commit c72e21d
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions pyscf/solvent/grad/smd.py
Original file line number Diff line number Diff line change
Expand Up @@ -52,7 +52,7 @@ def grad_solver(pcmobj, dm):

vK_1 = np.linalg.solve(K.T, v_grids)

dF, dA = pcm_grad.get_dD_dSt_dF_dA(pcmobj.surface)
dF, dA = pcm_grad.get_dF_dA(pcmobj.surface)
dD, dS, dSii = pcm_grad.get_dD_dS(pcmobj.surface, dF, with_D=True, with_S=True)
DA = D*A

Expand All @@ -71,7 +71,7 @@ def grad_solver(pcmobj, dm):
# dR = f_eps/(2*pi) * (dD*A + D*dA),
# dK = dS - f_eps/(2*pi) * (dD*A*S + D*dA*S + D*A*dS)
f_epsilon = (epsilon - 1.0)/(epsilon + 1.0)
fac = f_epsilon/(2.0*PI)
fac = f_epsilon/(2.0*np.pi)

Av = A*v_grids
de_dR = 0.5*fac * np.einsum('i,xij,j->ix', vK_1, dD, Av)
Expand Down

0 comments on commit c72e21d

Please sign in to comment.