LSST Applications g04e9c324dd+8c5ae1fdc5,g134cb467dc+b203dec576,g18429d2f64+358861cd2c,g199a45376c+0ba108daf9,g1fd858c14a+dd066899e3,g262e1987ae+ebfced1d55,g29ae962dfc+72fd90588e,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+b668f15bc5,g4595892280+3897dae354,g47891489e3+abcf9c3559,g4d44eb3520+fb4ddce128,g53246c7159+8c5ae1fdc5,g67b6fd64d1+abcf9c3559,g67fd3c3899+1f72b5a9f7,g74acd417e5+cb6b47f07b,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+abcf9c3559,g8d7436a09f+bcf525d20c,g8ea07a8fe4+9f5ccc88ac,g90f42f885a+6054cc57f1,g97be763408+06f794da49,g9dd6db0277+1f72b5a9f7,ga681d05dcb+7e36ad54cd,gabf8522325+735880ea63,gac2eed3f23+abcf9c3559,gb89ab40317+abcf9c3559,gbf99507273+8c5ae1fdc5,gd8ff7fe66e+1f72b5a9f7,gdab6d2f7ff+cb6b47f07b,gdc713202bf+1f72b5a9f7,gdfd2d52018+8225f2b331,ge365c994fd+375fc21c71,ge410e46f29+abcf9c3559,geaed405ab2+562b3308c0,gf9a733ac38+8c5ae1fdc5,w.2025.35
LSST Data Management Base Package
Loading...
Searching...
No Matches
matrixBuilder.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
25#include "ndarray/pybind11.h"
26
29
30namespace py = pybind11;
31using namespace pybind11::literals;
32
33namespace lsst {
34namespace shapelet {
35
36namespace {
37
38template <typename T>
39py::classh<MatrixBuilder<T>> declareMatrixBuilder(
40 lsst::cpputils::python::WrapperCollection &wrappers, std::string const &suffix) {
41 using Class = MatrixBuilder<T>;
42 using PyCLass = py::classh<Class>;
43 std::string name = "MatrixBuilder" + suffix;
44
45 return wrappers.wrapType(PyCLass(wrappers.module, name.c_str()), [](auto &mod, auto &cls) {
46 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &, int>(),
47 "x"_a, "y"_a, "order"_a);
48 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &, int,
49 ShapeletFunction const &>(),
50 "x"_a, "y"_a, "order"_a, "psf"_a);
51 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &,
52 MultiShapeletBasis const &>(),
53 "x"_a, "y"_a, "basis"_a);
54 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &,
55 MultiShapeletBasis const &, MultiShapeletFunction const &>(),
56 "x"_a, "y"_a, "basis"_a, "psf"_a);
57
58 cls.def("getDataSize", &Class::getDataSize);
59 cls.def("getBasisSize", &Class::getBasisSize);
60 cls.def("allocateOutput", &Class::allocateOutput);
61
62 cls.def("__call__",
63 (void (Class::*)(ndarray::Array<T, 2, -1> const &, afw::geom::ellipses::Ellipse const &) const) &
64 Class::operator());
65 cls.def("__call__", (ndarray::Array<T, 2, -2> (Class::*)(afw::geom::ellipses::Ellipse const &) const) &
66 Class::operator());
67 });
68}
69
70template <typename T>
71py::classh<MatrixBuilderWorkspace<T>>
72declareMatrixBuilderWorkspace(lsst::cpputils::python::WrapperCollection &wrappers, std::string const &suffix) {
73 using Class = MatrixBuilderWorkspace<T>;
74 using PyClass = py::classh<Class>;
75 std::string name = "MatrixBuilderWorkspace" + suffix;
76
77 return wrappers.wrapType(PyClass(wrappers.module, name.c_str()), [](auto &mod, auto &cls) {
78
79 cls.def(py::init<int>(), "size"_a);
80 cls.def(py::init<Class const &>(), "other"_a);
81
82 cls.def("getRemaining", &Class::getRemaining);
83 });
84}
85
86template <typename T>
87py::classh<MatrixBuilderFactory<T>> declareMatrixBuilderFactory(
88 lsst::cpputils::python::WrapperCollection &wrappers, std::string const &suffix) {
89 using Class = MatrixBuilderFactory<T>;
90 using PyClass = py::classh<Class>;
91 std::string name = "MatrixBuilderFactory" + suffix;
92
93 return wrappers.wrapType(PyClass(wrappers.module, name.c_str()), [](auto &mod, auto &cls) {
94 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &, int>(),
95 "x"_a, "y"_a, "order"_a);
96 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &, int,
97 ShapeletFunction const &>(),
98 "x"_a, "y"_a, "order"_a, "psf"_a);
99 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &,
100 MultiShapeletBasis const &>(),
101 "x"_a, "y"_a, "basis"_a);
102 cls.def(py::init<ndarray::Array<T const, 1, 1> const &, ndarray::Array<T const, 1, 1> const &,
103 MultiShapeletBasis const &, MultiShapeletFunction const &>(),
104 "x"_a, "y"_a, "basis"_a, "psf"_a);
105
106 cls.def("__call__", (MatrixBuilder<T> (Class::*)() const) &Class::operator());
107 cls.def("__call__", (MatrixBuilder<T> (Class::*)(typename Class::Workspace &) const) &Class::operator(),
108 "workspace"_a);
109
110 cls.def("getDataSize", &Class::getDataSize);
111 cls.def("getBasisSize", &Class::getBasisSize);
112 cls.def("computeWorkspace", &Class::computeWorkspace);
113 });
114}
115
116template <typename T>
117void declareMatrixBuilderTemplates(lsst::cpputils::python::WrapperCollection &wrappers, std::string const &suffix) {
118 auto &mod = wrappers.module;
119 auto clsMatrixBuilder = declareMatrixBuilder<T>(wrappers, suffix);
120 auto clsMatrixBuilderWorkspace = declareMatrixBuilderWorkspace<T>(wrappers, suffix);
121 auto clsMatrixBuilderFactory = declareMatrixBuilderFactory<T>(wrappers, suffix);
122
123 clsMatrixBuilder.attr("Workspace") = clsMatrixBuilderWorkspace;
124 clsMatrixBuilder.attr("Factory") = clsMatrixBuilderFactory;
125
126 clsMatrixBuilderFactory.attr("Workspace") = clsMatrixBuilderWorkspace;
127 clsMatrixBuilderFactory.attr("Builder") = clsMatrixBuilder;
128}
129
130} // <anonymous>
131
133 declareMatrixBuilderTemplates<float>(wrappers, "F");
134 declareMatrixBuilderTemplates<double>(wrappers, "D");
135}
136
137} // shapelet
138} // lsst
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
pybind11::module module
The module object passed to the PYBIND11_MODULE block that contains this WrapperCollection.
Definition python.h:448
A factory class for MatrixBuilder, providing more control over workspace memory.
Class that evaluates a (multi-)shapelet basis at predefined points.
Reusable, shareable workspace for MatrixBuilder.
py::classh< PixelAreaBoundedField, BoundedField > PyClass
void wrapMatrixBuilder(WrapperCollection &wrappers)