diff --git a/test/TransientFEsTests/ForwardEulerHeatEquationTests.jl b/test/TransientFEsTests/ForwardEulerHeatEquationTests.jl index 63427bc..af78fe1 100644 --- a/test/TransientFEsTests/ForwardEulerHeatEquationTests.jl +++ b/test/TransientFEsTests/ForwardEulerHeatEquationTests.jl @@ -1,4 +1,4 @@ -module HeatEquationTests +module ForwardEulerHeatEquationTests using Gridap using ForwardDiff @@ -30,25 +30,29 @@ model = CartesianDiscreteModel(domain,partition) order = 2 +reffe = ReferenceFE(lagrangian,Float64,order) V0 = FESpace( - reffe=lagrangian, order=order, valuetype=Float64, - conformity=:H1, model=model, dirichlet_tags="boundary") + model, + reffe, + conformity=:H1, + dirichlet_tags="boundary" +) U = TransientTrialFESpace(V0,u) -trian = Triangulation(model) +Ω = Triangulation(model) degree = 2*order -quad = CellQuadrature(trian,degree) +dΩ = Measure(Ω,degree) # a(u,v) = ∇(v)⋅∇(u) b(v,t) = v*f(t) -res(t,(u,ut),v) = a(u,v) + ut*v - b(v,t) +res(t,(u,ut),v) = ∫( a(u,v) + ut*v - b(v,t) )dΩ jac(t,(u,ut),du,v) = a(du,v) -jac_t(t,(u,ut),dut,v) = dut*v +jac_t(t,(u,ut),dut,v) = ∫( dut*v )dΩ -t_Ω = FETerm(res,jac,jac_t,trian,quad) -op = TransientFEOperator(U,V0,t_Ω) +#t_Ω = FETerm(res,jac,jac_t,trian,quad) +op = TransientFEOperator(res,jac,jac_t,U,V0) t0 = 0.0 tF = 1.0