Arbitrary Precision LAPACK

Open discussion regarding features, bugs, issues, vendors, etc.

Arbitrary Precision LAPACK

Postby kdallen » Tue Dec 06, 2011 4:45 pm

Is there any plan to release an arbitrary precision version (Fortran) of LAPACK?
kdallen
 
Posts: 39
Joined: Sun Jan 28, 2007 12:48 am

Re: Arbitrary Precision LAPACK

Postby yurii_ » Thu Dec 08, 2011 4:44 am

And what precision is necessary to you?
yurii_
 
Posts: 4
Joined: Tue Dec 06, 2011 10:40 am

Re: Arbitrary Precision LAPACK

Postby kdallen » Mon Dec 19, 2011 3:01 pm

I found the following dated from 2008 and was wondering if the project ever completed.

We are currently working on a multiple precision version of LAPACK,
and hope to release it within a year. You are welcome to take LAPACK
and modify it for your own purposes, under the guidelines of the
mBSD software license.

Jim Demmel
kdallen
 
Posts: 39
Joined: Sun Jan 28, 2007 12:48 am

Re: Arbitrary Precision LAPACK

Postby Julien Langou » Mon Dec 19, 2011 3:30 pm

Hello Ken,

I think the effort Jim was speaking about is now stopped unfortunately.

Two things:
1) we still have this in our mind and are seriously thinking about it, but there is nothing concrete yet.
2) A relevant software package is: http://sourceforge.net/projects/mplapack/

Julien.
Julien Langou
 
Posts: 727
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA

Re: Arbitrary Precision LAPACK

Postby kdallen » Mon Dec 19, 2011 6:59 pm

Julien,

Thanks for your reply. I found MPACK in searches but would prefer a 100% Fortran solution. I've actually modified the routines for DGELS to be compatible with

http://myweb.lmu.edu/dmsmith/fmlib.html

and seem to have gotten it to work, but it's quite a bit of effort.
kdallen
 
Posts: 39
Joined: Sun Jan 28, 2007 12:48 am

Re: Arbitrary Precision LAPACK

Postby yurii_ » Thu Dec 22, 2011 12:52 am

kdallen wrote:Is there any plan to release an arbitrary precision version (Fortran) of LAPACK?

It is possible to receive fourfold accuracy easily. The compiler of Intel supports it. But it is necessary to change code LAPACK a little. The program with fourfold accuracy will slowly work.
yurii_
 
Posts: 4
Joined: Tue Dec 06, 2011 10:40 am

Re: Arbitrary Precision LAPACK

Postby CyLith » Thu Jan 05, 2012 7:18 am

I have (had) been working on translating LAPACK into templated C++ code, templating on the number type. I actually converted the entire computation chain for unblocked ZGEEV (but not entirely templated), and you can get it here: https://github.com/victorliu/RNP. I would be incredibly interested in having the entirety of LAPACK translated this way since we could then plug in whatever number types we want, as well as being able to more easily interface it with existing C/C++ code.
CyLith
 
Posts: 35
Joined: Sun Feb 08, 2009 7:23 am
Location: Stanford, CA


Return to User Discussion

Who is online

Users browsing this forum: Bing [Bot] and 1 guest