#include "tutorials/cor-2to5/src/PoissonOperators.hpp"
Definition at line 729 of file PoissonOperators.hpp.
◆ OpRes_g()
PoissonExample::OpRes_g::OpRes_g |
( |
FVal | f_value, |
|
|
boost::shared_ptr< VectorDouble > & | field_vals ) |
|
inline |
◆ iNtegrate()
Integrate local constrains vector.
- Examples
- PoissonOperators.hpp.
Definition at line 740 of file PoissonOperators.hpp.
740 {
742
744
746
747 const double area = getArea() *
bEta;
748
749 auto t_w = getFTensor0IntegrationWeight();
750
752
754
755 auto t_coords = getFTensor1CoordsAtGaussPts();
756
758
759
760 double alpha = area * t_w;
761
764 for (
int rr = 0; rr !=
nbRows; rr++) {
765 t_a += alpha * t_l *
766 (t_u -
fValue(t_coords(
NX), t_coords(
NY), t_coords(
NZ)));
767 ++t_a;
768 ++t_l;
769 }
770 ++t_w;
771 ++t_u;
772 ++t_coords;
773 }
775 }
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
FTensor::Number< 2 > NZ
z-direction index
FTensor::Number< 0 > NX
x-direction index
FTensor::Number< 1 > NY
y-direction index
FVal fValue
Function pointer evaluating values of "U" at the boundary.
◆ fieldVals
boost::shared_ptr<VectorDouble> PoissonExample::OpRes_g::fieldVals |
|
protected |
The documentation for this struct was generated from the following file: