def test_signed_add():
A = exprvars('A', 8)
B = exprvars('B', 8)
for adder in (rca, ksa, bka):
S, C = adder(A, B)
# 0 + 0 = 0
assert sadd(S, A, B, 0, 0) == 0
# -64 + -64 = -128
assert sadd(S, A, B, -64, -64) == -128
# -1 + 1 = 0
assert sadd(S, A, B, -1, 1) == 0
# -64 + 64 = 0
assert sadd(S, A, B, -64, 64) == 0
# signed random vectors
for i in range(NVECS):
ra = random.randint(-2**6, 2**6-1) # -64..63
rb = random.randint(-2**6, 2**6) # -64..64
assert sadd(S, A, B, ra, rb) == ra + rb
# 64 + 64, overflow
R = C.vrestrict({A: int2exprs(64, 8), B: int2exprs(64, 8)})
assert R[7] != R[6]
# -65 + -64, overflow
R = C.vrestrict({A: int2exprs(-65, 8), B: int2exprs(-64, 8)})
assert R[7] != R[6]
def test_signed_add():
A = exprvars('A', 8)
B = exprvars('B', 8)
for adder in (rca, ksa, bka):
S, C = adder(A, B)
# 0 + 0 = 0
assert sadd(S, A, B, 0, 0) == 0
# -64 + -64 = -128
assert sadd(S, A, B, -64, -64) == -128
# -1 + 1 = 0
assert sadd(S, A, B, -1, 1) == 0
# -64 + 64 = 0
assert sadd(S, A, B, -64, 64) == 0
# signed random vectors
for i in range(NVECS):
ra = random.randint(-2**6, 2**6 - 1) # -64..63
rb = random.randint(-2**6, 2**6) # -64..64
assert sadd(S, A, B, ra, rb) == ra + rb
# 64 + 64, overflow
R = C.vrestrict({A: int2exprs(64, 8), B: int2exprs(64, 8)})
assert R[7] != R[6]
# -65 + -64, overflow
R = C.vrestrict({A: int2exprs(-65, 8), B: int2exprs(-64, 8)})
assert R[7] != R[6]
def test_vars():
xs = bddvars('x', 4)
assert all(x.name == 'x' and x.indices == (i, ) for i, x in enumerate(xs))
xs = exprvars('x', 4)
assert all(x.name == 'x' and x.indices == (i, ) for i, x in enumerate(xs))
xs = ttvars('x', 4)
assert all(x.name == 'x' and x.indices == (i, ) for i, x in enumerate(xs))
def __init__(self):
self.X = exprvars('x', (1, 10), (1, 10), (1, 10))
V = And(*[
And(*[
OneHot(*[self.X[r, c, v] for v in range(1, 10)])
for c in range(1, 10)
]) for r in range(1, 10)
])
R = And(*[
And(*[
OneHot(*[self.X[r, c, v] for c in range(1, 10)])
for v in range(1, 10)
]) for r in range(1, 10)
])
C = And(*[
And(*[
OneHot(*[self.X[r, c, v] for r in range(1, 10)])
for v in range(1, 10)
]) for c in range(1, 10)
])
B = And(*[
And(*[
OneHot(*[
self.X[3 * br + r, 3 * bc + c, v] for r in range(1, 4)
for c in range(1, 4)
]) for v in range(1, 10)
]) for br in range(3) for bc in range(3)
])
self.litmap, self.S = expr2dimacscnf(And(V, R, C, B))
def test_espresso():
assert espresso.FTYPE == 1
assert espresso.DTYPE == 2
assert espresso.RTYPE == 4
A = exprvars('a', 16)
B = exprvars('b', 16)
S, C = ripple_carry_add(A, B)
s0, s1, s2, s3 = espresso_exprs(S[0].to_dnf(), S[1].to_dnf(),
S[2].to_dnf(), S[3].to_dnf())
assert s0.equivalent(S[0])
assert s1.equivalent(S[1])
assert s2.equivalent(S[2])
assert s3.equivalent(S[3])
X = exprvars('x', 4)
f1 = truthtable(X, "0000011111------")
f2 = truthtable(X, "0001111100------")
f1m, f2m = espresso_tts(f1, f2)
truthtable2expr(f1).equivalent(f1m)
truthtable2expr(f2).equivalent(f2m)
def test_unsigned_add():
N = 9
A = exprvars('A', N)
B = exprvars('B', N)
for adder in (rca, ksa, bka):
S, C = adder(A, B)
S = fcat(S, C[-1])
# 0 + 0 = 0
assert uadd(S, A, B, 0, 0) == 0
# 255 + 255 = 510
assert uadd(S, A, B, 2**N - 1, 2**N - 1) == (2**(N + 1) - 2)
# 255 + 1 = 256
assert uadd(S, A, B, 2**N - 1, 1) == 2**N
# unsigned random vectors
for i in range(NVECS):
ra = random.randint(0, 2**N - 1)
rb = random.randint(0, 2**N - 1)
assert uadd(S, A, B, ra, rb) == ra + rb
def test_unsigned_add():
N = 9
A = exprvars('A', N)
B = exprvars('B', N)
for adder in (rca, ksa, bka):
S, C = adder(A, B)
S = fcat(S, C[-1])
# 0 + 0 = 0
assert uadd(S, A, B, 0, 0) == 0
# 255 + 255 = 510
assert uadd(S, A, B, 2**N-1, 2**N-1) == (2**(N+1)-2)
# 255 + 1 = 256
assert uadd(S, A, B, 2**N-1, 1) == 2**N
# unsigned random vectors
for i in range(NVECS):
ra = random.randint(0, 2**N-1)
rb = random.randint(0, 2**N-1)
assert uadd(S, A, B, ra, rb) == ra + rb
def bitvec(name, *dims):
"""Return a new array of given dimensions, filled with Expressions.
Parameters
----------
name : str
dims : (int or (int, int))
An int N means a slice from [0:N]
A tuple (M, N) means a slice from [M:N]
"""
if dims:
return bfarray.exprvars(name, *dims)
else:
return expr.exprvar(name)
def bitvec(name, *slices):
"""Return a BitVector with an arbitrary number of slices.
Parameters
----------
name : str
slices : (int or (int, int))
An int N means a slice from [0:N]
A tuple (M, N) means a slice from [M:N]
"""
if slices:
return bfarray.exprvars(name, *slices)
else:
return expr.exprvar(name)
def __init__(self):
self.X = exprvars('x', (1, 10), (1, 10), (1, 10))
V = And(*[And(*[OneHot(*[self.X[r, c, v] for v in range(1, 10)])
for c in range(1, 10)])
for r in range(1, 10)])
R = And(*[And(*[OneHot(*[self.X[r, c, v] for c in range(1, 10)])
for v in range(1, 10)])
for r in range(1, 10)])
C = And(*[And(*[OneHot(*[self.X[r, c, v] for r in range(1, 10)])
for v in range(1, 10)])
for c in range(1, 10)])
B = And(*[And(*[OneHot(*[self.X[3*br+r, 3*bc+c, v]
for r in range(1, 4)
for c in range(1, 4)])
for v in range(1, 10)])
for br in range(3) for bc in range(3)])
self.litmap, self.S = expr2dimacscnf(And(V, R, C, B))
from pyeda.boolalg.expr import (
Zero, One,
exprvar, expr,
#expr2dimacscnf, expr2dimacssat,
Expression,
Not, Or, And, Xor, Equal, Implies, ITE,
Nor, Nand, Xnor, Unequal,
OneHot0, OneHot, Majority, AchillesHeel, Mux,
)
# Common variables
a, b, c, d, e, p, q, s, w, x, y, z = map(exprvar, 'abcdepqswxyz')
d1, d0 = map(exprvar, ('d1', 'd0'))
xs = exprvars('x', 16)
ys = exprvars('y', 16, 16, 16)
def test_exprnode_constants():
"""Test exprnode constants"""
assert exprnode.ZERO == 0x0
assert exprnode.ONE == 0x1
assert exprnode.COMP == 0x4
assert exprnode.VAR == 0x5
assert exprnode.OP_OR == 0x8
assert exprnode.OP_AND == 0x9
assert exprnode.OP_XOR == 0xA
assert exprnode.OP_EQ == 0xB
def test_farray():
# expected shape volume to match items
assert_raises(ValueError, farray, [X[0], X[1]], shape=((0, 42), ))
# could not determine ftype parameter
assert_raises(ValueError, farray, [])
# expected ftype to be a type
assert_raises(TypeError, farray, [X[0], X[1]], ftype=42)
# expected ftype to match items
assert_raises(ValueError, farray, [X[0], X[1]], ftype=BinaryDecisionDiagram)
# expected ftype to be a property subclass of Function
assert_raises(TypeError, farray, [], ftype=int)
# expected a sequence of Function
assert_raises(TypeError, farray, 42)
assert_raises(TypeError, farray, [1, 2, 3, 4])
# expected uniform dimensions
assert_raises(ValueError, farray, [[a, b], [w, x, y, z], 42])
assert_raises(ValueError, farray, [[a, b], [w, x, y, z]])
# expected uniform types
assert_raises(ValueError, farray, [[a, b], [c, bddvar('d')]])
assert_raises(ValueError, farray, [[a, b], [bddvar('c'), bddvar('d')]])
# _check_shape errors
assert_raises(ValueError, farray, [a, b, c, d], shape=((-1, 3), ))
assert_raises(ValueError, farray, [a, b, c, d], shape=((3, -1), ))
assert_raises(ValueError, farray, [a, b, c, d], shape=((5, 1), ))
assert_raises(TypeError, farray, [a, b, c, d], shape=(('foo', 'bar'), ))
assert_raises(TypeError, farray, [a, b, c, d], shape=42)
temp = farray([[a, b], [c, d]])
assert str(temp) == """\
farray([[a, b],
[c, d]])\
"""
# __str__
Z = exprvars('z', 2, 2, 2)
assert str(Z) == """\
farray([[[z[0,0,0], z[0,0,1]],
[z[0,1,0], z[0,1,1]]],
[[z[1,0,0], z[1,0,1]],
[z[1,1,0], z[1,1,1]]]])\
"""
assert str(farray([], ftype=Expression)) == "farray([])"
# __getitem__
# expected <= M slice dimensions, got N
assert_raises(ValueError, X.__getitem__, (2, 2))
sel = exprvars('s', 2)
assert str(X[sel]) == "Or(And(~s[0], ~s[1], x[0]), And(s[0], ~s[1], x[1]), And(~s[0], s[1], x[2]), And(s[0], s[1], x[3]))"
assert str(X[:2][sel[0]]) == "Or(And(~s[0], x[0]), And(s[0], x[1]))"
# expected clog2(N) bits
assert_raises(ValueError, X.__getitem__, sel[0])
# slice step not supported
assert_raises(ValueError, X.__getitem__, slice(None, None, 2))
# type error
assert_raises(TypeError, X.__getitem__, 'foo')
# norm_index
assert X[-1] is X[3]
assert_raises(IndexError, X.__getitem__, 42)
# norm_indices
assert X[-3:-1]._items == [X[-3], X[-2]]
assert not X[-8:-10]._items
assert not X[-10:-8]._items
assert not X[8:10]._items
assert not X[10:8]._items
assert not X[3:1]._items
# __setitem__
Z = exprzeros(4, 4)
Z[0,0] = X[0]
assert Z._items[0] is X[0]
# expected item to be a Function
assert_raises(TypeError, Z.__setitem__, (0, 0), 42)
Z[0,:] = X[:4]
assert Z._items[0:4] == [X[0], X[1], X[2], X[3]]
# expected item to be an farray
assert_raises(TypeError, Z.__setitem__, (0, slice(None, None, None)), 42)
# expected item.size = ...
assert_raises(ValueError, Z.__setitem__, ..., X[:2])
# slice step not supported
assert_raises(ValueError, X.__setitem__, slice(None, None, 2), 42)
# type error
assert_raises(TypeError, X.__setitem__, 'foo', 42)
# __add__
assert (0 + X)._items[0].is_zero()
assert (X + 0)._items[4].is_zero()
assert (Y[0] + X)._items[0] is Y[0]
assert (X + Y[0])._items[4] is Y[0]
assert (X[:2] + Y[2:])._items == [X[0], X[1], Y[2], Y[3]]
# expected Function or farray
assert_raises(TypeError, X.__add__, 42)
assert_raises(TypeError, X.__radd__, 42)
A = exprvars('a', 2, 5, 6)
B = exprvars('b', 2, 5, 6)
C = exprvars('c', (1, 3), 5, 6)
# regular MDA will retain shape
assert (A+B).shape == ((0, 4), (0, 5), (0, 6))
# irregular MDA will not
assert (A+C).shape == ((0, 4*5*6), )
# regular MDA will retain shape
assert (A*2).shape == ((0, 4), (0, 5), (0, 6))
# irregular MDA will not
assert (C*2).shape == ((0, 4*5*6), )
# __mul__
# expected multiplier to be an int
assert_raises(TypeError, X.__mul__, 'foo')
# expected multiplier to be non-negative
assert_raises(ValueError, X.__mul__, -2)
assert (X[:2] * 2)._items == [X[0], X[1], X[0], X[1]]
assert (2 * X[:2])._items == [X[0], X[1], X[0], X[1]]
# offsets
Z = exprzeros((1, 5), (17, 21))
assert Z.offsets == (1, 17)
# reshape
assert Z.reshape(4, 4).shape == ((0, 4), (0, 4))
# expected shape with equal volume
assert_raises(ValueError, Z.reshape, 42, 42)
# restrict
assert str(X.vrestrict({X: '0101'})) == "farray([0, 1, 0, 1])"
# compose
assert X.compose({X[0]: Y[0]})._items[0] == Y[0]
# to_uint / to_int
assert uint2exprs(42).to_uint() == 42
assert uint2exprs(42, 8).to_uint() == 42
# expected all functions to be a constant (0 or 1) form
assert_raises(ValueError, X.to_uint)
# expected num >= 0
assert_raises(ValueError, uint2exprs, -1)
# overflow
assert_raises(ValueError, uint2exprs, 42, 2)
assert_raises(ValueError, int2exprs, 42, 2)
assert int2exprs(-42).to_int() == -42
assert int2exprs(-42, 8).to_int() == -42
assert int2exprs(42).to_int() == 42
assert int2exprs(42, 8).to_int() == 42
# zext, sext
assert X.zext(1)[4].is_zero()
assert X.sext(1)[4] is X[3]
# __invert__, __or__, __and__, __xor__
assert str(~X) == "farray([~x[0], ~x[1], ~x[2], ~x[3]])"
assert str(X | Y) == "farray([Or(x[0], y[0]), Or(x[1], y[1]), Or(x[2], y[2]), Or(x[3], y[3])])"
assert str(X & Y) == "farray([And(x[0], y[0]), And(x[1], y[1]), And(x[2], y[2]), And(x[3], y[3])])"
assert str(X ^ Y) == "farray([Xor(x[0], y[0]), Xor(x[1], y[1]), Xor(x[2], y[2]), Xor(x[3], y[3])])"
# _op_shape
# expected farray input
assert_raises(TypeError, X.__or__, 42)
Z = exprvars('z', 2, 2)
assert str(X | Z) == "farray([Or(x[0], z[0,0]), Or(x[1], z[0,1]), Or(x[2], z[1,0]), Or(x[3], z[1,1])])"
Z = exprvars('z', 2, 3)
# expected operand sizes to match
assert_raises(ValueError, X.__or__, Z)
# lsh, rsh
assert str(X.lsh(0)) == "(farray([x[0], x[1], x[2], x[3]]), farray([]))"
assert str(X << 0) == "farray([x[0], x[1], x[2], x[3]])"
assert str(X.lsh(2)) == "(farray([0, 0, x[0], x[1]]), farray([x[2], x[3]]))"
assert str(X << 2) == "farray([0, 0, x[0], x[1]])"
assert str(X << (2, Y[:2])) == "farray([y[0], y[1], x[0], x[1]])"
assert str(X.rsh(0)) == "(farray([x[0], x[1], x[2], x[3]]), farray([]))"
assert str(X >> 0) == "farray([x[0], x[1], x[2], x[3]])"
assert str(X.rsh(2)) == "(farray([x[2], x[3], 0, 0]), farray([x[0], x[1]]))"
assert str(X >> 2) == "farray([x[2], x[3], 0, 0])"
assert str(X >> (2, Y[:2])) == "farray([x[2], x[3], y[0], y[1]])"
assert_raises(TypeError, X.__lshift__, 'foo')
assert_raises(ValueError, X.__lshift__, -1)
assert_raises(ValueError, X.__lshift__, (2, Y))
assert_raises(TypeError, X.__rshift__, 'foo')
assert_raises(ValueError, X.__rshift__, -1)
assert_raises(ValueError, X.__rshift__, (2, Y))
# arsh
assert str(X.arsh(0)) == "(farray([x[0], x[1], x[2], x[3]]), farray([]))"
assert str(X.arsh(2)) == "(farray([x[2], x[3], x[3], x[3]]), farray([x[0], x[1]]))"
assert_raises(ValueError, X.arsh, -1)
# unary ops
assert str(X.uor()) == "Or(x[0], x[1], x[2], x[3])"
assert str(X.unor()) == "Not(Or(x[0], x[1], x[2], x[3]))"
assert str(X.uand()) == "And(x[0], x[1], x[2], x[3])"
assert str(X.unand()) == "Not(And(x[0], x[1], x[2], x[3]))"
assert str(X.uxor()) == "Xor(x[0], x[1], x[2], x[3])"
assert str(X.uxnor()) == "Not(Xor(x[0], x[1], x[2], x[3]))"
# decode
assert str(farray([], ftype=Expression).decode()) == "farray([1])"
assert str(X[:2].decode()) == "farray([And(~x[0], ~x[1]), And(x[0], ~x[1]), And(~x[0], x[1]), And(x[0], x[1])])"
from pyeda.boolalg import boolfunc
from pyeda.boolalg.expr import (
exprvar, expr,
expr2dimacscnf, expr2dimacssat,
Expression,
Not, Or, And, Nor, Nand, Xor, Xnor, Equal, Unequal, Implies, ITE,
OneHot0, OneHot, Majority, AchillesHeel, Mux,
EXPRZERO, EXPRONE,
)
from pyeda.boolalg.bfarray import exprvars
from nose.tools import assert_raises
a, b, c, d, e, p, q, s = map(exprvar, 'abcdepqs')
X = exprvars('x', 16)
Y = exprvars('y', 16, 16, 16)
def test_misc():
f = a & b | a & c | b & c
assert f.smoothing(a).equivalent(b | c)
assert f.consensus(a).equivalent(b & c)
assert f.derivative(a).equivalent(b & ~c | ~b & c)
assert_raises(TypeError, expr2dimacscnf, 'foo')
def test_issue81():
# Or(x) = x
assert str(Or(Or(a, b))) == "Or(a, b)"
assert str(Or(And(a, b))) == "And(a, b)"
def test_bin2gray():
B = exprvars('B', 4)
G = bin2gray(B)
gnums = [G.vrestrict({B: uint2exprs(i, 4)}).to_uint() for i in range(16)]
assert gnums == [0, 1, 3, 2, 6, 7, 5, 4, 12, 13, 15, 14, 10, 11, 9, 8]
"""
from nose.tools import assert_raises
from pyeda.boolalg.bdd import bddvar, BinaryDecisionDiagram
from pyeda.boolalg.bfarray import (
exprzeros, exprvars,
fcat, farray,
uint2exprs, int2exprs,
)
from pyeda.boolalg.boolfunc import Function
from pyeda.boolalg.expr import exprvar, Expression
from pyeda.boolalg.table import ttvar
X = exprvars('x', 4)
Y = exprvars('y', 4)
a, b, c, d, w, x, y, z = map(exprvar, 'abcdwxyz')
def test_fcat():
# expected Function or farray
assert_raises(TypeError, fcat, X, Y, 0)
assert str(fcat(X[0], X[2:], Y[3], Y[:-2])) == "farray([x[0], x[2], x[3], y[3], y[0], y[1]])"
def test_farray():
# expected shape volume to match items
assert_raises(ValueError, farray, [X[0], X[1]], shape=((0, 42), ))
# could not determine ftype parameter
assert_raises(ValueError, farray, [])
# expected ftype to be a type
assert_raises(TypeError, farray, [X[0], X[1]], ftype=42)
def test_gray2bin():
G = exprvars('G', 4)
B = gray2bin(G)
gnums = [0, 1, 3, 2, 6, 7, 5, 4, 12, 13, 15, 14, 10, 11, 9, 8]
bnums = [B.vrestrict({G: uint2exprs(i, 4)}).to_uint() for i in gnums]
assert bnums == list(range(16))
Equal,
Unequal,
Implies,
ITE,
OneHot0,
OneHot,
Majority,
AchillesHeel,
Mux,
EXPRZERO,
EXPRONE,
)
a, b, c, d, e, p, q, s = map(exprvar, 'abcdepqs')
X = exprvars('x', 16)
Y = exprvars('y', 16, 16, 16)
def test_misc():
f = a & b | a & c | b & c
assert f.smoothing(a).equivalent(b | c)
assert f.consensus(a).equivalent(b & c)
assert f.derivative(a).equivalent(b & ~c | ~b & c)
assert_raises(TypeError, expr2dimacscnf, 'foo')
def test_issue81():
# Or(x) = x
def test_errors():
A = exprvars('A', 7)
B = exprvars('B', 9)
for adder in (rca, ksa, bka):
assert_raises(ValueError, adder, A, B)
def test_errors():
A = exprvars('A', 7)
B = exprvars('B', 9)
for adder in (rca, ksa, bka):
assert_raises(ValueError, adder, A, B)
def test_gray2bin():
G = exprvars('G', 4)
B = gray2bin(G)
gnums = [0, 1, 3, 2, 6, 7, 5, 4, 12, 13, 15, 14, 10, 11, 9, 8]
bnums = [B.vrestrict({G: uint2exprs(i, 4)}).to_uint() for i in gnums]
assert bnums == list(range(16))
def test_bin2gray():
B = exprvars('B', 4)
G = bin2gray(B)
gnums = [G.vrestrict({B: uint2exprs(i, 4)}).to_uint() for i in range(16)]
assert gnums == [0, 1, 3, 2, 6, 7, 5, 4, 12, 13, 15, 14, 10, 11, 9, 8]
COLORS = ["red", "blue", "green", "yellow", "white", "gray", "cyan", "orange", "purple", "magenta"]
if "SHAPEWORLD_N_COLORS" in os.environ:
N_COLORS = int(os.environ['SHAPEWORLD_N_COLORS'])
assert N_COLORS < len(COLORS), f"{N_COLORS} requested, have {len(COLORS)}"
else:
N_COLORS = 6 # original setting
COLORS = COLORS[:N_COLORS]
AGGDRAW_COLORS = {
c: color_factory(hsv_c, {"brush": aggdraw.Brush, "pen": aggdraw.Pen})
for c, hsv_c in COLORS2HSV.items()
}
COLOR_VARS = exprvars('c', len(COLORS))
C2V = dict(zip(COLORS, COLOR_VARS))
V2C = dict(zip(COLOR_VARS, COLORS))
ONEHOT_VAR = expr.OneHot(*COLOR_VARS)
# FIXME - this should be OOP?
def random(colors=None):
if colors is None:
return np.random.choice(COLORS)
else:
return np.random.choice(colors)
def new(existing_color, colors=None):
if colors is None:
from pyeda.boolalg.expr import (
Zero, One,
exprvar, expr,
#expr2dimacscnf, expr2dimacssat,
Expression,
Not, Or, And, Xor, Equal, Implies, ITE,
Nor, Nand, Xnor, Unequal,
OneHot0, OneHot, Majority, AchillesHeel, Mux,
)
# Common variables
a, b, c, d, e, p, q, s, w, x, y, z = map(exprvar, 'abcdepqswxyz')
d1, d0 = map(exprvar, ('d1', 'd0'))
xs = exprvars('x', 16)
ys = exprvars('y', 16, 16, 16)
def test_exprnode_constants():
"""Test exprnode constants"""
assert exprnode.ZERO == 0x0
assert exprnode.ONE == 0x1
assert exprnode.COMP == 0x4
assert exprnode.VAR == 0x5
assert exprnode.OP_OR == 0x8
assert exprnode.OP_AND == 0x9
assert exprnode.OP_XOR == 0xA
assert exprnode.OP_EQ == 0xB
SHAPE_IMPLS = {
"circle": Circle,
"ellipse": Ellipse,
"square": Square,
"rectangle": Rectangle,
"triangle": Triangle,
"semicircle": Semicircle,
"pentagon": Pentagon,
"cross": Cross,
"hexagon": Hexagon,
}
SHAPES = [
"circle", "square", "rectangle", "ellipse", "triangle", "semicircle",
"pentagon", "cross", "hexagon"
]
if "SHAPEWORLD_N_SHAPES" in os.environ:
N_SHAPES = int(os.environ['SHAPEWORLD_N_SHAPES'])
assert N_SHAPES < len(SHAPES), f"{N_SHAPES} requested, have {len(SHAPES)}"
else:
N_SHAPES = 5 # original setting
SHAPES = SHAPES[:N_SHAPES]
SHAPE_VARS = exprvars('s', len(SHAPES))
S2V = dict(zip(SHAPES, SHAPE_VARS))
V2S = dict(zip(SHAPE_VARS, SHAPES))
ONEHOT_VAR = expr.OneHot(*SHAPE_VARS)