124 const size_t nb_integration_pts = row_data.getN().size1();
125 const size_t nb_row_base_functions = row_data.getN().size2();
126 auto t_w = getFTensor0IntegrationWeight();
128 auto t_row_base = row_data.getFTensor0N();
139 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->plasticFlow);
141 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
commonDataPtr->mDPtr);
143 auto t_D_Op = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
mDPtr);
152 for (
size_t gg = 0; gg != nb_integration_pts; ++gg) {
153 double alpha = getMeasure() * t_w;
160 t_diff_constrain_dgrad(
k,
l) =
161 t_diff_constrain_dstrain(
i,
j) * t_diff_grad_symmetrise(
i,
j,
k,
l);
167 auto t_col_diff_base = col_data.getFTensor1DiffN<
SPACE_DIM>(gg, 0);
170 t_mat(
i) +=
alpha * t_row_base * (t_diff_constrain_dgrad(
i,
j)) *
179 for (; rr != nb_row_base_functions; ++rr)
#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()
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
auto diff_constrain_dstress(double &&diff_constrain_df, FTensor::Tensor2_symmetric< T, SPACE_DIM > &t_plastic_flow)
FTensor::Index< 'j', SPACE_DIM > j
[Common data]
auto diff_constrain_dstrain(T1 &t_D, T2 &&t_diff_constrain_dstress)
static auto get_mat_scalar_dvector(MatrixDouble &mat, FTensor::Number< 2 >)
FTensor::Index< 'l', SPACE_DIM > l
FTensor::Index< 'k', SPACE_DIM > k
FTensor::Index< 'i', SPACE_DIM > i
auto diff_constrain_df(long double dot_tau, long double f, long double sigma_y)
int nbRows
number of dofs on rows
MatrixDouble locMat
local entity block matrix
int nbCols
number if dof on column