OperatorSpace_GPUtest.cu File Reference

Functions
template<typename Scalar >
__global__ void	createOperatorSpace_kernel (OperatorSpace< Scalar > **ptrptr, int dim)
int	main (int argc, char **argv)

Function Documentation

createOperatorSpace_kernel()

template<typename Scalar >

__global__ void createOperatorSpace_kernel	(	OperatorSpace< Scalar > **	ptrptr,
		int	dim
	)

                                                                                    {
    int const idx = blockIdx.x * blockDim.x + threadIdx.x;
    if(idx > 1) return;
    *ptrptr = new OperatorSpace<Scalar>(HilbertSpace<int>(dim));
    debug_printf("%s:\n\t dim=%d\n\t Created an OperatorSpace object at ptr=%p\n ",
                 __PRETTY_FUNCTION__, dim, *ptrptr);
}

main()

int main	(	int	argc,
		char **	argv
	)

                                {
    if(argc != 2) {
        std::cerr << "Usage: 1.(This) 2.(dim)\n"
                  << "       argc=" << argc << std::endl;
        std::exit(EXIT_FAILURE);
    }
    int const dim = std::atoi(argv[1]);
 
    OperatorSpace< Complex_t<double> >** ptrTo_dPtr;
    cuCHECK(cudaMalloc((void***)&ptrTo_dPtr, sizeof(OperatorSpace< Complex_t<double> >*)));
    createOperatorSpace_kernel<<<1, 1>>>(ptrTo_dPtr, dim);
    cuCHECK(cudaGetLastError());
    cuCHECK(cudaDeviceSynchronize());
 
    OperatorSpace< Complex_t<double>>* dPtr;
    cuCHECK(cudaMemcpy(&dPtr, ptrTo_dPtr, sizeof(OperatorSpace< Complex_t<double> >*),
                       cudaMemcpyDeviceToHost));
    std::cout << "dPtr=" << dPtr << std::endl;
    cuCHECK(cudaFree(ptrTo_dPtr));
 
    testOperatorSpaceBase_kernel<<<1, 1>>>(dPtr);
    cuCHECK(cudaGetLastError());
    cuCHECK(cudaDeviceSynchronize());
    std::cout << "Passed \"testOperatorSpaceBase_kernel\"\n" << std::endl;
 
    return EXIT_SUCCESS;
}