my code and CMakeLists are here
#include <mkl.h>
#include <iostream>
#include <fstream>
#include <istream>
#include <complex>
using namespace std;
int main(){
int inputsize=41*41;
ifstream input("input.svd");
complex<double>* svndinput=(complex<double>*)malloc(inputsize*inputsize*sizeof(complex<double>));
input.read((char*)svndinput, inputsize*inputsize*sizeof(complex<double>));
input.close();
int m=inputsize;
int n=inputsize;
complex<double>* a=svndinput;
complex<double>* u=new complex<double>[inputsize*inputsize]();
complex<double>* v=new complex<double>[inputsize*inputsize]();
double* s=new double[inputsize]();
typedef lapack_complex_double lcd;
auto result=LAPACKE_zgesdd(LAPACK_ROW_MAJOR, 'A', m, n, (lcd*)a, n, (double*)s, (lcd*)u, m, (lcd*)v, n);
cout<<"result: "<<result<<endl;
ofstream saves("s.svd");
saves.write((char*)s, inputsize*sizeof(s[0]));
saves.flush();
saves.close();
ofstream savev("v.svd");
savev.write((char*)v, inputsize*inputsize*sizeof(v[0]));
savev.flush();
savev.close();
return 0;
}
CMakeLists
project (Demo2)
add_compile_options(-fPIC)
SET(INTEL_MKL_INSTALL_DIR "/opt/intel/oneapi/mkl/latest")
SET(MKL_LINK static)
find_package(MKL PATHS /opt/intel/oneapi/mkl/latest NO_DEFAULT_PATH)
INCLUDE_DIRECTORIES(/home/grp4data/anaconda3/envs/smitopc20221111/include)
INCLUDE_DIRECTORIES(${MKL_INCLUDE})
add_executable(testsvd main.cpp)
target_link_libraries(testsvd $<LINK_ONLY:MKL::MKL>)
