dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Nov 25, 2010 at 05:07:40PM +0000, Garth N. Wells wrote:
> There is a dependency on
>
>    dolfin/ode/SORSolver.h
>
> in dolfin/common/real.cpp. This seems odd to me. Is there a reason for it?
>
> Garth

real.h/cpp provide some basic algorithms that work on high precision
floats. One of these is the matrix exponential which is computed by a
rational approximation, hence the need for solving a linear system.

The SORSolver is the only solver in DOLFIN that handles high precision
floats (uBLAS etc can't be used). It would be more natural to put it
under dolfin/la, although it would still be a bit non-optimal to have
a dependency from common to la. Any other suggestion?

Benjamin has quite a few improvements to the ODE solvers and the
multiprecision in a separate branch (on Launchpad). I will merge it
back soon, but it would be good to sort this out first.

--
Anders