LSST Applications  21.0.0+c4f5df5339,21.0.0+e70536a077,21.0.0-1-ga51b5d4+7c60f8a6ea,21.0.0-10-gcf60f90+74aa0801fd,21.0.0-12-g63909ac9+643a1044a5,21.0.0-15-gedb9d5423+1041c3824f,21.0.0-2-g103fe59+a356b2badb,21.0.0-2-g1367e85+6d3f3f41db,21.0.0-2-g45278ab+e70536a077,21.0.0-2-g5242d73+6d3f3f41db,21.0.0-2-g7f82c8f+8d7c04eab9,21.0.0-2-g8f08a60+9c9a9cfcc8,21.0.0-2-ga326454+8d7c04eab9,21.0.0-2-gde069b7+bedfc5e1fb,21.0.0-2-gecfae73+6cb6558258,21.0.0-2-gfc62afb+6d3f3f41db,21.0.0-3-g21c7a62+f6e98b25aa,21.0.0-3-g357aad2+bd62456bea,21.0.0-3-g4be5c26+6d3f3f41db,21.0.0-3-g65f322c+03a4076c01,21.0.0-3-g7d9da8d+c4f5df5339,21.0.0-3-gaa929c8+c6b98066dc,21.0.0-3-gc44e71e+a26d5c1aea,21.0.0-3-ge02ed75+04b527a9d5,21.0.0-35-g64f566ff+b27e5ef93e,21.0.0-4-g591bb35+04b527a9d5,21.0.0-4-g88306b8+8773047b2e,21.0.0-4-gc004bbf+80a0b7acb7,21.0.0-4-gccdca77+a5c54364a0,21.0.0-4-ge8fba5a+ccfc328107,21.0.0-5-gdf36809+87b8d260e6,21.0.0-6-g00874e7+7eeda2b6ba,21.0.0-6-g2d4f3f3+e70536a077,21.0.0-6-g4e60332+04b527a9d5,21.0.0-6-g5ef7dad+f53629abd8,21.0.0-7-gc8ca178+b63e69433b,21.0.0-8-gfbe0b4b+c6b98066dc,w.2021.06
LSST Data Management Base Package
interpolate.cc
Go to the documentation of this file.
1 /*
2  * LSST Data Management System
3  * Copyright 2008-2016 AURA/LSST.
4  *
5  * This product includes software developed by the
6  * LSST Project (http://www.lsst.org/).
7  *
8  * This program is free software: you can redistribute it and/or modify
9  * it under the terms of the GNU General Public License as published by
10  * the Free Software Foundation, either version 3 of the License, or
11  * (at your option) any later version.
12  *
13  * This program is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16  * GNU General Public License for more details.
17  *
18  * You should have received a copy of the LSST License Statement and
19  * the GNU General Public License along with this program. If not,
20  * see <https://www.lsstcorp.org/LegalNotices/>.
21  */
22 
23 #include <pybind11/pybind11.h>
24 //#include <pybind11/operators.h>
25 #include <pybind11/stl.h>
26 
27 #include "ndarray/pybind11.h"
28 
30 
31 namespace py = pybind11;
32 using namespace pybind11::literals;
33 
34 using namespace lsst::afw::math;
35 
36 PYBIND11_MODULE(interpolate, mod) {
37  py::class_<Interpolate, std::shared_ptr<Interpolate>> clsInterpolate(mod, "Interpolate");
38  py::enum_<Interpolate::Style>(clsInterpolate, "Style")
39  .value("UNKNOWN", Interpolate::Style::UNKNOWN)
40  .value("CONSTANT", Interpolate::Style::CONSTANT)
41  .value("LINEAR", Interpolate::Style::LINEAR)
42  .value("NATURAL_SPLINE", Interpolate::Style::NATURAL_SPLINE)
43  .value("CUBIC_SPLINE", Interpolate::Style::CUBIC_SPLINE)
44  .value("CUBIC_SPLINE_PERIODIC", Interpolate::Style::CUBIC_SPLINE_PERIODIC)
45  .value("AKIMA_SPLINE", Interpolate::Style::AKIMA_SPLINE)
46  .value("AKIMA_SPLINE_PERIODIC", Interpolate::Style::AKIMA_SPLINE_PERIODIC)
47  .value("NUM_STYLES", Interpolate::Style::NUM_STYLES)
48  .export_values();
49 
50  clsInterpolate.def("interpolate", [](Interpolate &t, double const x) {
51  /*
52  We use a lambda function here because interpolate (with a double) is a virtual function and therefor
53  cannot be wrapped directly.
54  */
55  return t.interpolate(x);
56  });
57 
58  clsInterpolate.def("interpolate",
60  Interpolate::interpolate);
61  clsInterpolate.def(
62  "interpolate",
63  (ndarray::Array<double, 1> (Interpolate::*)(ndarray::Array<double const, 1> const &) const) &
64  Interpolate::interpolate);
65 
66  mod.def("makeInterpolate",
69  "x"_a, "y"_a, "style"_a = Interpolate::AKIMA_SPLINE);
70  mod.def("makeInterpolate", (std::shared_ptr<Interpolate>(*)(ndarray::Array<double const, 1> const &,
71  ndarray::Array<double const, 1> const &y,
73  "x"_a, "y"_a, "style"_a = Interpolate::AKIMA_SPLINE);
74 
75  mod.def("stringToInterpStyle", stringToInterpStyle, "style"_a);
76  mod.def("lookupMaxInterpStyle", lookupMaxInterpStyle, "n"_a);
77  mod.def("lookupMinInterpPoints", lookupMinInterpPoints, "style"_a);
78 }
virtual double interpolate(double const x) const =0
Interpolate::Style lookupMaxInterpStyle(int const n)
Get the highest order Interpolation::Style available for 'n' points.
Definition: Interpolate.cc:275
int lookupMinInterpPoints(Interpolate::Style const style)
Get the minimum number of points needed to use the requested interpolation style.
Definition: Interpolate.cc:314
Interpolate::Style stringToInterpStyle(std::string const &style)
Conversion function to switch a string to an Interpolate::Style.
Definition: Interpolate.cc:257
std::shared_ptr< Interpolate > makeInterpolate(std::vector< double > const &x, std::vector< double > const &y, Interpolate::Style const style=Interpolate::AKIMA_SPLINE)
A factory function to make Interpolate objects.
Definition: Interpolate.cc:343
PYBIND11_MODULE(interpolate, mod)
Definition: interpolate.cc:36
int y
Definition: SpanSet.cc:49
double x