void swanBindToTexture1DEx( const char *modname, const char *texname, size_t width, void *ptr, size_t typesize, int flags ) {
CUresult err;
CUtexref cu_texref;
int mode, channels;
// get the module
CUmodule mod = swanGetModule( modname );
// get the texture
err = cuModuleGetTexRef(&cu_texref, mod, texname );
if( err != CUDA_SUCCESS) { error( "swanBindToTexture1D failed -- texture not found" ); }
// now bind
err = cuTexRefSetAddress( NULL, cu_texref, PTR_TO_CUDEVPTR(ptr), width * typesize );
if( err != CUDA_SUCCESS) {
printf("EEERRR = %d\n", err );
error( "swanBindToTexture1D failed -- bind failed" );
}
// does not work for linear memory
/*
if( (flags & TEXTURE_INTERPOLATE) == TEXTURE_INTERPOLATE ) {
err = cuTexRefSetFilterMode( cu_texref, CU_TR_FILTER_MODE_LINEAR );
}
else {
err = cuTexRefSetFilterMode( cu_texref, CU_TR_FILTER_MODE_POINT );
}
if( err != CUDA_SUCCESS) { error( "swanBindToTexture1D failed -- setfiltermode failed" ); }
*/
mode = flags & TEXTURE_TYPE_MASK;
channels = typesize / sizeof(float);
switch( mode ) {
case TEXTURE_FLOAT:
err = cuTexRefSetFormat( cu_texref, CU_AD_FORMAT_FLOAT, channels );
break;
case TEXTURE_INT:
err = cuTexRefSetFormat( cu_texref, CU_AD_FORMAT_SIGNED_INT32, channels );
break;
case TEXTURE_UINT:
err = cuTexRefSetFormat( cu_texref, CU_AD_FORMAT_UNSIGNED_INT32, channels );
break;
default:
error( "swanBinToTexture1D failed -- invalid format" );
}
if( err != CUDA_SUCCESS) {
error( "swanBinToTexture1D failed -- setformat failed" );
}
return;
}
void CudaModule::setTexRef( const std::string& name,
CUdeviceptr ptr,
S64 size,
CUarray_format format,
int numComponents)
{
CUtexref texRef = getTexRef(name);
checkError("cuTexRefSetFormat", cuTexRefSetFormat(texRef, format, numComponents));
checkError("cuTexRefSetAddress", cuTexRefSetAddress(NULL, texRef, ptr, (U32)size));
}
void CuModule::SetFormat(CuModule::TexBind* texBind, CUarray_format format,
int numChannels) {
if((texBind->sampler.fmt != format) ||
(texBind->sampler.numPackedComponents != numChannels)) {
cuTexRefSetFormat(texBind->texRef, format, numChannels);
texBind->sampler.fmt = format;
texBind->sampler.numPackedComponents = numChannels;
}
}
sparseStatus_t sparseEngine_d::LoadKernel(sparsePrec_t prec,
sparseEngine_d::Kernel** ppKernel) {
// First attempt to load the finalize module if it is not yet loaded.
CUresult result = CUDA_SUCCESS;
// Check if the requested kernel is available, and if not, load it.
int p = (int)prec;
if(!multiply[p].get()) {
std::auto_ptr<Kernel> k(new Kernel);
std::string filename = kernelPath + "spmxv_" + PrecNames[p] +
".cubin";
result = context->LoadModuleFilename(filename, &k->module);
if(CUDA_SUCCESS != result) return SPARSE_STATUS_KERNEL_NOT_FOUND;
// Load the five SpMxV kernels for different valuesPerThread counts.
for(int i(0); i < NumVT; ++i) {
std::ostringstream oss;
oss<< "SpMxV_"<< ValuesPerThread[i];
result = k->module->GetFunction(oss.str(),
make_int3(BlockSize, 1,1), &k->func[i]);
if(CUDA_SUCCESS != result) return SPARSE_STATUS_KERNEL_ERROR;
}
// Load the finalize function.
result = k->module->GetFunction("Finalize", make_int3(BlockSize, 1, 1),
&k->finalize);
if(CUDA_SUCCESS != result) return SPARSE_STATUS_KERNEL_ERROR;
// Cache the texture reference
result = cuModuleGetTexRef(&k->xVec_texture, k->module->Handle(),
"xVec_texture");
if(CUDA_SUCCESS != result) return SPARSE_STATUS_KERNEL_ERROR;
result = cuTexRefSetFlags(k->xVec_texture, CU_TRSF_READ_AS_INTEGER);
if(CUDA_SUCCESS != result) return SPARSE_STATUS_KERNEL_ERROR;
result = cuTexRefSetFormat(k->xVec_texture, PrecTerms[p].vecFormat,
PrecTerms[p].vecChannels);
if(CUDA_SUCCESS != result) return SPARSE_STATUS_KERNEL_ERROR;
multiply[p] = k;
}
*ppKernel = multiply[p].get();
return SPARSE_STATUS_SUCCESS;
}
void swanMakeTexture1DEx( const char *modname, const char *texname, size_t width, void *ptr, size_t typesize, int flags ) {
int err;
// get the texture
CUtexref cu_texref;
int mode, channels;
CUarray array;
CUDA_MEMCPY2D copyParam;
CUDA_ARRAY_DESCRIPTOR p;
// get the module
CUmodule mod = swanGetModule( modname );
err = cuModuleGetTexRef(&cu_texref, mod, texname );
if( err != CUDA_SUCCESS) { error( "swanMakeTexture1D failed -- texture not found" ); }
p.Width = width;
p.Height= 1;
mode = flags & TEXTURE_TYPE_MASK;
channels = typesize / sizeof(float);
switch( mode ) {
case TEXTURE_FLOAT:
p.Format = CU_AD_FORMAT_FLOAT;
p.NumChannels = channels;
break;
case TEXTURE_INT:
p.Format = CU_AD_FORMAT_SIGNED_INT32;
p.NumChannels = channels;
break;
case TEXTURE_UINT:
p.Format = CU_AD_FORMAT_UNSIGNED_INT32;
p.NumChannels = channels;
break;
default:
error( "swanMakeTexture1D failed -- invalid format" );
}
err = cuArrayCreate( &array , &p);
if( err != CUDA_SUCCESS) { error( "swanMakeTexture1D failed -- array create failed" ); }
memset(©Param, 0, sizeof(copyParam));
copyParam.dstMemoryType = CU_MEMORYTYPE_ARRAY;
copyParam.dstArray = array;
copyParam.srcMemoryType = CU_MEMORYTYPE_HOST;
copyParam.srcHost = ptr;
copyParam.srcPitch = width * sizeof(float);
copyParam.WidthInBytes = copyParam.srcPitch;
copyParam.Height = 1;
// err = cuMemcpy2D(©Param);
err = cuMemcpyHtoA( array, 0, ptr, typesize * width );
if( err != CUDA_SUCCESS) { error( "swanMakeTexture1D failed -- memcpy failed" ); }
err = cuTexRefSetArray ( cu_texref, array, CU_TRSA_OVERRIDE_FORMAT );
if( err != CUDA_SUCCESS) { error( "swanMakeTexture1D failed -- setarray failed" ); }
if( (flags & TEXTURE_INTERPOLATE) == TEXTURE_INTERPOLATE ) {
err = cuTexRefSetFilterMode( cu_texref, CU_TR_FILTER_MODE_LINEAR );
}
else {
err = cuTexRefSetFilterMode( cu_texref, CU_TR_FILTER_MODE_POINT );
}
if( err != CUDA_SUCCESS) { error( "swanBindToTexture1D failed -- setfiltermode failed" ); }
if( (flags & TEXTURE_NORMALISE ) == TEXTURE_NORMALISE ) {
err = cuTexRefSetFlags(cu_texref, CU_TRSF_NORMALIZED_COORDINATES);
err |= cuTexRefSetAddressMode(cu_texref, 0, CU_TR_ADDRESS_MODE_CLAMP);
err |= cuTexRefSetAddressMode(cu_texref, 1, CU_TR_ADDRESS_MODE_CLAMP);
if( err != CUDA_SUCCESS) { error( "swanBindToTexture1D failed -- setflags 1 failed" ); }
}
err = cuTexRefSetFormat( cu_texref, CU_AD_FORMAT_FLOAT, channels );
if( err != CUDA_SUCCESS) { error( "swanBindToTexture1D failed -- setformat failed" ); }
//printf("TEX BIND DONE\n");
}
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) {
CUresult cudastatus = CUDA_SUCCESS;
if (nrhs != 2)
mexErrMsgTxt("Wrong number of arguments");
if (init == 0) {
// Initialize function
//mexLock();
// load GPUmat
gm = gmGetGPUmat();
// load module
CUmodule *drvmod = gmGetModule("examples_texture");
// load float GPU function
CUresult status = cuModuleGetFunction(&drvfunf, *drvmod, "LININTERF");
if (CUDA_SUCCESS != status) {
mexErrMsgTxt("Unable to load user function.");
}
// load double GPU function
status = cuModuleGetFunction(&drvfund, *drvmod, "LININTERD");
if (CUDA_SUCCESS != status) {
mexErrMsgTxt("Unable to load user function.");
}
// load textures defined in module
status = cuModuleGetTexRef(&texf, *drvmod, "texref_f1_a");
if (CUDA_SUCCESS != status) {
mexErrMsgTxt("Unable to load texture.");
}
status = cuModuleGetTexRef(&texd, *drvmod, "texref_d1_a");
if (CUDA_SUCCESS != status) {
mexErrMsgTxt("Unable to load texture.");
}
// no complex function support
init = 1;
}
// mex parameters are:
// 1. IN1. Input array
// 2. IN2. Input indexes array
//IN1 is the input GPU array
GPUtype IN1 = gm->gputype.getGPUtype(prhs[0]);
//IN2 is the input GPU array
GPUtype IN2 = gm->gputype.getGPUtype(prhs[1]);
//OUT is the output GPU array (result)
// Create of the same size of IN1
gpuTYPE_t in1_t = gm->gputype.getType(IN1);
int in1_d = gm->gputype.getNdims(IN1);
const int * in1_s = gm->gputype.getSize(IN1);
int in1_n = gm->gputype.getNumel(IN1);
int in1_b = gm->gputype.getDataSize(IN1);
gpuTYPE_t in2_t = gm->gputype.getType(IN2);
int in2_d = gm->gputype.getNdims(IN2);
const int * in2_s = gm->gputype.getSize(IN2);
int in2_n = gm->gputype.getNumel(IN2);
if ((in1_t==gpuCFLOAT) || (in1_t==gpuCDOUBLE)) {
mexErrMsgTxt("Complex TYPE not supported");
}
if (in1_t != in2_t) {
mexErrMsgTxt("Input arguments must be of the same type");
}
if (in1_n != in2_n) {
mexErrMsgTxt("Input arguments must have the same number of elements");
}
//OUT is the output GPU array (result)
// Create of the same size of IN1
GPUtype OUT = gm->gputype.create(in1_t, in1_d, in1_s, NULL);
// I need the pointers to GPU memory
CUdeviceptr d_IN1 = (CUdeviceptr) (UINTPTR gm->gputype.getGPUptr(IN1));
CUdeviceptr d_IN2 = (CUdeviceptr) (UINTPTR gm->gputype.getGPUptr(IN2));
CUdeviceptr d_OUT = (CUdeviceptr) (UINTPTR gm->gputype.getGPUptr(OUT));
// The GPU kernel depends on the type of input/output
CUfunction drvfun;
CUtexref drvtex;
CUarray_format_enum drvtexformat;
int drvtexnum;
if (in1_t == gpuFLOAT) {
drvfun = drvfunf;
drvtex = texf;
drvtexformat = CU_AD_FORMAT_FLOAT;
drvtexnum = 1;
} else if (in1_t == gpuDOUBLE) {
drvfun = drvfund;
drvtex = texd;
drvtexformat = CU_AD_FORMAT_SIGNED_INT32;
drvtexnum = 2;
}
if (CUDA_SUCCESS != cuTexRefSetFormat(drvtex, drvtexformat, drvtexnum)) {
mexErrMsgTxt("Execution error (texture).");
}
if (CUDA_SUCCESS != cuTexRefSetAddress(NULL, drvtex, UINTPTR d_IN1, in1_n*in1_b)) {
mexErrMsgTxt("Execution error (texture).");
}
if (CUDA_SUCCESS != cuParamSetTexRef(drvfun, CU_PARAM_TR_DEFAULT, drvtex)) {
mexErrMsgTxt("Execution error (texture1).");
}
hostdrv_pars_t gpuprhs[2];
int gpunrhs = 2;
gpuprhs[0] = hostdrv_pars(&d_IN2,sizeof(d_IN2),__alignof(d_IN2));
gpuprhs[1] = hostdrv_pars(&d_OUT,sizeof(d_OUT),__alignof(d_OUT));
int N = in1_n;
hostGPUDRV(drvfun, N, gpunrhs, gpuprhs);
// return result
plhs[0] = gm->gputype.createMxArray(OUT);
}
