Skip to content

Commit

Permalink
fix
Browse files Browse the repository at this point in the history
  • Loading branch information
MartinAchondo committed Jul 2, 2024
1 parent 4603d59 commit 12fa6a4
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions xppbe/Model/Equations.py
Original file line number Diff line number Diff line change
Expand Up @@ -195,15 +195,15 @@ def __init__(self,*args,**kwargs):
def get_phi(self,X,flag,model,value='phi'):

if flag=='molecule':
slp = self.bempp.api.operators.potential.laplace.single_layer(self.neumann_space, X.numpy().transpose())
dlp = self.bempp.api.operators.potential.laplace.double_layer(self.dirichl_space, X.numpy().transpose())
slp = self.bempp.operators.potential.laplace.single_layer(self.neumann_space, X.numpy().transpose())
dlp = self.bempp.operators.potential.laplace.double_layer(self.dirichl_space, X.numpy().transpose())
phi_c,dphi_c = self.get_grid_coefficients_faces(model)
phi = slp * dphi_c - dlp * phi_c
phi = phi.reshape(-1,1) + self.G(X)

elif flag=='solvent':
slp = self.bempp.api.operators.potential.helmholtz_modified.single_layer(self.neumann_space, X.numpy().transpose(),self.kappa)
dlp = self.bempp.api.operators.potential.helmholtz_modified.double_layer(self.dirichl_space, X.numpy().transpose(),self.kappa)
slp = self.bempp.operators.potential.helmholtz_modified.single_layer(self.neumann_space, X.numpy().transpose(),self.kappa)
dlp = self.bempp.operators.potential.helmholtz_modified.double_layer(self.dirichl_space, X.numpy().transpose(),self.kappa)
phi_c,dphi_c = self.get_grid_coefficients_faces(model)
phi = slp * dphi_c - dlp * phi_c
phi = phi.reshape(-1,1)
Expand Down

0 comments on commit 12fa6a4

Please sign in to comment.