def products_layer(inputs,
num_inputs,
num_input_cols,
num_prods,
inf_type,
indices=None,
log=False,
output=None):
products_inputs = []
num_or_size_prods = []
if isinstance(inputs,
list): # Is a list of RawLeaf inputs - Multiple inputs
for inps, n_inp_cols, n_prods in zip(inputs, num_input_cols,
num_prods):
num_inputs = len(inps)
# Create permuted indices based on number and size of inputs
inds = map(int, np.arange(n_inp_cols))
permuted_inds = list(product(inds, repeat=num_inputs))
permuted_inds_list = [list(elem) for elem in permuted_inds]
permuted_inds_list_of_list = []
for elem in permuted_inds_list:
permuted_inds_list_of_list.append(
[elem[i:i + 1] for i in range(0, len(elem), 1)])
# Create inputs list by combining inputs and indices
permuted_inputs = []
for indices in permuted_inds_list_of_list:
permuted_inputs.append(
[tuple(i) for i in zip(inps, indices)])
products_inputs += list(chain.from_iterable(permuted_inputs))
# Create products-size list
num_or_size_prods += [num_inputs] * n_prods
else: # Is a single input of type RawLeaf - A single input
outer_offset = 0
permuted_inds_list = []
for n_inps, n_inp_cols in zip(num_inputs, num_input_cols):
# Create permuted indices based on number and size of inputs
inds = map(int, np.arange(n_inp_cols))
permuted_inds = list(product(inds, repeat=n_inps))
offsets = \
np.array(list(range(0, (n_inps * n_inp_cols), n_inp_cols))) \
+ outer_offset
outer_offset += n_inps * n_inp_cols
for perm_inds in permuted_inds:
permuted_inds_list.append([
p_ind + off
for p_ind, off in zip(list(perm_inds), offsets)
])
# Content of list object 'perm_inds' needs to be of type int, if not
# input_parser in Input class complains
products_inputs = [(inputs, list(map(int, perm_inds)))
for perm_inds in permuted_inds_list]
num_or_size_prods = [
len(perm_inds) for perm_inds in permuted_inds_list
]
# Generate a single ProductsLayer node, modeling 'sum(num_prods)' products
# within, connecting it to inputs
p = spn.ProductsLayer(*products_inputs,
num_or_size_prods=num_or_size_prods)
# Connect all product nodes to a single root Sum node and generate its
# weights
root = spn.Sum(p)
root.generate_weights()
if log:
mpe_path_gen = spn.MPEPath(value_inference_type=inf_type, log=True)
else:
mpe_path_gen = spn.MPEPath(value_inference_type=inf_type,
log=False)
mpe_path_gen.get_mpe_path(root)
if isinstance(inputs, list):
path_ops = [
mpe_path_gen.counts[inp]
for inp in list(chain.from_iterable(inputs))
]
else:
path_ops = mpe_path_gen.counts[inputs]
return spn.initialize_weights(root), path_ops
def test(input_sizes,
num_or_size_prods,
batch_size=7,
single_input=False,
indexed=True):
self.tearDown()
with self.subTest(input_sizes=input_sizes,
batch_size=batch_size,
num_or_size_prods=num_or_size_prods,
single_input=single_input,
indexed=indexed):
# Calculate num-prods or prod-input-size, based on type and
# values of the parameter 'num_or_size_prods'
if isinstance(num_or_size_prods, int):
num_prods = num_or_size_prods
prod_input_sizes = [sum(input_sizes) // num_prods
] * num_prods
elif isinstance(num_or_size_prods, list):
num_prods = len(num_or_size_prods)
prod_input_sizes = num_or_size_prods
inputs = []
indices = []
true_outputs = []
if single_input:
# Size of the only input variable - would be = size of the
# biggest product modelled
num_vars = max(prod_input_sizes)
# Randomly choose to create either an IndicatorLeaf or a RawLeaf node
contvar = random.choice([True, False])
# Create an input node
input_var = spn.RawLeaf(num_vars=num_vars) if contvar else \
spn.IndicatorLeaf(num_vars=num_vars, num_vals=num_vals)
inputs = list(
itertools.repeat(input_var, len(prod_input_sizes)))
# Generate random input-indices for each product modelled,
# based on product-input-sizes
indices = [
random.sample(range(num_vars if contvar else num_vars *
num_vals),
k=p_size) for p_size in prod_input_sizes
]
# Create Zeros matrices, with the same size as the created
# input node, one per product modelled
true_outputs = [
np.zeros(
(batch_size,
num_vars if contvar else num_vars * num_vals))
for _ in prod_input_sizes
]
else:
# Create random (IndicatorLeaf or ContVar) inputs of random size each
for i_size in input_sizes:
if indexed:
# Choose a random size for an input node
num_vars = np.random.randint(low=i_size,
high=i_size + 5)
else:
num_vars = i_size
# Randomly choose to create either a IndicatorLeaf or a ContVar input
if random.choice([True, False]):
# Create an IndicatorLeaf input node
inputs.append(
spn.IndicatorLeaf(num_vars=num_vars,
num_vals=num_vals))
# Generate random indices for the created input
indices.append(
random.sample(range(num_vars * num_vals),
k=i_size))
# Create a Zeros matrix, with the same size as the
# created input node
true_outputs.append(
np.zeros((batch_size, num_vars * num_vals)))
else:
# Create anRawLeaf input node
inputs.append(spn.RawLeaf(num_vars=num_vars))
# Generate random indices for the created input
indices.append(
random.sample(range(num_vars), k=i_size))
# Create a Zeros matrix, with the same size as the
# created input node
true_outputs.append(
np.zeros((batch_size, num_vars)))
# Create the ProductsLayer node
p = spn.ProductsLayer(*list(zip(inputs, indices)),
num_or_size_prods=num_or_size_prods)
# Create counts as a placeholder
counts = tf.placeholder(tf.float32, shape=(None, num_prods))
# Create mpe path op
ops = p._compute_log_mpe_path(
tf.identity(counts), *[inp.get_value() for inp in inputs])
ops = [op for op in ops if op is not None]
# Generate random counts feed
coutnts_feed = np.random.randint(100,
size=(batch_size, num_prods))
# Create a session and execute the generated op
with self.test_session() as sess:
outputs = sess.run(ops, feed_dict={counts: coutnts_feed})
# Calculate true-output
# (1) Split counts into num_prod sub-matrices
split_counts = np.split(coutnts_feed, num_prods, axis=1)
# (2) Tile each counts matrix to the size of prod-input
tiled_counts_list = [
np.tile(sc, [1, p_inp_size])
for sc, p_inp_size in zip(split_counts, prod_input_sizes)
]
# (3) Concat all tiled counts to a single matrix
concatenated_counts = np.concatenate(tiled_counts_list, axis=1)
# (4) Split the concatenated counts matrix into input-sizes
if single_input:
value_counts = np.split(concatenated_counts,
np.cumsum(prod_input_sizes),
axis=1)
else:
value_counts = np.split(concatenated_counts,
np.cumsum(input_sizes),
axis=1)
# (5) Scatter each count matrix to inputs
for t_out, i_indices, v_counts in zip(true_outputs, indices,
value_counts):
t_out[:, i_indices] = v_counts
# (6) Sum scattered counts together, for the single-input case
if single_input:
true_outputs = np.sum(np.array(true_outputs),
axis=0,
keepdims=True)
# Assert outputs with true_outputs
for out, t_out in zip(outputs, true_outputs):
np.testing.assert_array_almost_equal(out, t_out)
def test_comput_scope(self):
"""Calculating scope of ProductsLayer"""
# Create graph
v12 = spn.IndicatorLeaf(num_vars=2, num_vals=4, name="V12")
v34 = spn.RawLeaf(num_vars=2, name="V34")
s1 = spn.Sum((v12, [0, 1, 2, 3]), name="S1")
s1.generate_latent_indicators()
s2 = spn.Sum((v12, [4, 5, 6, 7]), name="S2")
pl1 = spn.ProductsLayer((v12, [0, 5, 6, 7]), (v12, [3, 4]),
v34,
num_or_size_prods=[4, 3, 1],
name="PL1")
n1 = spn.Concat(s1, s2, (pl1, [2]), name="N1")
n2 = spn.Concat((pl1, [0]), (pl1, [1]), name="N2")
s3 = spn.Sum(pl1, name="S3")
s3.generate_latent_indicators()
pl2 = spn.ProductsLayer((n1, [0, 1]), (n1, 2), (n2, 0), (pl1, [1]),
n2,
s3, (n2, 1),
s3,
pl1,
num_or_size_prods=[2, 3, 3, 5],
name="PL2")
s4 = spn.Sum((pl2, 0), n2, name="S4")
s5 = spn.Sum(pl2, name="S5")
s6 = spn.Sum((pl2, [1, 3]), name="S6")
s6.generate_latent_indicators()
pl3 = spn.ProductsLayer(s4, (n1, 2), num_or_size_prods=1, name="PL3")
pl4 = spn.ProductsLayer(s4,
s5,
s6,
s4,
s5,
s6,
num_or_size_prods=2,
name="PL4")
# Test
self.assertListEqual(v12.get_scope(), [
spn.Scope(v12, 0),
spn.Scope(v12, 0),
spn.Scope(v12, 0),
spn.Scope(v12, 0),
spn.Scope(v12, 1),
spn.Scope(v12, 1),
spn.Scope(v12, 1),
spn.Scope(v12, 1)
])
self.assertListEqual(
v34.get_scope(),
[spn.Scope(v34, 0), spn.Scope(v34, 1)])
self.assertListEqual(
s1.get_scope(),
[spn.Scope(v12, 0) | spn.Scope(s1.latent_indicators.node, 0)])
self.assertListEqual(s2.get_scope(), [spn.Scope(v12, 1)])
self.assertListEqual(pl1.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v12, 1)
| spn.Scope(v12, 1),
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0),
spn.Scope(v34, 1)
])
self.assertListEqual(n1.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(s1.latent_indicators.node, 0),
spn.Scope(v12, 1),
spn.Scope(v34, 1)
])
self.assertListEqual(n2.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v12, 1)
| spn.Scope(v12, 1),
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
])
self.assertListEqual(s3.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s3.latent_indicators.node, 0)
])
self.assertListEqual(pl2.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1)
| spn.Scope(s1.latent_indicators.node, 0),
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1),
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s3.latent_indicators.node, 0),
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s3.latent_indicators.node, 0)
])
self.assertListEqual(s4.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1)
| spn.Scope(s1.latent_indicators.node, 0) | spn.Scope(v34, 0)
])
self.assertListEqual(s5.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1)
| spn.Scope(s1.latent_indicators.node, 0) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s3.latent_indicators.node, 0)
])
self.assertListEqual(s6.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s3.latent_indicators.node, 0)
| spn.Scope(s6.latent_indicators.node, 0)
])
self.assertListEqual(pl3.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1)
| spn.Scope(s1.latent_indicators.node, 0) | spn.Scope(v34, 0)
| spn.Scope(v34, 1)
])
self.assertListEqual(pl4.get_scope(), [
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s1.latent_indicators.node, 0)
| spn.Scope(s3.latent_indicators.node, 0)
| spn.Scope(s6.latent_indicators.node, 0),
spn.Scope(v12, 0) | spn.Scope(v12, 1) | spn.Scope(v34, 0)
| spn.Scope(v34, 1) | spn.Scope(s1.latent_indicators.node, 0)
| spn.Scope(s3.latent_indicators.node, 0)
| spn.Scope(s6.latent_indicators.node, 0)
])
def test_compute_valid(self):
"""Calculating validity of ProductsLayer"""
v12 = spn.IndicatorLeaf(num_vars=2, num_vals=3)
v345 = spn.IndicatorLeaf(num_vars=3, num_vals=3)
v678 = spn.RawLeaf(num_vars=3)
v910 = spn.RawLeaf(num_vars=2)
p1 = spn.ProductsLayer((v12, [0, 3]), (v345, [1, 4, 7]),
(v678, [0, 1, 2]), (v910, [0]), (v910, 1),
num_or_size_prods=1)
p2 = spn.ProductsLayer((v12, [0, 3]), (v345, [1, 4, 7]),
v678,
v910,
num_or_size_prods=10)
p3 = spn.ProductsLayer(
(v345, 0), (v345, 3), (v345, 6), (v345, 0), (v345, 3), (v345, 7),
(v345, 0), (v345, 3), (v345, 8), (v345, 0), (v345, 4), (v345, 6),
(v345, 0), (v345, 4), (v345, 7), (v345, 0), (v345, 4), (v345, 8),
(v345, 0), (v345, 5), (v345, 6), (v345, 0), (v345, 5), (v345, 7),
(v345, 0), (v345, 5), (v345, 8), (v345, 1), (v345, 3), (v345, 6),
(v345, 1), (v345, 3), (v345, 7), (v345, 1), (v345, 3), (v345, 8),
(v345, 1), (v345, 4), (v345, 6), (v345, 1), (v345, 4), (v345, 7),
(v345, 1), (v345, 4), (v345, 8), (v345, 1), (v345, 5), (v345, 6),
(v345, 1), (v345, 5), (v345, 7), (v345, 1), (v345, 5), (v345, 8),
(v345, 2), (v345, 3), (v345, 6), (v345, 2), (v345, 3), (v345, 7),
(v345, 2), (v345, 3), (v345, 8), (v345, 2), (v345, 4), (v345, 6),
(v345, 2), (v345, 4), (v345, 7), (v345, 2), (v345, 4), (v345, 8),
(v345, 2), (v345, 5), (v345, 6), (v345, 2), (v345, 5), (v345, 7),
(v345, 2), (v345, 5), (v345, 8),
num_or_size_prods=27)
p4 = spn.ProductsLayer((v12, 0), (v12, [3]), (v345, [1, 4, 6]),
v678, (v910, 1),
num_or_size_prods=3)
p5 = spn.ProductsLayer((v12, 0), (v345, 4), (v678, [0, 1]), (v12, 0),
(v345, 4), (v678, [0, 1]),
num_or_size_prods=4)
p6 = spn.ProductsLayer((v678, 1), num_or_size_prods=[1])
p7 = spn.ProductsLayer((v12, 1), (v12, [5]), (v345, [0, 4, 8]),
v678,
v910,
num_or_size_prods=[10])
p8 = spn.ProductsLayer(
(v12, 1), (v12, [5]), (v345, [0, 4, 8]),
v678,
v910,
num_or_size_prods=[1, 1, 1, 1, 1, 1, 1, 1, 1, 1])
p9 = spn.ProductsLayer((v12, [1, 4]), (v345, [2, 5, 8]),
v678,
v910,
num_or_size_prods=[2, 2, 2, 2, 2])
p10 = spn.ProductsLayer((v12, [0]), (v12, 5), (v345, [0, 3, 6]),
v678, (v910, [1, 0]),
num_or_size_prods=[2, 3, 3, 2])
p11 = spn.ProductsLayer(v678, (v12, [2, 5]), (v345, [0, 3, 6]),
v910,
num_or_size_prods=[1, 2, 3, 4])
p12 = spn.ProductsLayer((v12, [0, 3, 4, 5]), (v345, [0, 3, 6, 7, 8]),
v678,
v910,
num_or_size_prods=[2, 1, 4, 1, 6])
self.assertTrue(p1.is_valid())
self.assertTrue(p2.is_valid())
self.assertTrue(p3.is_valid())
self.assertTrue(p4.is_valid())
self.assertTrue(p5.is_valid())
self.assertTrue(p6.is_valid())
self.assertTrue(p7.is_valid())
self.assertTrue(p8.is_valid())
self.assertTrue(p9.is_valid())
self.assertTrue(p10.is_valid())
self.assertTrue(p11.is_valid())
self.assertTrue(p12.is_valid())
p13 = spn.ProductsLayer(v12, num_or_size_prods=1)
p14 = spn.ProductsLayer((v910, 1), (v910, 1), num_or_size_prods=1)
p15 = spn.ProductsLayer((v12, 0), (v345, 1), (v345, 3), (v12, 0),
(v345, 1), (v345, 2),
num_or_size_prods=2)
p16 = spn.ProductsLayer((v345, [3, 4]), num_or_size_prods=[2])
p17 = spn.ProductsLayer((v910, 1), (v910, 1), num_or_size_prods=[2])
p18 = spn.ProductsLayer((v12, [1, 3, 4]), (v345, [3, 4]),
num_or_size_prods=[2, 3])
self.assertFalse(p13.is_valid())
self.assertFalse(p14.is_valid())
self.assertFalse(p15.is_valid())
self.assertFalse(p16.is_valid())
self.assertFalse(p17.is_valid())
self.assertFalse(p18.is_valid())
def test(input_sizes,
num_or_size_prods,
batch_size=7,
single_input=False,
indexed=True):
self.tearDown()
with self.subTest(input_sizes=input_sizes,
batch_size=batch_size,
num_or_size_prods=num_or_size_prods,
single_input=single_input,
indexed=indexed):
# Calculate num-prods or prod-input-size, based on type and
# values of the parameter 'num_or_size_prods'
if isinstance(num_or_size_prods, int):
num_prods = num_or_size_prods
prod_input_sizes = [sum(input_sizes) // num_prods
] * num_prods
elif isinstance(num_or_size_prods, list):
num_prods = len(num_or_size_prods)
prod_input_sizes = num_or_size_prods
max_prod_input_sizes = max(prod_input_sizes)
inputs = []
indices = []
inputs_feed = collections.OrderedDict()
if single_input:
# Size of the only input variable - would be = size of the
# biggest product modelled
num_vars = max(prod_input_sizes)
# Create a RawLeaf input node
input_var = spn.RawLeaf(num_vars=num_vars)
inputs = list(
itertools.repeat(input_var, len(prod_input_sizes)))
# Generate random input-indices for each product modelled,
# based on product-input-sizes
indices = [
random.sample(range(num_vars), k=p_size)
for p_size in prod_input_sizes
]
# Generate random values for the input feed
inputs_feed[input_var] = np.random.random(
(batch_size, num_vars))
else:
# Create ContVar inputs of random size each
for i_size in input_sizes:
if indexed:
# Choose a random size for an input node
num_vars = np.random.randint(low=i_size,
high=i_size + 5)
else:
num_vars = i_size
# Create a RawLeaf input node
inputs.append(spn.RawLeaf(num_vars=num_vars))
if indexed:
# Generate random indices for the created input
indices.append(
random.sample(range(num_vars), k=i_size))
else:
indices.append(None)
# Generate random values for the inputs feed
inputs_feed[inputs[-1]] = \
np.random.random((batch_size, num_vars))
# Create the ProductsLayer node
p = spn.ProductsLayer(*list(zip(inputs, indices)),
num_or_size_prods=num_or_size_prods)
# Create value ops
op_marg = p.get_value(spn.InferenceType.MARGINAL)
op_log_marg = p.get_log_value(spn.InferenceType.MARGINAL)
op_mpe = p.get_value(spn.InferenceType.MPE)
op_log_mpe = p.get_log_value(spn.InferenceType.MPE)
# Create an empty outputs ordered-dictionary
outputs = collections.OrderedDict()
# Create a session and execute the generated op
with self.test_session() as sess:
outputs[op_marg] = sess.run(op_marg, feed_dict=inputs_feed)
outputs[op_log_marg] = sess.run(tf.exp(op_log_marg),
feed_dict=inputs_feed)
outputs[op_mpe] = sess.run(op_mpe, feed_dict=inputs_feed)
outputs[op_log_marg] = sess.run(tf.exp(op_log_mpe),
feed_dict=inputs_feed)
# Calculate true-output
# Concat all inputs_feed into a single matrix
if single_input:
concatenated_inputs = np.hstack(
[inputs_feed[input_var][:, ind] for ind in indices])
else:
concatenated_inputs = np.hstack([
(i_feed[:, ind] if indexed else i_feed)
for (key,
i_feed), ind in zip(inputs_feed.items(), indices)
])
# Split concatenated-inputs based on product-size
inputs_per_prod = np.hsplit(concatenated_inputs,
np.cumsum(prod_input_sizes))
del inputs_per_prod[-1]
# Pad those inputs, with value '1.', whose size < max-prod-input-size
padded_values = \
np.stack([np.pad(inp, ((0, 0), (0, max_prod_input_sizes-inp.shape[1])),
'constant', constant_values=(1.0)) for
inp in inputs_per_prod])
# Reduce-prod and transpose to calculate true-output
true_output = np.transpose(np.prod(padded_values, axis=-1))
# Assert outputs with true_outputs
for _, out in outputs.items():
np.testing.assert_array_almost_equal(
out,
np.array(true_output,
dtype=spn.conf.dtype.as_numpy_dtype()))
