topic dgesvd under gcc in Intel® oneAPI Math Kernel Library

dgesvd under gcc

Patrick38 — Thu, 03 Oct 2024 10:29:25 GMT

Hello,

I downloaded the following file

and I compiled as follow

gcc example_svd.c -m64 -Wl,--start-group ${MKLROOT}/lib/libmkl_intel_ilp64.a ${MKLROOT}/lib/libmkl_gnu_thread.a ${MKLROOT}/lib/libmkl_core.a -Wl,--end-group -lgomp -lpthread -lm -ldl

running, I get

DGESVD Example Program Results

Intel oneMKL ERROR: Parameter 6 was incorrect on entry to DGESVD.

Could you clarify the point ? It seems to be an issue fortran/c

---------------------------------------------------------

#include <stdlib.h>
#include <stdio.h>

/* DGESVD prototype */
extern void dgesvd( char* jobu, char* jobvt, int* m, int* n, double* a,
int* lda, double* s, double* u, int* ldu, double* vt, int* ldvt,
double* work, int* lwork, int* info );
/* Auxiliary routines prototypes */
extern void print_matrix( char* desc, int m, int n, double* a, int lda );

/* Parameters */
#define M 6
#define N 5
#define LDA M
#define LDU M
#define LDVT N

/* Main program */
int main() {
/* Locals */
int m = M, n = N, lda = LDA, ldu = LDU, ldvt = LDVT, info, lwork;
double wkopt;
double* work;
/* Local arrays */
double s[N], u[LDU*M], vt[LDVT*N];
double a[LDA*N] = {
8.79, 6.11, -9.15, 9.57, -3.49, 9.84,
9.93, 6.91, -7.93, 1.64, 4.02, 0.15,
9.83, 5.04, 4.86, 8.83, 9.80, -8.99,
5.45, -0.27, 4.85, 0.74, 10.00, -6.02,
3.16, 7.98, 3.01, 5.80, 4.27, -5.31
};
/* Executable statements */
printf( " DGESVD Example Program Results\n" );
/* Query and allocate the optimal workspace */
lwork = -1;
dgesvd( "All", "All", &m, &n, a, &lda, s, u, &ldu, vt, &ldvt, &wkopt, &lwork,
&info );
lwork = (int)wkopt;
work = (double*)malloc( lwork*sizeof(double) );
/* Compute SVD */
dgesvd( "All", "All", &m, &n, a, &lda, s, u, &ldu, vt, &ldvt, work, &lwork,
&info );
/* Check for convergence */
if( info > 0 ) {
printf( "The algorithm computing SVD failed to converge.\n" );
exit( 1 );
}
/* Print singular values */
print_matrix( "Singular values", 1, n, s, 1 );
/* Print left singular vectors */
print_matrix( "Left singular vectors (stored columnwise)", m, n, u, ldu );
/* Print right singular vectors */
print_matrix( "Right singular vectors (stored rowwise)", n, n, vt, ldvt );
/* Free workspace */
free( (void*)work );
exit( 0 );
} /* End of DGESVD Example */

/* Auxiliary routine: printing a matrix */
void print_matrix( char* desc, int m, int n, double* a, int lda ) {
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+j*lda] );
printf( "\n" );
}
}

Re: dgesvd under gcc

Mark_L_Intel — Wed, 16 Oct 2024 02:02:47 GMT

Hello @Patrick38,

Using compiler options from the documentation page "Compiling and Building C Examples", I built example as

gcc -fopenmp -m64 -I${MKLROOT}/include dgesvd.c -Wl,--no-as-needed -L${MKLROOT}/lib/intel64 -lmkl_intel_lp64 -lmkl_core -lmkl_gnu_thread -lpthread -lm -ldl

and running example produces

./a.out
DGESVD Example Program Results

Singular values
27.47 22.64 8.56 5.99 2.01

Left singular vectors (stored columnwise)
-0.59 0.26 0.36 0.31 0.23
-0.40 0.24 -0.22 -0.75 -0.36
-0.03 -0.60 -0.45 0.23 -0.31
-0.43 0.24 -0.69 0.33 0.16
-0.47 -0.35 0.39 0.16 -0.52
0.29 0.58 -0.02 0.38 -0.65

Right singular vectors (stored rowwise)
-0.25 -0.40 -0.69 -0.37 -0.41
0.81 0.36 -0.25 -0.37 -0.10
-0.26 0.70 -0.22 0.39 -0.49
0.40 -0.45 0.25 0.43 -0.62
-0.22 0.14 0.59 -0.63 -0.44

this experiment was done with

>gcc --version
gcc (Ubuntu 13.2.0-23ubuntu4) 13.2.0

echo $MKLROOT
/opt/intel/oneapi/mkl/2024.2

Re: dgesvd under gcc

Mark_L_Intel — Wed, 30 Oct 2024 19:17:13 GMT

Hello @Patrick38,

Have you tried the the options outlined in my previous post?

Re: dgesvd under gcc

Mark_L_Intel — Thu, 07 Nov 2024 19:13:31 GMT

Hello @Patrick38,

With no response from you, the support will be limited to community support.