headercheck fixes

[[Imported from SVN: r5559]]

src/localgeodesicfestiffness.hh

@@ -8,6 +8,8 @@
 #include <dune/disc/operators/localstiffness.hh>
 #include<dune/disc/operators/boundaryconditions.hh>
 
+#include "rigidbodymotion.hh"
+
 template<class GridView, class TargetSpace>
 class LocalGeodesicFEStiffness 
     : public Dune::LocalStiffness<GridView,double,TargetSpace::TangentVector::size>

src/makestraightrod.hh

@@ -5,7 +5,7 @@
 #include <dune/common/fvector.hh>
 #include <dune/ag-common/crossproduct.hh>
 
-#include "rotation.hh"
+#include "rigidbodymotion.hh"
 
 /** \brief Make a straight rod from two given endpoints
 

src/maxnormtrustregion.hh

 #ifndef MAX_NORM_TRUST_REGION_HH
 #define MAX_NORM_TRUST_REGION_HH
 
+#include <vector>
+
+#include <dune/solvers/boxconstraint.hh>
+
 template <int blocksize>
 class MaxNormTrustRegion
 {

src/planarrodassembler.hh

@@ -7,6 +7,8 @@
 #include <dune/istl/matrixindexset.hh>
 #include <dune/istl/matrix.hh>
 
+#include "rigidbodymotion.hh"
+
 namespace Dune 
 {
 

src/riemanniantrsolver.hh

@@ -10,7 +10,7 @@
 
 #include <dune/solvers/boxconstraint.hh>
 #include <dune/solvers/norms/h1seminorm.hh>
-#include <dune/solvers/solvers/solver.hh>
+#include <dune/solvers/solvers/iterativesolver.hh>
 
 #include "geodesicfeassembler.hh"
 

src/rodrefine.hh

 #ifndef ROD_REFINE_HH
 #define ROD_REFINE_HH
 
+#include <vector>
+#include <map>
+
 #include "rigidbodymotion.hh"
 
 template <class GridType>

src/rodwriter.hh

@@ -2,6 +2,7 @@
 #define ROD_WRITER_HH
 
 #include <fstream>
+#include <vector>
 
 #include <dune/common/exceptions.hh>
 #include <dune/istl/bvector.hh>

src/svd.hh

@@ -4,7 +4,7 @@
 #ifndef SVD_HH
 #define SVD_HH
 
-#include <math.h>
+#include <cmath.h>
 
 
 template <class T>
@@ -19,10 +19,10 @@ template <class T>
 T pythag(T a, T b)
 {
     T absa,absb;
-    absa=std::abs(a);
-    absb=std::abs(b);
+    absa=std::fabs(a);
+    absb=std::fabs(b);
     if (absa > absb) 
-        return absa*sqrt(1.0+SQR(absb/absa));
+        return absa*std::sqrt(1.0+SQR(absb/absa));
     else 
         return (absb == 0.0 ? 0.0 : absb*sqrt(1.0+SQR(absa/absb)));
 }