diff options
Diffstat (limited to 'biology/gcta/files/patch-main_mkl.cpp')
-rw-r--r-- | biology/gcta/files/patch-main_mkl.cpp | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/biology/gcta/files/patch-main_mkl.cpp b/biology/gcta/files/patch-main_mkl.cpp new file mode 100644 index 000000000000..a74cbecfa6fd --- /dev/null +++ b/biology/gcta/files/patch-main_mkl.cpp @@ -0,0 +1,24 @@ +--- main/mkl.cpp.orig 2023-11-20 14:09:49 UTC ++++ main/mkl.cpp +@@ -365,7 +365,9 @@ bool gcta::comput_inverse_logdet_LDLT_mkl(eigenMatrix + #if GCTA_CPU_x86 + dpotrf(&uplo, &int_n, Vi_mkl, &int_n, &info); + #else +- dpotrf_(&uplo, &int_n, Vi_mkl, &int_n, &info); ++ // OpenBLAS 0.3.25 requires a length argument. Is 1 the right value? ++ // Should all dpotrf_() calls be patched? ++ dpotrf_(&uplo, &int_n, Vi_mkl, &int_n, &info, 1); + #endif + //LOGGER << "Finished decompose" << endl; + //spotrf( &uplo, &n, Vi_mkl, &n, &info ); +@@ -386,7 +388,9 @@ bool gcta::comput_inverse_logdet_LDLT_mkl(eigenMatrix + #if GCTA_CPU_x86 + dpotri(&uplo, &int_n, Vi_mkl, &int_n, &info); + #else +- dpotri_(&uplo, &int_n, Vi_mkl, &int_n, &info); ++ // OpenBLAS 0.3.25 requires a length argument. Is 1 the right value? ++ // Should all dpotri_() calls be patched? ++ dpotri_(&uplo, &int_n, Vi_mkl, &int_n, &info, 1); + #endif + //LOGGER << "Inverse finished" << endl; + //spotri( &uplo, &n, Vi_mkl, &n, &info ); |