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
Functions
lsst.geom.testUtils Namespace Reference

Functions

 extraMsg (msg)
 
 assertAnglesAlmostEqual (testCase, ang0, ang1, maxDiff=0.001 *arcseconds, ignoreWrap=True, msg="Angles differ")
 
 assertPairsAlmostEqual (testCase, pair0, pair1, maxDiff=1e-7, msg="Pairs differ")
 
 assertPairListsAlmostEqual (testCase, list0, list1, maxDiff=1e-7, msg=None)
 
 assertSpherePointsAlmostEqual (testCase, sp0, sp1, maxSep=0.001 *arcseconds, msg="")
 
 assertSpherePointListsAlmostEqual (testCase, splist0, splist1, maxSep=0.001 *arcseconds, msg=None)
 
 assertBoxesAlmostEqual (testCase, box0, box1, maxDiff=1e-7, msg="Boxes differ")
 

Detailed Description

Utilities that should be imported into the lsst.geom namespace when
lsst.geom is used

In the case of the assert functions, importing them makes them available in lsst.utils.tests.TestCase

Function Documentation

◆ assertAnglesAlmostEqual()

lsst.geom.testUtils.assertAnglesAlmostEqual ( testCase,
ang0,
ang1,
maxDiff = 0.001*arcseconds,
ignoreWrap = True,
msg = "Angles differ" )
Assert that two `~lsst.geom.Angle`\ s are almost equal, ignoring
wrap differences by default.

If both arguments are NaN the assert will pass.  If one of the arguments
is NaN but the other is not the assert will fail.

Parameters
----------
testCase : `unittest.TestCase`
    test case the test is part of; an object supporting one method:
    fail(self, msgStr)
ang0 : `lsst.geom.Angle`
    angle 0
ang1 : `lsst.geom.Angle`
    angle 1
maxDiff : `lsst.geom.Angle`
    maximum difference between the two angles
ignoreWrap : `bool`
    ignore wrap when comparing the angles?

    - if True then wrap is ignored, e.g. 0 and 360 degrees are considered
      equal
    - if False then wrap matters, e.g. 0 and 360 degrees are considered
      different

msg : `str`
    exception message prefix; details of the error are appended after ": "

Raises
------
AssertionError
    Raised if the difference is greater than ``maxDiff``

Definition at line 46 of file testUtils.py.

47 ignoreWrap=True, msg="Angles differ"):
48 r"""Assert that two `~lsst.geom.Angle`\ s are almost equal, ignoring
49 wrap differences by default.
50
51 If both arguments are NaN the assert will pass. If one of the arguments
52 is NaN but the other is not the assert will fail.
53
54 Parameters
55 ----------
56 testCase : `unittest.TestCase`
57 test case the test is part of; an object supporting one method:
58 fail(self, msgStr)
59 ang0 : `lsst.geom.Angle`
60 angle 0
61 ang1 : `lsst.geom.Angle`
62 angle 1
63 maxDiff : `lsst.geom.Angle`
64 maximum difference between the two angles
65 ignoreWrap : `bool`
66 ignore wrap when comparing the angles?
67
68 - if True then wrap is ignored, e.g. 0 and 360 degrees are considered
69 equal
70 - if False then wrap matters, e.g. 0 and 360 degrees are considered
71 different
72
73 msg : `str`
74 exception message prefix; details of the error are appended after ": "
75
76 Raises
77 ------
78 AssertionError
79 Raised if the difference is greater than ``maxDiff``
80 """
81 isNan0 = math.isnan(ang0.asRadians())
82 isNan1 = math.isnan(ang1.asRadians())
83 if isNan0 and isNan1:
84 return
85 if isNan0:
86 testCase.fail("ang0 is NaN")
87 if isNan1:
88 testCase.fail("ang1 is NaN")
89 measDiff = ang1 - ang0
90 if ignoreWrap:
91 measDiff = measDiff.wrapCtr()
92 if abs(measDiff) > maxDiff:
93 testCase.fail("%s: measured difference %s arcsec > max allowed %s arcsec" %
94 (msg, measDiff.asArcseconds(), maxDiff.asArcseconds()))
95
96
97@lsst.utils.tests.inTestCase

◆ assertBoxesAlmostEqual()

lsst.geom.testUtils.assertBoxesAlmostEqual ( testCase,
box0,
box1,
maxDiff = 1e-7,
msg = "Boxes differ" )
Assert that two boxes (`~lsst.geom.Box2D` or `~lsst.geom.Box2I`) are
almost equal

Parameters
----------
testCase : `unittest.TestCase`
    test case the test is part of; an object supporting one method:
    fail(self, msgStr)
box0 : `lsst.geom.Box2D` or `lsst.geom.Box2I`
    box 0
box1 : `lsst.geom.Box2D` or `lsst.geom.Box2I`
    box 1
maxDiff : `float`
    maximum radial separation between the min points and max points
msg : `str`
    exception message prefix; details of the error are appended after ": "

Raises
------
AssertionError
    Raised if the radial difference of the min points or max points is
    greater than maxDiff

Notes
-----
.. warning::

   Does not compare types, just compares values.

Definition at line 250 of file testUtils.py.

250def assertBoxesAlmostEqual(testCase, box0, box1, maxDiff=1e-7, msg="Boxes differ"):
251 """Assert that two boxes (`~lsst.geom.Box2D` or `~lsst.geom.Box2I`) are
252 almost equal
253
254 Parameters
255 ----------
256 testCase : `unittest.TestCase`
257 test case the test is part of; an object supporting one method:
258 fail(self, msgStr)
259 box0 : `lsst.geom.Box2D` or `lsst.geom.Box2I`
260 box 0
261 box1 : `lsst.geom.Box2D` or `lsst.geom.Box2I`
262 box 1
263 maxDiff : `float`
264 maximum radial separation between the min points and max points
265 msg : `str`
266 exception message prefix; details of the error are appended after ": "
267
268 Raises
269 ------
270 AssertionError
271 Raised if the radial difference of the min points or max points is
272 greater than maxDiff
273
274 Notes
275 -----
276 .. warning::
277
278 Does not compare types, just compares values.
279 """
280 assertPairsAlmostEqual(testCase, box0.getMin(),
281 box1.getMin(), maxDiff=maxDiff, msg=msg + ": min")
282 assertPairsAlmostEqual(testCase, box0.getMax(),
283 box1.getMax(), maxDiff=maxDiff, msg=msg + ": max")

◆ assertPairListsAlmostEqual()

lsst.geom.testUtils.assertPairListsAlmostEqual ( testCase,
list0,
list1,
maxDiff = 1e-7,
msg = None )
Assert that two lists of Cartesian points are almost equal

Each point can be any indexable pair of two floats, including
Point2D or Extent2D, a list or a tuple.

Parameters
----------
testCase : `unittest.TestCase`
    test case the test is part of; an object supporting one method: fail(self, msgStr)
list0 : `list` of pairs of `float`
    list of pairs 0
list1 : `list` of pairs of `float`
    list of pairs 1
maxDiff : `float`
    maximum radial separation between the two points
msg : `str`
    additional information for the error message; appended after ": "

Raises
------
AssertionError
    Raised if the radial difference is greater than ``maxDiff``

Notes
-----
.. warning::

   Does not compare types, just values.

Definition at line 141 of file testUtils.py.

141def assertPairListsAlmostEqual(testCase, list0, list1, maxDiff=1e-7, msg=None):
142 """Assert that two lists of Cartesian points are almost equal
143
144 Each point can be any indexable pair of two floats, including
145 Point2D or Extent2D, a list or a tuple.
146
147 Parameters
148 ----------
149 testCase : `unittest.TestCase`
150 test case the test is part of; an object supporting one method: fail(self, msgStr)
151 list0 : `list` of pairs of `float`
152 list of pairs 0
153 list1 : `list` of pairs of `float`
154 list of pairs 1
155 maxDiff : `float`
156 maximum radial separation between the two points
157 msg : `str`
158 additional information for the error message; appended after ": "
159
160 Raises
161 ------
162 AssertionError
163 Raised if the radial difference is greater than ``maxDiff``
164
165 Notes
166 -----
167 .. warning::
168
169 Does not compare types, just values.
170 """
171 testCase.assertEqual(len(list0), len(list1))
172 lenList1 = np.array([len(val) for val in list0])
173 lenList2 = np.array([len(val) for val in list1])
174 testCase.assertTrue(np.all(lenList1 == 2))
175 testCase.assertTrue(np.all(lenList2 == 2))
176
177 diffArr = np.array([(val0[0] - val1[0], val0[1] - val1[1])
178 for val0, val1 in zip(list0, list1)], dtype=float)
179 sepArr = np.hypot(diffArr[:, 0], diffArr[:, 1])
180 badArr = sepArr > maxDiff
181 if np.any(badArr):
182 maxInd = np.argmax(sepArr)
183 testCase.fail("PairLists differ in %s places; max separation is at %s: %s > %s%s" %
184 (np.sum(badArr), maxInd, sepArr[maxInd], maxDiff, extraMsg(msg)))
185
186
187@lsst.utils.tests.inTestCase

◆ assertPairsAlmostEqual()

lsst.geom.testUtils.assertPairsAlmostEqual ( testCase,
pair0,
pair1,
maxDiff = 1e-7,
msg = "Pairs differ" )
Assert that two Cartesian points are almost equal.

Each point can be any indexable pair of two floats, including
Point2D or Extent2D, a list or a tuple.

Parameters
----------
testCase : `unittest.TestCase`
    test case the test is part of; an object supporting one method: fail(self, msgStr)
pair0 : pair of `float`
    pair 0
pair1 : pair of `floats`
    pair 1
maxDiff : `float`
    maximum radial separation between the two points
msg : `str`
    exception message prefix; details of the error are appended after ": "

Raises
------
AssertionError
    Raised if the radial difference is greater than ``maxDiff``

Notes
-----
.. warning::

   Does not compare types, just compares values.

Definition at line 98 of file testUtils.py.

98def assertPairsAlmostEqual(testCase, pair0, pair1, maxDiff=1e-7, msg="Pairs differ"):
99 """Assert that two Cartesian points are almost equal.
100
101 Each point can be any indexable pair of two floats, including
102 Point2D or Extent2D, a list or a tuple.
103
104 Parameters
105 ----------
106 testCase : `unittest.TestCase`
107 test case the test is part of; an object supporting one method: fail(self, msgStr)
108 pair0 : pair of `float`
109 pair 0
110 pair1 : pair of `floats`
111 pair 1
112 maxDiff : `float`
113 maximum radial separation between the two points
114 msg : `str`
115 exception message prefix; details of the error are appended after ": "
116
117 Raises
118 ------
119 AssertionError
120 Raised if the radial difference is greater than ``maxDiff``
121
122 Notes
123 -----
124 .. warning::
125
126 Does not compare types, just compares values.
127 """
128 if len(pair0) != 2:
129 raise RuntimeError("len(pair0)=%s != 2" % (len(pair0),))
130 if len(pair1) != 2:
131 raise RuntimeError("len(pair1)=%s != 2" % (len(pair1),))
132
133 pairDiff = [float(pair1[i] - pair0[i]) for i in range(2)]
134 measDiff = math.hypot(*pairDiff)
135 if measDiff > maxDiff:
136 testCase.fail("%s: measured radial distance = %s > maxDiff = %s, pair0=(%r, %r), pair1=(%r, %r)" %
137 (msg, measDiff, maxDiff, pair0[0], pair0[1], pair1[0], pair1[1]))
138
139
140@lsst.utils.tests.inTestCase

◆ assertSpherePointListsAlmostEqual()

lsst.geom.testUtils.assertSpherePointListsAlmostEqual ( testCase,
splist0,
splist1,
maxSep = 0.001*arcseconds,
msg = None )
Assert that two lists of `~lsst.geom.SpherePoint`\ s are almost equal

Parameters
----------
testCase : `unittest.TestCase`
test case the test is part of; an object supporting one method:
fail(self, msgStr)
splist0 : `list` of `lsst.geom.SpherePoint`
list of SpherePoints 0
splist1 : `list` of `lsst.geom.SpherePoint`
list of SpherePoints 1
maxSep : `lsst.geom.Angle`
maximum separation
msg : `str`
exception message prefix; details of the error are appended after ": "

Raises
------
AssertionError
The SpherePoint lists are not equal.

Definition at line 216 of file testUtils.py.

216def assertSpherePointListsAlmostEqual(testCase, splist0, splist1, maxSep=0.001*arcseconds, msg=None):
217 r"""Assert that two lists of `~lsst.geom.SpherePoint`\ s are almost equal
218
219 Parameters
220 ----------
221 testCase : `unittest.TestCase`
222 test case the test is part of; an object supporting one method:
223 fail(self, msgStr)
224 splist0 : `list` of `lsst.geom.SpherePoint`
225 list of SpherePoints 0
226 splist1 : `list` of `lsst.geom.SpherePoint`
227 list of SpherePoints 1
228 maxSep : `lsst.geom.Angle`
229 maximum separation
230 msg : `str`
231 exception message prefix; details of the error are appended after ": "
232
233 Raises
234 ------
235 AssertionError
236 The SpherePoint lists are not equal.
237 """
238 testCase.assertEqual(len(splist0), len(splist1), msg=msg)
239 sepArr = np.array([sp0.separation(sp1)
240 for sp0, sp1 in zip(splist0, splist1)])
241 badArr = sepArr > maxSep
242 if np.any(badArr):
243 maxInd = np.argmax(sepArr)
244 testCase.fail("SpherePointLists differ in %s places; max separation is at %s: %s\" > %s\"%s" %
245 (np.sum(badArr), maxInd, sepArr[maxInd].asArcseconds(),
246 maxSep.asArcseconds(), extraMsg(msg)))
247
248
249@lsst.utils.tests.inTestCase

◆ assertSpherePointsAlmostEqual()

lsst.geom.testUtils.assertSpherePointsAlmostEqual ( testCase,
sp0,
sp1,
maxSep = 0.001*arcseconds,
msg = "" )
Assert that two `~lsst.geom.SpherePoint`\ s are almost equal

Parameters
----------
testCase : `unittest.TestCase`
test case the test is part of; an object supporting one method:
fail(self, msgStr)
sp0 : `lsst.geom.SpherePoint`
SpherePoint 0
sp1 : `lsst.geom.SpherePoint`
SpherePoint 1
maxSep : `lsst.geom.Angle`
maximum separation
msg : `str`
extra information to be printed with any error message

Raises
------
AssertionError
The SpherePoints are not equal.

Definition at line 188 of file testUtils.py.

188def assertSpherePointsAlmostEqual(testCase, sp0, sp1, maxSep=0.001*arcseconds, msg=""):
189 r"""Assert that two `~lsst.geom.SpherePoint`\ s are almost equal
190
191 Parameters
192 ----------
193 testCase : `unittest.TestCase`
194 test case the test is part of; an object supporting one method:
195 fail(self, msgStr)
196 sp0 : `lsst.geom.SpherePoint`
197 SpherePoint 0
198 sp1 : `lsst.geom.SpherePoint`
199 SpherePoint 1
200 maxSep : `lsst.geom.Angle`
201 maximum separation
202 msg : `str`
203 extra information to be printed with any error message
204
205 Raises
206 ------
207 AssertionError
208 The SpherePoints are not equal.
209 """
210 if sp0.separation(sp1) > maxSep:
211 testCase.fail("Angular separation between %s and %s = %s\" > maxSep = %s\"%s" %
212 (sp0, sp1, sp0.separation(sp1).asArcseconds(), maxSep.asArcseconds(), extraMsg(msg)))
213
214
215@lsst.utils.tests.inTestCase

◆ extraMsg()

lsst.geom.testUtils.extraMsg ( msg)
Format extra error message, if any

Definition at line 37 of file testUtils.py.

37def extraMsg(msg):
38 """Format extra error message, if any
39 """
40 if msg:
41 return ": " + msg
42 return ""
43
44
45@lsst.utils.tests.inTestCase