Question about dgecon routine

[Mr. President]

As the subject name suggets I have some problems running the dgecon routine to calculate the reciprocal condition number.
I make a LU-decomposition with the dgetrf-routine.
No problems there.

But in dgecon, the 5 parameter (standard ANORM) is an input estimate of the norm of the original matrix A. Dgecon then proceeds to calculate the condition number by calculating the norm of the inverse of A and multiplying this with your input.

However I can't find any routine to calculate the norm of A.
Does anybody have any suggestions.?
It would help me out a lot.

[ravi]

u can try with the follwing routine

ANORM = PDLANGE( 'I', N, N, A, 1, 1, DESCA, WORK )

ravi

[Julien Langou]

Hello,
just one comment regarding the first post:
DGECON does not compute the reciprocal condition number but does provide an estimate for it (in the same range of idea, during the computation, it does not compute the norm of the inverse of A, it only computes an estimation of it).
To get the 1-norm (Inf, or Fro), you can use the LAPACK routine DLANGE.
Julien

[Mr. President]

I just need to know if the matrix is well conditioned or not.
So an estimate will do just fine.

Thanks for the help