2561 {
2563
2564 int row_nb_dofs = row_data.
getIndices().size();
2565 int col_nb_dofs = col_data.
getIndices().size();
2566 auto &locMat = OP::locMat;
2567 locMat.resize(row_nb_dofs, col_nb_dofs, false);
2568 locMat.clear();
2569
2570
2571 EntityHandle fe_ent = OP::getFEEntityHandle();
2572
2574
2575 if (bc.faces.find(fe_ent) != bc.faces.end()) {
2576
2577 auto t_normal = OP::getFTensor1NormalsAtGaussPts();
2578 auto t_w = OP::getFTensor0IntegrationWeight();
2579
2583
2584 int nb_integration_pts = OP::getGaussPts().size2();
2585 int nb_base_functions = row_data.
getN().size2();
2587
2589
2590 for (int gg = 0; gg != nb_integration_pts; ++gg) {
2591
2593 t_N(
i) = t_normal(
i);
2595
2596 int rr = 0;
2597 for (; rr != row_nb_dofs /
SPACE_DIM; ++rr) {
2598 auto t_mat = getFTensor2FromArray<SPACE_DIM, SPACE_DIM, SPACE_DIM>(
2601 for (
auto cc = 0; cc != col_nb_dofs /
SPACE_DIM; ++cc) {
2602 t_mat(
i,
j) += ((t_w * t_row_base) * (t_N(
k) * t_col_base(
k))) *
2604 ++t_mat;
2605 ++t_col_base;
2606 }
2607 ++t_row_base;
2608 }
2609
2610 for (; rr != nb_base_functions; ++rr)
2611 ++t_row_base;
2612
2613 ++t_w;
2614 ++t_normal;
2615 }
2616
2617 locMat *= OP::getMeasure();
2618 }
2619 }
2621}
#define FTENSOR_INDEX(DIM, I)
Tensor1< T, Tensor_Dim > normalize()
#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()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
auto getFTensor1N(FieldApproximationBase base)
Get base functions for Hdiv/Hcurl spaces.
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.