[examples.elliptic.continuousgalerkin] update

examples/elliptic/continuousgalerkin.cc

@@ -214,8 +214,8 @@ int main(int argc, char** argv)
     info << "initializing matrix (of size " << testSpace.map().size() << "x" << ansatzSpace.map().size()
          << ") and vectors... " << std::flush;
     timer.reset();
-    typedef ContainerFactory::SparseMatrixType MatrixType;
-    Dune::shared_ptr<MatrixType> systemMatrix   = ContainerFactory::createSparseMatrix(testSpace, ansatzSpace);
+    typedef ContainerFactory::RowMajorSparseMatrixType MatrixType;
+    Dune::shared_ptr<MatrixType> systemMatrix   = ContainerFactory::createRowMajorSparseMatrix(testSpace, ansatzSpace);
     Dune::shared_ptr<VectorType> forceVector    = ContainerFactory::createDenseVector(testSpace);
     Dune::shared_ptr<VectorType> neumannVector  = ContainerFactory::createDenseVector(testSpace);
     Dune::shared_ptr<VectorType> rhsVector      = ContainerFactory::createDenseVector(testSpace);
@@ -264,8 +264,8 @@ int main(int argc, char** argv)
     const double solverPrecision     = paramTree.get("solver.precision", 1e-12);
     info << "  using '" << solverType << "'... " << std::flush;
     timer.reset();
-    typedef typename Dune::Stuff::LA::Solver::Eigen::Interface<MatrixType> SolverType;
-    SolverType* solver = Dune::Stuff::LA::Solver::Eigen::create<MatrixType>(solverType);
+    typedef typename Dune::Stuff::LA::Solver::Eigen::Interface<MatrixType, VectorType> SolverType;
+    SolverType* solver = Dune::Stuff::LA::Solver::Eigen::create<MatrixType, VectorType>(solverType);
     solver->init(*systemMatrix);
     const bool success = solver->apply(*rhsVector, *solutionVector, solverMaxIter, solverPrecision);
     if (!success)