|
- #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( " %6" LAPACK_IFMT, vec[j] );
- printf( "\n" );
- }
|