GGEMSDoseRecording.hh

Go to the documentation of this file.

#ifndef GUARD_GGEMS_NAVIGATORS_GGEMSDOSERECORDING_HH
#define GUARD_GGEMS_NAVIGATORS_GGEMSDOSERECORDING_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/>.      *
// *                                                                      *
// ************************************************************************
 
#ifndef __OPENCL_C_VERSION__
 
#include <memory>
 
typedef struct GGEMSDoseRecording_t
{
  cl::Buffer** edep_; 
  cl::Buffer** edep_squared_; 
  cl::Buffer** hit_; 
  cl::Buffer** photon_tracking_; 
  cl::Buffer** dose_; 
  cl::Buffer** uncertainty_dose_; 
} GGEMSDoseRecording; 
#endif
 
#ifdef __OPENCL_C_VERSION__
 
#include "GGEMS/navigators/GGEMSDoseParams.hh"
#include "GGEMS/geometries/GGEMSGeometryConstants.hh"
 
 
inline void dose_photon_tracking(global GGEMSDoseParams* dose_params, global GGint* photon_tracking, GGfloat3 const* position)
{
  // Check position of photon inside dosemap limits
  if (position->x < dose_params->border_min_xyz_.x + EPSILON6 || position->x > dose_params->border_max_xyz_.x - EPSILON6) return;
  if (position->y < dose_params->border_min_xyz_.y + EPSILON6 || position->y > dose_params->border_max_xyz_.y - EPSILON6) return;
  if (position->z < dose_params->border_min_xyz_.z + EPSILON6 || position->z > dose_params->border_max_xyz_.z - EPSILON6) return;
 
  // Get index in dose map
  GGint3 dosel_id = convert_int3((*position - dose_params->border_min_xyz_) * dose_params->inv_size_of_dosels_);
 
  GGint global_dosel_id = dosel_id.x + dosel_id.y * dose_params->number_of_dosels_.x + dosel_id.z * dose_params->number_of_dosels_.x * dose_params->number_of_dosels_.y;
 
  if (dosel_id.x < 0 || dosel_id.x >= dose_params->number_of_dosels_.x) return;
  if (dosel_id.y < 0 || dosel_id.y >= dose_params->number_of_dosels_.y) return;
  if (dosel_id.z < 0 || dosel_id.z >= dose_params->number_of_dosels_.z) return;
 
  atomic_add(&photon_tracking[global_dosel_id], 1);
}
 
 
inline void dose_record_standard(global GGEMSDoseParams* dose_params, global GGDosiType* edep_tracking, global GGDosiType* edep_squared_tracking, global GGint* hit_tracking, GGfloat edep, GGfloat3 const* position)
{
  // Check position of photon inside dosemap limits
  if (position->x < dose_params->border_min_xyz_.x + EPSILON6 || position->x > dose_params->border_max_xyz_.x - EPSILON6) return;
  if (position->y < dose_params->border_min_xyz_.y + EPSILON6 || position->y > dose_params->border_max_xyz_.y - EPSILON6) return;
  if (position->z < dose_params->border_min_xyz_.z + EPSILON6 || position->z > dose_params->border_max_xyz_.z - EPSILON6) return;
 
  // Get index in dose map
  GGint3 dosel_id = convert_int3((*position - dose_params->border_min_xyz_) * dose_params->inv_size_of_dosels_);
 
  GGint global_dosel_id = dosel_id.x + dosel_id.y * dose_params->number_of_dosels_.x + dosel_id.z * dose_params->number_of_dosels_.x * dose_params->number_of_dosels_.y;
 
  if (dosel_id.x < 0 || dosel_id.x >= dose_params->number_of_dosels_.x) return;
  if (dosel_id.y < 0 || dosel_id.y >= dose_params->number_of_dosels_.y) return;
  if (dosel_id.z < 0 || dosel_id.z >= dose_params->number_of_dosels_.z) return;
 
  if (hit_tracking) atomic_add(&hit_tracking[global_dosel_id], 1);
  #ifdef DOSIMETRY_DOUBLE_PRECISION
  AtomicAddDouble(&edep_tracking[global_dosel_id], (GGDosiType)edep);
  if (edep_squared_tracking) AtomicAddDouble(&edep_squared_tracking[global_dosel_id], (GGDosiType)edep*(GGDosiType)edep);
  #else
  AtomicAddFloat(&edep_tracking[global_dosel_id], (GGDosiType)edep);
  if (edep_squared_tracking) AtomicAddFloat(&edep_squared_tracking[global_dosel_id], (GGDosiType)edep*(GGDosiType)edep);
  #endif
}
 
#endif
 
#endif // End of GUARD_GGEMS_NAVIGATORS_GGEMSDOSERECORDING_HH