dgebd2 routine: 3rd parameter passed to dlarfdg is an array?

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

dgebd2 routine: 3rd parameter passed to dlarfdg is an array?

Postby marcoi » Thu May 04, 2017 7:19 am

Hi,
I'm trying to understand the inner functioning of dgebd2 routine
What I do not understand is the third parameter passed to dlarfdg routine : https://github.com/Reference-LAPACK/lapack/blob/master/SRC/dgebd2.f#L245
The third parameter, as far as I understood, should be an array: https://github.com/Reference-LAPACK/lapack/blob/master/SRC/dlarfg.f#L107
but....is A( MIN( I+1, M ), I ) an array?

Looking forward to your kind explanations and hints.
Marco
marcoi
 
Posts: 3
Joined: Thu May 04, 2017 7:10 am

Re: dgebd2 routine: 3rd parameter passed to dlarfdg is an ar

Postby Julien Langou » Mon May 15, 2017 6:15 am

Hi Marco.

(1) Before looking at DGEBD2, please look at DGEQR2 then. This will help you understand the behavior of DLARFG. This is easier to understand, I believe.

(2) A( MIN( I+1, M ), I ) points to the start of the array. Yes it represents an array. Its start. Please look at the interface of DLARFG as well.

If you are more comfortable with C++, please have a look at:
https://github.com/langou/latl/blob/master/include/geqr2.h

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


Return to Algorithm / Data

Who is online

Users browsing this forum: No registered users and 3 guests