The LAPACK forum has moved to https://github.com/Reference-LAPACK/lapack/discussions.

In-place extra precision ZGESV solution ?

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

In-place extra precision ZGESV solution ?

Postby mocsaitamas » Thu Mar 11, 2010 5:15 am

Hi,

Generally I use ZGESV to solve complex linear algebraic equations, which are ill-conditioned however. So I need extra precision to get precise solution. For this I started to use the new ZGESVXX routine, which uses the iterative refinement based on extended precision arithmetics. This of course not an in-place solver anymore, as both the LU factorized system matrix and the original system matrix are kept in the memory. Because of this I am running out of memory. Could anybody give me a hint how to overcome this, except hardware upgrade ?
Thx
mocsaitamas
 
Posts: 4
Joined: Thu Mar 11, 2010 5:03 am

Re: In-place extra precision ZGESV solution ?

Postby Julien Langou » Thu Mar 11, 2010 11:04 am

Unfortunately, this is an intrinsic limitation of all iterative refinement methods.
You need the original problem (to refine) and the factorization (to solve).
Best wishes,
Julien.
Julien Langou
 
Posts: 835
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA

Re: In-place extra precision ZGESV solution ?

Postby mocsaitamas » Thu Mar 11, 2010 1:15 pm

Are there extra precision in place solvers currently available ? What is the current status of XBLAS for example ?
mocsaitamas
 
Posts: 4
Joined: Thu Mar 11, 2010 5:03 am


Return to User Discussion

Who is online

Users browsing this forum: No registered users and 5 guests