LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
multiShapeletFunction.cc
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 *
4 * This product includes software developed by the
5 * LSST Project (http://www.lsst.org/).
6 * See the COPYRIGHT file
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#include "pybind11/pybind11.h"
24#include "pybind11/stl.h"
25
26#include "ndarray/pybind11.h"
27
29
30namespace py = pybind11;
31using namespace pybind11::literals;
32
33namespace lsst {
34namespace shapelet {
35
36namespace {
37
38template <typename PyClass>
39void declareMultiShapeletFunctionMembers(PyClass &cls) {
40 using Class = MultiShapeletFunction;
41
42 cls.def(py::init<>());
43 cls.def(py::init<Class const &>());
44 cls.def(py::init<typename Class::ComponentList const &>());
45 cls.def(py::init<ShapeletFunction const &>());
46
47 cls.def("getComponents", [](Class const &self) {
48 auto &components = self.getComponents();
49
50 py::tuple t(components.size());
51 for (size_t i = 0; i < components.size(); ++i) {
52 t[i] = py::cast(components[i]);
53 }
54
55 return t;
56 });
57 cls.def("addComponent",
58 [](Class &self, typename Class::Component const &c) { self.getComponents().push_back(c); });
59 cls.def("normalize", &Class::normalize, "value"_a = 1.0);
60 cls.def("shiftInPlace", &Class::shiftInPlace);
61 cls.def("transformInPlace", &Class::transformInPlace);
62 cls.def("convolve", (Class (Class::*)(ShapeletFunction const &) const) & Class::convolve);
63 cls.def("convolve", (Class (Class::*)(Class const &) const) & Class::convolve);
64 cls.def("evaluate", &Class::evaluate);
65}
66
67template <typename PyClass>
68void declareMultiShapeletFunctionEvaluatorMembers(PyClass &cls) {
69 using Class = MultiShapeletFunctionEvaluator;
70
71 cls.def(py::init<MultiShapeletFunction const &>());
72
73 cls.def("__call__", (double (Class::*)(double, double) const) & Class::operator());
74 cls.def("__call__", (double (Class::*)(geom::Point2D const &) const) & Class::operator());
75 cls.def("__call__", (double (Class::*)(geom::Extent2D const &) const) & Class::operator());
76 cls.def("__call__",
77 (ndarray::Array<double, 1, 1> (Class::*)(ndarray::Array<double const, 1> const &,
78 ndarray::Array<double const, 1> const &) const) &
79 Class::operator());
80
81 cls.def("addToImage",
82 (void (Class::*)(ndarray::Array<double, 2, 1> const &, geom::Point2I const &) const) &
83 Class::addToImage,
84 "array"_a, "xy0"_a = geom::Point2I());
85 cls.def("addToImage", (void (Class::*)(afw::image::Image<double> &) const) & Class::addToImage,
86 "image"_a);
87
88 cls.def("integrate", &Class::integrate);
89 cls.def("computeMoments", &Class::computeMoments);
90 cls.def("update", &Class::update);
91}
92
93} // <anonymous>
94
96 using PyMultiShapeletFunction = py::class_<MultiShapeletFunction, std::shared_ptr<MultiShapeletFunction>>;
97
98 wrappers.wrapType(PyMultiShapeletFunction(wrappers.module, "MultiShapeletFunction"), [](auto &mod, auto &cls) {
99 declareMultiShapeletFunctionMembers(cls);
100
101 });
102
103 using PyMultiShapeletFunctionEvaluator =
104 py::class_<MultiShapeletFunctionEvaluator, std::shared_ptr<MultiShapeletFunctionEvaluator>>;
105
106 wrappers.wrapType(PyMultiShapeletFunctionEvaluator(wrappers.module, "MultiShapeletFunctionEvaluator"), [](auto &mod, auto &cls) {
107 declareMultiShapeletFunctionEvaluatorMembers(cls);
108 });
109}
110
111} // shapelet
112} // lsst
table::Key< table::Array< int > > components
A helper class for subdividing pybind11 module across multiple translation units (i....
Definition python.h:242
PyType wrapType(PyType cls, ClassWrapperCallback function, bool setModuleName=true)
Add a type (class or enum) wrapper, deferring method and other attribute definitions until finish() i...
Definition python.h:391
Extent< double, 2 > Extent2D
Definition Extent.h:400
Point< double, 2 > Point2D
Definition Point.h:324
Point< int, 2 > Point2I
Definition Point.h:321
void wrapMultiShapeletFunction(WrapperCollection &wrappers)