Resolved conflicts:
	dune/istl/paamg/fastamg.hh

Replace it with	typename FieldTraits<field_type>::real_type

This fixes two problems:
- norms don't return complex number types when instantiated
  for complex coefficients
- You can now instantiate BlockVector for field_types that do not
  auto-convert to 'double'

A reference is an lvalue. Assigning to it later changes the state
of the rhe referencee (and not the reference). Therefore the code
actually stored references to the right hand and left hand side.
In the case where we did not want to reuse vectors we created two
dense SuperMatrices to used in this apply call and then assigned
them to the references. In fact we actually overwrote the right
and left hand side by doing this.

This patch uses pointers to accompish the intended behaviour, i.e.
point to the correct dense matrices, the ones stored in SuperLU
or the ones stored in SuperLU::apply.

Arne Morten Kvarving authored 11 years ago and Markus Blatt committed 11 years ago

A reference is an lvalue. Assigning to it later changes the state
of the rhe referencee (and not the reference). Therefore the code
actually stored references to the right hand and left hand side.
In the case where we did not want to reuse vectors we created two
dense SuperMatrices to used in this apply call and then assigned
them to the references. In fact we actually overwrote the right
and left hand side by doing this.

This patch uses pointers to accompish the intended behaviour, i.e.
point to the correct dense matrices, the ones stored in SuperLU
or the ones stored in SuperLU::apply.

Otherwise I get name clashes: DUNE_SOLVER_HH is used in dune-istl
and dune-solvers.

Andreas Dedner authored 11 years ago and Markus Blatt committed 11 years ago

removed again - sorry about that.

Unfortunately, SuperLU always overwrites the right hand side with
the solution. Therefore we now copy it if there is only one level
in the matrix hierarchy of AMG.

Christian Engwer authored 11 years ago and Markus Blatt committed 11 years ago

modify constructor with capacity to catch at least some of the
wrong invocations (see doxygen comment).

Oliver Sander authored 11 years ago and Markus Blatt committed 11 years ago

This fixes FS 1292, but with my own slightly different patch.

removed again - sorry about that.

Unfortunately, SuperLU always overwrites the right hand side with
the solution. Therefore we now copy it if there is only one level
in the matrix hierarchy of AMG.

modify constructor with capacity to catch at least some of the
wrong invocations (see doxygen comment).

This fixes FS 1292, but with my own slightly different patch.

The last commits caused compile errors for e.g. clang as I used
incomplete (only forward declared) types in template classes.
This patch moves the definition of InverseOperator and Preconditioners
from the headers solvers.hh and preconditioners into new separate
headers solver.hh and preconditioner.hh. This allows to include these
definitions in the previous headers without cause circular include directives.

This fixes FS #1224

Previously if the solver did not converge the number of iterations
that was printed and stored was one higher than the actual value.

With OpenMPI there appeared cases where MPI_Comm_free
was called after MPI_Finalize. This was caused by the
destructor of OwnerOverlapCopyCommunication being called
after MPI_Finalize.
This patch moves the destruction out of the main method.
If MPI 2 functionality is available we also check whether
MPI_Finalize was called before freeing the communicator and
only call free if not.

This patch introduces shared_ptrs at various places,
namely for the matrix hierarchies and coarse solvers.
This allows copying an AMG and using the same matrix
hierarchies in both preconditioners. Thus one could
solve the same system for various right hand sides,
where each solve happens in a different thread.