The point of the matter is that posix_memalign() is not allowed in this
context according to the relevant standards, so we use mmap() instad.

`mprotect` can only be used for `mmap`'ed memory

Makes testdebugallocator pass on OS X, too.

See merge request !43
(cherry picked from commit 6467775c)

Signed-off-by: Steffen Müthing <muething@dune-project.org>

This addresses warnings such as

  'const' type qualifier on return type has no effect

(cherry picked from commit 1de4cb89)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

feature/use MPI_Type_create_hindexed instead of MPI_Type_hindexed which is deprecated for MPI>=2

:)

See merge request !29
(cherry picked from commit dea4e7d0)

Signed-off-by: Steffen Müthing <muething@dune-project.org>

Fix FieldMatrix::rightmultiply()

This addresses an issue mentioned in #5 and additionally turns runtime checks into compile-time checks whenever possible.

See merge request !22
(cherry picked from commit 71e2e4b8)

Signed-off-by: Steffen Müthing <muething@dune-project.org>

Feature/fix shadow warnings

If you have better name suggestions, please tell me and I change the names.

See merge request !21
(cherry picked from commit 394c5cde)

Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 2abb3bd4)

Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 47830f03)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit c232686d)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 2125b283)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

Test complex<T> for floating point T only

The standard leaves (the effect of instantiating) complex<T> for more
general T unspecified.

While we're at it, kill duplicate tests and add long double to the list.

This fixes bug flyspray/FS#1726: Now also with clang 3.7 and libstdc++ 4.9.2 on Debian 8.2, the test fvectortest can be compiled again.

See merge request !11
(cherry picked from commit 627fc0b0)

Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 49879f16)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 1906105f)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 28762f7b)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

The section is empty since the code was split into the dune-geometry
module.

(cherry picked from commit c37cd9cb)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 60459f97)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

floating comparisons

(cherry picked from commit 94a327ac)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

(cherry picked from commit 7369e156)
Signed-off-by: Steffen Müthing <muething@dune-project.org>

by testing whether we work on a end iterator. In this case the distance
function will fail. Found that out by gcc's stdlib with debugging turned on.

(cherry picked from commit c494c58b)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 9588afdf)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 3bef4991)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 6121bbd9)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

... and don't take absolute value of comparison result.
Thanks to Max Kahnt for reporting this.

(cherry picked from commit 57a43efa)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

It's more CMake style.

(cherry picked from commit 635500eb)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit a46de680)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

IteratorRanges are convenient for storing ranges of iterators into a
larger container, e.g. for the per-codim GeometryTypes in an
IndexSet. Unfortunately, you can't just create an array of
IteratorRanges right now because they can neither be copied, nor are
they default-constructible.

This patch fixes that problem and also adds a second typedef
const_iterator to make IteratorRange work in contexts where the
iterating code isn't using range-based for and assumes a const object.

(cherry picked from commit 2c9596f9)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

With the recent changes to bigunsignedint, the index cannot
be initialized with a negative value anymore. We change this
to 1 to overcome the problem.

This commit together with the recent bugfixes by Ansgar
fixes the AMG code on master.

(cherry picked from commit bf023947)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 6100b778)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

Since 0c697d8d messages of size zero
are no longer sent. This can reduce the number of requests we have to
wait for.

(cherry picked from commit f35437e9)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

The storage of bigunsignedint was changed to std::uint16_t in commit
541b3dfc. This implements the same
change for the MPITraits<bigunsignedint<k>>

Note that this assumes that
  MPITraits<std::uint16_t>
is defined. This is the case as long as uint16_t maps to one of the
standard integer types (unsigned char/short/int/long).

Thanks to Markus for pointing this out.

Reference: https://dune-project.org/flyspray/index.php?do=details&task_id=1657#comment6393


(cherry picked from commit 586375b4)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 648dff2b)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

This bug could only be found outside Germany :D

(cherry picked from commit ecdb7f7b)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

the first problem was that I accidentially enabled the ScalarOrderingTest, but this forced me
to investigate further and I found out how this works, or not works.

We have to add the comparison operators to FieldVector, as the cast operator does not always
match, e.g. in the case of GMPField.
(cherry picked from commit e0ac8d9a)

Signed-off-by: Carsten Gräser <graeser@dune-project.org>

Before the test was not active since '#if HAVE_GMP' fails
if the gmp-flags where not added.

(cherry picked from commit 045c6d93)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

in order to make GMPField work in more situations, we add a constrctor, which matches
all scalar data types, which can me converted to K.
We deliberately exclude K from the template match.

(cherry picked from commit 9d46150e)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 7135cbb0)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

we first import function like std::max into the local scope and the call max(a,b) instead of
std::max(a,b). This allows to match specializations for other data types via ADL.

(cherry picked from commit 515f312a)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

(cherry picked from commit 49200b29)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>

The templated constructor of GMPField fatched everything. We reduce the constructor
to those scalar values, which are convertible to mpf_class. As this removes the constructor
from char* we explicitly add constructor from char* and std::string.

This is the main part to fix FS#1665

(cherry picked from commit 3fe52527)
Signed-off-by: Carsten Gräser <graeser@dune-project.org>