Very strange DGECON problem

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

Very strange DGECON problem

Postby GIBBOGLE » Sun Oct 12, 2014 10:26 pm

I am using Intel Fortran on Windows, and I have installed the 32-bit BLAS and LAPACK libraries. I am exploring the condition number of some matrices that I generate. A simple example matrix is 18x18. The '1' norm is calculated by DLANGE as 0.773E-01. Factorisation with DGETRF is successful, then DGECON returns condition number = Infinity. So it seems that this matrix is singular. BUT if I insert a statement to write out the array, write(*,'(6e12.3)') A, before the DLANGE statement, the condition number is determined to be 0.658E-03.

I have turned on the usual Fortran Runtime Error Checking, but it shows no problems.


Posts: 2
Joined: Fri Oct 10, 2014 6:11 pm

Return to User Discussion

Who is online

Users browsing this forum: No registered users and 6 guests