LAPACK Archives

[Lapack] WORK matrix in DGEQRT routine

Dear Sir/Madam,

Could you please tell me the WORK matrix elements usage? What will store in
WORK?

I initialize WORK matrix elements In this way:

double WORK[NB*N];
for(int a=0;a<NB*N; ++a)
WORK[a]=0;
and then run dgeqrt routine:

dgeqrt_(&M, &N, &NB, H, &LDA, T, &LDT, WORK, &INFO );

but the elements of WORK did not change. I pasted the result of running
dgeqrt below:

AMatrix
35 13 65 14
51 98 88 23
37 58 77 97
60 57 94 8

INFO 0
AMatrix
-76 0.12 0.58 0.13
-1.2e+02 -79 -0.073 0.01
-1.1e+02 -43 78 -0.96
-1.2e+02 -37 -15 -5.7

T
1.5 0 0 0
-0.22 2 0 0
-0.72 0.17 1 0
-0 -0 0 0

WORK
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

Best regards,
Nima
-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
<http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20150409/8efb6638/attachment.html>

<Prev in Thread] Current Thread [Next in Thread>
  • [Lapack] WORK matrix in DGEQRT routine, Nima Sahraneshin-Samani <=


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or