[Dune] [#571] MonomLocalInterpolation.interpolate() is extremly unstable

[Dune]

THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

The following task has a new comment added:

FS#571 - MonomLocalInterpolation.interpolate() is extremly unstable
User who did this - Carsten Gräser (Carsten)

----------
Unfortunately, I have no alternative implementation for dim>1. I just testet this with matlab in 1d.

There are a few papers on the 'stable' inversion of the Vandermonde matrix (see english wikipedia entry). I think the basic idea is to compute the monomial coefficients from a stable interpolation using a different basis. Perhaps one could embed a simplex in the geometry and use Lagrange interpolation first. However, I'm not sure if the efford (esp. for arbitrary dim) is worth it.

If you're happy with the current accuracy one could also document the instability and perhaps include a warning for e.g. order>5.
----------

More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=571#comment1000

You are receiving this message because you have requested it from the Flyspray bugtracking system.  If you did not expect this message or don't want to receive mails in future, you can change your notification settings at the URL shown above.