1205 {
1207
1210 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
1212
1213 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
1214 const auto nb_row_base_functions = row_data.getN().size2();
1215
1216 auto t_c_dstrain =
1217 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->resCdStrain);
1218 auto t_c_dplastic_strain =
1219 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->resCdPlasticStrain);
1220
1221 auto next = [&]() {
1222 ++t_c_dstrain;
1223 ++t_c_dplastic_strain;
1224 };
1225
1227
1228 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
1229 auto t_row_base = row_data.getFTensor0N();
1230 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
1231 const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
1232 ++t_w;
1233
1235 t_res_vec(L) =
1236 t_L(
i,
j, L) * (t_c_dplastic_strain(
i,
j) - t_c_dstrain(
i,
j));
1237 next();
1238
1241 size_t rr = 0;
1242 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
1243 const auto row_base = alpha * t_row_base;
1244 auto t_col_base = col_data.getFTensor0N(gg, 0);
1245 for (
size_t cc = 0; cc != AssemblyDomainEleOp::nbCols /
size_symm; cc++) {
1246 t_mat(L) += (row_base * t_col_base) * t_res_vec(L);
1247 ++t_mat;
1248 ++t_col_base;
1249 }
1250 ++t_row_base;
1251 }
1252 for (; rr != nb_row_base_functions; ++rr)
1253 ++t_row_base;
1254 }
1255
1257}
#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
auto get_mat_scalar_dtensor_sym(MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)