diff --git a/scripts/Embedded/Bugs/change_domain_bug.jl b/scripts/Embedded/Bugs/change_domain_bug.jl index a2de91a..9977b7e 100644 --- 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Γ \ No newline at end of file +∫(2vh)dΓ \ No newline at end of file