diff --git a/ProcessLib/GroundwaterFlowProcess.h b/ProcessLib/GroundwaterFlowProcess.h index 91a7fc14433f15a584a9770b4aba4680edec8667..f42a16a2bddb7c909c820fb9f6d7f72f3b9e03d7 100644 --- a/ProcessLib/GroundwaterFlowProcess.h +++ b/ProcessLib/GroundwaterFlowProcess.h @@ -68,6 +68,26 @@ public: } } + //! \name ODESystem interface + //! @{ + + bool isLinear() const override + { + return true; + } + + //! @} + +private: + Parameter<double, MeshLib::Element const&> const& _hydraulic_conductivity; + + using LocalAssembler = GroundwaterFlow::LocalAssemblerDataInterface< + typename GlobalSetup::MatrixType, typename GlobalSetup::VectorType>; + + using GlobalAssembler = AssemblerLib::VectorMatrixAssembler< + GlobalMatrix, GlobalVector, LocalAssembler, + NumLib::ODESystemTag::FirstOrderImplicitQuasilinear>; + template <unsigned GlobalDim> void createLocalAssemblers(AssemblerLib::LocalToGlobalIndexMap const& dof_table, MeshLib::Mesh const& mesh, @@ -121,27 +141,6 @@ public: assert(false); } - //! \name ODESystem interface - //! @{ - - bool isLinear() const override - { - return true; - } - - //! @} - -private: - Parameter<double, MeshLib::Element const&> const& _hydraulic_conductivity; - - using LocalAssembler = GroundwaterFlow::LocalAssemblerDataInterface< - typename GlobalSetup::MatrixType, typename GlobalSetup::VectorType>; - - using GlobalAssembler = AssemblerLib::VectorMatrixAssembler< - GlobalMatrix, GlobalVector, LocalAssembler, - NumLib::ODESystemTag::FirstOrderImplicitQuasilinear>; - - void assembleConcreteProcess(const double t, GlobalVector const& x, GlobalMatrix& M, GlobalMatrix& K, GlobalVector& b) override {