22 {
24
26
27 const size_t nb_integration_pts = row_data.getN().size1();
28 const size_t nb_row_base_functions = row_data.getN().size2();
29
33 false);
34 auto t_res_diff =
35 getFTensor3FromMat<SPACE_DIM, SPACE_DIM, size_symm>(
resDiff);
36 auto t_D = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
mDPtr);
37 auto t_logC_dC = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM>(
39 auto t_grad = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(
43 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
47 t_DLogC_dC;
48 t_DLogC_dC(
i,
j,
k,
l) = t_D(
m,
n,
k,
l) * t_logC_dC(
m,
n,
i,
j);
50 t_FDLogC_dC;
51 t_FDLogC_dC(
i,
j,
k,
l) = t_F(
i,
m) * t_DLogC_dC(
m,
j,
k,
l);
53 t_res_diff(
i,
j,
L) = t_FDLogC_dC(
i,
j,
k,
l) * t_L(
k,
l,
L);
54 ++t_logC_dC;
55 ++t_grad;
56 ++t_res_diff;
57 }
58 }
59
60 auto t_w = getFTensor0IntegrationWeight();
61 auto t_row_diff_base = row_data.getFTensor1DiffN<
SPACE_DIM>();
62 auto t_res_diff =
63 getFTensor3FromMat<SPACE_DIM, SPACE_DIM, size_symm>(
resDiff);
64
65 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
66 double alpha = getMeasure() * t_w;
67
68 size_t rr = 0;
70
71 auto t_mat =
73
75 t_tmp(
i,
L) = (t_res_diff(
i,
j,
L) * (
alpha * t_row_diff_base(
j)));
76
77 auto t_col_base = col_data.getFTensor0N(gg, 0);
79
80 t_mat(
i,
L) -= (t_col_base * t_tmp(
i,
L));
81
82 ++t_mat;
83 ++t_col_base;
84 }
85
86 ++t_row_diff_base;
87 }
88
89 for (; rr < nb_row_base_functions; ++rr)
90 ++t_row_diff_base;
91
92 ++t_w;
93 ++t_res_diff;
94 }
95
97}
#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< 'j', SPACE_DIM > j
FTensor::Index< 'L', size_symm > L
FTensor::Index< 'l', SPACE_DIM > l
FTensor::Index< 'k', SPACE_DIM > k
FTensor::Index< 'i', SPACE_DIM > i
[Common data]
FTensor::Index< 'm', SPACE_DIM > m
static FTensor::Tensor2< FTensor::PackPtr< double *, 3 >, 2, 3 > get_mat_vector_dtensor_sym(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
FTensor::Index< 'n', SPACE_DIM > n
int rowSide
row side number
int nbRows
number of dofs on rows
MatrixDouble locMat
local entity block matrix
int nbCols
number if dof on column
EntityType rowType
row type