291 {
293
298
300
301 AssemblyDomainEleOp::locMat.resize(AssemblyDomainEleOp::nbRows, AssemblyDomainEleOp::nbCols,
302 false);
303 AssemblyDomainEleOp::locMat.clear();
304
305 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
306 const auto nb_row_base_functions = row_data.getN().size2();
307
308 if (AssemblyDomainEleOp::colType == MBVERTEX &&
309 AssemblyDomainEleOp::colSide == 0) {
310
311 resDiff.resize(nb_integration_pts, DIM * DIM,
false);
312 auto t_res_diff = getFTensor2FromMat<DIM, DIM>(
resDiff);
313 auto t_grad = getFTensor2FromMat<DIM, DIM>(*(
commonDataPtr->mGradPtr));
314 auto t_logC_dC =
316 auto t_c_dstrain = getFTensor2SymmetricFromMat<DIM>(
commonDataPtr->resCdStrain);
317
318 auto next = [&]() {
319 ++t_grad;
320 ++t_logC_dC;
321 ++t_c_dstrain;
322 ++t_res_diff;
323 };
324
325 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
329
330 t_diff_ls_dlog_c(
k,
l) =
331 (t_c_dstrain(
i,
j)) * (t_logC_dC(
i,
j,
k,
l) / 2);
334
336 t_res_diff(
i,
j) = (t_diff_ls_dlog_c(
k,
l) * t_dC_dF(
k,
l,
i,
j));
337 next();
338 }
339 }
340
341 auto t_res_diff = getFTensor2FromMat<DIM, DIM>(
resDiff);
342
343 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
344 auto t_row_base = row_data.getFTensor0N();
345 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
346 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
347 ++t_w;
348
349 auto t_mat =
350 getFTensor1FromPtr<DIM>(AssemblyDomainEleOp::locMat.data().data());
351 size_t rr = 0;
352 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
353 const auto row_base = alpha * t_row_base;
354 auto t_col_diff_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
355 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / DIM; cc++) {
356 t_mat(
i) += row_base * (t_res_diff(
i,
j) * t_col_diff_base(
j));
357 ++t_mat;
358 ++t_col_diff_base;
359 }
360 ++t_row_base;
361 }
362 for (; rr != nb_row_base_functions; ++rr)
363 ++t_row_base;
364
365 ++t_res_diff;
366 }
367
369}
#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< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
