ProjMatrixByDensel.inl

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2000- 2002, IRSL
    See STIR/LICENSE.txt for details
*/
 
#include "stir/recon_buildblock/SymmetryOperation.h"
#include "stir/Succeeded.h"
 
START_NAMESPACE_STIR
 
void
ProjMatrixByDensel::get_proj_matrix_elems_for_one_densel(ProjMatrixElemsForOneDensel& probabilities, const Densel& densel) const
{
  // set to empty
  probabilities.erase();
 
  // find basic densel
  Densel basic_densel = densel;
 
  unique_ptr<SymmetryOperation> symm_ptr = get_symmetries_ptr()->find_symmetry_operation_from_basic_densel(basic_densel);
 
  probabilities.set_densel(basic_densel);
 
  // check if in cache
  if (get_cached_proj_matrix_elems_for_one_densel(probabilities) == Succeeded::no)
 
    {
      // call 'calculate' just for the basic densel
      calculate_proj_matrix_elems_for_one_densel(probabilities);
#ifndef NDEBUG
      probabilities.check_state();
#endif
      cache_proj_matrix_elems_for_one_densel(probabilities);
    }
 
  // now transform to original densel
  symm_ptr->transform_proj_matrix_elems_for_one_densel(probabilities);
}

ProjMatrixByDensel::CacheKey
ProjMatrixByDensel::cache_key(const Densel& densel)
{
  assert(densel[1] >= 0);
  assert(densel[1] < (1 << 10));
  assert(densel[2] >= 0);
  assert(densel[2] < (1 << 10));
  assert(densel[3] >= 0);
  assert(densel[3] < (1 << 10));
  return (CacheKey)((static_cast<unsigned int>(densel[1]) << 20) | (static_cast<unsigned int>(densel[2]) << 10)
                    | (static_cast<unsigned int>(densel[3])));
}

void
ProjMatrixByDensel::cache_proj_matrix_elems_for_one_densel(const ProjMatrixElemsForOneDensel& probabilities) const
{
  if (cache_disabled)
    return;
 
  // insert probabilities into the collection
  cache_collection.insert(MapProjMatrixElemsForOneDensel::value_type(cache_key(probabilities.get_densel()), probabilities));
}
 
END_NAMESPACE_STIR