ProjMatrixElemsForOneDensel.h

Go to the documentation of this file.

//
//

#ifndef __stir_recon_buildblock_ProjMatrixElemsForOneDensel__
#define __stir_recon_buildblock_ProjMatrixElemsForOneDensel__

/*
    Copyright (C) 2000- 2009, Hammersmith Imanet Ltd
    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0

    See STIR/LICENSE.txt for details
*/

#include "stir/recon_buildblock/ProjMatrixElemsForOneDenselValue.h"
#include "stir/Densel.h"
#include <vector>

START_NAMESPACE_STIR

class RelatedDensels;
template <int num_dimensions, typename elemT>
class DiscretisedDensity;
class Succeeded;

/*
  it might be a bit faster to derive this (privately) from
  std::vector<value_type> as opposed to having a member of
  that type.
  TODO: check
*/
class ProjMatrixElemsForOneDensel
{
public:
  typedef ProjMatrixElemsForOneDenselValue value_type;

private:
  typedef std::vector<value_type> Element_vector;

public:
  typedef Element_vector::iterator iterator;
  typedef Element_vector::const_iterator const_iterator;
  typedef Element_vector::size_type size_type;
  typedef Element_vector::difference_type difference_type;
  typedef std::random_access_iterator_tag iterator_category;

  typedef value_type& reference;
  typedef const value_type& const_reference;

  ProjMatrixElemsForOneDensel();
  explicit ProjMatrixElemsForOneDensel(const Densel& Densel, const int default_capacity = 300);

  /* rely on compiler-generated versions
  ProjMatrixElemsForOneDensel( const ProjMatrixElemsForOneDensel&);
  ProjMatrixElemsForOneDensel& operator=(const ProjMatrixElemsForOneDensel&) ;
  */

  Succeeded check_state() const;

  inline const Densel& get_densel() const;
  inline void set_densel(const Densel&);

  inline iterator begin();
  inline const_iterator begin() const;
  inline iterator end();
  inline const_iterator end() const;

  void erase();

  inline void push_back(const value_type&);
  void reserve(size_type max_number);
  inline size_type size() const;

  ProjMatrixElemsForOneDensel& operator*=(const float d);
  ProjMatrixElemsForOneDensel& operator/=(const float d);


  void sort();


  // TODO make sure we can have a const argument
  void merge(ProjMatrixElemsForOneDensel& lor);

#if 0  
  void write(fstream&fst) const;        
  void read(fstream&fst );
#endif


  float square_sum() const;

  //******************** projection operations ********************//
#if 0
  void back_project(DiscretisedDensity<3,float>&,
                    const Densel&) const;

  void forward_project(Densel&,
                      const DiscretisedDensity<3,float>&) const;
  void back_project(DiscretisedDensity<3,float>&,
                    const RelatedDensels&) const; 
  void forward_project(RelatedDensels&,
                       const DiscretisedDensity<3,float>&) const;

#endif

private:
  std::vector<value_type> elements;
  Densel densel;

  inline iterator erase(iterator it);
};

END_NAMESPACE_STIR

#include "stir/recon_buildblock/ProjMatrixElemsForOneDensel.inl"

#endif