On Thu, Nov 25, 2010 at 05:55:36PM +0000, Garth N. Wells wrote: > > > On 25/11/10 17:34, Anders Logg 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). > > What is the typical system size? > > Garth
We have used it for (dense) matrices of size up to 100 x 100 for computing quadrature weights etc for the cG(100) method. -- Anders _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : dolfin@lists.launchpad.net Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
