Newer
Older
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
# imports
# standard
import random
# local, external
# local, internal
import src.topupopt.problems.esipp.signal as signal
# ******************************************************************************
# ******************************************************************************
def examples():
# **************************************************************************
# **************************************************************************
# test creating fixed signals
example_fixed_signals()
# test creating free signals
example_free_signals()
# test creating bounded signals
example_bounded_signals()
# test setting a signal
example_set_signal()
# test non-negative reals
example_nnr_signals()
# test binary signals
example_binary_signals()
# trigger errors that can only happen by messing with private/interval vars
example_peculiar_errors()
# test amplitude constrained signals
example_amplitude_constrained_signals()
# test amplitude constrained non-negative real signals
example_amplitude_constrained_nnr_signals()
# **************************************************************************
# **************************************************************************
# ******************************************************************************
# ******************************************************************************
def example_amplitude_constrained_nnr_signals():
# number of time intervals
number_intervals = 3
# **************************************************************************
# error-free examples
# create signal with max positive amplitude limit
sig = signal.AmplitudeConstrainedNNRSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=None,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert not sig.has_min_pos_amp_limit
assert sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# create signal with min positive amplitude limit
sig = signal.AmplitudeConstrainedNNRSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=3,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert not sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# create signal with positive constraints only
sig = signal.AmplitudeConstrainedNNRSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# **************************************************************************
# trigger errors
# by providing a non-numeric nr. of samples without specific lower bounds
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedNNRSignal(
number_samples=(number_intervals,), max_pos_amp_limit=10
)
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedNNRSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
lower_bounds=[-1 for i in range(number_intervals)],
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
# ******************************************************************************
# ******************************************************************************
def example_amplitude_constrained_signals():
# number of time intervals
number_intervals = 3
# **************************************************************************
# error-free examples
# create signal with max positive amplitude limit
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=None,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert not sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
sig.set_positive_amplitude(positive_amplitude=5)
sig.validate_positive_amplitude()
# use the tolerances and validate an otherwise invalid amplitude
sig.set_positive_amplitude(positive_amplitude=12)
sig.validate_positive_amplitude(tolerance=2)
# create signal with min positive amplitude limit
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert not sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
sig.set_positive_amplitude(positive_amplitude=5)
sig.validate_positive_amplitude()
# use the tolerances and validate an otherwise invalid amplitude
sig.set_positive_amplitude(positive_amplitude=1)
sig.validate_positive_amplitude(tolerance=2)
# create signal with max negative amplitude limit
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=None,
)
assert sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert not sig.has_max_pos_amp_limit
assert not sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
sig.set_negative_amplitude(negative_amplitude=5)
sig.validate_negative_amplitude()
# use the tolerances and validate an otherwise invalid amplitude
sig.set_negative_amplitude(negative_amplitude=12)
sig.validate_negative_amplitude(tolerance=2)
# create signal with min negative amplitude limit
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=None,
min_neg_amp_limit=3,
)
assert not sig.has_max_neg_amp_limit
assert sig.has_min_neg_amp_limit
assert not sig.has_max_pos_amp_limit
assert not sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
sig.set_negative_amplitude(negative_amplitude=5)
sig.validate_negative_amplitude()
# use the tolerances and validate an otherwise invalid amplitude
sig.set_negative_amplitude(negative_amplitude=1)
sig.validate_negative_amplitude(tolerance=2)
# create signal with positive constraints only
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# create signal with negative constraints only
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
assert sig.has_max_neg_amp_limit
assert sig.has_min_neg_amp_limit
assert not sig.has_max_pos_amp_limit
assert not sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# create amplitude constrained signal with all limits but without bounds
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=11,
min_neg_amp_limit=4,
)
assert sig.has_max_neg_amp_limit
assert sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# create amplitude constrained signal with all limits and with bounds
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=11,
min_neg_amp_limit=4,
lower_bounds=[-7 for i in range(number_intervals)],
upper_bounds=[15 for i in range(number_intervals)],
)
assert sig.has_max_neg_amp_limit
assert sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert sig.is_signal_bounded()
assert not sig.is_signal_fixed()
# set the signal using samples that do not violate the limits nor the bounds
sig.set_signal([5 for i in range(number_intervals)])
assert not sig.violates_amplitude_limits()
assert not sig.violates_bounds()
# set the signal using samples that violate the positive limits and bounds
sig.set_signal([25 for i in range(number_intervals)])
assert sig.violates_amplitude_limits()
assert sig.violates_bounds()
# set the signal using samples that violate the positive limits and bounds
sig.set_signal([-25 for i in range(number_intervals)])
assert sig.violates_amplitude_limits()
assert sig.violates_bounds()
# test external samples that do not violate the limits
assert not sig.violates_amplitude_limits(
samples=[0 for i in range(number_intervals)]
)
# test external samples that violate the positive limits
assert sig.violates_amplitude_limits(samples=[15 for i in range(number_intervals)])
# test external samples that violate the negative limits
assert sig.violates_amplitude_limits(samples=[-15 for i in range(number_intervals)])
# create amplitude constrained signal without limits or bounds
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert not sig.has_max_pos_amp_limit
assert not sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed()
assert not sig.violates_amplitude_limits() # because it has none
# **************************************************************************
# trigger errors
# by providing negative 'positive' amplitude limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=-10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=-3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing negative 'negative' amplitude limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=-11,
min_neg_amp_limit=4,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=11,
min_neg_amp_limit=-4,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing non-numeric or not None amplitude limits (e.g. tuple)
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=(10,),
min_pos_amp_limit=None,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=(3,),
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=(10,),
min_neg_amp_limit=None,
)
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=None,
min_neg_amp_limit=(3,),
)
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing bounds incompatible with positive limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=9,
min_pos_amp_limit=None,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
upper_bounds=None,
lower_bounds=[10 for i in range(number_intervals)],
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing bounds incompatible with negative limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=6,
min_neg_amp_limit=None,
upper_bounds=[-10 for i in range(number_intervals)],
lower_bounds=None,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing incompatible maximum and minimum positive limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=5,
min_pos_amp_limit=10,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing incompatible maximum and minimum negative limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=6,
min_neg_amp_limit=11,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing non-numeric or not None amplitude limits (e.g. tuple)
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
sig.set_positive_amplitude(positive_amplitude=(5,))
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing non-numeric or not None amplitude limits (e.g. tuple)
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
sig.set_negative_amplitude(negative_amplitude=(5,))
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by checking if bounds have been violated without there being samples
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
assert not sig.has_max_neg_amp_limit
assert not sig.has_min_neg_amp_limit
assert sig.has_max_pos_amp_limit
assert sig.has_min_pos_amp_limit
assert not sig.is_signal_bounded()
assert not sig.is_signal_fixed() # signal is not set
assert not sig.violates_amplitude_limits() # since the sig is not set
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a positive amplitude when there are no positive
# amplitude limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
sig.validate_negative_amplitude()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a negative amplitude when there are no negative
# amplitude limits
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
sig.validate_positive_amplitude()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a positive amplitude that exceeds its tolerated
# maximum, using the internal positive amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
sig.set_positive_amplitude(12)
sig.validate_positive_amplitude()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a positive amplitude that exceeds its tolerated
# maximum, using an externally supplied amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
sig.validate_positive_amplitude(12)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a positive amplitude that is below its tolerated
# minimum, using the internal positive amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
sig.set_positive_amplitude(2)
sig.validate_positive_amplitude()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a positive amplitude that is below its tolerated
# minimum, using an externally supplied amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=10,
min_pos_amp_limit=3,
max_neg_amp_limit=None,
min_neg_amp_limit=None,
)
sig.validate_positive_amplitude(2)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a negative amplitude that exceeds its tolerated
# maximum, using the internal negative amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
sig.set_negative_amplitude(12)
sig.validate_negative_amplitude()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a negative amplitude that exceeds its tolerated
# maximum, using an externally supplied amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
sig.validate_negative_amplitude(12)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a negative amplitude that is below its tolerated
# minimum, using the internal negative amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
sig.set_negative_amplitude(2)
sig.validate_negative_amplitude()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by seeking to validate a negative amplitude that is below its tolerated
# minimum, using an externally supplied amplitude
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.AmplitudeConstrainedSignal(
number_samples=number_intervals,
max_pos_amp_limit=None,
min_pos_amp_limit=None,
max_neg_amp_limit=10,
min_neg_amp_limit=3,
)
sig.validate_negative_amplitude(2)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# ******************************************************************************
# ******************************************************************************
def example_peculiar_errors():
# number of time intervals
number_intervals = 3
# **************************************************************************
# by providing samples as something other than a list, e.g. tuples
Pedro L. Magalhães
committed
error_was_raised = False
try:
_ = signal.Signal(
number_samples=number_intervals,
samples=(random.random() for i in range(number_intervals)),
lower_bounds=None,
upper_bounds=None,
)
except TypeError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing an incorrect number of samples
Pedro L. Magalhães
committed
error_was_raised = False
try:
_ = signal.Signal(
number_samples=number_intervals,
samples=[random.random() for i in range(number_intervals + 1)],
lower_bounds=None,
upper_bounds=None,
)
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# **************************************************************************
# **************************************************************************
# the tests below require messing with the internals
# by providing an incorrect number of lower bounds
lower_bounds = [5 for i in range(number_intervals)]
upper_bounds = [7 for i in range(number_intervals)]
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.Signal(
number_samples=number_intervals,
samples=None,
lower_bounds=lower_bounds,
upper_bounds=upper_bounds,
)
sig.lower_bounds = [random.random() for i in range(number_intervals + 1)]
sig.has_lower_bounds()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing an incorrect number of upper bounds
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.Signal(
number_samples=number_intervals,
samples=None,
lower_bounds=lower_bounds,
upper_bounds=upper_bounds,
)
sig.upper_bounds = [random.random() for i in range(number_intervals - 1)]
sig.has_upper_bounds()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing an incorrect number of samples
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.Signal(
number_samples=number_intervals,
samples=[random.random() for i in range(number_intervals)],
lower_bounds=None,
upper_bounds=None,
)
sig.samples = [random.random() for i in range(number_intervals - 1)]
sig.is_signal_fixed()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by deleting the lower bounds after creating the object
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.NonNegativeRealSignal(number_samples=number_intervals)
sig.lower_bounds = None
sig.is_lower_bounded = False
if not sig.are_bounds_nnr():
raise ValueError()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
# by providing negative upper bounds (requires even lower lower bounds)
Pedro L. Magalhães
committed
error_was_raised = False
try:
sig = signal.NonNegativeRealSignal(number_samples=number_intervals)
sig.is_upper_bounded = True
sig.upper_bounds = [-1 for i in range(number_intervals)]
if not sig.are_bounds_nnr():
raise ValueError()
except ValueError:
Pedro L. Magalhães
committed
error_was_raised = True
assert error_was_raised
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
# ******************************************************************************
# ******************************************************************************
def example_binary_signals():
# number of time intervals
number_intervals = 3
# **************************************************************************
# error-free examples
# create a binary signal and insert non-binary integer numbers in it
sig = signal.BinarySignal(number_samples=number_intervals)
assert sig.is_signal_fixed() == False
assert sig.is_signal_bounded() == True # it has upper and lower bounds
assert sig.violates_bounds() == False # because it is not fixed
assert sig.is_nnr() == False # because it is not fixed
assert sig.are_bounds_nnr() == True # because it is within the range [0,1]
sig.set_signal(samples=[3 for _ in range(number_intervals)])
assert sig.is_signal_binary_only(integrality_tolerance=0.0) == False
assert sig.is_signal_integer_only(integrality_tolerance=0.0) == True
assert sig.is_signal_binary_only(integrality_tolerance=None) == False
# create a binary signal and insert a noisy binary numbers in it
sig = signal.BinarySignal(number_samples=number_intervals)
amplitude = 0.1
deviation = [amplitude * (random.random() - 0.5) for _ in range(number_intervals)]
samples = [random.randint(0, 1) + deviation[i] for i in range(number_intervals)]
sig.set_signal(samples=samples)
assert sig.is_signal_fixed() == True
assert sig.is_signal_bounded() == True # it has upper and lower bounds
assert sig.violates_bounds(tolerance=amplitude) == False # it should not
assert sig.is_nnr(tolerance=amplitude) == True # it should be
assert sig.are_bounds_nnr() == True # because it is within the range [0,1]
assert (
sig.is_signal_binary_only(integrality_tolerance=amplitude) == True
) # since the tol. is the ampli.
assert (
sig.is_signal_integer_only(integrality_tolerance=amplitude) == True
) # since the tol. is the ampli.
assert (
sig.is_signal_binary_only(
integrality_tolerance=max([abs(max(deviation)), abs(min(deviation))])
* (1 - 0.1)
)
== False
) # because the tolerance was set below the maximum deviation
assert (
sig.is_signal_integer_only(
integrality_tolerance=max([abs(max(deviation)), abs(min(deviation))])
* (1 - 0.1)
)
== False
) # because the tolerance was set below the maximum deviation
assert sig.is_signal_binary_only(integrality_tolerance=None) == True
# create a binary signal
sig = signal.BinarySignal(number_samples=number_intervals)
assert sig.is_signal_binary_only() == False # because it is not fixed yet
assert sig.is_signal_integer_only() == False # because it is not fixed yet
assert sig.is_signal_fixed() == False
assert sig.is_signal_bounded() == True # it has upper and lower bounds