debugging Forward Euler

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