update

PaulOlyslager · Feb 19, 2025 · 85ddc13 · 85ddc13
1 parent 5194151
commit 85ddc13
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 2 deletions.
diff --git a/src/kernels.jl b/src/kernels.jl
@@ -49,6 +49,16 @@ end
 struct HH3DInt1{T} <: BEAST.Kernel{T}
     gamma::T
 end
+function (op::HH3DInt1)(x,y)
+    gamma = op.gamma
+
+    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
+end
+
 struct HH3DInt2{T} <: BEAST.Kernel{T}
     gamma::T
 end

diff --git a/src/potentials/math_expressions.jl b/src/potentials/math_expressions.jl
@@ -8,9 +8,9 @@ function BEAST.potential(op::BEAST.PotentialIntegralOperator{D,U, <: _TupleFunct
     type=SVector{3,ComplexF64},
 	quadstrat=nothing) where {D,U,W}
     if quadstrat === nothing
-        return potential(PotentialIntegralOperatorKern(op.kernel,operator(op.op2)),points,coeffs,op.bfunc(basis);type)
+        return potential(BEAST.PotentialIntegralOperatorKern(op.kernel,operator(op.op2)),points,coeffs,op.bfunc(basis);type)
     else
-        return potential(PotentialIntegralOperatorKern(op.kernel,operator(op.op2)),points,coeffs,op.bfunc(basis);type,quadstrat)
+        return potential(BEAST.PotentialIntegralOperatorKern(op.kernel,operator(op.op2)),points,coeffs,op.bfunc(basis);type,quadstrat)
 
     end
 end