Gibt es eine Möglichkeit, die Größe der BLAS IDAMAX-Funktion zu bestimmen?

Question

Ich möchte die Größe der Integer-Funktion IDAMAX in jeder BLAS Fortran-Bibliothek bestimmen. Es kann entweder 4-Byte (i32lp64-Modell) oder 8-Byte (ilp64) sein.

Wenn Sie diese Größe kennen, können Sie bestimmen, welche Gesamt-Integer-Deklaration (Ganzzahl * 4 oder Ganzzahl * 8) in einer vorgefertigten BLAS-Bibliothek enthalten ist.

Das Problem ist, dass sizeof (* idamax_ (& n, dx, & INCX)) das Programm unter immer 4 zurückkehrt, wenn I 8 mit MKL-integer * 8 blas erwarten würde.

Irgendwelche Kommentare, bitte?

#include <stdio.h> 
void main() { 
extern * idamax_(); // external Fortran BLAS function IDAMAX 
int n=2; int incx=1; 
//long n=2; long incx=1; 
double dx[2]; dx[0]=1.0; dx[1]=2.0; 
printf("sizeof(n)=%i\n",sizeof(n)); 
printf("sizeof(*idamax_(&n, dx, &incx))=%i\n",sizeof(*idamax_(&n, dx,&incx))); // still returns four !!! 
//printf("sizeof(idamax_(&n, dx, &incx))=%i\n",sizeof(idamax_(&n, dx, &incx))); 
// idamax call sometimes crashes with wrong integer sizes - with MKL, but with GNU ibblas.a ! The same with Fortran 
idamax_(&n, dx, &incx); 
} 

Answer 1

Ihr Code wird jede Funktion in jeder BLAS Fortran Bibliothek nicht bestimmen. Es ist eine Funktion deklariert, indem Sie sich mit

extern * idamax_(); 

Bestimmen Sie nicht den Rückgabetyp dieser Funktion angegeben haben, dann ist es standardmäßig int, mit einer konstanten Größe von 4 auf den meisten Maschinen.

Wenn Sie den Rückgabetyp von MKL BLAS-Funktion bestimmen möchten, ist das einzige, was Sie tun müssen, überprüfen Sie die Header-Datei. In ${MKL_ROOT}/include/mkl_blas.h können Sie

MKL_INT IDAMAX(const MKL_INT *n, const double *x, const MKL_INT *incx); 

Der Rückgabetyp finden MKL_INT ist, und die Größe ist sizeof(MKL_INT), die auf der Kompilierung bestimmt wird. Sie können sich auf etwas Ähnliches wie dies in ${MKL_ROOT}/include/mkl_types.h

#ifdef MKL_ILP64 
    #ifndef MKL_INT 
    #define MKL_INT MKL_INT64 
    #endif 
#else 
    #ifndef MKL_INT 
    #define MKL_INT int 
    #endif 
#endif 

So ist die Größe von MKL_INT hängt nur, ob Sie MKL_ILP64 oder nicht an der Kompilierung definieren.