Skip to content
Snippets Groups Projects
Commit 71250bf8 authored by akbib's avatar akbib Committed by akbib
Browse files

Don't hardwire double and use more methods from SymmetricTensor

[[Imported from SVN: r12661]]
parent a020e36d
No related branches found
No related tags found
No related merge requests found
...@@ -192,10 +192,9 @@ protected: ...@@ -192,10 +192,9 @@ protected:
SymmetricTensor<dim> stress = strain; SymmetricTensor<dim> stress = strain;
stress *= E_/(1+nu_); stress *= E_/(1+nu_);
double f = E_*nu_ / ((1+nu_)*(1-2*nu_)) * strain.trace(); ctype f = E_*nu_ / ((1+nu_)*(1-2*nu_)) * strain.trace();
for (int i=0; i<dim; i++) stress.addToDiag(f);
stress[i] += f;
return stress; return stress;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment