ManyBodySpaceBase_GPUtest_kernel.cuh File Reference

Functions
template<class Derived >
__global__ void	testManyBodySpaceBase_kernel (ManyBodySpaceBase< Derived > *const ptr)

Function Documentation

testManyBodySpaceBase_kernel()

template<class Derived >

__global__ void testManyBodySpaceBase_kernel

(

ManyBodySpaceBase< Derived > *const

ptr

)

                                                                                    {
    int const idx = blockIdx.x * blockDim.x + threadIdx.x;
    if(idx > 1) return;
 
    ManyBodySpaceBase<Derived>& hSpace = *ptr;
    debug_printf("\t hSpace.dim()=%d, hSpace.sysSize()=%d, hSpace.dimLoc()=%d\n", hSpace.dim(),
                 hSpace.sysSize(), hSpace.dimLoc());
    {
        debug_printf("\t Calling the copy constructor\n");
        ManyBodySpaceBase<Derived> hSpace2(hSpace);
        debug_printf("\t (END) Calling the copy constructor\n\n");
    }
 
    debug_printf("\t Computing translation equivalent classes\n");
    hSpace.computeTransEqClass();
    debug_printf("\t hSpace.transEqDim()=%d\n", hSpace.transEqDim());
    printf("(GPU)\t");
    for(int j = 0; j != hSpace.transEqDim(); ++j) { printf("%d ", hSpace.transEqClassRep(j)); }
    printf("\n");
    printf("(GPU)\t");
    for(int j = 0; j != hSpace.transEqDim(); ++j) { printf("%d ", hSpace.transPeriod()[j]); }
    printf("\n");
    debug_printf("\t (END) Computing translation equivalent classes\n\n");
    debug_printf("\t hSpace.transEqDim()=%d, hSpace.transEqClassRep().data()=%p, "
                 "hSpace.transPeriod().data()=%p\n",
                 hSpace.transEqDim(), hSpace.transEqClassRep().data(), hSpace.transPeriod().data());
    {
        debug_printf("\t Calling the copy constructor\n");
        ManyBodySpaceBase<Derived> hSpace2(hSpace);
        debug_printf("\t (END) Calling the copy constructor\n");
        debug_printf("\t hSpace2.transEqDim()=%d, hSpace2.transEqClassRep().data()=%p, "
                     "hSpace2.transPeriod().data()=%p\n",
                     hSpace2.transEqDim(), hSpace2.transEqClassRep().data(),
                     hSpace2.transPeriod().data());
        printf("(GPU)\t");
        for(int j = 0; j != hSpace2.transEqDim(); ++j) {
            printf("%d ", hSpace2.transEqClassRep(j));
        }
        printf("\n");
        printf("(GPU)\t");
        for(int j = 0; j != hSpace2.transEqDim(); ++j) { printf("%d ", hSpace2.transPeriod()[j]); }
        printf("\n");
 
        debug_printf("\t Calling the move constructor\n");
        ManyBodySpaceBase<Derived> hSpace3(std::move(hSpace2));
        debug_printf("\t (END) Calling the move constructor\n");
        debug_printf("\t hSpace3.transEqDim()=%d, hSpace3.transEqClassRep().data()=%p, "
                     "hSpace3.transPeriod().data()=%p\n",
                     hSpace3.transEqDim(), hSpace3.transEqClassRep().data(),
                     hSpace3.transPeriod().data());
        printf("(GPU)\t");
        for(int j = 0; j != hSpace3.transEqDim(); ++j) {
            printf("%d ", hSpace3.transEqClassRep(j));
        }
        printf("\n");
        printf("(GPU)\t");
        for(int j = 0; j != hSpace3.transEqDim(); ++j) { printf("%d ", hSpace3.transPeriod()[j]); }
        printf("\n");
    }
 
    debug_printf("\t Testing translation operation\n");
    for(int j = 0; j != hSpace.transEqDim(); ++j) {
        for(int trans = 0; trans != hSpace.transPeriod(j); ++trans) {
            printf("(GPU)\t");
            for(int pos = 0; pos != hSpace.sysSize(); ++pos)
                printf("%d ", hSpace.locState(hSpace.transEqClassRep(j, trans), pos));
            printf("\n");
        }
        printf("\n");
    }
}