message: Removed duplicate calculation in gradient
diff: @@ -596,6 +596,7 @@ void IncompressibleNeoHookeanIntegrator::AssembleElementGrad(
  double pres = Sh_p * *elfun[1];
  double mu = c_mu->Eval(Tr, ip);
  double dJ = F.Det();
  + double dJ_FinvT_DS;
  CalcInverseTranspose(F, FinvT);
  @@ -647,10 +648,10 @@ void IncompressibleNeoHookeanIntegrator::AssembleElementGrad(
  {
  for (int l=0; l<dim; l++)
  {
  - (*elmats(1,0))(i_p, j_u + dof_u * dim_u) += dJ * FinvT(dim_u,l) * DS_u(j_u,
  - l) * Sh_p(i_p) * ip.weight * Tr.Weight();
  - (*elmats(0,1))(j_u + dof_u * dim_u, i_p) -= dJ * FinvT(dim_u,l) * DS_u(j_u,
  - l) * Sh_p(i_p) * ip.weight * Tr.Weight();
  + dJ_FinvT_DS = dJ * FinvT(dim_u,l) * DS_u(j_u,l) * Sh_p(i_p) * ip.weight * Tr.Weight();
  + (*elmats(1,0))(i_p, j_u + dof_u * dim_u) += dJ_FinvT_DS;
  + (*elmats(0,1))(j_u + dof_u * dim_u, i_p) -= dJ_FinvT_DS;
  +
  }
  }
  }
  