Skip to content
Snippets Groups Projects
Commit ac02c1be authored by Jonathan Youett's avatar Jonathan Youett
Browse files

Initialize using the correct type

parent 045d0095
No related branches found
No related tags found
No related merge requests found
......@@ -18,7 +18,7 @@ namespace Dune {
*/
template <int dim, class field_type=double>
void computeNonlinearStrain(const Dune::FieldMatrix<field_type, dim, dim>& grad, SymmetricTensor<dim,field_type>& strain) {
strain = 0;
strain = field_type(0);
for (int i=0; i<dim ; ++i) {
strain(i,i) +=grad[i][i];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment