These constructors and assignments from base class are
all using a downcast which may result in undefined behaviour.
Since this is potentially dangerous and not used anywhere
in the core modules, this patch removes those methods.

These constructors and assignments are all implemented using
a downcast to derived class and may lead to undefined behaviour.
This commit removes the checks for these methods in preparation
for their removal.

The doc now correctly describes
* that the overflow area is used during compress and not during assembly
* the meaning of avg and overflowsize
* the numbering dependence of this build mode
* the magic number 4
* the late failure during compress()

Stephan Hilb authored 8 years ago and Christoph Grüninger committed 8 years ago

This fixes the clang compiler warnings about tautological compare
when idxtype is unsigned. For clang behaviour see [1].

[1]: https://llvm.org/bugs/show_bug.cgi?id=8682

This reverts ddf49c23 and removes the check if
the average parameter behaves as expected.

similar to Matlab matrix writer

Remove stamp-vc file.

Since the istl view is, that a matrix is a container of rows,
it should have a size method.

This is in fact an _unsafe_ downcast to a derived
class. Accessing members via this pointer is undefined
if it is not guaranteed that the passed pointer
points to an object of the _derived_ class - even
if the memory layout is the same. In this case
this static_cast acts like a reinterpret_cast.

If the correct average is provided, compress() may still
fail if more than the average numebr of entries are added
to an early row (e.g. the first one)

Due to a confusion between the number of lines N and the number of
columns M, Matrix::transpose() gave wrong results for matrices with
N != M.

Due to a bug in Matrix::transpose() it will give wrong results for
non-quadratic matrices. Thus the new test case fails for the moment.

Remove output that printed the matrix for parallel runs for small grids n<20.

According to git blame the last time these constructors were changed was
in 2011.

Simplifies the memory management.

Reimplement this using Hybrid::forEach. Along the way this
also fixes the bug that operator<< ignores the passed stream
and uses std::cout instead.

The new implementation is based on Hybrid::forEach

Test with non-zero matrix and print results

Reimplement this using Hybrid::forEach. This also adds some
comments and uses proper variable names to make the nan-tracking
easier to understand.

Reimplement this using Hybrid::forEach.

Reimplement this using Hybrid::forEach.

Reimplement this using Hybrid::forEach.

These separate impelmentations could be replaced by a single one,
once we have IsNumber. Currently we could use std::is_arithmetic
instead, but that's not the topic of this change.

Reimplement this using Hybrid::forEach.

Reimplement this using Hybrid::forEach. Along the way this
also fixes the bug that operator<< ignores the passed stream
and uses std::cout instead.

The implementation in MultiTypeBlockMatrix would be even
simpler, if it would provide a static size() method, because
we could then use a simple Hybrid::forEach over its entries
instead of one over all indices.

Except for those parts that do still rely on the enable trick to
test a certain feature set.

* Use correct macro
* adjust changes for other types then double, too

Arne Morten Kvarving authored 8 years ago and Christoph Grüninger committed 8 years ago

API for solver routine changed. Since there is no versioning
in SuperLU we have to do a probe to detect the new library.

Here we try to solve A*x=b for x, with
```
      / 1  1 \       / 1 \
  A = |      |,  b = |   |,
      \ 1  1 /       \ 2 /
```
which has no solution.  In the CGSolver, the this leads to a NaN defect after
46 iterations.  In BiCGSTABSolver this leads to "abs(h) < EPSILON" after 1.5
iterations.

BiCGSTABSolver and RestartedGMResSolver also throw SolverAbort when they
detect "breakdown", checking for that is left as a todo.

Fixes dune-common#23