Skip to content
Snippets Groups Projects
Commit d485546b authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

use ctype instead of hardcoded 'double'

[[Imported from SVN: r9299]]
parent b82398ae
No related branches found
No related tags found
No related merge requests found
...@@ -232,7 +232,7 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local, const TargetSpace ...@@ -232,7 +232,7 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local, const TargetSpace
// the matrix that turns coordinates on the reference simplex into coordinates on the standard simplex // the matrix that turns coordinates on the reference simplex into coordinates on the standard simplex
std::vector<Dune::FieldMatrix<ctype,1,dim> > BNested(coefficients_.size()); std::vector<Dune::FieldMatrix<ctype,1,dim> > BNested(coefficients_.size());
localFiniteElement_.localBasis().evaluateJacobian(local, BNested); localFiniteElement_.localBasis().evaluateJacobian(local, BNested);
Dune::Matrix<Dune::FieldMatrix<double,1,1> > B(coefficients_.size(), dim); Dune::Matrix<Dune::FieldMatrix<ctype,1,1> > B(coefficients_.size(), dim);
for (size_t i=0; i<coefficients_.size(); i++) for (size_t i=0; i<coefficients_.size(); i++)
for (size_t j=0; j<dim; j++) for (size_t j=0; j<dim; j++)
B[i][j] = BNested[i][0][j]; B[i][j] = BNested[i][0][j];
......
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