message: Modified spectralonesteptest to include a case with Cuda projectors
diff: @@ -316,7 +316,7 @@ int main(int, char** )
  CheckVectorImageQuality<MaterialVolumeType>(mechlemOneStep->GetOutput(), composeVols->GetOutput(), 0.08, 23, 2.0);
  std::cout << "\n\nTest PASSED! " << std::endl;
  - std::cout << "\n\n****** Case 3: Voxel-based Backprojector, 4 subsets ******" << std::endl;
  + std::cout << "\n\n****** Case 3: Voxel-based Backprojector, 4 subsets, with regularization ******" << std::endl;
  typedef rtk::ReorderProjectionsImageFilter<PhotonCountsType> ReorderProjectionsFilterType;
  ReorderProjectionsFilterType::Pointer reorder = ReorderProjectionsFilterType::New();
  @@ -330,13 +330,6 @@ int main(int, char** )
  mechlemOneStep->SetBackProjectionFilter( MechlemType::BP_VOXELBASED );
  mechlemOneStep->SetNumberOfSubsets(4);
  mechlemOneStep->SetNumberOfIterations( 5 );
  - TRY_AND_EXIT_ON_ITK_EXCEPTION( mechlemOneStep->Update() );
  -
  - CheckVectorImageQuality<MaterialVolumeType>(mechlemOneStep->GetOutput(), composeVols->GetOutput(), 0.08, 23, 2.0);
  - std::cout << "\n\nTest PASSED! " << std::endl;
  -
  - std::cout << "\n\n****** Case 4: Voxel-based Backprojector, 4 subsets, with regularization ******" << std::endl;
  -
  MaterialVolumeType::RegionType::SizeType radius;
  radius.Fill(1);
  MaterialVolumeType::PixelType weights;
  @@ -348,7 +341,15 @@ int main(int, char** )
  CheckVectorImageQuality<MaterialVolumeType>(mechlemOneStep->GetOutput(), composeVols->GetOutput(), 0.08, 23, 2.0);
  std::cout << "\n\nTest PASSED! " << std::endl;
  -#ifdef USE_CUDA
  +#ifdef RTK_USE_CUDA
  + std::cout << "\n\n****** Case 4: CUDA voxel-based Backprojector, 4 subsets, with regularization ******" << std::endl;
  +
  + mechlemOneStep->SetForwardProjectionFilter( MechlemType::FP_CUDARAYCAST );
  + mechlemOneStep->SetBackProjectionFilter( MechlemType::BP_CUDAVOXELBASED );
  + TRY_AND_EXIT_ON_ITK_EXCEPTION( mechlemOneStep->Update() );
  +
  + CheckVectorImageQuality<MaterialVolumeType>(mechlemOneStep->GetOutput(), composeVols->GetOutput(), 0.08, 23, 2.0);
  + std::cout << "\n\nTest PASSED! " << std::endl;
  #endif
  return EXIT_SUCCESS;
  