#include <MoFEM.hpp>
#include <boost/math/quadrature/gauss_kronrod.hpp>

Classes
struct	ElementsAndOps< DIM >

struct	ElementsAndOps< 2 >

struct	OpURhs

struct	OpULhs_dU

struct	OpULhs_dH

struct	Example
	[Example] More...

struct	Monitor
	[Push operators to pipeline] More...

Typedefs
using	EntData = DataForcesAndSourcesCore::EntData

using	DomainEle = ElementsAndOps< FE_DIM >::DomainEle

using	DomainEleOp = ElementsAndOps< FE_DIM >::DomainEleOp

using	PostProcEle = ElementsAndOps< FE_DIM >::PostProcEle

using	AssemblyDomainEleOp = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::OpBase

using	OpBaseDivU = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMixScalarTimesDiv< 3 >

using	OpBaseGradH = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMixVectorTimesGrad< 1, 3, 3 >

using	OpMassUU = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, 3 >

using	OpMassHH = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, 1 >

using	OpBaseTimesDotH = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpBaseTimesScalarField< 1 >

using	OpBaseTimesDivU = OpBaseTimesDotH

using	OpSourceU = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< 1, 3 >

using	OpSourceH = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< 1, 1 >

using	OpConvectiveH = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpConvectiveTermRhs< 1, 1, 3 >

using	OpConvectiveH_dU = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpConvectiveTermLhsDu< 1, 1, 3 >

using	OpConvectiveH_dGradH = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpConvectiveTermLhsDy< 1, 1, 3 >

Functions
int	main (int argc, char *argv[])

Variables
constexpr int	FE_DIM = 2

constexpr double	omega = 7.292 * 1e-5

constexpr double	g = 9.80616

constexpr double	mu = 1e4

constexpr double	h0 = 1e4

constexpr double	h_hat = 120

constexpr double	u_max = 80

constexpr double	phi_0 = M_PI / 7

constexpr double	phi_1 = M_PI / 2 - phi_0

constexpr double	phi_2 = M_PI / 4

constexpr double	alpha_montain = 1. / 3.

constexpr double	beta_montain = 1. / 15.

constexpr double	penalty = 1

FTensor::Index< 'i', 3 >	i

FTensor::Index< 'j', 3 >	j

FTensor::Index< 'l', 3 >	l

FTensor::Index< 'm', 3 >	m

static char	help [] = "...\n\n"
	[Check] More...

Typedef Documentation

◆ AssemblyDomainEleOp

using AssemblyDomainEleOp = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::OpBase

Definition at line 46 of file shallow_wave.cpp.

◆ DomainEle

using DomainEle = ElementsAndOps<FE_DIM>::DomainEle

Definition at line 42 of file shallow_wave.cpp.

◆ DomainEleOp

using DomainEleOp = ElementsAndOps<FE_DIM>::DomainEleOp

Definition at line 43 of file shallow_wave.cpp.

◆ EntData

using EntData = DataForcesAndSourcesCore::EntData

Examples: shallow_wave.cpp.

Definition at line 41 of file shallow_wave.cpp.

◆ OpBaseDivU

using OpBaseDivU = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpMixScalarTimesDiv<3>

Examples: shallow_wave.cpp.

Definition at line 50 of file shallow_wave.cpp.

◆ OpBaseGradH

using OpBaseGradH = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpMixVectorTimesGrad<1, 3, 3>

Definition at line 52 of file shallow_wave.cpp.

◆ OpBaseTimesDivU

using OpBaseTimesDivU = OpBaseTimesDotH

Examples: shallow_wave.cpp.

Definition at line 63 of file shallow_wave.cpp.

◆ OpBaseTimesDotH

using OpBaseTimesDotH = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpBaseTimesScalarField<1>

Examples: shallow_wave.cpp.

Definition at line 61 of file shallow_wave.cpp.

◆ OpConvectiveH

using OpConvectiveH = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpConvectiveTermRhs<1, 1, 3>

Examples: shallow_wave.cpp.

Definition at line 70 of file shallow_wave.cpp.

◆ OpConvectiveH_dGradH

using OpConvectiveH_dGradH = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpConvectiveTermLhsDy<1, 1, 3>

Examples: shallow_wave.cpp.

Definition at line 74 of file shallow_wave.cpp.

◆ OpConvectiveH_dU

using OpConvectiveH_dU = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpConvectiveTermLhsDu<1, 1, 3>

Examples: shallow_wave.cpp.

Definition at line 72 of file shallow_wave.cpp.

◆ OpMassHH

using OpMassHH = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpMass<1, 1>

Examples: shallow_wave.cpp.

Definition at line 58 of file shallow_wave.cpp.

◆ OpMassUU

using OpMassUU = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpMass<1, 3>

Examples: shallow_wave.cpp.

Definition at line 56 of file shallow_wave.cpp.

◆ OpSourceH

using OpSourceH = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm< GAUSS>::OpSource<1, 1>

Examples: shallow_wave.cpp.

Definition at line 67 of file shallow_wave.cpp.

◆ OpSourceU

using OpSourceU = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm< GAUSS>::OpSource<1, 3>

Examples: shallow_wave.cpp.

Definition at line 65 of file shallow_wave.cpp.

◆ PostProcEle

using PostProcEle = ElementsAndOps<FE_DIM>::PostProcEle

Definition at line 44 of file shallow_wave.cpp.

Function Documentation

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

Examples: shallow_wave.cpp.

Definition at line 862 of file shallow_wave.cpp.

                                  {
  
   // Initialisation of MoFEM/PETSc and MOAB data structures
   const char param_file[] = "param_file.petsc";
   MoFEM::Core::Initialize(&argc, &argv, param_file, help);
  
   // Add logging channel for example
   auto core_log = logging::core::get();
   core_log->add_sink(LogManager::createSink(LogManager::getStrmWorld(), "SW"));
   LogManager::setLog("SW");
   MOFEM_LOG_TAG("SW", "example");
  
   try {
  
     //! [Register MoFEM discrete manager in PETSc]
     DMType dm_name = "DMMOFEM";
     CHKERR DMRegister_MoFEM(dm_name);
     //! [Register MoFEM discrete manager in PETSc
  
     //! [Create MoAB]
     moab::Core mb_instance;              ///< mesh database
     moab::Interface &moab = mb_instance; ///< mesh database interface
     //! [Create MoAB]
  
     //! [Create MoFEM]
     MoFEM::Core core(moab);           ///< finite element database
     MoFEM::Interface &m_field = core; ///< finite element database insterface
     //! [Create MoFEM]
  
     //! [Example]
     Example ex(m_field);
     CHKERR ex.runProblem();
     //! [Example]
   }
   CATCH_ERRORS;
  
   CHKERR MoFEM::Core::Finalize();
 }