#include <deal.II/lac/precondition.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/precondition.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
-// test the various VECTOR::extract_subvector_to functions
+// test the various VectorType::extract_subvector_to functions
#include "../tests.h"
#include <deal.II/lac/generic_linear_algebra.h>
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
{};
void
- local_mass_operator (const MatrixFree<dim,Number> &data,
- VECTOR &dst,
- const VECTOR &src,
+ local_mass_operator (const MatrixFree<dim,Number> &data,
+ VectorType &dst,
+ const VectorType &src,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,1,Number> fe_eval (data);
}
void
- local_inverse_mass_operator (const MatrixFree<dim,Number> &data,
- VECTOR &dst,
- const VECTOR &src,
+ local_inverse_mass_operator (const MatrixFree<dim,Number> &data,
+ VectorType &dst,
+ const VectorType &src,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,1,Number> fe_eval (data);
}
}
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_mass_operator,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_mass_operator,
this, dst, src);
};
- void apply_inverse (VECTOR &dst,
- const VECTOR &src) const
+ void apply_inverse (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_inverse_mass_operator,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_inverse_mass_operator,
this, dst, src);
};
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
{};
void
- local_mass_operator (const MatrixFree<dim,Number> &data,
- VECTOR &dst,
- const VECTOR &src,
+ local_mass_operator (const MatrixFree<dim,Number> &data,
+ VectorType &dst,
+ const VectorType &src,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,dim,Number> fe_eval (data);
}
void
- local_inverse_mass_operator (const MatrixFree<dim,Number> &data,
- VECTOR &dst,
- const VECTOR &src,
+ local_inverse_mass_operator (const MatrixFree<dim,Number> &data,
+ VectorType &dst,
+ const VectorType &src,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,dim,Number> fe_eval (data);
}
}
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_mass_operator,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_mass_operator,
this, dst, src);
};
- void apply_inverse (VECTOR &dst,
- const VECTOR &src) const
+ void apply_inverse (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_inverse_mass_operator,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_inverse_mass_operator,
this, dst, src);
};
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
{};
void
- local_mass_operator (const MatrixFree<dim,Number> &data,
- VECTOR &dst,
- const VECTOR &src,
+ local_mass_operator (const MatrixFree<dim,Number> &data,
+ VectorType &dst,
+ const VectorType &src,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,3,Number> fe_eval (data);
}
void
- local_inverse_mass_operator (const MatrixFree<dim,Number> &data,
- VECTOR &dst,
- const VECTOR &src,
+ local_inverse_mass_operator (const MatrixFree<dim,Number> &data,
+ VectorType &dst,
+ const VectorType &src,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,3,Number> fe_eval (data);
}
}
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_mass_operator,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_mass_operator,
this, dst, src);
};
- void apply_inverse (VECTOR &dst,
- const VECTOR &src) const
+ void apply_inverse (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_inverse_mass_operator,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_inverse_mass_operator,
this, dst, src);
};
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
- MatrixFreeTest(const DoFHandler<dim> &dof_handler,
+ MatrixFreeTest(const DoFHandler<dim> &dof_handler,
const ConstraintMatrix &constraints)
:
dof_handler (dof_handler),
constraints (constraints)
{}
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
- VECTOR src_cpy = src;
+ VectorType src_cpy = src;
constraints.distribute(src_cpy);
FEEvaluation<dim,fe_degree,fe_degree+1,1,Number>
fe_eval(dof_handler.get_fe(), QGauss<1>(fe_degree+1),
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
- MatrixFreeTest(const DoFHandler<dim> &dof_handler,
- const ConstraintMatrix &constraints)
+ MatrixFreeTest (const DoFHandler<dim> &dof_handler,
+ const ConstraintMatrix &constraints)
:
dof_handler (dof_handler),
constraints (constraints)
{}
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
- VECTOR src_cpy = src;
+ VectorType src_cpy = src;
constraints.distribute(src_cpy);
FEEvaluation<dim,fe_degree,fe_degree+1,dim,Number> fe_eval
(dof_handler.get_fe(), QGauss<1>(fe_degree+1),
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
data (data_in)
{};
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VECTOR>::local_operation,
+ data.cell_loop (&MatrixFreeTest<dim,fe_degree,Number,VectorType>::local_operation,
this, dst, src);
};
private:
const MatrixFree<dim,Number> &data;
- void local_operation (const MatrixFree<dim,Number> &data,
- VECTOR &out,
- const VECTOR &in,
+ void local_operation (const MatrixFree<dim,Number> &data,
+ VectorType &out,
+ const VectorType &in,
const std::pair<unsigned int,unsigned int> &cell_range) const
{
FEEvaluation<dim,fe_degree,fe_degree+1,1,Number> fe_eval (data);
#include <deal.II/lac/parallel_vector.h>
-template <int dim, int fe_degree, typename VECTOR>
+template <int dim, int fe_degree, typename VectorType>
void
-helmholtz_operator (const MatrixFree<dim,typename VECTOR::value_type> &data,
- VECTOR &dst,
- const VECTOR &src,
- const std::pair<unsigned int,unsigned int> &cell_range)
+helmholtz_operator (const MatrixFree<dim,typename VectorType::value_type> &data,
+ VectorType &dst,
+ const VectorType &src,
+ const std::pair<unsigned int,unsigned int> &cell_range)
{
- typedef typename VECTOR::value_type Number;
+ typedef typename VectorType::value_type Number;
FEEvaluation<dim,fe_degree,fe_degree+1,1,Number> fe_eval (data);
const unsigned int n_q_points = fe_eval.n_q_points;
-template <int dim, int fe_degree, typename Number, typename VECTOR=Vector<Number> >
+template <int dim, int fe_degree, typename Number, typename VectorType=Vector<Number> >
class MatrixFreeTest
{
public:
data (data_in)
{};
- void vmult (VECTOR &dst,
- const VECTOR &src) const
+ void vmult (VectorType &dst,
+ const VectorType &src) const
{
dst = 0;
- const std_cxx11::function<void(const MatrixFree<dim,typename VECTOR::value_type> &,
- VECTOR &,
- const VECTOR &,
+ const std_cxx11::function<void(const MatrixFree<dim,typename VectorType::value_type> &,
+ VectorType &,
+ const VectorType &,
const std::pair<unsigned int,unsigned int> &)>
- wrap = helmholtz_operator<dim,fe_degree,VECTOR>;
+ wrap = helmholtz_operator<dim,fe_degree,VectorType>;
data.cell_loop (wrap, dst, src);
};
#include <fstream>
#define N 3
-typedef Vector<double> VECTOR;
+typedef Vector<double> VectorType;
class MGAll
:
- public MGSmootherBase<VECTOR>,
- public MGTransferBase<VECTOR>,
- public MGCoarseGridBase<VECTOR>
+ public MGSmootherBase<VectorType>,
+ public MGTransferBase<VectorType>,
+ public MGCoarseGridBase<VectorType>
{
public:
virtual ~MGAll()
{}
virtual void smooth (const unsigned int,
- VECTOR &, const VECTOR &) const
+ VectorType &, const VectorType &) const
{}
virtual void prolongate (const unsigned int,
- VECTOR &, const VECTOR &) const
+ VectorType &, const VectorType &) const
{}
virtual void restrict_and_add (const unsigned int,
- VECTOR &, const VECTOR &) const
+ VectorType &, const VectorType &) const
{}
virtual void clear ()
{}
virtual void operator() (const unsigned int,
- VECTOR &, const VECTOR &) const
+ VectorType &, const VectorType &) const
{}
};
MGLevelObject<FullMatrix<double> > level_matrices(0, maxlevel);
for (unsigned int i=0; i<=maxlevel; ++i)
level_matrices[i].reinit(N, N);
- mg::Matrix<VECTOR> mgmatrix(level_matrices);
+ mg::Matrix<VectorType> mgmatrix(level_matrices);
- Multigrid<VECTOR> mg1(minlevel, maxlevel, mgmatrix, all, all, all, all,
- Multigrid<VECTOR>::v_cycle);
+ Multigrid<VectorType> mg1(minlevel, maxlevel, mgmatrix, all, all, all, all,
+ Multigrid<VectorType>::v_cycle);
mg1.set_debug(3);
for (unsigned int i=minlevel; i<=maxlevel; ++i)
mg1.defect[i].reinit(N);
mg1.cycle();
deallog << std::endl;
- mg1.set_cycle(Multigrid<VECTOR>::w_cycle);
+ mg1.set_cycle(Multigrid<VectorType>::w_cycle);
mg1.cycle();
deallog << std::endl;
- mg1.set_cycle(Multigrid<VECTOR>::f_cycle);
+ mg1.set_cycle(Multigrid<VectorType>::f_cycle);
mg1.cycle();
}
/**
* Apply preconditioner.
*/
- template<class VECTOR>
- void vmult (VECTOR &, const VECTOR &) const;
+ template<typename VectorType>
+ void vmult (VectorType &, const VectorType &) const;
/**
* Apply transpose
* the same as
* vmult().
*/
- template<class VECTOR>
- void Tvmult (VECTOR &, const VECTOR &) const;
+ template<typename VectorType>
+ void Tvmult (VectorType &, const VectorType &) const;
/**
* Apply preconditioner, adding to the previous value.
*/
- template<class VECTOR>
- void vmult_add (VECTOR &, const VECTOR &) const;
+ template<typename VectorType>
+ void vmult_add (VectorType &, const VectorType &) const;
/**
* Apply transpose
* the same as
* vmult_add().
*/
- template<class VECTOR>
- void Tvmult_add (VECTOR &, const VECTOR &) const;
+ template<typename VectorType>
+ void Tvmult_add (VectorType &, const VectorType &) const;
private:
number factor;
template<typename number>
-template<class VECTOR>
+template<typename VectorType>
inline void
-ScalingMatrix<number>::vmult (VECTOR &dst, const VECTOR &src) const
+ScalingMatrix<number>::vmult (VectorType &dst, const VectorType &src) const
{
dst.equ(factor, src);
}
template<typename number>
-template<class VECTOR>
+template<typename VectorType>
inline void
-ScalingMatrix<number>::Tvmult (VECTOR &dst, const VECTOR &src) const
+ScalingMatrix<number>::Tvmult (VectorType &dst, const VectorType &src) const
{
dst.equ(factor, src);
}
template<typename number>
-template<class VECTOR>
+template<typename VectorType>
inline void
-ScalingMatrix<number>::vmult_add (VECTOR &dst, const VECTOR &src) const
+ScalingMatrix<number>::vmult_add (VectorType &dst, const VectorType &src) const
{
dst.add(factor, src);
}
template<typename number>
-template<class VECTOR>
+template<typename VectorType>
inline void
-ScalingMatrix<number>::Tvmult_add (VECTOR &dst, const VECTOR &src) const
+ScalingMatrix<number>::Tvmult_add (VectorType &dst, const VectorType &src) const
{
dst.add(factor, src);
}
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
+check_solve(SOLVER &solver,
const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;
#include <deal.II/lac/vector_memory.h>
#include <typeinfo>
-template<class SOLVER, class MATRIX, class VECTOR, class PRECONDITION>
+template<class SOLVER, class MATRIX, typename VectorType, class PRECONDITION>
void
-check_solve( SOLVER &solver,
- const SolverControl &solver_control,
- const MATRIX &A,
- VECTOR &u, VECTOR &f, const PRECONDITION &P)
+check_solve(SOLVER &solver,
+ const SolverControl &solver_control,
+ const MATRIX &A,
+ VectorType &u,
+ VectorType &f,
+ const PRECONDITION &P)
{
deallog << "Solver type: " << typeid(solver).name() << std::endl;