From 184975b32b442121790896d8e62901362c68f9e1 Mon Sep 17 00:00:00 2001 From: Umberto Zerbinati Date: Wed, 29 May 2024 15:28:44 +0100 Subject: [PATCH] Using absorption length Signed-off-by: Umberto Zerbinati --- ngsPETSc/utils/firedrake/pml.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/ngsPETSc/utils/firedrake/pml.py b/ngsPETSc/utils/firedrake/pml.py index c71b880..f348280 100644 --- a/ngsPETSc/utils/firedrake/pml.py +++ b/ngsPETSc/utils/firedrake/pml.py @@ -95,7 +95,8 @@ def adiabatic_layer(self, k, deg_absorb=2): absk = fd.assemble(fd.interpolate(fd.Constant(0), self.V)) for dalet in self.dalets: absk += fd.assemble(fd.interpolate( - self.alpha*k*sigma0*(fd.real(dalet)**deg_absorb)-(sigma0*fd.real(dalet))**2, + self.alpha*k*sigma0*(fd.real(d_absorb*dalet)**deg_absorb) + -(sigma0*fd.real(d_absorb*dalet))**2, self.V)) fd.File("absk.pvd").write(absk) return absk