MFEM
v4.3.0
Finite element discretization library
|
#include <bilininteg.hpp>
Public Member Functions | |
DGTraceIntegrator (VectorCoefficient &u_, double a) | |
Construct integrator with rho = 1, b = 0.5*a. More... | |
DGTraceIntegrator (VectorCoefficient &u_, double a, double b) | |
Construct integrator with rho = 1. More... | |
DGTraceIntegrator (Coefficient &rho_, VectorCoefficient &u_, double a, double b) | |
virtual void | AssembleFaceMatrix (const FiniteElement &el1, const FiniteElement &el2, FaceElementTransformations &Trans, DenseMatrix &elmat) |
virtual void | AssemblePAInteriorFaces (const FiniteElementSpace &fes) |
virtual void | AssemblePABoundaryFaces (const FiniteElementSpace &fes) |
virtual void | AddMultTransposePA (const Vector &x, Vector &y) const |
Method for partially assembled transposed action. More... | |
virtual void | AddMultPA (const Vector &, Vector &) const |
Method for partially assembled action. More... | |
virtual void | AssembleEAInteriorFaces (const FiniteElementSpace &fes, Vector &ea_data_int, Vector &ea_data_ext, const bool add) |
virtual void | AssembleEABoundaryFaces (const FiniteElementSpace &fes, Vector &ea_data_bdr, const bool add) |
Public Member Functions inherited from mfem::BilinearFormIntegrator | |
virtual void | AssemblePA (const FiniteElementSpace &fes) |
Method defining partial assembly. More... | |
virtual void | AssemblePA (const FiniteElementSpace &trial_fes, const FiniteElementSpace &test_fes) |
virtual void | AssembleDiagonalPA (Vector &diag) |
Assemble diagonal and add it to Vector diag. More... | |
virtual void | AssembleDiagonalPA_ADAt (const Vector &D, Vector &diag) |
Assemble diagonal of ADA^T (A is this integrator) and add it to diag. More... | |
virtual void | AssembleEA (const FiniteElementSpace &fes, Vector &emat, const bool add=true) |
Method defining element assembly. More... | |
virtual void | AssembleMF (const FiniteElementSpace &fes) |
Method defining matrix-free assembly. More... | |
virtual void | AddMultMF (const Vector &x, Vector &y) const |
virtual void | AddMultTransposeMF (const Vector &x, Vector &y) const |
virtual void | AssembleDiagonalMF (Vector &diag) |
Assemble diagonal and add it to Vector diag. More... | |
virtual void | AssembleElementMatrix (const FiniteElement &el, ElementTransformation &Trans, DenseMatrix &elmat) |
Given a particular Finite Element computes the element matrix elmat. More... | |
virtual void | AssembleElementMatrix2 (const FiniteElement &trial_fe, const FiniteElement &test_fe, ElementTransformation &Trans, DenseMatrix &elmat) |
virtual void | AssembleFaceMatrix (const FiniteElement &trial_face_fe, const FiniteElement &test_fe1, const FiniteElement &test_fe2, FaceElementTransformations &Trans, DenseMatrix &elmat) |
virtual void | AssembleElementVector (const FiniteElement &el, ElementTransformation &Tr, const Vector &elfun, Vector &elvect) |
Perform the local action of the BilinearFormIntegrator. Note that the default implementation in the base class is general but not efficient. More... | |
virtual void | AssembleFaceVector (const FiniteElement &el1, const FiniteElement &el2, FaceElementTransformations &Tr, const Vector &elfun, Vector &elvect) |
Perform the local action of the BilinearFormIntegrator resulting from a face integral term. Note that the default implementation in the base class is general but not efficient. More... | |
virtual void | AssembleElementGrad (const FiniteElement &el, ElementTransformation &Tr, const Vector &elfun, DenseMatrix &elmat) |
Assemble the local gradient matrix. More... | |
virtual void | AssembleFaceGrad (const FiniteElement &el1, const FiniteElement &el2, FaceElementTransformations &Tr, const Vector &elfun, DenseMatrix &elmat) |
Assemble the local action of the gradient of the NonlinearFormIntegrator resulting from a face integral term. More... | |
virtual void | ComputeElementFlux (const FiniteElement &el, ElementTransformation &Trans, Vector &u, const FiniteElement &fluxelem, Vector &flux, bool with_coef=true) |
Virtual method required for Zienkiewicz-Zhu type error estimators. More... | |
virtual double | ComputeFluxEnergy (const FiniteElement &fluxelem, ElementTransformation &Trans, Vector &flux, Vector *d_energy=NULL) |
Virtual method required for Zienkiewicz-Zhu type error estimators. More... | |
virtual | ~BilinearFormIntegrator () |
Public Member Functions inherited from mfem::NonlinearFormIntegrator | |
virtual void | SetIntRule (const IntegrationRule *ir) |
Prescribe a fixed IntegrationRule to use (when ir != NULL) or let the integrator choose (when ir == NULL). More... | |
void | SetIntegrationRule (const IntegrationRule &ir) |
Prescribe a fixed IntegrationRule to use. More... | |
void | SetPAMemoryType (MemoryType mt) |
const IntegrationRule * | GetIntegrationRule () const |
Get the integration rule of the integrator (possibly NULL). More... | |
virtual double | GetElementEnergy (const FiniteElement &el, ElementTransformation &Tr, const Vector &elfun) |
Compute the local energy. More... | |
virtual void | AssembleGradPA (const Vector &x, const FiniteElementSpace &fes) |
Prepare the integrator for partial assembly (PA) gradient evaluations on the given FE space fes at the state x. More... | |
virtual double | GetLocalStateEnergyPA (const Vector &x) const |
Compute the local (to the MPI rank) energy with partial assembly. More... | |
virtual void | AddMultGradPA (const Vector &x, Vector &y) const |
Method for partially assembled gradient action. More... | |
virtual void | AssembleGradDiagonalPA (Vector &diag) const |
Method for computing the diagonal of the gradient with partial assembly. More... | |
virtual bool | SupportsCeed () const |
Indicates whether this integrator can use a Ceed backend. More... | |
ceed::Operator & | GetCeedOp () |
virtual | ~NonlinearFormIntegrator () |
Static Public Member Functions | |
static const IntegrationRule & | GetRule (Geometry::Type geom, int order, FaceElementTransformations &T) |
Protected Attributes | |
Coefficient * | rho |
VectorCoefficient * | u |
double | alpha |
double | beta |
Vector | pa_data |
const DofToQuad * | maps |
Not owned. More... | |
const FaceGeometricFactors * | geom |
Not owned. More... | |
int | dim |
int | nf |
int | nq |
int | dofs1D |
int | quad1D |
Protected Attributes inherited from mfem::NonlinearFormIntegrator | |
const IntegrationRule * | IntRule |
ceed::Operator * | ceedOp |
MemoryType | pa_mt = MemoryType::DEFAULT |
Additional Inherited Members | |
Protected Member Functions inherited from mfem::BilinearFormIntegrator | |
BilinearFormIntegrator (const IntegrationRule *ir=NULL) | |
Protected Member Functions inherited from mfem::NonlinearFormIntegrator | |
NonlinearFormIntegrator (const IntegrationRule *ir=NULL) | |
Integrator for the DG form: alpha < rho_u (u.n) {v},[w] > + beta < rho_u |u.n| [v],[w] >, where v and w are the trial and test variables, respectively, and rho/u are given scalar/vector coefficients. {v} represents the average value of v on the face and [v] is the jump such that {v}=(v1+v2)/2 and [v]=(v1-v2) for the face between elements 1 and 2. For boundary elements, v2=0. The vector coefficient, u, is assumed to be continuous across the faces and when given the scalar coefficient, rho, is assumed to be discontinuous. The integrator uses the upwind value of rho, rho_u, which is value from the side into which the vector coefficient, u, points.
One use case for this integrator is to discretize the operator -u.grad(v) with a DG formulation. The resulting formulation uses the ConvectionIntegrator (with coefficient u, and parameter alpha = -1) and the transpose of the DGTraceIntegrator (with coefficient u, and parameters alpha = 1, beta = -1/2 to use the upwind face flux, see also NonconservativeDGTraceIntegrator). This discretization and the handling of the inflow and outflow boundaries is illustrated in Example 9/9p.
Another use case for this integrator is to discretize the operator -div(u v) with a DG formulation. The resulting formulation is conservative and consists of the ConservativeConvectionIntegrator (with coefficient u, and parameter alpha = -1) plus the DGTraceIntegrator (with coefficient u, and parameters alpha = -1, beta = -1/2 to use the upwind face flux).
Definition at line 2805 of file bilininteg.hpp.
|
inline |
Construct integrator with rho = 1, b = 0.5*a.
Definition at line 2822 of file bilininteg.hpp.
|
inline |
Construct integrator with rho = 1.
Definition at line 2826 of file bilininteg.hpp.
|
inline |
Definition at line 2829 of file bilininteg.hpp.
Method for partially assembled action.
Perform the action of integrator on the input x and add the result to the output y. Both x and y are E-vectors, i.e. they represent the element-wise discontinuous version of the FE space.
This method can be called only after the method AssemblePA() has been called.
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 1133 of file bilininteg_dgtrace_pa.cpp.
Method for partially assembled transposed action.
Perform the transpose action of integrator on the input x and add the result to the output y. Both x and y are E-vectors, i.e. they represent the element-wise discontinuous version of the FE space.
This method can be called only after the method AssemblePA() has been called.
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 1140 of file bilininteg_dgtrace_pa.cpp.
|
virtual |
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 427 of file bilininteg_dgtrace_ea.cpp.
|
virtual |
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 348 of file bilininteg_dgtrace_ea.cpp.
|
virtual |
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 2922 of file bilininteg.cpp.
|
virtual |
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 298 of file bilininteg_dgtrace_pa.cpp.
|
virtual |
Reimplemented from mfem::BilinearFormIntegrator.
Definition at line 293 of file bilininteg_dgtrace_pa.cpp.
|
static |
Definition at line 3056 of file bilininteg.cpp.
|
protected |
Definition at line 2810 of file bilininteg.hpp.
|
protected |
Definition at line 2810 of file bilininteg.hpp.
|
protected |
Definition at line 2815 of file bilininteg.hpp.
|
protected |
Definition at line 2815 of file bilininteg.hpp.
|
protected |
Not owned.
Definition at line 2814 of file bilininteg.hpp.
|
protected |
Not owned.
Definition at line 2813 of file bilininteg.hpp.
|
protected |
Definition at line 2815 of file bilininteg.hpp.
|
protected |
Definition at line 2815 of file bilininteg.hpp.
|
protected |
Definition at line 2812 of file bilininteg.hpp.
|
protected |
Definition at line 2815 of file bilininteg.hpp.
|
protected |
Definition at line 2808 of file bilininteg.hpp.
|
protected |
Definition at line 2809 of file bilininteg.hpp.