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

use the AverageDistanceAssembler to assemble the Hessian of the average distance functional 

[[Imported from SVN: r5596]]
parent a32673ec
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/localgeodesicfefunction.hh
+ 3
− 3
View file @ a9889474
...@@ -130,10 +130,10 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local) ...@@ -130,10 +130,10 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local)
// the actual system matrix // the actual system matrix
Dune::FieldVector<ctype,dim+1> w = barycentricCoordinates(local); Dune::FieldVector<ctype,dim+1> w = barycentricCoordinates(local);
AverageDistanceAssembler<TargetSpace> assembler(coefficients_, w);
Dune::FieldMatrix<ctype,dim+1,dim+1> dFdq(0); Dune::FieldMatrix<ctype,dim+1,dim+1> dFdq(0);
for (int i=0; i<dim+1; i++) assembler.assembleHessian(q,dFdq);
dFdq.axpy(w[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], q));
// //////////////////////////////////// // ////////////////////////////////////
// solve the system // solve the system
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment