BSplines1DRegularGrid.h

Go to the documentation of this file.

/*
  Copyright (C) 2005 - 2009-10-27, Hammersmith Imanet Ltd
  Copyright (C) 2013, University College London
  This file is part of STIR.

  SPDX-License-Identifier: Apache-2.0

  See STIR/LICENSE.txt for details
*/
#ifndef __stir_numerics_BSplines1DRegularGrid__H__
#define __stir_numerics_BSplines1DRegularGrid__H__

#include "stir/numerics/BSplines.h"
#include <vector>

START_NAMESPACE_STIR

namespace BSpline
{

template <typename out_elemT, typename in_elemT = out_elemT, typename constantsT = in_elemT>
class BSplines1DRegularGrid
{
private:
  typedef typename std::vector<out_elemT>::iterator RandIterOut;
  int input_size; // create in the constructor
  constantsT z1;
  constantsT z2;
  constantsT lambda;

  inline out_elemT BSplines_product(const int index, const pos_type relative_position, const bool if_deriv) const;

  inline void set_private_values(BSplineType);

  inline out_elemT compute_BSplines_value(const pos_type relative_position, const bool if_deriv) const;

  // sadly,VC6.0 needs definition of template members in the class definition
  template <class IterT>
  inline void set_coef(IterT input_begin_iterator, IterT input_end_iterator)
  {
    BSplines1DRegularGrid::input_size = static_cast<int>(input_end_iterator - input_begin_iterator);
    BSplines_coef_vector.resize(input_size);
    BSplines_coef(
        BSplines_coef_vector.begin(), BSplines_coef_vector.end(), input_begin_iterator, input_end_iterator, z1, z2, lambda);
  }

public:
  std::vector<out_elemT> BSplines_coef_vector;
  BSplineType spline_type;

  inline BSplines1DRegularGrid();

  inline explicit BSplines1DRegularGrid(const std::vector<in_elemT>& input_vector);

  template <class IterT>
  inline BSplines1DRegularGrid(const IterT input_begin_iterator, const IterT input_end_iterator)
  {
    set_private_values(cubic);
    set_coef(input_begin_iterator, input_end_iterator);
  }
  template <class IterT>
  inline BSplines1DRegularGrid(const IterT input_begin_iterator, const IterT input_end_iterator, const BSplineType this_type)
  {
    set_private_values(this_type);
    set_coef(input_begin_iterator, input_end_iterator);
  }

  inline BSplines1DRegularGrid(const std::vector<in_elemT>& input_vector, const BSplineType this_type);

  inline ~BSplines1DRegularGrid();

  inline out_elemT BSplines(const pos_type relative_position) const;

  inline out_elemT BSplines_1st_der(const pos_type relative_position) const;

  inline const out_elemT operator()(const pos_type relative_position) const;

  inline const std::vector<out_elemT>
  BSplines_output_sequence(RandIterOut output_relative_position_begin_iterator, // relative_position might be better float
                           RandIterOut output_relative_position_end_iterator);
  inline const std::vector<out_elemT> BSplines_output_sequence(std::vector<pos_type> output_relative_position);
};
} // namespace BSpline

END_NAMESPACE_STIR

#include "stir/numerics/BSplines1DRegularGrid.inl"

#endif