diff --git a/test/frameinvariancetest.cc b/test/frameinvariancetest.cc index 95180a69423e3c017acd89001dd6c8c8500f0c7f..f9a355a8aaaf7ca687124359f478d12474f8242d 100644 --- a/test/frameinvariancetest.cc +++ b/test/frameinvariancetest.cc @@ -23,8 +23,6 @@ using std::string; int main (int argc, char *argv[]) try { // Some types that I need - typedef BCRSMatrix<FieldMatrix<double, blocksize, blocksize> > MatrixType; - typedef BlockVector<FieldVector<double, blocksize> > CorrectionType; typedef std::vector<RigidBodyMotion<double,3> > SolutionType; // Problem settings @@ -35,8 +33,10 @@ int main (int argc, char *argv[]) try // /////////////////////////////////////// typedef OneDGrid GridType; GridType grid(numRodBaseElements, 0, 1); + using GridView = GridType::LeafGridView; + GridView gridView = grid.leafGridView(); - SolutionType x(grid.size(1)); + SolutionType x(gridView.size(1)); // ////////////////////////// // Initial solution @@ -57,10 +57,7 @@ int main (int argc, char *argv[]) try // Create a second, rotated copy of the configuration // ///////////////////////////////////////////////////////////////////// - FieldVector<double,3> displacement; - displacement[0] = 0; - displacement[1] = 0; - displacement[2] = 0; + FieldVector<double,3> displacement {0, 0, 0}; FieldVector<double,3> axis(0); axis[0]=1; Rotation<double,3> rotation(axis,M_PI/2); @@ -88,20 +85,20 @@ int main (int argc, char *argv[]) try writeRod(x,"rod"); writeRod(rotatedX, "rotated"); - RodLocalStiffness<GridType::LeafGridView,double> assembler(grid.leafGridView(), + RodLocalStiffness<GridView,double> assembler(gridView, 1,1,1,1e6,0.3); for (int i=1; i<2; i++) { double p = double(i)/2; - assembler.getStrain(x,*grid.lbegin<0>(0), p); - assembler.getStrain(rotatedX,*grid.lbegin<0>(0), p); + assembler.getStrain(x,*gridView.begin<0>(), p); + assembler.getStrain(rotatedX,*gridView.begin<0>(), p); } } catch (Exception e) { - std::cout << e << std::endl; + std::cout << e.what() << std::endl; } diff --git a/test/globalgfetestfunctionbasistest.cc b/test/globalgfetestfunctionbasistest.cc index beb1d7197b7c233dd603ebd2ff31243f93f39015..16e8cc5399c75cbbdcc4601652760d5a8a79cc0b 100644 --- a/test/globalgfetestfunctionbasistest.cc +++ b/test/globalgfetestfunctionbasistest.cc @@ -41,19 +41,17 @@ void test() typedef GlobalGFETestFunctionBasis<P1Basis,TargetSpace> GlobalBasis; GlobalBasis basis(p1Basis,testPoints); - typedef typename OneDGrid::Codim<0>::LeafIterator ElementIterator; - ElementIterator eIt = grid.leafbegin<0>(); - ElementIterator eEndIt = grid.leafend<0>(); - - for (; eIt != eEndIt; ++eIt) { - - const typename GlobalBasis::LocalFiniteElement& lfe = basis.getLocalFiniteElement(*eIt); + for (const auto element : elements(grid.leafGridView())) + { + const typename GlobalBasis::LocalFiniteElement& lfe = basis.getLocalFiniteElement(element); FieldVector<double,1> stupidTestPoint(0); std::vector<std::array<typename TargetSpace::EmbeddedTangentVector, TargetSpace::TangentVector::dimension> > values; lfe.localBasis().evaluateFunction(stupidTestPoint, values); for(size_t i=0;i<values.size();i++) { - std::cout<<values[i]<<std::endl; + for (auto v : values[i]) + std::cout << v << " "; + std::cout << std::endl; std::cout<<lfe.localCoefficients().localKey(i)<<std::endl; } //int i = basis.index(*eIt,1);