Tutorials icon indicating copy to clipboard operation
Tutorials copied to clipboard

Broken Tutorial 17

Open alex-fu27 opened this issue 1 year ago • 1 comments

In Gridap 0.18, Tutorial 17 (transient equations) is broken in multiple ways:

Most importantly,

u₀ = interpolate_everywhere(0.0,U(0.0))
t₀ = 0.0
T = 10.0
uₕₜ = solve(ode_solver,op,u₀,t₀,T)

yields

MethodError: no method matching solve(::ThetaMethod, ::Gridap.ODEs.TransientFEOpFromWeakForm, ::Gridap.FESpaces.SingleFieldFEFunction{Gridap.CellData.GenericCellField{ReferenceDomain}}, ::Float64, ::Float64)

Closest candidates are:
  solve(::Gridap.ODEs.ODESolver, ::TransientFEOperator, !Matched::Real, ::Real, !Matched::CellField)
   @ Gridap ~/.julia/packages/Gridap/Ort08/src/ODEs/TransientFESolutions.jl:130
  solve(::Gridap.ODEs.ODESolver, ::TransientFEOperator, !Matched::Real, ::Real, !Matched::Tuple{Vararg{CellField}})
   @ Gridap ~/.julia/packages/Gridap/Ort08/src/ODEs/TransientFESolutions.jl:123
  solve(::Gridap.ODEs.ODESolver, !Matched::Gridap.ODEs.ODEOperator, !Matched::Real, ::Real, !Matched::AbstractVector)
   @ Gridap ~/.julia/packages/Gridap/Ort08/src/ODEs/ODESolutions.jl:142
  ...

I will attach the Jupyter Notebook of the example, showing the error message t017_poisson_transient.ipynb.txt but Github does not allow uploads with the .ipynb filename extension, so you need to remove the .txt.

Also, from the alternative FEOperator constructors, TransientConstantMatrixFEOperator seems to be missing:

UndefVarError: `TransientConstantMatrixFEOperator` not defined

alex-fu27 avatar Aug 01 '24 14:08 alex-fu27

Ah I found this very handy pull request #178 :+1:

alex-fu27 avatar Aug 01 '24 14:08 alex-fu27