26 #ifndef __GPUImplHelper__
27 #define __GPUImplHelper__
30 #include "libhmsbeagle/config.h"
33 #include "libhmsbeagle/GPU/GPUImplDefs.h"
35 void checkHostMemory(
void* ptr);
40 template<
typename Real>
41 void transposeSquareMatrix(Real* mat,
43 for (
int i = 0; i < size - 1; i++) {
44 for (
int j = i + 1; j < size; j++) {
45 Real tmp = mat[i * size + j];
46 mat[i * size + j] = mat[j * size + i];
47 mat[j * size + i] = tmp;
52 template<
typename Real>
53 void printfVector(Real* ptr,
55 fprintf(stderr,
"[ %1.5e", ptr[0]);
57 for (i = 1; i < length; i++)
58 fprintf(stderr,
" %1.5e", ptr[i]);
59 fprintf(stderr,
" ]\n");
62 void printfInt(
int* ptr,
65 #endif // __GPUImplHelper__