Skip to content

Commit

Permalink
Update derivatives.cpp
Browse files Browse the repository at this point in the history
  • Loading branch information
chemiskyy committed Jan 4, 2023
1 parent 8f131ea commit b06f66a
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions src/Continuum_mechanics/Functions/derivatives.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -66,15 +66,15 @@ mat dinvSdSsym(const mat &S) {
mat invS = inv(S);

Tensor2<double,3,3> invS_ = mat_FTensor2(invS);
Tensor4<double,3,3,3,3> dinvSdS_;
Tensor4<double,3,3,3,3> dinvSdSsym_;

Index<'i', 3> i;
Index<'j', 3> j;
Index<'k', 3> k;
Index<'l', 3> l;

dinvSdS_(i,j,k,l) = invS_(i,k)*invS_(j,l)+invS_(i,l)*invS_(j,k);
return 0.5*FTensor4_mat(dinvSdS_);
dinvSdSsym_(i,j,k,l) = invS_(i,k)*invS_(j,l)+invS_(i,l)*invS_(j,k);
return 0.5*FTensor4_mat(dinvSdSsym_);
}

} //namespace simcoon

0 comments on commit b06f66a

Please sign in to comment.