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

remove some unused code

[[Imported from SVN: r7133]]
parent 8d4f32bc
No related branches found
No related tags found
No related merge requests found
......@@ -364,28 +364,9 @@ evaluateFDDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>
LocalGeodesicFEFunction<dim,double,TargetSpace> fPlus(cornersPlus);
LocalGeodesicFEFunction<dim,double,TargetSpace> fMinus(cornersMinus);
Dune::FieldMatrix<double,TargetSpace::EmbeddedTangentVector::size,dim> h = evaluateDerivative(local);
Dune::FieldMatrix<double,TargetSpace::EmbeddedTangentVector::size,dim> hPlus = fPlus.evaluateDerivative(local);
Dune::FieldMatrix<double,TargetSpace::EmbeddedTangentVector::size,dim> hMinus = fMinus.evaluateDerivative(local);
#if 0
double l = sqrt(h[0][0]*h[0][0] + h[1][0]*h[1][0]);
double lPlus = sqrt(hPlus[0][0]*hPlus[0][0] + hPlus[1][0]*hPlus[1][0]);
double lMinus = sqrt(hMinus[0][0]*hMinus[0][0] + hMinus[1][0]*hMinus[1][0]);
std::cout << "h:\n" << h << "length: " << l << std::endl<< std::endl;
std::cout << "hPlus:\n" << hPlus << "length: " << lPlus << std::endl<< std::endl;
std::cout << "hMinus:\n" << hMinus << "length: " << lMinus << std::endl << std::endl;
#endif
/* h /= l;
hPlus /= lPlus;
hMinus /= lMinus;*/
/* hPlus -= h;
h -= hMinus;*/
/* std::cout << "forward:\n" << hPlus << std::endl;
std::cout << "backward:\n" << h << std::endl;*/
result[j] = hPlus;
result[j] -= hMinus;
result[j] /= 2*eps;
......
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