Hi,
I modified one of the LAPACKE_examples to compute the matrix condition number
using LAPACKE_dlange and LAPACKE_dgecon. See changes to the example
highlighted in yellow, at end of the e-mail.
When I run the modified example, I get the following output:
Entry Matrix A
6.80 -6.05 -0.45 8.32 -9.67
-2.11 -3.30 2.58 2.71 -5.14
5.66 5.36 -2.70 4.35 -7.26
5.97 -4.44 0.27 -7.17 6.08
8.23 1.08 9.04 2.14 -6.87
Right Rand Side
4.02 -1.56 9.81
6.19 4.00 -4.09
-8.22 -8.67 -4.57
-7.57 1.75 -8.61
-3.03 2.86 8.99
LAPACKE_dgecon Example Program Results
Wrong parameter 1 in LAPACKE_dlange
** On entry to DGECON parameter number 5 had an illegal value
I've tried several things but I can't get past this error. Any ideas?
Thanks.
-Scott Lee
Modified example:
---------------------------
#include <stdlib.h>
#include <stdio.h>
#include "lapacke.h"
/* Auxiliary routines prototypes */
extern void print_matrix( char* desc, lapack_int m, lapack_int n, double* a,
lapack_int lda );
extern void print_int_vector( char* desc, lapack_int n, lapack_int* a );
/* Parameters */
#define N 5
#define NRHS 3
#define LDA N
#define LDB NRHS
/* Main program */
int main() {
/* Locals */
lapack_int n = N, nrhs = NRHS, lda = LDA, ldb = LDB, info;
/* Local arrays */
lapack_int ipiv[N];
double a[LDA*N] = {
6.80, -6.05, -0.45, 8.32, -9.67,
-2.11, -3.30, 2.58, 2.71, -5.14,
5.66, 5.36, -2.70, 4.35, -7.26,
5.97, -4.44, 0.27, -7.17, 6.08,
8.23, 1.08, 9.04, 2.14, -6.87
};
double b[LDB*N] = {
4.02, -1.56, 9.81,
6.19, 4.00, -4.09,
-8.22, -8.67, -4.57,
-7.57, 1.75, -8.61,
-3.03, 2.86, 8.99
};
double aNorm;
double rcond[1] = { 0.0 };
char ONE_NORM = '1';
int MATRIX_ORDER = n;
int NROWS = n;
int NCOLS = n;
int LEADING_DIMENSION_A = n;
double WORK[N];
/* Print Entry Matrix */
print_matrix( "Entry Matrix A", n, n, a, lda );
/* Print Right Rand Side */
print_matrix( "Right Rand Side", n, nrhs, b, ldb );
printf( "\n" );
/* Executable statements */
printf( "LAPACKE_dgecon Example Program Results\n" );
// DLANGE( NORM, M, N, A, LDA, WORK )
// LEFT OFF...When running this example, I get "Wrong parameter 1
in LAPACKE_dlange
aNorm = LAPACKE_dlange(MATRIX_ORDER, ONE_NORM, NROWS, NCOLS, a,
LEADING_DIMENSION_A); //, WORK);
info = LAPACKE_dgecon( LAPACK_ROW_MAJOR, ONE_NORM, MATRIX_ORDER, a,
LEADING_DIMENSION_A, aNorm, rcond); // aNorm should be 35.019999999999996
/* Check for the exact singularity */
if (info == 0)
{
printf("LAPACKE_dgecon completed SUCCESSFULLY...\n");
}
else if ( info < 0 )
{
printf( "Element %d of A had an illegal value\n", -info );
exit( 1 );
}
else
{
printf( "Unrecognized value of INFO = %d\n", info );
exit( 1 );
}
/* Print solution */
printf("LAPACKE_dlange / One-norm of A = %lf\n", aNorm);
printf("LAPACKE_dgecon / RCOND of A = %e", rcond[1] );
exit( 0 );
} /* End of LAPACKE_dgesv Example */
/* Auxiliary routine: printing a matrix */
void print_matrix( char* desc, lapack_int m, lapack_int n, double* a,
lapack_int lda ) {
lapack_int i, j;
printf( "\n %s\n", desc );
for( i = 0; i < m; i++ ) {
for( j = 0; j < n; j++ ) printf( " %6.2f", a[i*lda+j] );
printf( "\n" );
}
}
/* Auxiliary routine: printing a vector of integers */
void print_int_vector( char* desc, lapack_int n, lapack_int* a ) {
lapack_int j;
printf( "\n %s\n", desc );
for( j = 0; j < n; j++ ) printf( " %6i", a[j] );
printf( "\n" );
}
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20120308/01e03e28/attachment-0001.html
|