372 const int number_of_active_variables =
physicsPtr->activeVariables.size();
373 const int number_of_dependent_variables =
375 std::vector<double *> jac_ptr(number_of_dependent_variables);
376 for (
unsigned int n = 0;
n != number_of_dependent_variables; ++
n) {
379 ->dependentVariablesPiolaDirevatives[
n *
380 number_of_active_variables]);
383 const auto nb_integration_pts = getGaussPts().size2();
385 auto create_data_mat = [nb_integration_pts](
auto &
m) {
386 m.resize(9, nb_integration_pts,
false);
389 dataAtPts->P_du.resize(81, nb_integration_pts,
false);
391 auto iu = getFTensor2SymmetricFromMat<3>(
dataAtPts->stretchTensorAtPts);
392 auto t_grad_h1 = getFTensor2FromMat<3, 3>(
dataAtPts->wGradH1AtPts);
393 auto r_P_du = getFTensor4FromMat<3, 3, 3, 3>(
dataAtPts->P_du);
401 for (
unsigned int gg = 0; gg != nb_integration_pts; ++gg) {
421 int r = ::jacobian(
tAg, number_of_dependent_variables,
422 number_of_active_variables,
423 &
physicsPtr->activeVariables[0], &jac_ptr[0]);
426 "ADOL-C function evaluation with error");
441 (t_P_dh_tmp(
i,
j, o, p) * t_h1_du(o, p,
m,
n)) * t_h1_du(
i,
j,
k,
l);
444 r_P_du(
i,
j,
m,
n) = t_P_dh_tmp(
i,
j,
m,
n);