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

Reenable tests for Euclidean Hesse matrix

[[Imported from SVN: r9913]]
parent 0787ae4a
No related branches found
No related tags found
No related merge requests found
......@@ -528,7 +528,7 @@ int main (int argc, char *argv[]) try
Matrix<FieldMatrix<double,embeddedBlocksize,embeddedBlocksize> > localADHessian;
Matrix<FieldMatrix<double,embeddedBlocksize,embeddedBlocksize> > localADVMHessian; // VM: vector-mode
Matrix<FieldMatrix<double,embeddedBlocksize,embeddedBlocksize> > localFDHessian;
if (false) {
// Assemble Euclidean derivatives
localADOLCStiffness.assembleGradientAndHessian(*it,
feBasis.getLocalFiniteElement(*it),
......@@ -553,7 +553,7 @@ int main (int argc, char *argv[]) try
// compare
compareMatrices(localADHessian, "AD", localFDHessian, "FD");
compareMatrices(localADHessian, "AD scalar", localADVMHessian, "AD vector");
}
// Assemble Riemannian derivatives
std::vector<Dune::FieldVector<double,blocksize> > localRiemannianADGradient(numOfBaseFct);
std::vector<Dune::FieldVector<double,blocksize> > localRiemannianFDGradient(numOfBaseFct);
......
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