minor

zjwegert · Nov 15, 2024 · c74c61a · c74c61a
1 parent f9d848f
commit c74c61a
Showing 1 changed file with 9 additions and 4 deletions.
diff --git a/scripts/Embedded/Bugs/change_domain_bug.jl b/scripts/Embedded/Bugs/change_domain_bug.jl
@@ -12,13 +12,18 @@ model = ref_model.model
 Ω = Triangulation(model)
 dΩ = Measure(Ω,2*order)
 V_φ = TestFESpace(model,ReferenceFE(lagrangian,Float64,order))
-
 φh = interpolate(x->-cos(4π*x[1])*cos(4π*x[2])-0.2,V_φ)
+vh = zero(V_φ)
+
 geo = DiscreteGeometry(φh,model)
 cutgeo = cut(model,geo)
-Γ = DifferentiableTriangulation(EmbeddedBoundary(cutgeo))
+
+Γ = EmbeddedBoundary(cutgeo)
 dΓ = Measure(Γ,2order)
+∫(vh)dΓ
+∫(2vh)dΓ
 
-vh = zero(V_φ)
+Γ = DifferentiableTriangulation(EmbeddedBoundary(cutgeo))
+dΓ = Measure(Γ,2order)
 ∫(vh)dΓ
-∫(vh*vh)dΓ
+∫(2vh)dΓ