Skip to content
Snippets Groups Projects
Commit 3fa30a78 authored by Youett, Jonathan's avatar Youett, Jonathan Committed by akbib@FU-BERLIN.DE
Browse files

use new method coefficient(int i) to make the basis work for RigidBodyMotions

[[Imported from SVN: r8032]]
parent 16818f45
No related branches found
No related tags found
No related merge requests found
...@@ -160,7 +160,7 @@ void LocalGfeTestFunctionBasis<LocalFiniteElement,TargetSpace>::evaluateFunction ...@@ -160,7 +160,7 @@ void LocalGfeTestFunctionBasis<LocalFiniteElement,TargetSpace>::evaluateFunction
/** \todo This call internally keeps computing the value of the gfe function at 'local'. /** \todo This call internally keeps computing the value of the gfe function at 'local'.
* This is expensive. Eventually we should precompute it once and reused the result. */ * This is expensive. Eventually we should precompute it once and reused the result. */
localGFEFunction_.evaluateDerivativeOfValueWRTCoefficient (local, i, derivative); localGFEFunction_.evaluateDerivativeOfValueWRTCoefficient (local, i, derivative);
Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficients_[i].orthonormalFrame(); Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficient(i).orthonormalFrame();
for (int j=0; j<spaceDim; j++) for (int j=0; j<spaceDim; j++)
derivative.mv(basisVectors[j], out[i][j]); derivative.mv(basisVectors[j], out[i][j]);
...@@ -182,7 +182,7 @@ void LocalGfeTestFunctionBasis<LocalFiniteElement,TargetSpace>::evaluateJacobian ...@@ -182,7 +182,7 @@ void LocalGfeTestFunctionBasis<LocalFiniteElement,TargetSpace>::evaluateJacobian
Tensor3< double, embeddedDim, embeddedDim, dim > derivative; Tensor3< double, embeddedDim, embeddedDim, dim > derivative;
localGFEFunction_.evaluateDerivativeOfGradientWRTCoefficient (local, i, derivative); localGFEFunction_.evaluateDerivativeOfGradientWRTCoefficient (local, i, derivative);
Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficients_[i].orthonormalFrame(); Dune::FieldMatrix<ctype,spaceDim,embeddedDim> basisVectors = localGFEFunction_.coefficient(i).orthonormalFrame();
for (int j=0; j<spaceDim; j++) { for (int j=0; j<spaceDim; 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