Commit 936f0db3 authored by RoteKekse's avatar RoteKekse
Browse files

Update blasLapackWrapper.cpp, code breaking changes to Lapack lib

parent 82dbb250
Pipeline #2163 failed with stages
in 4 minutes and 43 seconds
......@@ -223,7 +223,7 @@ namespace xerus {
double work = 0;
lapack_int lwork = -1;
lapack_int min = std::min(m,n);
dgesdd_( &job, &n, &m, nullptr, &n, nullptr, nullptr, &n, nullptr, &min, &work, &lwork, nullptr, &info );
dgesdd_( &job, &n, &m, nullptr, &n, nullptr, nullptr, &n, nullptr, &min, &work, &lwork, nullptr, &info, 1 );
REQUIRE(info == 0, "work array size query of dgesdd returned " << info);
return lapack_int(work);
}
......@@ -236,7 +236,7 @@ namespace xerus {
lapack_int min = std::min(m,n);
// if A = U*S*V^T, then A^T = V^T*S*U^T, so instead of transposing all input and output matrices we can simply exchange the order of U and Vt
dgesdd_( &job, &n, &m, a, &n, s, vt, &n, u, &min, work, &lwork, iwork, &info );
dgesdd_( &job, &n, &m, a, &n, s, vt, &n, u, &min, work, &lwork, iwork, &info, 1 );
REQUIRE(info == 0, "dgesdd failed with info " << info);
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment