35 const int nb_row_dofs = row_data.
getIndices().size();
36 const int nb_col_dofs = col_data.
getIndices().size();
38 if (nb_row_dofs && nb_col_dofs) {
39 locLhs.resize(nb_row_dofs, nb_col_dofs,
false);
41 const int nb_integration_pts = getGaussPts().size2();
42 auto t_w = getFTensor0IntegrationWeight();
46 const double vol = getMeasure();
47 for (
int gg = 0; gg != nb_integration_pts; ++gg) {
48 const double a = t_w * vol;
49 for (
int rr = 0; rr != nb_row_dofs; ++rr) {
51 for (
int cc = 0; cc != nb_col_dofs; ++cc) {
52 locLhs(rr, cc) +=
a * t_row_diff_base(
i) * t_col_diff_base(
i);
60 &
locLhs(0, 0), ADD_VALUES);
61 if (row_side != col_side || row_type != col_type) {
62 transLocLhs.resize(nb_col_dofs, nb_row_dofs,
false);
#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()
#define CHKERR
Inline error check.
constexpr auto MatSetValues
FTensor::Index< 'i', 2 > i
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > getFTensor1DiffN(const FieldApproximationBase base)
Get derivatives of base functions.
const VectorInt & getIndices() const
Get global indices of dofs on entity.