TEST(TestGPUGpuNUFFTFFT,GPUTest_FactorTwoTest)
{
//oversampling ratio
float osr = DEFAULT_OVERSAMPLING_RATIO;
//kernel width
int kernel_width = 3;
//Image
int im_width = 16;
//Data
int data_entries = 5;
DType2* data = (DType2*) calloc(data_entries,sizeof(DType2)); //2* re + im
int data_cnt = 0;
data[data_cnt].x = 0.5f;
data[data_cnt++].y = 0.5f;
data[data_cnt].x = 0.5f;
data[data_cnt++].y = 0.5f;
data[data_cnt].x = 1;
data[data_cnt++].y = 1;
data[data_cnt].x = 1;
data[data_cnt++].y = 1;
data[data_cnt].x = 1;
data[data_cnt++].y = 1;
//Coords
//Scaled between -0.5 and 0.5
//in triplets (x,y,z)
DType* coords = (DType*) calloc(3*data_entries,sizeof(DType));//3* x,y,z
int coord_cnt = 0;
//7.Sektor
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0.5f;
coords[coord_cnt++] = 0.3f;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0.3f;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
coords[coord_cnt++] = 0;
//sectors of data, count and start indices
int sector_width = 8;
gpuNUFFT::Array<DType> kSpaceData;
kSpaceData.data = coords;
kSpaceData.dim.length = data_entries;
gpuNUFFT::Dimensions imgDims;
imgDims.width = im_width;
imgDims.height = im_width;
imgDims.depth = im_width;
gpuNUFFT::GpuNUFFTOperatorFactory factory(false,true,true);
gpuNUFFT::GpuNUFFTOperator *gpuNUFFTOp = factory.createGpuNUFFTOperator(kSpaceData,kernel_width,sector_width,osr,imgDims);
gpuNUFFT::Array<DType2> dataArray;
dataArray.data = data;
dataArray.dim.length = data_entries;
gpuNUFFT::Array<CufftType> gdataArray;
gdataArray = gpuNUFFTOp->performGpuNUFFTAdj(dataArray,gpuNUFFT::FFT);
//Output Grid
CufftType* gdata = gdataArray.data;
int index = get3DC2lin(5,5,5,im_width);
if (DEBUG) printf("index to test %d\n",index);
if (DEBUG)
for (int j=0; j<im_width; j++)
{
for (int i=0; i<im_width; i++)
printf("%.4f ",gdata[get3DC2lin(i,j,5,im_width)].x);
printf("\n");
}
EXPECT_NEAR(gdata[get3DC2lin(8,8,5,16)].x,0.0514f,epsilon);
EXPECT_NEAR(gdata[get3DC2lin(2,2,5,16)].x,4.0694f,epsilon);
EXPECT_NEAR(gdata[get3DC2lin(4,3,5,16)].x,3.22198f,epsilon);
free(data);
free(coords);
free(gdata);
delete gpuNUFFTOp;
}
TEST(TestGPUGpuNUFFTFFT,GPUTest_Kernel5w64)
{
//oversampling ratio
float osr = DEFAULT_OVERSAMPLING_RATIO;
//kernel width
int kernel_width = 5;
//Image
int im_width = 64;
//Data
int data_entries = 1;
DType2* data = (DType2*) calloc(data_entries,sizeof(DType2)); //2* re + im
data[0].x = 1;
data[0].y = 0;
//Coords
//Scaled between -0.5 and 0.5
//in triplets (x,y,z)
DType* coords = (DType*) calloc(3*data_entries,sizeof(DType));//3* x,y,z
coords[0] = 0; //should result in 7,7,7 center
coords[1] = 0;
coords[2] = 0;
//sectors of data, count and start indices
int sector_width = 8;
//gpuNUFFT_gpu_adj(data,data_entries,1,coords,&gdata,grid_size,dims_g[1],kern,kernel_entries, kernel_width,sectors,sector_count,sector_centers,sector_width, im_width,osr,false,NULL,FFT);
gpuNUFFT::Array<DType> kSpaceData;
kSpaceData.data = coords;
kSpaceData.dim.length = data_entries;
gpuNUFFT::Dimensions imgDims;
imgDims.width = im_width;
imgDims.height = im_width;
imgDims.depth = im_width;
gpuNUFFT::GpuNUFFTOperatorFactory factory(false,true,true);
gpuNUFFT::GpuNUFFTOperator *gpuNUFFTOp = factory.createGpuNUFFTOperator(kSpaceData,kernel_width,sector_width,osr,imgDims);
gpuNUFFT::Array<DType2> dataArray;
dataArray.data = data;
dataArray.dim.length = data_entries;
gpuNUFFT::Array<CufftType> gdataArray;
gdataArray = gpuNUFFTOp->performGpuNUFFTAdj(dataArray,gpuNUFFT::FFT);
//Output Grid
CufftType* gdata = gdataArray.data;
if (DEBUG) printf("test: 57,33,32 = %.4f\n",gdata[get3DC2lin(57,33,32,im_width)].x);
EXPECT_NEAR(0.0001,gdata[get3DC2lin(57,33,32,im_width)].x,epsilon*10.0f);
/*for (int j=0; j<im_width; j++)
{
for (int i=0; i<im_width; i++)
if (abs(gdata[get3DC2lin(i,im_width-j,32,im_width)].x)>0.0f)
printf("(%d,%d,%d):%.4f ",i,j,32,gdata[get3DC2lin(i,im_width-j,32,im_width)].x);
printf("\n");
}*/
free(data);
free(coords);
free(gdata);
delete gpuNUFFTOp;
//free(sectors);
//free(sector_centers);
}
TEST(TestGPUGpuNUFFTDeapo,KernelCall1Sector2Coils)
{
int kernel_width = 3;
//oversampling ratio
float osr = 2.0f;//DEFAULT_OVERSAMPLING_RATIO;
//Image
int im_width = 32;
//Data
int data_entries = 1;
int n_coils = 2;
DType2* data = (DType2*) calloc(data_entries*n_coils,sizeof(DType2)); //2* re + im
data[0].x = 1;//Re
data[0].y = 0;//Im
data[1].x = 1;//Re
data[1].y = 0;//Im
//Coords, same per coil
//Scaled between -0.5 and 0.5
//in triplets (x,y,z)
DType* coords = (DType*) calloc(3*data_entries,sizeof(DType));//3* x,y,z
coords[0] = 0;
coords[1] = 0;
coords[2] = 0;
//sectors of data, count and start indices
int sector_width = 8;
gpuNUFFT::Array<DType> kSpaceData;
kSpaceData.data = coords;
kSpaceData.dim.length = data_entries;
gpuNUFFT::Dimensions imgDims;
imgDims.width = im_width;
imgDims.height = im_width;
imgDims.depth = im_width;
gpuNUFFT::GpuNUFFTOperatorFactory factory; gpuNUFFT::GpuNUFFTOperator *gpuNUFFTOp = factory.createGpuNUFFTOperator(kSpaceData,kernel_width,sector_width,osr,imgDims);
gpuNUFFT::Array<DType2> dataArray;
dataArray.data = data;
dataArray.dim.length = data_entries;
dataArray.dim.channels = n_coils;
gpuNUFFT::Array<CufftType> gdataArray;
gdataArray = gpuNUFFTOp->performGpuNUFFTAdj(dataArray,gpuNUFFT::DEAPODIZATION);
//Output Grid
CufftType* gdata = gdataArray.data;
if (DEBUG) printf("test %f \n",gdata[4].x);
int index = get3DC2lin(5,5,5,im_width);
if (DEBUG) printf("index to test %d\n",index);
int coil_offset = 1 * gpuNUFFTOp->getImageDims().count();
if (DEBUG) printf("grid_size: %d\n",gpuNUFFTOp->getImageDims().count());
/*for (int j=0; j<im_width; j++)
{
for (int i=0; i<im_width; i++)
printf("A: %.2f B: %.2f",gdata[get3DC2lin(i,j,16,im_width)].x,gdata[coil_offset+get3DC2lin(i,j,16,im_width)].x);
printf("\n");
}*/
EXPECT_NEAR(gdata[get3DC2lin(0,0,16,im_width)].x,0.016846,epsilon);
//EXPECT_NEAR(gdata[coil_offset + get3DC2lin(0,0,16,im_width)].x,0.00928f,epsilon);
EXPECT_NEAR(gdata[get3DC2lin(10,0,16,im_width)].x,gdata[coil_offset + get3DC2lin(10,0,16,im_width)].x,epsilon);
EXPECT_NEAR(gdata[get3DC2lin(0,5,16,im_width)].x,gdata[coil_offset + get3DC2lin(0,5,16,im_width)].x,epsilon);
EXPECT_NEAR(gdata[get3DC2lin(4,4,16,im_width)].x,gdata[coil_offset + get3DC2lin(4,4,16,im_width)].x,epsilon);
EXPECT_NEAR(gdata[get3DC2lin(15,14,16,im_width)].x,gdata[coil_offset + get3DC2lin(15,14,16,im_width)].x,epsilon);
free(data);
free(coords);
free(gdata);
delete gpuNUFFTOp;
EXPECT_EQ(1, 1);
}
TEST(TestGPUGpuNUFFTFFT,KernelCall1Sector)
{
int kernel_width = 3;
//Image
int im_width = 10;
//Data
int data_entries = 1;
DType2* data = (DType2*) calloc(data_entries,sizeof(DType2)); //2* re + im
data[0].x = 1;
data[0].y = 1;
//Coords
//Scaled between -0.5 and 0.5
//in triplets (x,y,z)
DType* coords = (DType*) calloc(3*data_entries,sizeof(DType));//3* x,y,z
coords[0] = 0; //should result in 7,7,7 center
coords[1] = 0;
coords[2] = 0;
//oversampling ratio
float osr = DEFAULT_OVERSAMPLING_RATIO;
//sectors of data, count and start indices
int sector_width = 5;
gpuNUFFT::Array<DType> kSpaceData;
kSpaceData.data = coords;
kSpaceData.dim.length = data_entries;
gpuNUFFT::Dimensions imgDims;
imgDims.width = im_width;
imgDims.height = im_width;
imgDims.depth = im_width;
gpuNUFFT::GpuNUFFTOperatorFactory factory(false,true,true);
gpuNUFFT::GpuNUFFTOperator *gpuNUFFTOp = factory.createGpuNUFFTOperator(kSpaceData,kernel_width,sector_width,osr,imgDims);
gpuNUFFT::Array<DType2> dataArray;
dataArray.data = data;
dataArray.dim.length = data_entries;
gpuNUFFT::Array<CufftType> gdataArray;
gdataArray = gpuNUFFTOp->performGpuNUFFTAdj(dataArray,gpuNUFFT::FFT);
//Output Grid
CufftType* gdata = gdataArray.data;
// gpuNUFFT_gpu_adj(data,data_entries,1,coords,&gdata,grid_size,dims_g[1],kern,kernel_entries, kernel_width,sectors,sector_count,sector_centers,sector_width, im_width,osr,false,NULL,FFT);
if (DEBUG) printf("test %f \n",gdata[4].x);
int index = get3DC2lin(5,5,5,im_width);
if (DEBUG) printf("index to test %d\n",index);
EXPECT_EQ(index,555);
EXPECT_NEAR(0.00097f,gdata[index].x,epsilon);
EXPECT_NEAR(0.0027f,gdata[get3DC2lin(5,4,5,im_width)].x,epsilon*10.0f);
EXPECT_NEAR(0.0027f,gdata[get3DC2lin(4,5,5,im_width)].x,epsilon*10.0f);
EXPECT_NEAR(0.3262f,gdata[get3DC2lin(0,1,5,im_width)].x,epsilon*10.0f);
/*for (int j=0; j<im_width; j++)
{
for (int i=0; i<im_width; i++)
printf("%.4f ",gdata[get3DC2lin(i,j,5,im_width)].x);
printf("\n");
}*/
free(data);
free(coords);
free(gdata);
EXPECT_EQ(1, 1);
}
TEST(TestGPUGpuNUFFTDeapo,KernelCall1Sector)
{
int kernel_width = 3;
//oversampling ratio
float osr = 1.75;//DEFAULT_OVERSAMPLING_RATIO;
//Image
int im_width = 64;
//Data
int data_entries = 1;
DType2* data = (DType2*) calloc(data_entries,sizeof(DType2)); //2* re + im
data[0].x = 1;//Re
data[0].y = 0;//Im
//Coords
//Scaled between -0.5 and 0.5
//in triplets (x,y,z)
DType* coords = (DType*) calloc(3*data_entries,sizeof(DType));//3* x,y,z
coords[0] = 0; //should result in 7,7,7 center
coords[1] = 0;
coords[2] = 0;
//sectors of data, count and start indices
int sector_width = 8;
gpuNUFFT::Array<DType> kSpaceData;
kSpaceData.data = coords;
kSpaceData.dim.length = data_entries;
gpuNUFFT::Dimensions imgDims;
imgDims.width = im_width;
imgDims.height = im_width;
imgDims.depth = im_width;
gpuNUFFT::GpuNUFFTOperatorFactory factory; gpuNUFFT::GpuNUFFTOperator *gpuNUFFTOp = factory.createGpuNUFFTOperator(kSpaceData,kernel_width,sector_width,osr,imgDims);
gpuNUFFT::Array<DType2> dataArray;
dataArray.data = data;
dataArray.dim.length = data_entries;
gpuNUFFT::Array<CufftType> gdataArray;
gdataArray = gpuNUFFTOp->performGpuNUFFTAdj(dataArray,gpuNUFFT::DEAPODIZATION);
//Output Grid
CufftType* gdata = gdataArray.data;
if (DEBUG) printf("test %f \n",gdata[4].x);
/*for (int y = 0; y < im_width; y++)
{
for (int x = 0; x < im_width; x++)
printf(" %.3f ", gdata[get3DC2lin(x,y,16,im_width)].x,epsilon);
printf(" \n ");
}*/
int index = get3DC2lin(5,5,5,im_width);
if (DEBUG) printf("index to test %d\n",index);
EXPECT_NEAR(0.00637,gdata[get3DC2lin(10,0,16,im_width)].x,epsilon);
free(data);
free(coords);
free(gdata);
delete gpuNUFFTOp;
}
