On 25 November 2010 18:34, Anders Logg <[email protected]> wrote: > 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?
Should the implementation of the matrix exponential be put somewhere else? (I mean, if common/real.cpp is reserved for the most basic functionality on the high precision floats) In dolfin/la together with the SORSolver? > > 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 > _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
