lsst.jointcal g26f71f4252+3e455882a9
Loading...
Searching...
No Matches
Tripletlist.h
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2/*
3 * This file is part of jointcal.
4 *
5 * Developed for the LSST Data Management System.
6 * This product includes software developed by the LSST Project
7 * (https://www.lsst.org).
8 * See the COPYRIGHT file at the top-level directory of this distribution
9 * for details of code ownership.
10 *
11 * This program is free software: you can redistribute it and/or modify
12 * it under the terms of the GNU General Public License as published by
13 * the Free Software Foundation, either version 3 of the License, or
14 * (at your option) any later version.
15 *
16 * This program is distributed in the hope that it will be useful,
17 * but WITHOUT ANY WARRANTY; without even the implied warranty of
18 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19 * GNU General Public License for more details.
20 *
21 * You should have received a copy of the GNU General Public License
22 * along with this program. If not, see <https://www.gnu.org/licenses/>.
23 */
24
25#ifndef LSST_JOINTCAL_TRIPLETLIST_H
26#define LSST_JOINTCAL_TRIPLETLIST_H
27
28#include "Eigen/Sparse"
29
30#include <vector>
31
32namespace lsst {
33namespace jointcal {
34
35using Trip = Eigen::Triplet<double>;
36
37// at the moment this class implements the eigen format.
38// it would be wise to implement it differently if talking to cholmod
39class TripletList : public std::vector<Trip> {
40public:
41 TripletList(std::size_t count) : _nextFreeIndex(0) { reserve(count); };
42
43 void addTriplet(Eigen::Index i, Eigen::Index j, double val) {
44 push_back(Trip(i, j, val));
45 }
46
47 Eigen::Index getNextFreeIndex() const { return _nextFreeIndex; }
48
49 void setNextFreeIndex(Eigen::Index index) { _nextFreeIndex = index; }
50
51private:
52 Eigen::Index _nextFreeIndex;
53};
54} // namespace jointcal
55} // namespace lsst
56
57#endif // LSST_JOINTCAL_TRIPLETLIST_H
ImageT val
TripletList(std::size_t count)
Definition: Tripletlist.h:41
Eigen::Index getNextFreeIndex() const
Definition: Tripletlist.h:47
void setNextFreeIndex(Eigen::Index index)
Definition: Tripletlist.h:49
void addTriplet(Eigen::Index i, Eigen::Index j, double val)
Definition: Tripletlist.h:43
T count(T... args)
Eigen::Triplet< double > Trip
Definition: Tripletlist.h:35
Trip push_back(Trip ... args)
Trip reserve(Trip ... args)