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

about work space size

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

about work space size

Postby neodreamer » Thu Dec 27, 2007 5:45 pm

many lapack routines require a work space,
for the optimal work space lwork >= N*NB.

I dont want to use lwork=-1 to get the optimal size.
If I want to allocate sufficient work space to the routines
such as zgeqrf or zunmqr, Is lwork=N*N large enough?
I thought NB<=N, is this correct?
.
neodreamer
 
Posts: 17
Joined: Thu Apr 05, 2007 3:52 pm

Postby Julien Langou » Wed Jan 02, 2008 12:34 pm

You really should use LWORK = -1.That's definetely the best way to go.

Now taking LWORK = N*N should indeed give you a workspace bigger than
the optimal workspace (for very small N that might not work though). This means
that the LAPACK code will be able to block the operations (and so run fast), this also
means that you will waste space.
Julien Langou
 
Posts: 835
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA


Return to User Discussion

Who is online

Users browsing this forum: No registered users and 5 guests