CUBLAS Fortran issues

Hi,

I’m try to use the CUBLAS Fortran wrapper (fortran.c). The code runs OK on my desktop (Intel Core2 Duo, GTX285, Ubuntu10.04 64bit). But the same code will return “CUBLAS_STATUS_MAPPING_ERROR” at the cublas_set_vector, if run on one node of a GPU cluster (6-core Intel Xeon, Tesla M2050, Linux x86-64).

I also test the C version of the code and it runs OK on both machines.

And it seems that there are some similar issues posted before but didn’t get solved:

http://forums.nvidia.com/index.php?showtopic=61060

Here is the Fortran code:

program main

	implicit none

	integer, parameter :: N=8

	integer, parameter :: sizeof_real=4

	real(kind=4), dimension(:),pointer :: h_p

	integer :: d_p

	integer stat

	integer cublas_init, cublas_alloc, cublas_set_vector, cublas_free, cublas_shutdown

	allocate(h_p(N))

	h_p = 0.0d0

	stat = cublas_init()

	if (stat .NE. 0) then

		write(*,*) "cublas_int failed"

		write(*,*) "stat=",stat

		stop

	endif

	stat = cublas_alloc(N, sizeof_real, d_p)

	if (stat .NE. 0) then

		write(*,*) " cublas_alloc failed"

		write(*,*) "stat=",stat

		stop

	endif

	write(*,*) "d_p =", d_p

	stat = cublas_set_vector (N, sizeof_real, h_p, 1, d_p, 1)

	if (stat .NE. 0) then

		write(*,*) "cublas_set_vector failed"

		write(*,*) "stat=",stat

		stop

	endif

	deallocate(h_p)

	stat = cublas_free (d_p)

	stat = cublas_shutdown()

	write(*,*) "test pass..."

	stop

end

Here is C version code:

#include <stdio.h>

#include <stdint.h>

#include <cublas.h>

int main(int argc, char *argv[])

{

	int N=8;

        float* h_p;

	float* d_p;

	int devPtr;

        cublasStatus stat;

h_p = (float *)malloc (N * sizeof(float));

        for (int i = 0; i < N; i++)  h_p[i] = 0;

stat = cublasInit();

        if (stat != 0) {

                printf ("cublasInit failed\n");

		printf ("stat=%d\n",stat);

                return 1;

        }

stat = cublasAlloc (N, sizeof(float), (void**)&d_p);

        if (stat != 0) {

                printf ("cublasAlloc failed");

		printf ("stat=%d\n",stat);

                return 1;

        }

	devPtr = (int)(uintptr_t)d_p;

	printf("devPtr=%d\n",devPtr);

stat = cublasSetVector (N, sizeof(float), h_p, 1, d_p, 1);

        if (stat != 0) {

                printf ("cublasSetVector failed");

		printf ("stat=%d\n",stat);

                return 1;

        }

cublasFree (d_p);

        cublasShutdown();

	printf("test pass...\n");

        return 0;

}

The output of the Fortran code on GTX 285:

d_p =     1048576

 test pass...

The output of the C code on GTX 285:

devPtr=1048576

test pass...

The output of the Fortran code on M2050:

d_p =           0

 cublas_set_vector failed

 stat=          11

The output of the C code on M2050:

devPtr=0

test pass...

For all the cases, the compilers and options are:

gcc -O3

gfortran -cpp -O3 -ffree-form -ffree-line-length-0 -fno-second-underscore

nvcc -O3 -arch sm_13 -use_fast_math

Any clue?

Thanks,

Kong