update

PaulOlyslager · Feb 25, 2025 · 3fe6318 · 3fe6318
1 parent 846f703
commit 3fe6318
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/kernels.jl b/src/kernels.jl
@@ -55,7 +55,7 @@ function (op::HH3DInt1)(x,y)
     r = cartesian(x) - cartesian(y)
     R = norm(r)
     iR = 1/R
-    green = -1/gamma^2*(expm1(-gamma*R)+gamma*R)*(i4pi*iR^2)*r/R
+    green = -1/gamma^2*(expm1(-gamma*R)+gamma*R*exp(-gamma*R))*(i4pi*iR^2)*r/R
     green
 end
 

diff --git a/src/potentials/potentials.jl b/src/potentials/potentials.jl
@@ -22,7 +22,7 @@ end
 function potential(op::LinearCombinationOfPotentials,points, coeffs, basis; 
     type=SVector{3,ComplexF64},
 	quadstrat=nothing)
-    if quadstrat == nothing
+    if quadstrat === nothing
         return sum(ci*potential(opi,points,coeffs,basis,type=type) for (opi,ci) in zip(op.potentials,op.coefficients))
     else
         return sum(ci*potential(opi,points,coeffs,basis,type=type,quadstrat=quadstrat) for (opi,ci) in zip(op.potentials,op.coefficients))