Changeset 1968 for trunk/LMDZ.TITAN/libf/chimtitan/solve_lapack.c

Timestamp:

Jul 7, 2018, 7:48:18 PM (6 years ago)

Author:

jvatant

Message:

Minor follow-on of r1965 - to be able to call Fortran routines with gcc compiler
you need trailing underscores. Also add prototype for safety.
--JVO

File:

• trunk/LMDZ.TITAN/libf/chimtitan/solve_lapack.c (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/chimtitan/solve_lapack.c

@@ -1 @@
-/* Call to LAPACKE matrix inversion dgesv */
+/* Call to LAPACK matrix inversion dgesv */
 /* ~60 times quicker than the Numerical Recipes standard solver */
 /* Author : Jan Vatant d'Ollone 2018 */
@@ -5 +5 @@
 #include "titan.h"
 
+/* DGESV prototype */
+extern void dgesv_(int* n, int* nrhs, double* a, int* lda, int* ipiv,
+               double* b, int* ldb, int* info);
+
+/*Main*/
 void solve_lapack( double ***aa, int m, int n0, int n1 )
 {
@@ -22 +27 @@
    }
 
-   dgesv(&n,&n,a,&n,ipiv,b,&n,&ok); /* Lapack subroutine */
+   dgesv_(&n,&n,a,&n,ipiv,b,&n,&ok); /* Lapack subroutine */
    
    if( ok > 0 ) {