GGEMSSolid.hh

Go to the documentation of this file.

#ifndef GUARD_GGEMS_GEOMETRIES_GGEMSSOLID_HH
#define GUARD_GGEMS_GEOMETRIES_GGEMSSOLID_HH
 
// ************************************************************************
// * This file is part of GGEMS.                                          *
// *                                                                      *
// * GGEMS is free software: you can redistribute it and/or modify        *
// * it under the terms of the GNU General Public License as published by *
// * the Free Software Foundation, either version 3 of the License, or    *
// * (at your option) any later version.                                  *
// *                                                                      *
// * GGEMS is distributed in the hope that it will be useful,             *
// * but WITHOUT ANY WARRANTY; without even the implied warranty of       *
// * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
// * GNU General Public License for more details.                         *
// *                                                                      *
// * You should have received a copy of the GNU General Public License    *
// * along with GGEMS.  If not, see <https://www.gnu.org/licenses/>.      *
// *                                                                      *
// ************************************************************************
 
#include <limits>
 
#include "GGEMS/io/GGEMSTextReader.hh"
#include "GGEMS/io/GGEMSHistogramMode.hh"
#include "GGEMS/tools/GGEMSRAMManager.hh"
#include "GGEMS/navigators/GGEMSNavigatorManager.hh"
 
class GGEMSGeometryTransformation;
 
class GGEMS_EXPORT GGEMSSolid
{
  public:
    GGEMSSolid(void);
 
    virtual ~GGEMSSolid(void);
 
    GGEMSSolid(GGEMSSolid const& solid) = delete;
 
    GGEMSSolid& operator=(GGEMSSolid const& solid) = delete;
 
    GGEMSSolid(GGEMSSolid const&& solid) = delete;
 
    GGEMSSolid& operator=(GGEMSSolid const&& solid) = delete;
 
    void EnableTracking(void);
 
    inline cl::Buffer* GetSolidData(GGsize const& thread_index) const {return solid_data_[thread_index];};
 
    inline cl::Buffer* GetLabelData(GGsize const& thread_index) const {return label_data_[thread_index];};
 
    void SetRotation(GGfloat3 const& rotation_xyz);
 
    void SetPosition(GGfloat3 const& position_xyz);
 
    template<typename T>
    void SetSolidID(GGsize const& solid_id, GGsize const& thread_index);
 
    virtual void UpdateTransformationMatrix(GGsize const& thread_index) = 0;
 
    virtual void Initialize(GGEMSMaterials* materials) = 0;
 
    virtual void EnableScatter(void) = 0;
 
    virtual void PrintInfos(void) const = 0;
 
    inline std::string GetRegisteredDataType(void) const {return data_reg_type_;};
 
    inline cl::Kernel* GetKernelParticleSolidDistance(GGsize const& thread_index) const {return kernel_particle_solid_distance_[thread_index];}
 
    inline cl::Kernel* GetKernelProjectToSolid(GGsize const& thread_index) const {return kernel_project_to_solid_[thread_index];}
 
    inline cl::Kernel* GetKernelTrackThroughSolid(GGsize const& thread_index) const {return kernel_track_through_solid_[thread_index];}
 
    inline cl::Buffer* GetHistogram(GGsize const& thread_index) const {return histogram_.histogram_[thread_index];}
 
    inline cl::Buffer* GetScatterHistogram(GGsize const& thread_index) const {return histogram_.scatter_[thread_index];}
 
  protected:
    virtual void InitializeKernel(void) = 0;
 
  protected:
    // Solid data infos and label (for voxelized solid)
    cl::Buffer** solid_data_; 
    cl::Buffer** label_data_; 
    GGsize number_activated_devices_; 
    // Geometric transformation applyied to solid
    GGEMSGeometryTransformation* geometry_transformation_; 
    // OpenCL kernels and options for kernel
    cl::Kernel** kernel_particle_solid_distance_; 
    cl::Kernel** kernel_project_to_solid_; 
    cl::Kernel** kernel_track_through_solid_; 
    std::string kernel_option_; 
    // Output data
    std::string data_reg_type_; 
    GGEMSHistogramMode histogram_; 
    bool is_scatter_; 
};
 
 
template<typename T>
void GGEMSSolid::SetSolidID(GGsize const& solid_id, GGsize const& thread_index)
{
  // Get the OpenCL manager
  GGEMSOpenCLManager& opencl_manager = GGEMSOpenCLManager::GetInstance();
 
  // Get pointer on OpenCL device
  T* solid_data_device = opencl_manager.GetDeviceBuffer<T>(solid_data_[thread_index], sizeof(T), thread_index);
 
  solid_data_device->solid_id_ = static_cast<GGint>(solid_id);
 
  // Release the pointer
  opencl_manager.ReleaseDeviceBuffer(solid_data_[thread_index], solid_data_device, thread_index);
}
 
#endif // End of GUARD_GGEMS_GEOMETRIES_GGEMSSOLID_HH