LSSTApplications  10.0-2-g4f67435,11.0.rc2+1,11.0.rc2+12,11.0.rc2+3,11.0.rc2+4,11.0.rc2+5,11.0.rc2+6,11.0.rc2+7,11.0.rc2+8
LSSTDataManagementBasePackage
SincCoeffs.h
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2014 LSST Corporation.
5  *
6  * This product includes software developed by the
7  * LSST Project (http://www.lsst.org/).
8  *
9  * This program is free software: you can redistribute it and/or modify
10  * it under the terms of the GNU General Public License as published by
11  * the Free Software Foundation, either version 3 of the License, or
12  * (at your option) any later version.
13  *
14  * This program is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17  * GNU General Public License for more details.
18  *
19  * You should have received a copy of the LSST License Statement and
20  * the GNU General Public License along with this program. If not,
21  * see <http://www.lsstcorp.org/LegalNotices/>.
22  */
23 
24 #ifndef LSST_MEAS_BASE_SincCoeffs_h_INCLUDED
25 #define LSST_MEAS_BASE_SincCoeffs_h_INCLUDED
26 
27 #include <map>
28 
29 #include "lsst/afw/image/Image.h"
31 
32 namespace lsst { namespace meas { namespace base {
33 
42 template<typename PixelT>
43 class SincCoeffs {
44 public:
45 
47 
53  static void cache(float rInner, float rOuter);
54 
60  static PTR(CoeffT const)
61  get(afw::geom::ellipses::Axes const & outerEllipse, float const innerRadiusFactor=0.0);
62 
64  static PTR(CoeffT)
65  calculate(afw::geom::ellipses::Axes const& outerEllipse, double const innerFactor=0.0);
66 
67 private:
68 
69  // A comparison function that doesn't require equality closer than machine epsilon
70  template <typename T>
71  struct FuzzyCompare {
72  bool operator()(T x, T y) const {
73  if (isEqual(x, y)) {
74  return false;
75  }
76  return (x - y < 0) ? true : false;
77  }
78  bool isEqual(T x, T y) const {
79  return ::fabs(x - y) < std::numeric_limits<T>::epsilon();
80  }
81  };
82 
83  typedef std::map<float, PTR(CoeffT), FuzzyCompare<float> > CoeffMap;
84  typedef std::map<float, CoeffMap, FuzzyCompare<float> > CoeffMapMap;
85  SincCoeffs() : _cache() {};
86  SincCoeffs(SincCoeffs const&); // unimplemented: singleton
87  void operator=(SincCoeffs const&); // unimplemented: singleton
88 
89  static SincCoeffs& getInstance();
90 
91  /*
92  * Search the cache for coefficients for an aperture
93  *
94  * If the coefficients are not cached, a null shared_ptr will be returned.
95  */
96  PTR(CoeffT const)
97  _lookup(afw::geom::ellipses::Axes const & outerEllipse, double const innerRadiusFactor=0.0) const;
98 
99  CoeffMapMap _cache; //< Cache of coefficients
100 };
101 
102 }}} // namespace lsst::meas::base
103 
104 #endif // !LSST_MEAS_BASE_SincCoeffs_h_INCLUDED
int y
static boost::shared_ptr< CoeffT const > get(afw::geom::ellipses::Axes const &outerEllipse, float const innerRadiusFactor=0.0)
Definition: SincCoeffs.cc:532
#define PTR(...)
Definition: base.h:41
Definitions and inlines for Axes.
std::map< float, CoeffMap, FuzzyCompare< float > > CoeffMapMap
Definition: SincCoeffs.h:84
static SincCoeffs & getInstance()
Definition: SincCoeffs.cc:508
boost::shared_ptr< CoeffT const > _lookup(afw::geom::ellipses::Axes const &outerEllipse, double const innerRadiusFactor=0.0) const
Definition: SincCoeffs.cc:540
Support for 2-D images.
static void cache(float rInner, float rOuter)
Definition: SincCoeffs.cc:515
int x
std::map< float, boost::shared_ptr< CoeffT >, FuzzyCompare< float > > CoeffMap
Definition: SincCoeffs.h:83
static boost::shared_ptr< CoeffT > calculate(afw::geom::ellipses::Axes const &outerEllipse, double const innerFactor=0.0)
Calculate the coefficients for an aperture.
Definition: SincCoeffs.cc:563
void operator=(SincCoeffs const &)
afw::image::Image< PixelT > CoeffT
Definition: SincCoeffs.h:46