I am currently implementing a SDPA wrapper. I use the BLAS/LAPACK libraries used by Julia which can be obtained by
Libdl.dlpath(BLAS.libblas). The value is typically
/usr/lib/julia/libopenblas64_.so. As Julia adds a
to the BLAS calls, I need to redefine BLAS functions when compiling the library using, e.g.
-Ddgemm_=dgemm_64_ when compiling SDPA.
However, for some fortran code, redefinining
dgemm_64 gives me
Error: Function ‘alpuip21k’ at (1) has no IMPLICIT type dmumps_part6.F:3575:14: & root%SCHUR_POINTER(1),
I have never done any fortran but it seems to me that the fact that
gfortran is not able to infer this type may be due to the fact that it does not know any external function
dgemm_64. Therefore, it seems to me that this could be fixed by using an header with the suffixed functions. Is Julia shipped with the openblas header in addition to the openblas library ? Or would you have another idea to fix this ?