Hi Paul,
I have used ublas in the past with GNU multi precision library through the C++ class based interface and did not have any troubles, although I did not test it extensively. The only downside is the dramatic speed decrease that after all you cannot avoid with multi-prec. arithmetic.

On the other hand I don't find why uBlas should have trouble, since it just requires appropriate numeric semantics. If a type library is not implementing those properly you may be getting problems like you mention (although I believe it is difficult to be the case). 

Here is a long shot: there is the possibility that qq_real classes require some aligned memory access that afaik is not supported by ublas containers, but I think this may be something worth considering in the future having in mind the hype with GPUs and vectorization of algorithms.

I am not aware of any ready-implemented tests that would allow you to infiltrate through the problem you describe. I guess some exhaustive debuging on small matrices (like 2x2 or 3x3) could be the only way. If you could post an example I might be able to test it.

Best
Nasos


> From: paul.leopardi@iinet.net.au
> To: ublas@lists.boost.org
> Date: Wed, 24 Mar 2010 12:38:09 +1100
> Subject: [ublas] Does uBLAS work with dd_real and qd_real ?
>
> Hi all,
> Has anyone here worked with dd_real or qd_real (from qd-2.3.7
> http://www.cs.berkeley.edu/~yozo/ ) and uBLAS? My testing of these classes in
> the context of GluCat has resulted in peculiar errors in matrix results that
> look like they occur somewhere within uBLAS matrix multiplication (or possibly
> other operations) for matrix< dd_real, row_major > dense matrices.
>
> Is there a test suite I could use to narrow down the cause of the errors, or
> are there any known errors in the use of uBLAS with high and arbitrary
> precision floating point?
> Thanks, Paul
> _______________________________________________
> ublas mailing list
> ublas@lists.boost.org
> http://lists.boost.org/mailman/listinfo.cgi/ublas
> Sent to: nasos_i@hotmail.com


The New Busy is not the old busy. Search, chat and e-mail from your inbox. Get started.