Skip to content
Snippets Groups Projects
Commit f0fa230e authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

Make BasisType a template parameter

[[Imported from SVN: r9967]]
parent 2f73f105
No related branches found
No related tags found
No related merge requests found
...@@ -20,11 +20,13 @@ ...@@ -20,11 +20,13 @@
#include <dune/gfe/feassembler.hh> #include <dune/gfe/feassembler.hh>
/** \brief Trust-region solver */ /** \brief Trust-region solver */
template <class GridType, class VectorType> template <class BasisType, class VectorType>
class TrustRegionSolver class TrustRegionSolver
: public IterativeSolver<VectorType, : public IterativeSolver<VectorType,
Dune::BitSetVector<VectorType::value_type::dimension> > Dune::BitSetVector<VectorType::value_type::dimension> >
{ {
typedef typename BasisType::GridView::Grid GridType;
const static int blocksize = VectorType::value_type::dimension; const static int blocksize = VectorType::value_type::dimension;
const static int gridDim = GridType::dimension; const static int gridDim = GridType::dimension;
...@@ -37,14 +39,6 @@ class TrustRegionSolver ...@@ -37,14 +39,6 @@ class TrustRegionSolver
typedef Dune::BlockVector<Dune::FieldVector<field_type, blocksize> > CorrectionType; typedef Dune::BlockVector<Dune::FieldVector<field_type, blocksize> > CorrectionType;
typedef VectorType SolutionType; typedef VectorType SolutionType;
#ifdef THIRD_ORDER
typedef P3NodalBasis<typename GridType::LeafGridView,double> BasisType;
#elif defined SECOND_ORDER
typedef P2NodalBasis<typename GridType::LeafGridView,double> BasisType;
#else
typedef P1NodalBasis<typename GridType::LeafGridView,double> BasisType;
#endif
public: public:
TrustRegionSolver() TrustRegionSolver()
......
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