This patch removes the assumption that the field_type and / or real_type
can always be implicitly converted to double and can be implicitly
constructed from an int (in constructs like field_type d = 0;).

This patch is required for compatibility with the FLOP counter in PDELab.

It belongs to dune-common's config.h.cmake. Otherwise
it cannot be shared with other modules using ParMETIS.

It should be std::endl.
Courtesy of GCC 5-svn, it refuses to compile these.

While c++11 allows a trailing comma in enums,
we should avoid this because it provokes
gcc-4.4 warnings.

The check gets removed from dune-common, thus ISTL has to call it itself.

This solves warnings in UG for redefined macros.
Thanks to Christian for hinting the proper way
fixing this.

Jonathan Youett authored 10 years ago and Markus Blatt committed 10 years ago

When assigning a BCRSMatrix to another BCRSMatrix which was setup for
a different size/pattern, it segfaulted and I got the following
Exception:

"InvalidStateException [allocate:.../dune/istl/bcrsmatrix.hh:2108]:
Rows have already been allocated, cannot allocate a second time."

The problem was, that inside the assignment operator the rows (r) were
never deallocated. This commit fixes this and adds a testcase.

Testcase slightly modified by Markus.

This mode is not implemented and therefore throws
an assertion that makes the test fail everytime.

The elapsed wallclock time is readily available by UMFPack
and can be written to InverseOperatorResult.elapsed!

Namely we try to document the template parameters in more detail and
make sure that cross reference to OwnerOverlapCopyCommunication,
vectors, and matrix implemetation appear.

Fix broken handling of complex problems in MINRES and GMRes. With these
patches, both solvers now support complex matrices as well as complex
right-hand sides. Due to a couple of bugs, the convergence rate of
complex problems was off even before the current breakage.

* feature/fs1467-fix-gmres-and-minres-for-complex:
  [Documentation] Add file comment and formatting hints to complexdata.hh
  [complexmatrixtest] Rename assemble_random_matrix_vectors.hh -> complexdata.hh
  [CMake] Properly exclude header from headercheck.
  [Buildsystem] Ignore helper header for complexmatrixtest
  [Bugfix] GMRes now needs enable_if, so include common/typetraits.hh
  [Tests] Move tests to exercise MINRES and GMRes on complex matrices into separate executable
  [Cleanup][Tests] Rearranged the automatically parsed code for the assembling of the random matrices and random vectors in assemble_random_matrix_vectors.hh for better readability
  [Autotools] added new header file in the Makefile.am
  [Tests] added tests using complex nonsingular matrices and complex hermitian matrices to test the MinRes and GMRes solver on purely complex problems
  [Bugfix] Exploit analytical properties in MINRES to avoid special-casing for complex case
  [Bugfix] fixed order of the scalar product when applying the Arnoldi algorithm in GMRes solver, which is important in the complex case
  [Bugfix] implemented complex Givens rotations and fixed MinRes and GMRes for complex valued problems

As Christoph pointed out in FS1467, that's a rather unusual name for a
DUNE header file, so make it a little more consistent.

Add newline at and of file.

That header is not supposed to be used on its own, so it makes no sense
to check whether that is possible.

The recently added tests for complex matrices added to complexrhstest
made that file really large and unwieldy, so this patch moves those
tests to a new file complexmatrixtest.

[Cleanup][Tests] Rearranged the automatically parsed code for the assembling of the random matrices and random vectors in assemble_random_matrix_vectors.hh for better readability

[Tests] added tests using complex nonsingular matrices and complex hermitian matrices to test the MinRes and GMRes solver on purely complex problems

[Bugfix] fixed order of the scalar product when applying the Arnoldi algorithm in GMRes solver, which is important in the complex case

minor cleanups in the MinRes solver

[Bugfix] implemented complex Givens rotations and fixed MinRes and GMRes for complex valued problems

The type idxtype is used by header repartition even if parmetis is not
used. Therefore it has to be typedefed meaningful in this case,
too. This was not the case before and is accomplished by this commit.

This was introduced for debugging only and could have side effects
onto other headers. Removes the accidentally committed line.

* feature/fix-gmres-and-minres:
  [GMRes] Fix deprecation warning for old constructors
  made the distinction between field_type and real_type which is the same for real numbers but differs for std::complex
  added the old constructors called with the variable recalc_defect to RestartedGMResSolver and set them as deprecated, since they should not be used anymore
  changed the initialization of the arrays in the MinRes solver to remove the compiler warnings
  fixed MinRes with the following changes:
  fixed GMRes

Conflicts:
	dune/istl/solvers.hh

Function attribute have to be placed in front of the function name,
not behind it.

Marian Piatkowski authored 10 years ago and Steffen Müthing committed 10 years ago

made the distinction between field_type and real_type which is the same for real numbers but differs for std::complex

Marian Piatkowski authored 10 years ago and Steffen Müthing committed 10 years ago

added the old constructors called with the variable recalc_defect to RestartedGMResSolver and set them as deprecated, since they should not be used anymore

Marian Piatkowski authored 11 years ago and Steffen Müthing committed 10 years ago

replaced the usual arrays for the Givens-rotation coefficients and the search directions etc. with the arrays from Dune

implemented the BLAS-routine drotg for greater robustness which was a TODO note to compute the Givens rotation

Marian Piatkowski authored 11 years ago and Steffen Müthing committed 10 years ago

removed bool recalc_defect
reimplement the update function with respect to optimality
fixed strange behaviour when maxiter is reached and GMRes wants to restart infinitely many times
fixed strange behaviour when linear reduction is reached but GMRes wants to restart and continue calculation
fixed the printing of the results