lsst.sphgeom gd90a20c951+052faf71bd
Loading...
Searching...
No Matches
utils.h
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 * Copyright 2014-2015 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#ifndef LSST_SPHGEOM_UTILS_H_
24#define LSST_SPHGEOM_UTILS_H_
25
28
29#include "Angle.h"
30
31
32namespace lsst {
33namespace sphgeom {
34
35// Forward declarations
36class Vector3d;
37class UnitVector3d;
38
45double getMinSquaredChordLength(Vector3d const & v,
46 Vector3d const & a,
47 Vector3d const & b,
48 Vector3d const & n);
49
55double getMaxSquaredChordLength(Vector3d const & v,
56 Vector3d const & a,
57 Vector3d const & b,
58 Vector3d const & n);
59
63 return abs(x - c);
64}
65
69 Angle a = getMinAngleToCircle(x, c);
70 if (abs(x) <= abs(c)) {
71 return a + Angle(PI) - 2.0 * abs(c);
72 }
73 if (a < abs(x)) {
74 return Angle(PI) - 2.0 * abs(c) - a;
75 }
76 return Angle(PI) + 2.0 * abs(c) - a;
77}
78
82Vector3d getWeightedCentroid(UnitVector3d const & v0,
83 UnitVector3d const & v1,
84 UnitVector3d const & v2);
85
86}} // namespace lsst::sphgeom
87
88#endif // LSST_SPHGEOM_UTILS_H_
This file declares a class for representing angles.
Definition: Angle.h:43
Angle getMaxAngleToCircle(Angle x, Angle c)
Definition: utils.h:68
Angle getMinAngleToCircle(Angle x, Angle c)
Definition: utils.h:62
double getMaxSquaredChordLength(Vector3d const &v, Vector3d const &a, Vector3d const &b, Vector3d const &n)
Definition: utils.cc:58
double getMinSquaredChordLength(Vector3d const &v, Vector3d const &a, Vector3d const &b, Vector3d const &n)
Definition: utils.cc:36