MAGMA Forum

Code: Select all

#include<stdlib.h>
#include<stdio.h>
#include<cmath>
#include "magma.h"

int main ( int argc , char ** argv ){
magma_init (); // initialize Magma
magma_int_t m = 1024; // length of a
float *a; // a - m- vector on the host
float *d_a; // d_a - m- vector a on the device
magma_err_t err;
// allocate the vector on the host
err = magma_smalloc_cpu ( &a , m ); // host memory for a
// allocate the vector on the device
err = magma_smalloc ( &d_a , m ); // device memory for a
// a={ sin (0) , sin (1) ,... , sin (m -1)}
for(int j=0;j<m;j++) a[j]=sin((float)j);
// copy data from host to device
magma_ssetvector ( m, a, 1 , d_a , 1 ); // copy a -> d_a
// find the smallest index of the element of d_a with maximum
// absolute value
int i = magma_isamax( m, d_a, 1 );
printf ("max |a[i]|: %f\n", fabs(a[i -1]));
printf (" fortran index : %d\n",i);
free (a); // free host memory
magma_free (d_a ); // free device memory
magma_finalize (); // finalize Magma
return 0;
}
// max |a[i]|: 0.999990
//
// fortran index : 700

Code: Select all

#include<stdlib.h>
#include<stdio.h>
#include<cmath>
#include "magma.h"

int main ( int argc , char ** argv ){
magma_init (); // initialize Magma
magma_int_t m = 1024; // length of a
float *a; // a - m- vector on the host
float *d_a; // d_a - m- vector a on the device 

//magma_err_t err;
// allocate the vector on the host
 magma_smalloc_cpu ( &a , m ); // host memory for a
// allocate the vector on the device
magma_smalloc ( &d_a , m ); // device memory for a
// a={ sin (0) , sin (1) ,... , sin (m -1)}
for(int j=0;j<m;j++) a[j]=sin((float)j);
// copy data from host to device
magma_ssetvector ( m, a, 1 , d_a , 1 ); // copy a -> d_a
// find the smallest index of the element of d_a with maximum
// absolute value
int i = magma_isamax( m, d_a, 1 );
printf ("max |a[i]|: %f\n", fabs(a[i -1]));
printf (" fortran index : %d\n",i);
free (a); // free host memory
magma_free (d_a ); // free device memory
magma_finalize (); // finalize Magma
return 0;
}
// max |a[i]|: 0.999990
//
// fortran index : 700

Code: Select all

mint magma-trunk/testing> ./testing_auxiliary 
======================================================================
test_num_gpus
$MAGMA_NUM_GPUS     ngpu     expect
===================================
not set                   1       1

$MAGMA_NUM_GPUS='' is an invalid number; using 1 GPU.
                          1       1

$MAGMA_NUM_GPUS='-1' is an invalid number; using 1 GPU.
-1                        1       1

$MAGMA_NUM_GPUS='2junk' is an invalid number; using 1 GPU.
2junk                     1       1

$MAGMA_NUM_GPUS='0' is an invalid number; using 1 GPU.
0                         1       1

1                         1       1

$MAGMA_NUM_GPUS='2' exceeds MagmaMaxGPUs=8 or available GPUs=1; using 1 GPUs.
2                         1       2

$MAGMA_NUM_GPUS='4' exceeds MagmaMaxGPUs=8 or available GPUs=1; using 1 GPUs.
4                         1       4

$MAGMA_NUM_GPUS='8' exceeds MagmaMaxGPUs=8 or available GPUs=1; using 1 GPUs.
8                         1       8

$MAGMA_NUM_GPUS='16' exceeds MagmaMaxGPUs=8 or available GPUs=1; using 1 GPUs.
16                        1      16

$MAGMA_NUM_GPUS='1000' exceeds MagmaMaxGPUs=8 or available GPUs=1; using 1 GPUs.
1000                      1       1 (maxgpu)

======================================================================
test_num_threads
get;      parallel_numthread= 8, lapack_numthread= 1
set( 4);  parallel_numthread= 8, lapack_numthread= 1 (expect  4)
WARNING: testing_auxiliary.cpp:124: assertion l_nthread == 4 failed
set( 1);  parallel_numthread= 8, lapack_numthread= 1 (expect  1)
set( 8);  parallel_numthread= 8, lapack_numthread= 1 (expect  8)
WARNING: testing_auxiliary.cpp:140: assertion l_nthread == 8 failed
set( 1);  parallel_numthread= 8, lapack_numthread= 1 (expect  1)

using ncores=8, omp_num_threads=4

$MAGMA_NUM_THREADS  nthread  expect
===================================
not set                   8       4 (omp_threads)

WARNING: testing_auxiliary.cpp:172: assertion p_nthread == omp_threads failed
$MAGMA_NUM_THREADS='' is an invalid number; using 1 thread.
                          1       1

$MAGMA_NUM_THREADS='-1' is an invalid number; using 1 thread.
-1                        1       1

$MAGMA_NUM_THREADS='2junk' is an invalid number; using 1 thread.
2junk                     1       1

$MAGMA_NUM_THREADS='0' is an invalid number; using 1 thread.
0                         1       1

1                         1       1

2                         2       2

4                         4       4

8                         8       8

16                        8      16

1000                      8       8 (ncores)

Error in test_xerbla, function-specific error (info = 113)
Error in test_xerbla, function-specific error (info = 112)
Error in test_xerbla, function-specific error (info = 100)
Error in test_xerbla, function-specific error (info = 3)
Error in test_xerbla, function-specific error (info = 2)
Error in test_xerbla, function-specific error (info = 1)
No error, why is test_xerbla calling xerbla? (info = 0)
On entry to test_xerbla, parameter 1 had an illegal value (info = -1)
On entry to test_xerbla, parameter 2 had an illegal value (info = -2)
On entry to test_xerbla, parameter 3 had an illegal value (info = -3)
Error in test_xerbla, unknown error (info = -100)
Error in test_xerbla, cannot allocate memory on CPU host (info = -112)
Error in test_xerbla, cannot allocate memory on GPU device (info = -113)

3 tests failed.
mint magma-trunk/testing>