Skip to content
Snippets Groups Projects
Commit dc750bcd authored by Norihiro Watanabe's avatar Norihiro Watanabe
Browse files

[Mat/Frac] suggested optimization

parent 8d03f8f5
No related branches found
No related tags found
No related merge requests found
......@@ -102,10 +102,10 @@ void MohrCoulomb<DisplacementDim>::computeConstitutiveRelation(
else
{
for (int i=0; i<index_ns; i++)
dFs_dS[i] = boost::math::sign(tau) / tau * sigma[i];
dFs_dS[i] = sigma[i] / std::abs(tau);
for (int i=0; i<index_ns; i++)
dQs_dS[i] = boost::math::sign(tau) / tau * sigma[i];
dQs_dS[i] = sigma[i] / std::abs(tau);
}
dFs_dS[index_ns] = std::tan(mat.phi);
dQs_dS[index_ns] = std::tan(mat.psi);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment