Skip to content
This repository has been archived by the owner on Feb 28, 2022. It is now read-only.

Commit

Permalink
Adapted README.md snipped
Browse files Browse the repository at this point in the history
  • Loading branch information
oriolcg committed Dec 9, 2021
1 parent 3044f0f commit cdbec51
Showing 1 changed file with 3 additions and 6 deletions.
9 changes: 3 additions & 6 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -14,9 +14,6 @@ using ForwardDiff
using GridapODEs.ODETools
using GridapODEs.TransientFETools

import Gridap:
import GridapODEs.TransientFETools: ∂t

θ = 0.5

u(x,t) = (1.0-x[1])*x[1]*(1.0-x[2])*x[2]*t
Expand Down Expand Up @@ -48,9 +45,9 @@ a(u,v) = ∫( ∇(v)⋅∇(u) )dΩ
b(v,t) = ( v*f(t) )dΩ
m(u,v) = ( v*u )dΩ

res(t,(u,ut),v) = a(u,v) + m(ut,v) - b(v,t)
jac(t,(u,ut),du,v) = a(du,v)
jac_t(t,(u,ut),dut,v) = m(dut,v)
res(t,u,v) = a(u,v) + m(∂t(u),v) - b(v,t)
jac(t,u,du,v) = a(du,v)
jac_t(t,u,dut,v) = m(dut,v)

op = TransientFEOperator(res,jac,jac_t,U,V0)

Expand Down

0 comments on commit cdbec51

Please sign in to comment.