Skip to content

Broken Tutorial 17 #181

@alex-fu27

Description

@alex-fu27

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions