Class DeviceDnVecDouble

XM/include/Utils/memory.h:911–961

dense vector wrapper on device: double type

// dense vector wrapper on device: double type
class DeviceDnVecDouble {
    public:
        int gpu_id;
        int size;
        double* vals;
        cusparseDnVecDescr_t cusparse_descr;

        DeviceDnVecDouble(): gpu_id(0), size(0), vals(nullptr), cusparse_descr(NULL) {}
        DeviceDnVecDouble(const int gpu_id, const int size, bool as_byte = false):
            gpu_id(gpu_id), size(size), vals(nullptr), cusparse_descr(NULL) {
            this->allocate(this->gpu_id, this->size, as_byte);
        }

        inline void allocate(const int gpu_id, const int size, bool as_byte = false) {
            if (this->vals == nullptr) {
                this->gpu_id = gpu_id;
                this->size = size;
                CHECK_CUDA( cudaSetDevice(this->gpu_id) );
                // as_byte is used to allocate buffer size, which is usually given in terms of bytes
                if (as_byte) {
                    this->size = (size + sizeof(double) - 1) / sizeof(double);
                } else {
                    this->size = size;
                }
                CHECK_CUDA( cudaMalloc((void**) &this->vals, sizeof(double) * this->size) );
                CHECK_CUSPARSE( cusparseCreateDnVec(&this->cusparse_descr, this->size, this->vals, CUDA_R_64F) );
            }
            return;
        }
        inline double get_norm(const DeviceBlasHandle& cublas_H) {
            double norm;
            CHECK_CUDA( cudaSetDevice(this->gpu_id) );
            CHECK_CUBLAS( cublasDnrm2_v2(
                cublas_H.cublas_handle, this->size, this->vals, 1, &norm
            ) );
            return norm;
        }

        ~DeviceDnVecDouble() {
            CHECK_CUDA( cudaSetDevice(this->gpu_id) );
            if (this->vals != nullptr) {
                CHECK_CUDA( cudaFree(this->vals) );
                this->vals = nullptr;
            }
            if (this->cusparse_descr != NULL) {
                CHECK_CUSPARSE( cusparseDestroyDnVec(this->cusparse_descr) );
                this->cusparse_descr = NULL;
            }
            // std::cout << "DeviceDnVecDouble destructor called!" << std::endl;
        }
};

// dense vector wrapper on device: float type
class DeviceDnVecFloat {