Skip to content

Commit 0dc39c3

Browse files
committed
Adding reproduced of issue 614 to a test.
This is a safety measure since the cause of the bug is not yet well understood.
1 parent 6e1d5c2 commit 0dc39c3

File tree

2 files changed

+36
-0
lines changed

2 files changed

+36
-0
lines changed

test/GridapTests/issue_614.jl

+34
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,34 @@
1+
module Issue614
2+
3+
using Gridap
4+
using Gridap.CellData
5+
6+
domain = (0,2,0,1)
7+
partition = (2,2)
8+
model = CartesianDiscreteModel(domain,partition)
9+
const k = 2
10+
reffeᵤ = ReferenceFE(lagrangian,VectorValue{2,Float64},k)
11+
reffe_p = ReferenceFE(lagrangian,Float64,k-1)
12+
V₀ = TestFESpace(model, reffeᵤ, conformity=:H1)
13+
U = TrialFESpace(V₀)
14+
Q = TestFESpace(model, reffe_p, conformity=:C0)
15+
P = TrialFESpace(Q)
16+
Y = MultiFieldFESpace([V₀,Q])
17+
X = MultiFieldFESpace([U,P])
18+
19+
Ω = Triangulation(model)
20+
21+
degree = 1
22+
xh = interpolate_everywhere([VectorValue(0.0,1.0),x->x[1]],X)
23+
24+
Γc = BoundaryTriangulation(model)
25+
dΓc = Measure(Γc, degree)
26+
nΓc = get_normal_vector(Γc)
27+
28+
for i in 0:1000
29+
uh, ph = xh
30+
α = ph*nΓc
31+
CD, CL = ( ( α )dΓc )
32+
end
33+
34+
end # module

test/GridapTests/runtests.jl

+2
Original file line numberDiff line numberDiff line change
@@ -28,4 +28,6 @@ using Test
2828

2929
@time @testset "PoissonLagrangeMultiplier" begin include("PoissonLagrangeMultiplierTests.jl") end
3030

31+
@time @testset "Issue614" begin include("issue_614.jl") end
32+
3133
end # module

0 commit comments

Comments
 (0)