MFEM
v3.0
|
Class for grid function - Vector with associated FE space. More...
#include <gridfunc.hpp>
Public Member Functions | |
GridFunction () | |
GridFunction (FiniteElementSpace *f) | |
Creates grid function associated with *f. More... | |
GridFunction (Mesh *m, std::istream &input) | |
GridFunction (Mesh *m, GridFunction *gf_array[], int num_pieces) | |
void | MakeOwner (FiniteElementCollection *_fec) |
Make the GridFunction the owner of 'fec' and 'fes'. More... | |
FiniteElementCollection * | OwnFEC () |
int | VectorDim () const |
void | GetNodalValues (int i, Array< double > &nval, int vdim=1) const |
Returns the values in the vertices of i'th element for dimension vdim. More... | |
virtual double | GetValue (int i, const IntegrationPoint &ip, int vdim=1) const |
void | GetVectorValue (int i, const IntegrationPoint &ip, Vector &val) const |
void | GetValues (int i, const IntegrationRule &ir, Vector &vals, int vdim=1) const |
void | GetValues (int i, const IntegrationRule &ir, Vector &vals, DenseMatrix &tr, int vdim=1) const |
int | GetFaceValues (int i, int side, const IntegrationRule &ir, Vector &vals, DenseMatrix &tr, int vdim=1) const |
void | GetVectorValues (ElementTransformation &T, const IntegrationRule &ir, DenseMatrix &vals) const |
void | GetVectorValues (int i, const IntegrationRule &ir, DenseMatrix &vals, DenseMatrix &tr) const |
int | GetFaceVectorValues (int i, int side, const IntegrationRule &ir, DenseMatrix &vals, DenseMatrix &tr) const |
void | GetValuesFrom (GridFunction &) |
void | GetBdrValuesFrom (GridFunction &) |
void | GetVectorFieldValues (int i, const IntegrationRule &ir, DenseMatrix &vals, DenseMatrix &tr, int comp=0) const |
void | ReorderByNodes () |
For a vector grid function, makes sure that the ordering is byNODES. More... | |
void | GetNodalValues (Vector &nval, int vdim=1) const |
Return the values as a vector on mesh vertices for dimension vdim. More... | |
void | GetVectorFieldNodalValues (Vector &val, int comp) const |
void | ProjectVectorFieldOn (GridFunction &vec_field, int comp=0) |
void | GetDerivative (int comp, int der_comp, GridFunction &der) |
double | GetDivergence (ElementTransformation &tr) |
void | GetGradient (ElementTransformation &tr, Vector &grad) |
void | GetGradients (const int elem, const IntegrationRule &ir, DenseMatrix &grad) |
void | GetVectorGradient (ElementTransformation &tr, DenseMatrix &grad) |
void | GetElementAverages (GridFunction &avgs) |
void | ImposeBounds (int i, const Vector &weights, const Vector &_lo, const Vector &_hi) |
void | ImposeBounds (int i, const Vector &weights, double _min=0.0, double _max=std::numeric_limits< double >::infinity()) |
void | ProjectGridFunction (const GridFunction &src) |
void | ProjectCoefficient (Coefficient &coeff) |
void | ProjectCoefficient (Coefficient &coeff, Array< int > &dofs, int vd=0) |
void | ProjectCoefficient (VectorCoefficient &vcoeff) |
void | ProjectCoefficient (Coefficient *coeff[]) |
void | ProjectBdrCoefficient (Coefficient &coeff, Array< int > &attr) |
void | ProjectBdrCoefficient (Coefficient *coeff[], Array< int > &attr) |
void | ProjectBdrCoefficientNormal (VectorCoefficient &vcoeff, Array< int > &bdr_attr) |
void | ProjectBdrCoefficientTangent (VectorCoefficient &vcoeff, Array< int > &bdr_attr) |
double | ComputeL2Error (Coefficient &exsol, const IntegrationRule *irs[]=NULL) const |
double | ComputeL2Error (Coefficient *exsol[], const IntegrationRule *irs[]=NULL) const |
double | ComputeL2Error (VectorCoefficient &exsol, const IntegrationRule *irs[]=NULL, Array< int > *elems=NULL) const |
double | ComputeH1Error (Coefficient *exsol, VectorCoefficient *exgrad, Coefficient *ell_coef, double Nu, int norm_type) const |
double | ComputeMaxError (Coefficient &exsol, const IntegrationRule *irs[]=NULL) const |
double | ComputeMaxError (Coefficient *exsol[], const IntegrationRule *irs[]=NULL) const |
double | ComputeMaxError (VectorCoefficient &exsol, const IntegrationRule *irs[]=NULL) const |
double | ComputeL1Error (Coefficient &exsol, const IntegrationRule *irs[]=NULL) const |
double | ComputeW11Error (Coefficient *exsol, VectorCoefficient *exgrad, int norm_type, Array< int > *elems=NULL, const IntegrationRule *irs[]=NULL) const |
double | ComputeL1Error (VectorCoefficient &exsol, const IntegrationRule *irs[]=NULL) const |
double | ComputeLpError (const double p, Coefficient &exsol, Coefficient *weight=NULL, const IntegrationRule *irs[]=NULL) const |
double | ComputeLpError (const double p, VectorCoefficient &exsol, Coefficient *weight=NULL, VectorCoefficient *v_weight=NULL, const IntegrationRule *irs[]=NULL) const |
GridFunction & | operator= (double value) |
Redefine '=' for GridFunction = constant. More... | |
GridFunction & | operator= (const Vector &v) |
GridFunction & | operator= (const GridFunction &v) |
void | ConformingProlongate (const Vector &x) |
void | ConformingProlongate () |
void | ConformingProject (Vector &x) const |
void | ConformingProject () |
FiniteElementSpace * | FESpace () |
void | Update () |
void | Update (FiniteElementSpace *f) |
void | Update (FiniteElementSpace *f, Vector &v, int v_offset) |
virtual void | Save (std::ostream &out) const |
Save the GridFunction to an output stream. More... | |
void | SaveVTK (std::ostream &out, const std::string &field_name, int ref) |
void | SaveSTL (std::ostream &out, int TimesToRefine=1) |
virtual | ~GridFunction () |
Destroys grid function. More... | |
Public Member Functions inherited from mfem::Vector | |
Vector () | |
Default constructor for Vector. Sets size = 0 and data = NULL. More... | |
Vector (const Vector &) | |
Copy constructor. More... | |
Vector (int s) | |
Creates vector of size s. More... | |
Vector (double *_data, int _size) | |
void | Load (std::istream **in, int np, int *dim) |
Reads a vector from multiple files. More... | |
void | Load (std::istream &in, int Size) |
Load a vector from an input stream. More... | |
void | Load (std::istream &in) |
Load a vector from an input stream. More... | |
void | SetSize (int s) |
Resizes the vector if the new size is different. More... | |
void | SetData (double *d) |
void | SetDataAndSize (double *d, int s) |
void | NewDataAndSize (double *d, int s) |
void | MakeDataOwner () |
void | Destroy () |
Destroy a vector. More... | |
int | Size () const |
Returns the size of the vector. More... | |
double * | GetData () const |
operator double * () | |
operator const double * () const | |
bool | OwnsData () const |
void | StealData (double **p) |
Changes the ownership of the data; after the call the Vector is empty. More... | |
double * | StealData () |
Changes the ownership of the data; after the call the Vector is empty. More... | |
double & | Elem (int i) |
Sets value in vector. Index i = 0 .. size-1. More... | |
const double & | Elem (int i) const |
Sets value in vector. Index i = 0 .. size-1. More... | |
double & | operator() (int i) |
Sets value in vector. Index i = 0 .. size-1. More... | |
const double & | operator() (int i) const |
Sets value in vector. Index i = 0 .. size-1. More... | |
double | operator* (const double *) const |
double | operator* (const Vector &v) const |
Return the inner-product. More... | |
Vector & | operator= (const double *v) |
Vector & | operator= (const Vector &v) |
Redefine '=' for vector = vector. More... | |
Vector & | operator= (double value) |
Redefine '=' for vector = constant. More... | |
Vector & | operator*= (double c) |
Vector & | operator/= (double c) |
Vector & | operator-= (double c) |
Vector & | operator-= (const Vector &v) |
Vector & | operator+= (const Vector &v) |
Vector & | Add (const double a, const Vector &Va) |
(*this) += a * Va More... | |
Vector & | Set (const double a, const Vector &x) |
(*this) = a * x More... | |
void | SetVector (const Vector &v, int offset) |
void | Neg () |
(*this) = -(*this) More... | |
void | median (const Vector &lo, const Vector &hi) |
v = median(v,lo,hi) entrywise. Implementation assumes lo <= hi. More... | |
void | GetSubVector (const Array< int > &dofs, Vector &elemvect) const |
void | GetSubVector (const Array< int > &dofs, double *elem_data) const |
void | SetSubVector (const Array< int > &dofs, const Vector &elemvect) |
void | SetSubVector (const Array< int > &dofs, double *elem_data) |
void | AddElementVector (const Array< int > &dofs, const Vector &elemvect) |
Add (element) subvector to the vector. More... | |
void | AddElementVector (const Array< int > &dofs, double *elem_data) |
void | AddElementVector (const Array< int > &dofs, const double a, const Vector &elemvect) |
void | Print (std::ostream &out=std::cout, int width=8) const |
Prints vector to stream out. More... | |
void | Print_HYPRE (std::ostream &out) const |
Prints vector to stream out in HYPRE_Vector format. More... | |
void | Randomize (int seed=0) |
Set random values in the vector. More... | |
double | Norml2 () const |
Returns the l2 norm of the vector. More... | |
double | Normlinf () const |
Returns the l_infinity norm of the vector. More... | |
double | Norml1 () const |
Returns the l_1 norm of the vector. More... | |
double | Normlp (double p) const |
Returns the l_p norm of the vector. More... | |
double | Max () const |
Returns the maximal element of the vector. More... | |
double | Min () const |
Returns the minimal element of the vector. More... | |
double | Sum () const |
Return the sum of the vector entries. More... | |
double | DistanceTo (const double *p) const |
Compute the Euclidian distance to another vector. More... | |
int | CheckFinite () const |
~Vector () | |
Destroys vector. More... | |
Protected Member Functions | |
void | SaveSTLTri (std::ostream &out, double p1[], double p2[], double p3[]) |
void | GetVectorGradientHat (ElementTransformation &T, DenseMatrix &gh) |
void | ProjectDeltaCoefficient (DeltaCoefficient &delta_coeff, double &integral) |
Protected Attributes | |
FiniteElementSpace * | fes |
FE space on which grid function lives. More... | |
FiniteElementCollection * | fec |
Used when the grid function is read from a file. More... | |
Protected Attributes inherited from mfem::Vector | |
int | size |
int | allocsize |
double * | data |
Class for grid function - Vector with associated FE space.
Definition at line 26 of file gridfunc.hpp.
|
inline |
Definition at line 46 of file gridfunc.hpp.
|
inline |
Creates grid function associated with *f.
Definition at line 49 of file gridfunc.hpp.
mfem::GridFunction::GridFunction | ( | Mesh * | m, |
std::istream & | input | ||
) |
Definition at line 26 of file gridfunc.cpp.
mfem::GridFunction::GridFunction | ( | Mesh * | m, |
GridFunction * | gf_array[], | ||
int | num_pieces | ||
) |
Definition at line 52 of file gridfunc.cpp.
|
virtual |
Destroys grid function.
Definition at line 127 of file gridfunc.cpp.
double mfem::GridFunction::ComputeH1Error | ( | Coefficient * | exsol, |
VectorCoefficient * | exgrad, | ||
Coefficient * | ell_coef, | ||
double | Nu, | ||
int | norm_type | ||
) | const |
Definition at line 1324 of file gridfunc.cpp.
|
inline |
Definition at line 198 of file gridfunc.hpp.
|
inline |
Definition at line 206 of file gridfunc.hpp.
|
inline |
Definition at line 166 of file gridfunc.hpp.
double mfem::GridFunction::ComputeL2Error | ( | Coefficient * | exsol[], |
const IntegrationRule * | irs[] = NULL |
||
) | const |
Definition at line 1238 of file gridfunc.cpp.
double mfem::GridFunction::ComputeL2Error | ( | VectorCoefficient & | exsol, |
const IntegrationRule * | irs[] = NULL , |
||
Array< int > * | elems = NULL |
||
) | const |
Definition at line 1285 of file gridfunc.cpp.
double mfem::GridFunction::ComputeLpError | ( | const double | p, |
Coefficient & | exsol, | ||
Coefficient * | weight = NULL , |
||
const IntegrationRule * | irs[] = NULL |
||
) | const |
Definition at line 1593 of file gridfunc.cpp.
double mfem::GridFunction::ComputeLpError | ( | const double | p, |
VectorCoefficient & | exsol, | ||
Coefficient * | weight = NULL , |
||
VectorCoefficient * | v_weight = NULL , |
||
const IntegrationRule * | irs[] = NULL |
||
) | const |
When given a vector weight, compute the pointwise (scalar) error as the dot product of the vector error with the vector weight. Otherwise, the scalar error is the l_2 norm of the vector error.
Definition at line 1650 of file gridfunc.cpp.
|
inline |
Definition at line 181 of file gridfunc.hpp.
double mfem::GridFunction::ComputeMaxError | ( | Coefficient * | exsol[], |
const IntegrationRule * | irs[] = NULL |
||
) | const |
Definition at line 1456 of file gridfunc.cpp.
|
inline |
Definition at line 191 of file gridfunc.hpp.
double mfem::GridFunction::ComputeW11Error | ( | Coefficient * | exsol, |
VectorCoefficient * | exgrad, | ||
int | norm_type, | ||
Array< int > * | elems = NULL , |
||
const IntegrationRule * | irs[] = NULL |
||
) | const |
Definition at line 1503 of file gridfunc.cpp.
void mfem::GridFunction::ConformingProject | ( | Vector & | x | ) | const |
For partially conforming FE spaces, project this partially conforming GridFunction onto the conforming vector x.
Definition at line 1773 of file gridfunc.cpp.
void mfem::GridFunction::ConformingProject | ( | ) |
Same as above, but the destination conforming vector is 'this'. NOTE: the GridFunction's functionality is limited after this call, as the underlying vector shrinks to the number of conforming DOFs. Normal state is restored with ConformingProlongate.
Definition at line 1787 of file gridfunc.cpp.
void mfem::GridFunction::ConformingProlongate | ( | const Vector & | x | ) |
For partially conforming FE spaces, prolongate the conforming vector x to this partially conforming GridFunction.
Definition at line 1750 of file gridfunc.cpp.
void mfem::GridFunction::ConformingProlongate | ( | ) |
As above, but the destination is 'this'.
Definition at line 1764 of file gridfunc.cpp.
|
inline |
Definition at line 246 of file gridfunc.hpp.
void mfem::GridFunction::GetBdrValuesFrom | ( | GridFunction & | orig_func | ) |
Definition at line 454 of file gridfunc.cpp.
void mfem::GridFunction::GetDerivative | ( | int | comp, |
int | der_comp, | ||
GridFunction & | der | ||
) |
Definition at line 612 of file gridfunc.cpp.
double mfem::GridFunction::GetDivergence | ( | ElementTransformation & | tr | ) |
Definition at line 685 of file gridfunc.cpp.
void mfem::GridFunction::GetElementAverages | ( | GridFunction & | avgs | ) |
Compute , where are the basis functions for the FE space of avgs. Both FE spaces should be scalar and on the same mesh.
Definition at line 769 of file gridfunc.cpp.
int mfem::GridFunction::GetFaceValues | ( | int | i, |
int | side, | ||
const IntegrationRule & | ir, | ||
Vector & | vals, | ||
DenseMatrix & | tr, | ||
int | vdim = 1 |
||
) | const |
Definition at line 286 of file gridfunc.cpp.
int mfem::GridFunction::GetFaceVectorValues | ( | int | i, |
int | side, | ||
const IntegrationRule & | ir, | ||
DenseMatrix & | vals, | ||
DenseMatrix & | tr | ||
) | const |
Definition at line 380 of file gridfunc.cpp.
void mfem::GridFunction::GetGradient | ( | ElementTransformation & | tr, |
Vector & | grad | ||
) |
Definition at line 715 of file gridfunc.cpp.
void mfem::GridFunction::GetGradients | ( | const int | elem, |
const IntegrationRule & | ir, | ||
DenseMatrix & | grad | ||
) |
Definition at line 733 of file gridfunc.cpp.
void mfem::GridFunction::GetNodalValues | ( | int | i, |
Array< double > & | nval, | ||
int | vdim = 1 |
||
) | const |
Returns the values in the vertices of i'th element for dimension vdim.
Definition at line 172 of file gridfunc.cpp.
void mfem::GridFunction::GetNodalValues | ( | Vector & | nval, |
int | vdim = 1 |
||
) | const |
Return the values as a vector on mesh vertices for dimension vdim.
Definition at line 882 of file gridfunc.cpp.
|
virtual |
Reimplemented in mfem::ParGridFunction.
Definition at line 211 of file gridfunc.cpp.
void mfem::GridFunction::GetValues | ( | int | i, |
const IntegrationRule & | ir, | ||
Vector & | vals, | ||
int | vdim = 1 |
||
) | const |
Definition at line 256 of file gridfunc.cpp.
void mfem::GridFunction::GetValues | ( | int | i, |
const IntegrationRule & | ir, | ||
Vector & | vals, | ||
DenseMatrix & | tr, | ||
int | vdim = 1 |
||
) | const |
Definition at line 275 of file gridfunc.cpp.
void mfem::GridFunction::GetValuesFrom | ( | GridFunction & | orig_func | ) |
Definition at line 417 of file gridfunc.cpp.
void mfem::GridFunction::GetVectorFieldNodalValues | ( | Vector & | val, |
int | comp | ||
) | const |
Definition at line 550 of file gridfunc.cpp.
void mfem::GridFunction::GetVectorFieldValues | ( | int | i, |
const IntegrationRule & | ir, | ||
DenseMatrix & | vals, | ||
DenseMatrix & | tr, | ||
int | comp = 0 |
||
) | const |
Definition at line 491 of file gridfunc.cpp.
void mfem::GridFunction::GetVectorGradient | ( | ElementTransformation & | tr, |
DenseMatrix & | grad | ||
) |
Definition at line 757 of file gridfunc.cpp.
|
protected |
Definition at line 666 of file gridfunc.cpp.
void mfem::GridFunction::GetVectorValue | ( | int | i, |
const IntegrationPoint & | ip, | ||
Vector & | val | ||
) | const |
Definition at line 224 of file gridfunc.cpp.
void mfem::GridFunction::GetVectorValues | ( | ElementTransformation & | T, |
const IntegrationRule & | ir, | ||
DenseMatrix & | vals | ||
) | const |
Definition at line 328 of file gridfunc.cpp.
void mfem::GridFunction::GetVectorValues | ( | int | i, |
const IntegrationRule & | ir, | ||
DenseMatrix & | vals, | ||
DenseMatrix & | tr | ||
) | const |
Definition at line 371 of file gridfunc.cpp.
void mfem::GridFunction::ImposeBounds | ( | int | i, |
const Vector & | weights, | ||
const Vector & | _lo, | ||
const Vector & | _hi | ||
) |
Impose the given bounds on the function's DOFs while preserving its local integral (described in terms of the given weights) on the i'th element through SLBPQ optimization. Intended to be used for discontinuos FE functions.
Definition at line 826 of file gridfunc.cpp.
void mfem::GridFunction::ImposeBounds | ( | int | i, |
const Vector & | weights, | ||
double | _min = 0.0 , |
||
double | _max = std::numeric_limits<double>::infinity() |
||
) |
Definition at line 853 of file gridfunc.cpp.
|
inline |
Make the GridFunction the owner of 'fec' and 'fes'.
Definition at line 57 of file gridfunc.hpp.
GridFunction & mfem::GridFunction::operator= | ( | double | value | ) |
Redefine '=' for GridFunction = constant.
Definition at line 1731 of file gridfunc.cpp.
GridFunction & mfem::GridFunction::operator= | ( | const Vector & | v | ) |
Definition at line 1738 of file gridfunc.cpp.
GridFunction & mfem::GridFunction::operator= | ( | const GridFunction & | v | ) |
Definition at line 1745 of file gridfunc.cpp.
|
inline |
Definition at line 60 of file gridfunc.hpp.
|
inline |
Definition at line 146 of file gridfunc.hpp.
void mfem::GridFunction::ProjectBdrCoefficient | ( | Coefficient * | coeff[], |
Array< int > & | attr | ||
) |
Definition at line 1064 of file gridfunc.cpp.
void mfem::GridFunction::ProjectBdrCoefficientNormal | ( | VectorCoefficient & | vcoeff, |
Array< int > & | bdr_attr | ||
) |
Project the normal component of the given VectorCoefficient on the boundary. Only boundary attributes that are marked in 'bdr_attr' are projected. Assumes RT-type VectorFE GridFunction.
Definition at line 1150 of file gridfunc.cpp.
void mfem::GridFunction::ProjectBdrCoefficientTangent | ( | VectorCoefficient & | vcoeff, |
Array< int > & | bdr_attr | ||
) |
Project the tangential components of the given VectorCoefficient on the boundary. Only boundary attributes that are marked in 'bdr_attr' are projected. Assumes ND-type VectorFE GridFunction.
Definition at line 1217 of file gridfunc.cpp.
void mfem::GridFunction::ProjectCoefficient | ( | Coefficient & | coeff | ) |
Definition at line 967 of file gridfunc.cpp.
void mfem::GridFunction::ProjectCoefficient | ( | Coefficient & | coeff, |
Array< int > & | dofs, | ||
int | vd = 0 |
||
) |
Definition at line 994 of file gridfunc.cpp.
void mfem::GridFunction::ProjectCoefficient | ( | VectorCoefficient & | vcoeff | ) |
Definition at line 1018 of file gridfunc.cpp.
void mfem::GridFunction::ProjectCoefficient | ( | Coefficient * | coeff[] | ) |
Definition at line 1033 of file gridfunc.cpp.
|
protected |
Definition at line 906 of file gridfunc.cpp.
void mfem::GridFunction::ProjectGridFunction | ( | const GridFunction & | src | ) |
Project the given 'src' GridFunction to 'this' GridFunction, both of which must be on the same mesh. The current implementation assumes that all element use the same projection matrix.
Definition at line 797 of file gridfunc.cpp.
void mfem::GridFunction::ProjectVectorFieldOn | ( | GridFunction & | vec_field, |
int | comp = 0 |
||
) |
Definition at line 579 of file gridfunc.cpp.
void mfem::GridFunction::ReorderByNodes | ( | ) |
For a vector grid function, makes sure that the ordering is byNODES.
Definition at line 529 of file gridfunc.cpp.
|
virtual |
Save the GridFunction to an output stream.
Reimplemented in mfem::ParGridFunction.
Definition at line 1797 of file gridfunc.cpp.
void mfem::GridFunction::SaveSTL | ( | std::ostream & | out, |
int | TimesToRefine = 1 |
||
) |
Definition at line 1897 of file gridfunc.cpp.
|
protected |
Definition at line 1878 of file gridfunc.cpp.
void mfem::GridFunction::SaveVTK | ( | std::ostream & | out, |
const std::string & | field_name, | ||
int | ref | ||
) |
Write the GridFunction in VTK format. Note that Mesh::PrintVTK must be called first. The parameter ref must match the one used in Mesh::PrintVTK.
Definition at line 1807 of file gridfunc.cpp.
|
inline |
Definition at line 248 of file gridfunc.hpp.
void mfem::GridFunction::Update | ( | FiniteElementSpace * | f | ) |
Definition at line 136 of file gridfunc.cpp.
void mfem::GridFunction::Update | ( | FiniteElementSpace * | f, |
Vector & | v, | ||
int | v_offset | ||
) |
Definition at line 148 of file gridfunc.cpp.
int mfem::GridFunction::VectorDim | ( | ) | const |
Definition at line 160 of file gridfunc.cpp.
|
protected |
Used when the grid function is read from a file.
Definition at line 33 of file gridfunc.hpp.
|
protected |
FE space on which grid function lives.
Definition at line 30 of file gridfunc.hpp.