GridapROMs.ParamODEs

const JointODEParamOperator{O<:ODEParamOperatorType} = ODEParamOperator{O,JointDomains}

const JointTransientParamFEOperator{O<:ODEParamOperatorType} = TransientParamFEOperator{O,JointDomains}

struct LinearNonlinearParamODE <: ODEParamOperatorType end

struct LinearParamODE <: ODEParamOperatorType end

struct NonlinearParamODE <: ODEParamOperatorType end

const ODEParamOperator{T<:ODEParamOperatorType,T<:TriangulationStyle} <: ParamOperator{O,T}

Transient extension of the type ParamOperator.

abstract type ODEParamOperatorType <: UnEvalOperatorType end

Parametric extension of the type ODEOperatorType in Gridap.

Subtypes:

• NonlinearParamODE
• LinearParamODE
• LinearNonlinearParamODE

struct ODEParamSolution{V} <: ODESolution
  solver::ODESolver
  odeop::ODEParamOperator
  r::TransientRealization
  us0::Tuple{Vararg{V}}
end

struct ParamStageOperator{O} <: NonlinearParamOperator
  op::ODEParamOperator{O}
  r::TransientRealization
  state_update::Function
  ws::Tuple{Vararg{Real}}
  paramcache::AbstractParamCache
end

Stage operator to solve a parametric ODE with a time marching scheme

const SplitODEParamOperator{O<:ODEParamOperatorType} = ODEParamOperator{O,SplitDomains}

const SplitTransientParamFEOperator{O<:ODEParamOperatorType} = TransientParamFEOperator{O,SplitDomains}

const TransientMultiFieldParamFESpace = MultiFieldFESpace

struct TransientParamFEOpFromWeakForm{T} <: TransientParamFEOperator{NonlinearParamODE,T}
  res::Function
  jacs::Tuple{Vararg{Function}}
  tpspace::TransientParamSpace
  assem::Assembler
  trial::FESpace
  test::FESpace
  domains::FEDomains
  order::Integer
end

Instance of TransientParamFEOperator, to be used when the transient problem is nonlinear

const TransientParamFEOperator{O<:ODEParamOperatorType,T<:TriangulationStyle} = ParamFEOperator{O,T}

Parametric extension of a TransientFEOperator in Gridap. Compared to a standard TransientFEOperator, there are the following novelties:

• a TransientParamSpace is provided, so that parametric realizations can be extracted directly from the TransientParamFEOperator
• a function representing a norm matrix is provided, so that errors in the desired norm can be automatically computed

Subtypes:

• TransientParamFEOpFromWeakForm
• TransientParamLinearFEOpFromWeakForm

struct TransientParamFESolution{V} <: TransientFESolution
  odesol::ODEParamSolution{V}
  trial
end

Wrapper around a TransientParamFEOperator and ODESolver that represents the parametric solution at a set of time steps. It is an iterator that computes the solution at each time step in a lazy fashion when accessing the solution.

struct TransientParamLinearFEOpFromWeakForm{T} <: TransientParamFEOperator{LinearParamODE,T}
  res::Function
  jacs::Tuple{Vararg{Function}}
  constant_forms::Tuple{Vararg{Bool}}
  tpspace::TransientParamSpace
  assem::Assembler
  trial::FESpace
  test::FESpace
  domains::FEDomains
  order::Integer
end

Instance of TransientParamFEOperator, to be used when the transient problem is linear

const TransientTrialParamFESpace = UnEvalTrialFESpace

TransientParamLinearFEOperator(forms::Tuple{Vararg{Function}},res::Function, tpspace,trial,test;kwargs...) -> TransientParamLinearFEOpFromWeakForm{TriangulationStyle}

Returns a linear parametric FE operator