def test_with_string4(self):
"""
Correlations with finite DoF
"""
context._context = Context()
x,y = multiple_ureal([1,2],[1,1],4)
r = 0.5
set_correlation(r,x,y)
z = result( x + y )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
db = persistence.dumps_json(ar)
context._context = Context()
ar = persistence.loads_json(db)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
self.assertEqual( get_correlation(x,y), r )
def test_with_file5(self):
"""
Complex
"""
wdir = os.getcwd()
fname = 'test_file.json'
path = os.path.join(wdir,fname)
context._context = Context()
x = ucomplex(1,[10,2,2,10],5)
y = ucomplex(1-6j,[10,2,2,10],7)
z = result( log( x * y ) )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
with open(path,'w') as f:
persistence.dump_json(f,ar)
context._context = Context()
with open(path,'r') as f:
ar = persistence.load_json(f)
os.remove(path)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
def test_with_string5(self):
"""
Complex
"""
context._context = Context()
x = ucomplex(1,[10,2,2,10],5)
y = ucomplex(1-6j,[10,2,2,10],7)
z = result( log( x * y ) )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
db = persistence.dumps_json(ar)
context._context = Context()
ar = persistence.loads_json(db)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
def test_with_string3(self):
"""
Dependent elementary UNs
"""
context._context = Context()
x = ureal(1,1,independent=False)
y = ureal(2,1,independent=False)
r = 0.5
set_correlation(r,x,y)
z = result( x + y )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
db = persistence.dumps_json(ar)
context._context = Context()
ar = persistence.loads_json(db)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
self.assertEqual( get_correlation(x,y), r )
def test_with_string2(self):
"""
Save to a file and then restore several times
to test the effectiveness of GTC's uid system.
"""
context._context = Context()
x = ureal(1,1,3,label='x')
y = ureal(2,1,4,label='y')
z = result( x + y )
ar = persistence.Archive()
# Saving only `z` means that when the archive
# is restored `x` and `y` are not recreated as UNs.
# However, Leaf nodes are created. We need to make sure
# that only one Leaf node gets created.
ar.add(z=z)
s = persistence.dumps_json(ar)
context._context = Context()
ar1 = persistence.loads_json(s)
z1 = ar1.extract('z')
# The attempt to create a new uncertain number
# is allowed but a new node is not created
ar2 = persistence.loads_json(s)
z2 = ar2.extract('z')
self.assertTrue(z2.is_intermediate)
self.assertTrue( z2._node is z1._node)
def test_with_file(self):
"""
Save to a file and then restore by reading
"""
wdir = os.getcwd()
fname = 'test_file.pkl'
path = os.path.join(wdir, fname)
context._context = Context()
x = ureal(1, 1)
y = ureal(2, 1)
z = result(x + y)
ar = persistence.Archive()
ar.add(x=x, y=y, z=z)
f = open(path, 'wb')
persistence.dump(f, ar)
f.close()
context._context = Context()
f = open(path, 'rb')
ar = persistence.load(f)
f.close()
os.remove(path)
x1, y1, z1 = ar.extract('x', 'y', 'z')
self.assertEqual(repr(x1), repr(x))
self.assertEqual(repr(y1), repr(y))
self.assertEqual(repr(z1), repr(z))
def test_multiple_names(self):
"""
Different names may apply to the same UN
"""
context._context = Context()
x1 = ureal(1, 1)
x2 = ureal(1, 1)
x3 = ureal(1, 1)
x4 = result(x1 + x2)
x5 = result(x2 + x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar.add(x1=x1, x2=x2, x3=x3, x4=x4, x5=x5, x6=x6, z1=x1, z4=x4, z6=x6)
ar._freeze()
# unpack and test
context._context = Context()
ar._thaw()
y1, y2, y3 = ar.extract('x1', 'x2', 'x3')
yz1 = ar['z1']
# This used to be the case, but thawing now allows multiple UNs
# self.assertTrue( yz1 is y1 )
self.assertTrue(yz1._node.uid is y1._node.uid)
self.assertEqual(repr(x1), repr(y1))
self.assertEqual(repr(yz1), repr(y1))
# NB most of these tests are redundant
# now that multiple names are restored
# to the same object.
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
yz4 = ar['z4']
yz6 = ar['z6']
# self.assertTrue( yz4 is y4 )
# self.assertTrue( yz6 is y6 )
self.assertTrue(yz4._node.uid is y4._node.uid)
self.assertTrue(yz6._node.uid is y6._node.uid)
self.assertEqual(repr(x4), repr(y4))
self.assertEqual(repr(yz4), repr(y4))
self.assertEqual(repr(x6), repr(y6))
self.assertEqual(repr(yz6), repr(y6))
a = component(x6, x4)
b = component(y6, y4)
self.assertTrue(equivalent(a, b, TOL))
b = component(yz6, y4)
self.assertTrue(equivalent(a, b, TOL))
b = component(yz6, yz4)
self.assertTrue(equivalent(a, b, TOL))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(yz6._u_components))
def test_complex_4(self):
"""
Same as 'complex_3' but using [] indexing
"""
context._context = Context()
x1 = ucomplex(1, [3, 2], independent=False)
x2 = ucomplex(1, [1, 1], independent=False)
x3 = ucomplex(1, [4, 5])
R = [.1, .2, .3, .4]
set_correlation(R, x1, x2)
# Make sure we get what we expect
self.assertTrue(equivalent_sequence(R, get_correlation(x1, x2)))
x4 = result(x1 * x2)
x5 = result(x2 * x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar['x1'] = x1
ar['x2'] = x2
ar['x3'] = x3
ar['x4'] = x4
ar['x5'] = x5
ar['x6'] = x6
ar._freeze()
context._context = Context()
ar._thaw()
y4 = ar['x4']
y5 = ar['x5']
y6 = ar['x6']
self.assertEqual(repr(x4), repr(y4))
self.assertEqual(repr(x5), repr(y5))
self.assertEqual(repr(x6), repr(y6))
self.assertTrue(
equivalent_sequence(get_correlation(x6, x5),
get_correlation(y6, y5)))
self.assertTrue(
equivalent_sequence(get_correlation(x6, x4),
get_correlation(y6, y4)))
self.assertTrue(equivalent_sequence(variance(x6), variance(y6)))
# Now look at the elementary influences
y1, y2 = ar.extract('x1', 'x2')
self.assertTrue(equivalent_sequence(R, get_correlation(y1, y2)))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(x6.real._u_components))
self.assertTrue(is_ordered(x6.imag._u_components))
def testGUMH2_b(self):
"""
GUM H2 as an ensemble. Only save the results this time.
"""
context._context = Context()
x = [4.999, 0.019661, 1.04446]
u = [0.0032, 0.0000095, 0.00075]
v, i, phi = multiple_ureal(x, u, 5)
set_correlation(-0.36, v, i)
set_correlation(0.86, v, phi)
set_correlation(-0.65, i, phi)
r = result(v * cos(phi) / i)
x = result(v * sin(phi) / i)
z = result(v / i)
ar = persistence.Archive()
ar.add(r=r, x=x, z=z)
ar._freeze()
# unpack without elementary uns and test
context._context = Context()
ar._thaw()
r1, x1, z1 = ar.extract('r', 'x', 'z')
self.assertEqual(repr(r1), repr(r))
self.assertEqual(repr(x1), repr(x))
self.assertEqual(repr(z1), repr(z))
# The DoF calculation would fail if the inputs
# are not part of the same ensemble.
self.assertTrue(equivalent(dof(r1), 5, TOL))
self.assertTrue(equivalent(dof(x1), 5, TOL))
self.assertTrue(equivalent(dof(z1), 5, TOL))
self.assertTrue(equivalent(value(r1), value(r), TOL))
self.assertTrue(equivalent(value(x1), value(x), TOL))
self.assertTrue(equivalent(value(z1), value(z), TOL))
self.assertTrue(equivalent(uncertainty(r1), uncertainty(r), TOL))
self.assertTrue(equivalent(uncertainty(x1), uncertainty(x), TOL))
self.assertTrue(equivalent(uncertainty(z1), uncertainty(z), TOL))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(r1._u_components))
self.assertTrue(is_ordered(x1._u_components))
self.assertTrue(is_ordered(z1._u_components))
def test_complex_2(self):
"""
Complex logical correlation
"""
context._context = Context()
x1 = ucomplex(1, [3, 2], 4)
x2 = ucomplex(1, [1, 1], 5)
x3 = ucomplex(1, [4, 5], 6)
x4 = result(x1 * x2)
x5 = result(x2 * x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar.add(x4=x4, x5=x5, x6=x6)
ar._freeze()
context._context = Context()
ar._thaw()
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
self.assertEqual(repr(x4), repr(y4))
self.assertEqual(repr(x5), repr(y5))
self.assertEqual(repr(x6), repr(y6))
self.assertTrue(equivalent_complex(x4.x, y4.x))
self.assertTrue(equivalent_sequence(x4.u, y4.u))
self.assertTrue(equivalent(x4.df, y4.df))
self.assertTrue(equivalent_complex(x5.x, y5.x))
self.assertTrue(equivalent_sequence(x5.u, y5.u))
self.assertTrue(equivalent(x5.df, y5.df))
self.assertTrue(equivalent_complex(x6.x, y6.x))
self.assertTrue(equivalent_sequence(x6.u, y6.u))
self.assertTrue(equivalent(x6.df, y6.df))
self.assertTrue(
equivalent_sequence(get_correlation(x6, x5),
get_correlation(y6, y5)))
self.assertTrue(
equivalent_sequence(get_correlation(x6, x4),
get_correlation(y6, y4)))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(x6.real._u_components))
self.assertTrue(is_ordered(x6.imag._u_components))
def test2(self):
"""
Logical correlation works - I
with all objects archived
"""
context._context = Context()
x1 = ureal(1, 1)
x2 = ureal(1, 1)
x3 = ureal(1, 1)
x4 = result(x1 + x2)
x5 = result(x2 + x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar.add(x1=x1, x2=x2, x3=x3, x4=x4, x5=x5, x6=x6)
ar._freeze()
# unpack and test
context._context = Context()
ar._thaw()
y1, y2, y3 = ar.extract('x1', 'x2', 'x3')
self.assertEqual(str(x1), str(y1))
self.assertEqual(str(x2), str(y2))
self.assertEqual(str(x3), str(y3))
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
self.assertEqual(str(x4), str(y4))
self.assertEqual(str(x5), str(y5))
self.assertEqual(str(x6), str(y6))
a = component(x6, x4)
b = component(y6, y4)
self.assertTrue(equivalent(a, b, TOL))
a = component(x6, x2)
b = component(y6, y2)
self.assertTrue(equivalent(a, b, TOL))
a = get_correlation(x4, x5)
b = get_correlation(y4, y5)
self.assertTrue(equivalent(a, b, TOL))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(y6._u_components))
self.assertTrue(is_ordered(y6._d_components))
self.assertTrue(is_ordered(y6._i_components))
def test4(self):
"""
Explicit correlation works - I
with all objects archived
"""
context._context = Context()
x1 = ureal(1, 1, independent=False)
x2 = ureal(1, 1, independent=False)
x3 = ureal(1, 1, independent=False)
x4 = result(x1 + x2)
x5 = result(x2 + x3)
x6 = result(x4 + x5)
set_correlation(0.5, x1, x2)
set_correlation(-0.25, x3, x2)
ar = persistence.Archive()
ar.add(x1=x1, x2=x2, x3=x3, x4=x4, x5=x5, x6=x6)
ar._freeze()
# unpack without elementary uns and test
context._context = Context()
ar._thaw()
y1, y2, y3 = ar.extract('x1', 'x2', 'x3')
a = get_correlation(x1, x2)
b = get_correlation(y1, y2)
self.assertTrue(equivalent(a, b, TOL))
a = get_correlation(x3, x2)
b = get_correlation(y3, y2)
self.assertTrue(equivalent(a, b, TOL))
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
self.assertEqual(str(x4), str(y4))
self.assertEqual(str(x5), str(y5))
self.assertEqual(str(x6), str(y6))
a = get_correlation(x4, x5)
b = get_correlation(y4, y5)
self.assertTrue(equivalent(a, b, TOL))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(y6._u_components))
self.assertTrue(is_ordered(y6._d_components))
self.assertTrue(is_ordered(y6._i_components))
def test_complex_3(self):
"""
Complex logical and external correlation
"""
context._context = Context()
x1 = ucomplex(1, [3, 2], independent=False)
x2 = ucomplex(1, [1, 1], independent=False)
x3 = ucomplex(1, [4, 5])
R = [.1, .2, .3, .4]
set_correlation(R, x1, x2)
# Make sure we get what we expect
self.assertTrue(equivalent_sequence(R, get_correlation(x1, x2)))
x4 = result(x1 * x2)
x5 = result(x2 * x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar.add(x1=x1, x2=x2, x3=x3, x4=x4, x5=x5, x6=x6)
ar._freeze()
context._context = Context()
ar._thaw()
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
self.assertEqual(repr(x4), repr(y4))
self.assertEqual(repr(x5), repr(y5))
self.assertEqual(repr(x6), repr(y6))
self.assertTrue(
equivalent_sequence(get_correlation(x6, x5),
get_correlation(y6, y5)))
self.assertTrue(
equivalent_sequence(get_correlation(x6, x4),
get_correlation(y6, y4)))
self.assertTrue(equivalent_sequence(variance(x6), variance(y6)))
# Now look at the elementary influences
y1, y2 = ar.extract('x1', 'x2')
self.assertTrue(equivalent_sequence(R, get_correlation(y1, y2)))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(x6.real._u_components))
self.assertTrue(is_ordered(x6.imag._u_components))
def test_attributes(self):
"""
Dict-like attributes: keys(), values() items()
and the 'iter' variants, also len()
"""
context._context = Context()
x = ureal(1, 1)
y = ureal(2, 1)
z = result(x + y)
ar = persistence.Archive()
self.assertEqual(len(ar), 0)
ar.add(x=x, z=z)
self.assertEqual(len(ar), 2)
names = sorted(['x', 'z'])
objs = sorted([x, z])
self.assertEqual(names, sorted(ar.keys()))
self.assertEqual(objs, sorted(ar.values()))
self.assertEqual(sorted(zip(names, objs)), sorted(ar.items()))
for i, k in enumerate(sorted(ar.iterkeys())):
self.assertEqual(names[i], k)
for i, k in enumerate(sorted(ar.itervalues())):
self.assertEqual(objs[i], k)
items = sorted(zip(names, objs))
for i, k in enumerate(sorted(ar.iteritems())):
self.assertEqual(items[i], k)
def test_intermediates_multipath(self):
"""
Check the sensitivities when the tree structure
has intermediates that depend on intermediates
In this case y is directly connected to x1, x4 and x5
but x4 and x5 also depend on x1 and x5 depends on x4
"""
context._context = Context()
x1 = ureal(1.2, 0.6, label='x1')
x2 = ureal(2.5, 1.6, label='x2')
x3 = ureal(-5.3, .77, label='x3')
x4 = result(2.0 * (x1 + x2))
x5 = result((x4 + x3)**2)
y = result(sqrt(x5 + x4 + x1))
ar = persistence.Archive()
ar.add(y=y, x4=x4, x5=x5, x1=x1)
ar._freeze()
# unpack without elementary uns and test
context._context = Context()
ar._thaw()
yy, xx4, xx5, xx1 = ar.extract('y', 'x4', 'x5', 'x1')
self.assertTrue(
equivalent(rp.u_component(y, x4), rp.u_component(yy, xx4), TOL))
self.assertTrue(
equivalent(rp.u_component(y, x5), rp.u_component(yy, xx5), TOL))
self.assertTrue(
equivalent(rp.u_component(y, x1), rp.u_component(yy, xx1), TOL))
self.assertTrue(
equivalent(rp.u_component(x5, x1), rp.u_component(xx5, xx1), TOL))
self.assertTrue(
equivalent(rp.u_component(x4, x1), rp.u_component(xx4, xx1), TOL))
self.assertTrue(
equivalent(rp.u_component(x5, x4), rp.u_component(xx5, xx4), TOL))
self.assertTrue(
equivalent(rp.u_component(x4, x5), rp.u_component(xx4, xx5), TOL))
def test_complex_1(self):
"""
Simple arithmetic - all objects stored
and restored.
"""
context._context = Context()
x = ucomplex(1, [1, 2], 4)
y = ucomplex(1, [3, 12], 3)
z = result(x * y)
ar = persistence.Archive()
ar.add(x=x)
ar.add(y=y)
ar.add(z=z)
ar._freeze()
# unpack
context._context = Context()
ar._thaw()
x1, y1, z1 = ar.extract('x', 'y', 'z')
self.assertEqual(repr(y1), repr(y))
self.assertEqual(repr(x1), repr(x))
self.assertTrue(equivalent_complex(x1.x, x.x))
self.assertTrue(equivalent_sequence(x1.u, x.u))
self.assertTrue(equivalent(x1.df, x.df))
self.assertTrue(equivalent_complex(y1.x, y.x))
self.assertTrue(equivalent_sequence(y1.u, y.u))
self.assertTrue(equivalent(y1.df, y.df))
self.assertTrue(equivalent_complex(z1.x, z.x))
self.assertTrue(equivalent_sequence(z1.u, z.u))
self.assertTrue(equivalent(z1.df, z.df))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(y1.real._u_components))
self.assertTrue(is_ordered(x1.real._u_components))
self.assertTrue(is_ordered(z1.real._u_components))
self.assertTrue(is_ordered(y1.imag._u_components))
self.assertTrue(is_ordered(x1.imag._u_components))
self.assertTrue(is_ordered(z1.imag._u_components))
def test6(self):
"""
Same as test 5 but with [] indexing
"""
context._context = Context()
x1 = ureal(1, 1, independent=False)
x2 = ureal(1, 1, independent=False)
x3 = ureal(1, 1, independent=False)
x4 = result(x1 + x2)
x5 = result(x2 + x3)
x6 = result(x4 + x5)
set_correlation(0.5, x1, x2)
set_correlation(-0.25, x3, x2)
ar = persistence.Archive()
ar['x1'] = x1
ar['x2'] = x2
ar['x3'] = x3
ar['x4'] = x4
ar['x5'] = x5
ar['x6'] = x6
ar._freeze()
# unpack without elementary uns and test
context._context = Context()
ar._thaw()
y4 = ar['x4']
y5 = ar['x5']
y6 = ar['x6']
self.assertEqual(str(x4), str(y4))
self.assertEqual(str(x5), str(y5))
self.assertEqual(str(x6), str(y6))
a = get_correlation(x4, x5)
b = get_correlation(y4, y5)
self.assertTrue(equivalent(a, b, TOL))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(y6._u_components))
self.assertTrue(is_ordered(y6._d_components))
self.assertTrue(is_ordered(y6._i_components))
def test_with_file2(self):
"""
Save to a file and then restore several times
to test the effectiveness of GTC's uid system.
"""
wdir = os.getcwd()
fname = 'test_file.json'
path = os.path.join(wdir,fname)
context._context = Context()
x = ureal(1,1,3,label='x')
y = ureal(2,1,4,label='y')
z = result( x + y )
ar = persistence.Archive()
# Saving only `z` means that when the archive
# is restored `x` and `y` are not recreated as UNs.
# However, Leaf nodes are created. We need to make sure
# that only one Leaf node gets created.
ar.add(z=z)
with open(path,'w') as f:
persistence.dump_json(f,ar)
context._context = Context()
with open(path,'r') as f:
ar1 = persistence.load_json(f)
z1 = ar1.extract('z')
self.assertEqual( repr(z1), repr(z) )
with open(path,'r') as f:
# The attempt to create the uncertain number again is allowed
# but should not create a new node object
ar2 = persistence.load_json(f)
z2 = ar2.extract('z')
self.assertTrue( z2.is_intermediate )
self.assertTrue( z1._node is z2._node )
os.remove(path)
def test_context_management(self):
"""
Context registers should be updated when all dependence
on an influence quantity disappears.
"""
context._context = Context()
x_seq = [1, 1, 1]
u_seq = [1, 1, 1]
label_seq = ['x1', 'x2', 'x3']
x1, x2, x3 = multiple_ureal(x_seq, u_seq, df=5, label_seq=label_seq)
set_correlation(.5, x1, x2)
set_correlation(.5, x2, x3)
x4 = result(x1 + x2)
x5 = result(x2 + x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar.add(x1=x1, x2=x2, x3=x3, x4=x4, x5=x5, x6=x6)
ar._freeze()
# unpack and test
context._context = Context()
ar._thaw()
y1, y2, y3 = ar.extract('x1', 'x2', 'x3')
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
self.assertTrue(
equivalent(get_correlation(x1, x2), get_correlation(y1, y2), TOL))
self.assertEqual(repr(x1), repr(y1))
self.assertEqual(repr(x2), repr(y2))
self.assertEqual(repr(x3), repr(y3))
self.assertEqual(repr(x4), repr(y4))
self.assertEqual(repr(x5), repr(y5))
self.assertEqual(repr(x6), repr(y6))
def test(self):
"""
Simple x,y z problem, but don't save
one of the elementary uns.
"""
context._context = Context()
x = ureal(1, 1)
y = ureal(2, 1)
z = result(x + y)
self.assertTrue(z.is_intermediate)
ar = persistence.Archive()
ar.add(x=x, z=z)
ar._freeze()
# Unpack and check
context._context = Context()
ar._thaw()
# Single argument unpacks to an object
# not a sequence
x1 = ar.extract('x')
z1 = ar.extract('z')
self.assertTrue(z1.is_intermediate)
self.assertTrue(
all(z_i.uid == z1_i.uid
for (z_i, z1_i) in izip(z._i_components, z1._i_components)))
a = component(z, x)
b = component(z1, x1)
self.assertTrue(equivalent(a, b, TOL))
# Make sure the vectors are well-formed
self.assertTrue(is_ordered(z1._u_components))
self.assertTrue(is_ordered(z1._d_components))
self.assertTrue(is_ordered(z1._i_components))
def test_with_file2(self):
"""
Save to a file and then restore several times
to test the effectiveness of GTC's uid system.
"""
wdir = os.getcwd()
fname = 'test_file.pkl'
path = os.path.join(wdir, fname)
context._context = Context()
x = ureal(1, 1, 3, label='x')
y = ureal(2, 1, 4, label='y')
z = result(x + y)
ar = persistence.Archive()
# Saving only `z` means that when the archive
# is restored `x` and `y` are not recreated as UNs.
# However, Leaf nodes are created. We need to make sure
# that only one Leaf node gets created.
ar.add(z=z)
with open(path, 'wb') as f:
persistence.dump(f, ar)
context._context = Context()
with open(path, 'rb') as f:
ar1 = persistence.load(f)
z1 = ar1.extract('z')
with open(path, 'rb') as f:
# The attempt to create the node again is caught
self.assertRaises(RuntimeError, persistence.load, f)
os.remove(path)
def test_errors(self):
"""
Error conditions:
## - trying to archive a complex and one of its components, or vice versa
## - trying to archive the same object with different names
- trying to extract with the wrong name
- adding to an archive that has been frozen
- extracting from an archive that has not been thawed
- restoring an archive to the same context
"""
context._context = Context()
x = ureal(1, 1)
y = ureal(2, 1)
z = result(x + y)
ar = persistence.Archive()
ar.add(x=x, z=z)
self.assertRaises(RuntimeError, ar.add, x=y) # same name
# These used to be illegal but are now allowed
zc = ucomplex(1, 1)
ar.add(zc=zc)
ar.add(zr=zc.real) # a component
zz = ucomplex(1, 1)
ar.add(zr2=zz.real)
ar.add(zz=zz) # a component was archived
ar._freeze()
self.assertRaises(RuntimeError, ar.extract, 'x') # cannot extract yet
context._context = Context()
ar._thaw()
self.assertRaises(RuntimeError, ar.add, y=y) # cannot add now
def test_with_file3(self):
"""
Dependent elementary UNs
"""
wdir = os.getcwd()
fname = 'test_file.json'
path = os.path.join(wdir,fname)
context._context = Context()
x = ureal(1,1,independent=False)
y = ureal(2,1,independent=False)
r = 0.5
set_correlation(r,x,y)
z = result( x + y )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
with open(path,'w') as f:
persistence.dump_json(f,ar)
context._context = Context()
with open(path,'r') as f:
ar = persistence.load_json(f)
os.remove(path)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
self.assertEqual( get_correlation(x,y), r )
def test1(self):
"""
Simple tests: elementary and intermediate
uncertain numbers restored OK? Saving
all objects.
"""
context._context = Context()
x = ureal(1, 1)
y = ureal(2, 1)
z = result(x + y)
ar = persistence.Archive()
ar.add(x=x, y=y, z=z)
ar._freeze()
# Unpack and check
context._context = Context()
ar._thaw()
x1, y1, z1 = ar.extract('x', 'y', 'z')
self.assertEqual(x.x, x1.x)
self.assertEqual(x.u, x1.u)
self.assertEqual(x.df, x1.df)
self.assertEqual(y.x, y1.x)
self.assertEqual(y.u, y1.u)
self.assertEqual(y.df, y1.df)
self.assertEqual(z.x, z1.x)
self.assertEqual(z.u, z1.u)
self.assertEqual(z.df, z1.df)
a = component(z, x)
b = component(z1, x1)
self.assertTrue(equivalent(a, b, TOL))
def test3b(self):
"""
Logical correlation works - III
with not all objects archived or restored
"""
context._context = Context()
x1 = ureal(1, 1)
x2 = ureal(1, 1)
x3 = ureal(1, 1)
x4 = result(x1 + x2)
x5 = result(x2 + x3)
x6 = result(x4 + x5)
ar = persistence.Archive()
ar.add(x4=x4, x5=x5, x6=x6)
ar._freeze()
# unpack without elementary uns and test
context._context = Context()
ar._thaw()
y4, y5, y6 = ar.extract('x4', 'x5', 'x6')
self.assertEqual(str(x4), str(y4))
self.assertEqual(str(x5), str(y5))
self.assertEqual(str(x6), str(y6))
a = component(x6, x4)
b = component(y6, y4)
self.assertTrue(equivalent(a, b, TOL))
# Make sure the vector of components is well-formed
self.assertTrue(is_ordered(y6._u_components))
self.assertTrue(is_ordered(y6._d_components))
self.assertTrue(is_ordered(y6._i_components))
def test_with_file4(self):
"""
Correlations with finite DoF
"""
wdir = os.getcwd()
fname = 'test_file.json'
path = os.path.join(wdir,fname)
context._context = Context()
x,y = multiple_ureal([1,2],[1,1],4)
r = 0.5
set_correlation(r,x,y)
z = result( x + y )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
with open(path,'w') as f:
persistence.dump_json(f,ar)
context._context = Context()
with open(path,'r') as f:
ar = persistence.load_json(f)
os.remove(path)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
self.assertEqual( get_correlation(x1,y1), r )
def test_with_string1(self):
"""
Simple save with intermediate
"""
context._context = Context()
x = ureal(1,1)
y = ureal(2,1)
z = result( x + y )
ar = persistence.Archive()
ar.add(x=x,y=y,z=z)
db = persistence.dumps_json(ar)
context._context = Context()
ar = persistence.loads_json(db)
x1, y1, z1 = ar.extract('x','y','z')
self.assertEqual( repr(x1), repr(x) )
self.assertEqual( repr(y1), repr(y) )
self.assertEqual( repr(z1), repr(z) )
def test_with_string(self):
"""
Save to a file and then restore by reading
"""
context._context = Context()
x = ureal(1, 1)
y = ureal(2, 1)
z = result(x + y)
ar = persistence.Archive()
ar.add(x=x, y=y, z=z)
db = persistence.dumps(ar)
context._context = Context()
ar = persistence.loads(db)
x1, y1, z1 = ar.extract('x', 'y', 'z')
self.assertEqual(repr(x1), repr(x))
self.assertEqual(repr(y1), repr(y))
self.assertEqual(repr(z1), repr(z))
def test_nop_intermediates(self):
"""
Archiving irrelevant intermediate components
"""
context._context = Context()
x1 = ucomplex(1, [3, 2])
x2 = ureal(1, 1)
x3 = constant(1.3)
x4 = result(x1 * x2)
x5 = result(x2 * x3)
x6 = result(x4 + x5)
x7 = result(mag_squared(x6))
ar = persistence.Archive()
ar.add(x1=x1, x2=x2, x3=x3, x4=x4.real, x5=x5, x6=x6)
ar.add(x7=x7)
ar._freeze()
context._context = Context()
ar._thaw()
y1 = ar['x1']
y2 = ar['x2']
y3 = ar['x3']
y6 = ar['x6']
y7 = result(ar['x7'])
self.assertEqual(repr(x1), repr(y1))
self.assertEqual(repr(x2), repr(y2))
self.assertEqual(repr(x3), repr(y3))
self.assertEqual(repr(x6), repr(y6))
self.assertEqual(repr(x7), repr(y7))
# Archived as a complex influence of a complex quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x6, x1), rp.u_component(y6,
y1)))
# Archived as a real influence of a complex quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x6, x2), rp.u_component(y6,
y2)))
# Archived as a complex influence of a real quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x7, x1), rp.u_component(y7,
y1)))
# Archived as a real influence of a real quantity
self.assertTrue(
equivalent(rp.u_component(x7, x2), rp.u_component(y7, y2)))
# Archived as a constant influence of a real quantity
self.assertTrue(
equivalent(rp.u_component(x7, x3), rp.u_component(y7, y3)))
# Archived as a constant influence of a complex quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x6, x3), rp.u_component(y6,
y3)))
def test_complex_5(self):
"""
Archive intermediate complex components
"""
context._context = Context()
x1 = ucomplex(1, [3, 2], independent=False)
x2 = ucomplex(1, [1, 1], independent=False)
x3 = ucomplex(1, [4, 5])
R = [.1, .2, .3, .4]
set_correlation(R, x1, x2)
# Make sure we get what we expect
self.assertTrue(equivalent_sequence(R, get_correlation(x1, x2)))
x4 = result(x1 * x2)
x5 = result(x2 * x3)
x6 = result(x4 + x5)
x7 = result(mag_squared(x6) + mag_squared(x4))
ar = persistence.Archive()
ar['x1'] = x1
ar['x2'] = x2
ar['x3'] = x3
ar['x4'] = x4.real # just the real component
ar['x5'] = x5
ar['x6'] = x6
ar['x7'] = x7
ar._freeze()
context._context = Context()
ar._thaw()
y4 = ar['x4']
y5 = ar['x5']
y6 = ar['x6']
y7 = ar['x7']
self.assertEqual(repr(x4.real), repr(y4))
self.assertEqual(repr(x5), repr(y5))
self.assertEqual(repr(x6), repr(y6))
# Archived as a complex influence of a complex quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x6, x5), rp.u_component(y6,
y5)))
# Archived as a real influence of a complex quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x6, x4.real),
rp.u_component(y6, y4)))
# Archived as a complex influence of a real quantity
self.assertTrue(
equivalent_sequence(rp.u_component(x7, x5), rp.u_component(y7,
y5)))
# Archived as a real influence of a real quantity
self.assertTrue(
equivalent(rp.u_component(x7, x4.real), rp.u_component(y7, y4)))