34 lines
1022 B
C
34 lines
1022 B
C
#include <lapacke.h>
|
|
#include <stdio.h>
|
|
|
|
/* Auxiliary routine: printing a matrix */
|
|
void print_matrix_rowmajor( char* desc, lapack_int m, lapack_int n, double* mat, lapack_int ldm ) {
|
|
lapack_int i, j;
|
|
printf( "\n %s\n", desc );
|
|
|
|
for( i = 0; i < m; i++ ) {
|
|
for( j = 0; j < n; j++ ) printf( " %6.2f", mat[i*ldm+j] );
|
|
printf( "\n" );
|
|
}
|
|
}
|
|
|
|
|
|
/* Auxiliary routine: printing a matrix */
|
|
void print_matrix_colmajor( char* desc, lapack_int m, lapack_int n, double* mat, lapack_int ldm ) {
|
|
lapack_int i, j;
|
|
printf( "\n %s\n", desc );
|
|
|
|
for( i = 0; i < m; i++ ) {
|
|
for( j = 0; j < n; j++ ) printf( " %6.2f", mat[i+j*ldm] );
|
|
printf( "\n" );
|
|
}
|
|
}
|
|
|
|
/* Auxiliary routine: printing a vector of integers */
|
|
void print_vector( char* desc, lapack_int n, lapack_int* vec ) {
|
|
lapack_int j;
|
|
printf( "\n %s\n", desc );
|
|
for( j = 0; j < n; j++ ) printf( " %6i", vec[j] );
|
|
printf( "\n" );
|
|
}
|