diff --git a/scripts/Embedded/Bugs/change_domain_bug.jl b/scripts/Embedded/Bugs/change_domain_bug.jl index 11aea81f..a2de91a6 100644 --- a/scripts/Embedded/Bugs/change_domain_bug.jl +++ b/scripts/Embedded/Bugs/change_domain_bug.jl @@ -12,23 +12,13 @@ model = ref_model.model Ω = Triangulation(model) dΩ = Measure(Ω,2*order) V_φ = TestFESpace(model,ReferenceFE(lagrangian,Float64,order)) -Λ = SkeletonTriangulation(Ω) -dΛ = Measure(Λ,2order) φh = interpolate(x->-cos(4π*x[1])*cos(4π*x[2])-0.2,V_φ) -x_φ = get_free_dof_values(φh) -idx = findall(isapprox(0.0;atol=10^-10),x_φ) -@assert isempty(idx) - -Ωs = EmbeddedCollection(model,φh) do cutgeo - Γ = DifferentiableTriangulation(EmbeddedBoundary(cutgeo)) - (; - :Γ => Γ, - :dΓ => Measure(Γ,2*order) - ) -end +geo = DiscreteGeometry(φh,model) +cutgeo = cut(model,geo) +Γ = DifferentiableTriangulation(EmbeddedBoundary(cutgeo)) +dΓ = Measure(Γ,2order) vh = zero(V_φ) - -∫(vh)Ωs.dΓ -∫(vh*vh)Ωs.dΓ \ No newline at end of file +∫(vh)dΓ +∫(vh*vh)dΓ \ No newline at end of file