Skip to content
Snippets Groups Projects
Commit c5788af4 authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

fix another place where simplex grids where hardcoded

[[Imported from SVN: r7848]]
parent 99bfcfb5
No related branches found
No related tags found
No related merge requests found
......@@ -598,21 +598,12 @@ public:
Dune::FieldMatrix<ctype, embeddedDim, dim> result(0);
// get translation part
for (int i=0; i<translationCoefficients_.size(); i++) {
// get derivative of shape function
Dune::FieldVector<ctype,dim> sfDer;
if (i==0)
sfDer = -1;
else
for (int j=0; j<dim; j++)
sfDer[j] = (i-1)==j;
for (int j=0; j<3; j++)
result[j].axpy(translationCoefficients_[i][j], sfDer);
}
std::vector<Dune::FieldMatrix<ctype,1,dim> > sfDer(translationCoefficients_.size());
localFiniteElement_.localBasis().evaluateJacobian(local, sfDer);
for (int i=0; i<translationCoefficients_.size(); i++)
for (int j=0; j<3; j++)
result[j].axpy(translationCoefficients_[i][j], sfDer[i][0]);
// get orientation part
Dune::FieldMatrix<ctype,4,dim> qResult = orientationFEFunction_->evaluateDerivative(local);
......
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