Using absorption length

Signed-off-by: Umberto Zerbinati <[email protected]>
NGSolve · May 29, 2024 · 184975b · 184975b
1 parent 239247a
commit 184975b
Showing 1 changed file with 2 additions and 1 deletion.
diff --git 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