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
Public Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >