1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
|
# -*- coding: utf-8 -*-
from __future__ import division, unicode_literals, print_function, absolute_import
import os
import logging
from pint.compat import ndarray, unittest, np
from pint import logger, UnitRegistry
from pint.quantity import _Quantity
from logging.handlers import BufferingHandler
h = logging.StreamHandler()
f = logging.Formatter("%(levelname)s %(asctime)s %(funcName)s %(lineno)d %(message)s")
h.setLevel(logging.DEBUG)
h.setFormatter(f)
logger.addHandler(h)
logger.setLevel(logging.DEBUG)
class TestHandler(BufferingHandler):
def __init__(self):
# BufferingHandler takes a "capacity" argument
# so as to know when to flush. As we're overriding
# shouldFlush anyway, we can set a capacity of zero.
# You can call flush() manually to clear out the
# buffer.
BufferingHandler.__init__(self, 0)
def shouldFlush(self):
return False
def emit(self, record):
self.buffer.append(record.__dict__)
class TestCase(unittest.TestCase):
FORCE_NDARRAY = False
@classmethod
def setUpClass(cls):
cls.ureg = UnitRegistry(force_ndarray=cls.FORCE_NDARRAY)
cls.Q_ = cls.ureg.Quantity
def _get_comparable_magnitudes(self, first, second, msg):
if isinstance(first, _Quantity) and isinstance(second, _Quantity):
second = second.to(first)
self.assertEqual(first.units, second.units, msg=msg + 'Units are not equal.')
m1, m2 = first.magnitude, second.magnitude
elif isinstance(first, _Quantity):
self.assertTrue(first.dimensionless, msg=msg + 'The first is not dimensionless.')
first = first.to('')
m1, m2 = first.magnitude, second
elif isinstance(second, _Quantity):
self.assertTrue(second.dimensionless, msg=msg + 'The second is not dimensionless.')
second = second.to('')
m1, m2 = first, second.magnitude
else:
m1, m2 = first, second
return m1, m2
def assertSequenceEqual(self, seq1, seq2, msg=None, seq_type=None):
if isinstance(seq1, _Quantity):
if isinstance(seq2, _Quantity):
seq2 = seq2.to(seq1).magnitude
seq1 = seq1.magnitude
else:
seq1 = seq1.to('').magnitude
if isinstance(seq2, _Quantity):
if isinstance(seq1, _Quantity):
seq1 = seq1.to(seq2).magnitude
seq2 = seq2.magnitude
else:
seq2 = seq2.to('').magnitude
if isinstance(seq1, ndarray):
seq1 = seq1.tolist()
if isinstance(seq2, ndarray):
seq2 = seq2.tolist()
unittest.TestCase.assertSequenceEqual(self, seq1, seq2, msg, seq_type)
def assertQuantityEqual(self, first, second, msg=None):
if msg is None:
msg = 'Comparing %r and %r. ' % (first, second)
m1, m2 = self._get_comparable_magnitudes(first, second, msg)
if isinstance(m1, ndarray) or isinstance(m2, ndarray):
np.testing.assert_array_equal(m1, m2, err_msg=msg)
else:
unittest.TestCase.assertEqual(self, m1, m2, msg, )
def assertQuantityAlmostEqual(self, first, second, places=None, msg=None, delta=None):
if msg is None:
msg = 'Comparing %r and %r. ' % (first, second)
m1, m2 = self._get_comparable_magnitudes(first, second, msg)
if isinstance(m1, ndarray) or isinstance(m2, ndarray):
if delta is not None:
rtol, atol = 0, delta
else:
places = places or 7
rtol, atol = 10 ** (-places), 0
np.testing.assert_allclose(m1, m2, rtol=rtol, atol=atol, err_msg=msg)
else:
unittest.TestCase.assertAlmostEqual(self, m1, m2, places, msg, delta)
def assertAlmostEqual(self, first, second, places=None, msg=None, delta=None):
if isinstance(first, _Quantity) and isinstance(second, _Quantity):
second = second.to(first)
unittest.TestCase.assertAlmostEqual(self, first.magnitude, second.magnitude, places, msg, delta)
self.assertEqual(first.units, second.units)
elif isinstance(first, _Quantity):
self.assertTrue(first.dimensionless)
first = first.to('')
unittest.TestCase.assertAlmostEqual(self, first.magnitude, second, places, msg, delta)
elif isinstance(second, _Quantity):
self.assertTrue(second.dimensionless)
second = second.to('')
unittest.TestCase.assertAlmostEqual(self, first, second.magnitude, places, msg, delta)
else:
unittest.TestCase.assertAlmostEqual(self, first, second, places, msg, delta)
def assertAlmostEqualRelError(self, first, second, rel, msg=None):
if isinstance(first, _Quantity) and isinstance(second, _Quantity):
second = second.to(first)
val = abs((second - first) / (second + first))
elif isinstance(first, _Quantity):
self.assertTrue(first.dimensionless)
first = first.to('')
val = abs((second - first) / (second + first))
elif isinstance(second, _Quantity):
self.assertTrue(second.dimensionless)
second = second.to('')
val = abs((second - first) / (second + first))
else:
val = abs((second - first) / (second + first))
self.assertLess(val, rel, msg=msg)
def testsuite():
"""A testsuite that has all the pint tests.
"""
return unittest.TestLoader().discover(os.path.dirname(__file__))
def main():
"""Runs the testsuite as command line application."""
try:
unittest.main()
except Exception as e:
print('Error: %s' % e)
def run():
"""Run all tests i
"""
test_runner = unittest.TextTestRunner()
test_runner.run(testsuite())
|