Skip to content
Snippets Groups Projects
Commit 45e5e725 authored by Markus Blatt's avatar Markus Blatt
Browse files

Forgot to check my changes of the test. Hopefully this fixes

Christoph's issues.

[[Imported from SVN: r1562]]
parent 294525e2
No related branches found
No related tags found
No related merge requests found
......@@ -8,6 +8,30 @@
#include <laplacian.hh>
#include <dune/common/timer.hh>
#include <dune/istl/superlu.hh>
#include <complex>
#ifndef SUPERLU_NTYPE
#define SUPERLU_NTYPE 1
#endif
#if SUPERLU_NTYPE==1
typedef double FIELD_TYPE;
#endif
#if SUPERLU_NTYPE==0
typedef float FIELD_TYPE;
#endif
#if SUPERLU_NTYPE==2
typedef std::complex<float> FIELD_TYPE;
#endif
#if SUPERLU_NTYPE>=3
typedef std::complex<double> FIELD_TYPE;
#endif
int main(int argc, char** argv)
{
......@@ -19,9 +43,9 @@ int main(int argc, char** argv)
std::cout<<"testing for N="<<N<<" BS="<<1<<std::endl;
typedef Dune::FieldMatrix<double,BS,BS> MatrixBlock;
typedef Dune::FieldMatrix<FIELD_TYPE,BS,BS> MatrixBlock;
typedef Dune::BCRSMatrix<MatrixBlock> BCRSMat;
typedef Dune::FieldVector<double,BS> VectorBlock;
typedef Dune::FieldVector<FIELD_TYPE,BS> VectorBlock;
typedef Dune::BlockVector<VectorBlock> Vector;
typedef Dune::MatrixAdapter<BCRSMat,Vector,Vector> Operator;
......
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