106 const int number_of_active_variables =
physicsPtr->activeVariables.size();
107 const int number_of_dependent_variables =
109 std::vector<double *> jac_ptr(number_of_dependent_variables);
110 for (
unsigned int n = 0;
n != number_of_dependent_variables; ++
n) {
113 ->dependentVariablesPiolaDirevatives[
n *
114 number_of_active_variables]);
117 const auto nb_integration_pts = getGaussPts().size2();
119 auto create_data_mat = [nb_integration_pts](
auto &
m) {
120 m.resize(9, nb_integration_pts,
false);
123 dataAtPts->P_du.resize(81, nb_integration_pts,
false);
125 auto iu = getFTensor2SymmetricFromMat<3>(
dataAtPts->stretchTensorAtPts);
126 auto t_grad_h1 = getFTensor2FromMat<3, 3>(
dataAtPts->wGradH1AtPts);
127 auto r_P_du = getFTensor4FromMat<3, 3, 3, 3>(
dataAtPts->P_du);
135 for (
unsigned int gg = 0; gg != nb_integration_pts; ++gg) {
151 int r = ::jacobian(
tAg, number_of_dependent_variables,
152 number_of_active_variables,
153 &
physicsPtr->activeVariables[0], &jac_ptr[0]);
156 "ADOL-C function evaluation with error");
170 (t_P_dh_tmp(
i,
j, o, p) * t_h1_du(o, p,
m,
n)) * t_h1_du(
i,
j,
k,
l);
173 r_P_du(
i,
j,
m,
n) = t_P_dh_tmp(
i,
j,
m,
n);