passagemath-brial 10.6.33__cp311-cp311-manylinux_2_24_aarch64.manylinux_2_28_aarch64.whl

sage/rings/polynomial/pbori/blocks.py

@@ -0,0 +1,443 @@
+# sage_setup: distribution = sagemath-brial
+# sage.doctest: needs sage.rings.polynomial.pbori
+import sys
+from itertools import chain, islice
+
+from .pbori import VariableBlock
+from .PyPolyBoRi import (Ring, Polynomial, VariableFactory, Variable)
+
+
+class Block:
+    r"""
+    The block class represents a block of variables
+    <var_name>(start_index,...,start_index+size-1), it is the preferred
+    block type for simple one-dimensional variable sets
+    """
+    def __init__(self, var_name, size, start_index=0, reverse=False):
+        indices = range(start_index, start_index + size)
+        if reverse:
+            indices = reversed(indices)
+        self.names = [var_name + "(" + str(i) + ")" for i in indices]
+        self.var_name = var_name
+        self.start_index = start_index
+        self.reverse = reverse
+        self.size = size
+
+    def __iter__(self):
+        return iter(self.names)
+
+    def __getitem__(self, i):
+        return self.names[i]
+
+    def __len__(self):
+        return self.size
+
+    def register(self, start, context):
+        ring_context = context
+        while isinstance(ring_context, PrefixedDictProxy):
+            ring_context = ring_context.wrapped
+        ring = ring_context['r']
+
+        var_func = VariableBlock(self.size, self.start_index, start, self.
+                                 reverse, ring)
+        var_func.__name__ = self.var_name
+        context[self.var_name] = var_func
+
+
+class AlternatingBlock:
+    r"""
+    The Alternating Block class is used for doing tricky variable
+    schemes,where base names vary, e.g.
+    a(0),b(0),a(1),b(1),a(2),b(2)
+    """
+    def __init__(self, var_names, size_per_variable, start_index=0,
+                 reverse=False):
+        self.var_names = var_names
+        self.size_per_variable = size_per_variable
+        self.reverse = reverse
+        indices = range(start_index, start_index + size_per_variable)
+
+        if reverse:
+            indices = reversed(indices)
+        names = [f"{n}({str(i)})" for i in indices for n in var_names]
+        self.indices = indices
+        self.index2pos = {v: k for k, v in enumerate(indices)}
+        self.names = names
+
+    def __len__(self):
+        return self.size_per_variable * len(self.var_names)
+
+    def __iter__(self):
+        return iter(self.names)
+
+    def __getitem__(self, i):
+        return self.names[i]
+
+    def register(self, start, context):
+        def gen_var_func(var_pos):
+
+            class var_factory:
+                def __init__(self, ring, index2pos, size):
+                    self.ring = ring
+                    self.index2pos = index2pos
+                    self.size = size
+
+                def __call__(self, idx):
+                    return self.ring.variable(self.index2pos[idx] * self.size +
+                                              var_pos + start)
+            ring_context = context
+            while isinstance(ring_context, PrefixedDictProxy):
+                ring_context = ring_context.wrapped
+            ring = ring_context['r']
+
+            return var_factory(ring, self.index2pos, len(self.var_names))
+
+        for (var_pos, n) in enumerate(self.var_names):
+            var_func = gen_var_func(var_pos)
+            var_func.__name__ = n
+            context[n] = var_func
+
+
+def shift(f, i):
+    def g(j):
+        return f(i + j)
+    g.__name__ = f.__name__
+    return g
+
+
+class AdderBlock(AlternatingBlock):
+    def __init__(self, adder_bits, sums='s', carries='c', input1='a',
+                 input2='b', start_index=0):
+        AlternatingBlock.__init__(self, (sums, carries, input1, input2),
+                                  adder_bits, start_index=start_index,
+                                  reverse=True)
+        self.input1 = input1
+        self.input2 = input2
+        self.sums = sums
+        self.carries = carries
+        self.start_index = start_index
+        self.adder_bits = adder_bits
+
+    def register(self, start, context):
+        super().register(start, context)
+        a = context[self.input1]
+        b = context[self.input2]
+        self.s = shift(context[self.sums], self.start_index)
+        self.c = shift(context[self.carries], self.start_index)
+        a = shift(a, self.start_index)
+        b = shift(b, self.start_index)
+        carries = [Polynomial(a(0).ring().zero())]
+        last = carries[0]
+        for i in range(self.adder_bits):
+            c = 1 + (1 + a(i) * b(i)) * (1 + last * a(i)) * (1 + last * b(i))
+            carries.append(c)
+            last = c
+
+        self.add_results = [a(i) + b(i) + carries[i]
+                            for i in range(self.adder_bits)]
+        self.carries_polys = carries[1:]
+
+    # def s(i):
+    #     return self.add_results[i-self.start_index]
+    # def c(i):
+    #     return self.carries_polys[i-self.start_index]
+    # context[self.sums]=s
+    # context[self.carries]=c
+
+    def implement(self, equations):
+        for i in range(self.adder_bits):
+            equations.append(self.s(i) + self.add_results[i])
+            equations.append(self.c(i) + self.carries_polys[i])
+
+
+class HigherOrderBlock:
+    r"""
+    HigherOrderBlocks are multidimensional blocks of variables.
+
+    For each dimension a separate start_index and size can be specified.
+
+    var_name : variables will be called <var_name>(multiindex), where
+    multiindex is a tuple of the size <size_tuple>
+
+    size_tuple : specifies the sizes of the ranges of each component of the
+    multi-indices
+
+    start_index_tuple : the multi-indices will be of the form
+    start_index_tuple + a, where a is a multi-index with nonnegative components
+    """
+    def __init__(self, var_name, size_tuple, start_index_tuple=None,
+                 reverse=False):
+        if start_index_tuple is None:
+            start_index_tuple = len(size_tuple) * (0, )
+        cart = [()]
+        assert len(size_tuple) == len(start_index_tuple)
+        outer_indices = reversed(range(len(size_tuple)))
+        for i in outer_indices:
+            s_i = start_index_tuple[i]
+            s = size_tuple[i]
+            cart = [(j, ) + c for j in range(s_i, s_i + s) for c in cart]
+        if reverse:
+            cart.reverse()
+        self.cart = cart
+        self.cart2index = {v: k for k, v in enumerate(cart)}
+        self.var_name = var_name
+        self.names = [var_name + str(c) for c in cart]
+
+    def __getitem__(self, i):
+        return self.names[i]
+
+    def __iter__(self):
+        return iter(self.names)
+
+    def __len__(self):
+        return len(self.names)
+
+    def register(self, start, context):
+        def var_func(*indices):
+            return Variable(self.cart2index[indices] + start)
+        var_func.__name__ = self.var_name
+        context[self.var_name] = var_func
+
+
+class InOutBlock:
+    def __init__(self, out_size, in_size, output='out', input='in',
+                 in_start_index=0, out_start_index=0,
+                 out_reverse=False, in_reverse=False):
+        self.output = Block(var_name=output, start_index=out_start_index,
+                            size=out_size, reverse=out_reverse)
+        self.input = Block(var_name=input, start_index=in_start_index,
+                           size=in_size, reverse=in_reverse)
+        self.out_start_index = out_start_index
+
+        self.in_start_index = in_start_index
+
+    def __iter__(self):
+        return chain(self.output, self.input)
+
+    def __getitem__(self, i):
+        if i < len(self.output):
+            return self.output[i]
+        return self.input[i - len(self.output)]
+
+    def __len__(self):
+        return len(self.output) + len(self.input)
+
+    def register(self, start, context):
+        self.output.register(start, context)
+        self.input.register(start + len(self.output), context)
+        self.out_vars = shift(context[self.output.var_name], self.
+                              out_start_index)
+        self.in_vars = shift(context[self.input.var_name], self.in_start_index)
+
+
+class MultiBlock:
+    def __init__(self, sizes=None, var_names=["v"],
+                 start_indices=[], reverses=None):
+        if reverses is None:
+            reverses = []
+        if sizes is None:
+            sizes = []
+        self.start_indices = start_indices + [0] * (len(var_names) -
+                                                    len(start_indices))
+        reverses += [False] * (len(var_names) - len(reverses))
+        sizes += [1] * (len(var_names) - len(sizes))
+
+        self.blocks = [Block(var_name=var_names[idx], size=sizes[idx],
+                             start_index=self.start_indices[idx],
+                             reverse=reverses[idx])
+                       for idx in range(len(var_names))]
+
+    def __iter__(self):
+        return chain(*self.blocks)
+
+    def __getitem__(self, i):
+        return next(islice(chain(*self.blocks), i, i + 1))
+        # sum([bl.names for bl in self.blocks])[i]
+
+    def __len__(self):
+        return sum(len(bl) for bl in self.blocks)
+
+    def register(self, start, context):
+        offset = 0
+        for bl in self.blocks:
+            bl.register(start + offset, context)
+            offset += len(bl)
+
+        self.vars = [shift(context[self.blocks[idx].var_name],
+                           self.start_indices[idx])
+                     for idx in range(len(self.blocks))]
+
+
+class PrefixedDictProxy:
+    """docstring for PrefixedDictProxy"""
+
+    def __init__(self, wrapped, prefix):
+        super().__init__()
+        self.wrapped = wrapped
+        self.prefix = prefix
+
+    def __getitem__(self, k):
+        try:
+            return self.wrapped[self.prefix + k]
+        except KeyError:
+            print(self.prefix, k, list(self.wrapped))
+            raise KeyError
+
+    def __setitem__(self, k, v):
+        self.wrapped[self.prefix + k] = v
+
+
+class MacroBlock:
+    def __init__(self, prefix):
+
+        self.prefix = prefix
+        self.blocks = []
+        self.combinations = []
+        self.connections = []
+
+    def declare(self, blocks):
+        self.blocks = blocks
+
+    def connect(self, combinations):
+        self.combinations = combinations
+
+    def __iter__(self):
+        return (self.prefix + "_" + n for n in chain(*self.blocks))
+
+    def __getitem__(self, i):
+        return self.prefix + "_" + next(islice(chain(*self.blocks), i, i + 1))
+
+    def __len__(self):
+        return sum(len(bl) for bl in self.blocks)
+
+    def resolve(self, localname):
+        return self.prefix + "_" + localname
+
+    def register(self, start, context):
+        context = PrefixedDictProxy(context, self.prefix + "_")
+        offset = 0
+        for bl in self.blocks:
+            bl.register(start + offset, context)
+            offset += len(bl)
+
+        for ((con1, indices1), (con2, indices2)) in self.combinations:
+            for idx in range(min(len(indices1), len(indices2))):
+                self.connections += [context[con1](indices1[idx]) + context[
+                    con2](indices2[idx])]
+
+    def implement(self, equations):
+        for bl in self.blocks:
+            if hasattr(bl, "implement"):
+                bl.implement(equations)
+
+        equations += self.connections
+
+
+class IfThen:
+    def __init__(self, ifpart, thenpart, supposed_to_be_valid=True):
+        self.ifpart = [Polynomial(p) for p in ifpart]
+        self.thenpart = [Polynomial(p) for p in thenpart]
+        self.supposedToBeValid = supposed_to_be_valid
+
+    def __str__(self):
+        return ("If(AND(" + ", ".join(f"{p} == 0" for p in self.ifpart) +
+                ")), THEN " + ", ".join(f"{p} == 0" for p in self.thenpart))
+
+
+def if_then(i, t, supposed_to_be_valid=True):
+    return IfThen(i, t, supposed_to_be_valid)
+
+
+def declare_ring(blocks, context=None):
+    r"""
+    Declare Ring is the preferred function to create a ring and declare a variable scheme,
+    the number of variables is automatically determined, usually you pass globals() as context
+    argument to store the ring and the variable mapping.
+
+    EXAMPLES::
+
+        sage: from sage.rings.polynomial.pbori import *
+        sage: declare_ring([Block("x",10),Block("y",5)],globals())
+        Boolean PolynomialRing in x0, x1, x2, x3, x4, x5, x6, x7, x8, x9, y0, y1, y2, y3, y4
+
+    gives  a ring with x(0..9),y(0..4) and registers the ring as r, and the variable
+    blocks x and y in the context dictionary globals(), which consists of the global
+    variables of the python module
+    """
+    if context is None:
+        context = sys.modules['__main__'].__dict__
+
+    def canonicalize(blocks):
+        for elt in blocks:
+            if isinstance(elt, str):
+                yield elt
+            else:
+                yield from elt
+
+    blocks = list(blocks)
+    n = 0
+
+    for b in blocks:
+        n = n + 1 if isinstance(b, str) else n + len(b)
+
+    r = Ring(n, names=canonicalize(blocks))
+
+    context["internalVariable"] = VariableFactory(r)
+#  context["Monomial"] = MonomialFactory(r)
+    context["r"] = r
+    declare_block_scheme(blocks, context)
+    return r
+
+
+def declare_block_scheme(blocks, context):
+    start = 0
+    block_starts = []
+    for b in blocks:
+        if start:
+            block_starts.append(start)
+        if isinstance(b, str):
+            context[b] = context["internalVariable"](start)
+            start = start + 1
+        else:
+            b.register(start, context)
+            start = start + len(b)
+    context["block_start_hints"] = block_starts
+    context["number_of_declared_vars"] = start
+
+
+def main():
+    r = Ring(1000)
+    ablock = AlternatingBlock(["a", "b", "c"], 100)
+    declare_block_scheme([ablock], globals())
+    for i in range(10):
+        print(r.variable(i))
+
+    print(list(ablock))
+    declare_block_scheme([Block(var_name="x", size=100),
+                          HigherOrderBlock("y", (3, 4, 11, 2)),
+                          AlternatingBlock(["a", "b", "c"], 100)],
+                         globals())
+    for i in range(10):
+        print(x(i))
+    print(y(0, 0, 0, 0))
+    print(y(0, 0, 0, 1))
+    print(y(0, 0, 1, 0))
+    print(y(0, 0, 1, 1))
+    print(a(0), a(1), a(2), b(0), b(1), c(0))
+    declare_block_scheme([Block(var_name="x", size=100, reverse=True),
+                          HigherOrderBlock("y", (3, 4, 11, 2), reverse=True),
+                          AlternatingBlock(["a", "b", "c"], 100, reverse=True)],
+                         globals())
+    for i in range(10):
+        print(x(i))
+    print(y(0, 0, 0, 0))
+    print(y(0, 0, 0, 1))
+    print(y(0, 0, 1, 0))
+    print(y(0, 0, 1, 1))
+    print(a(0), a(1), a(2), b(0), b(1), c(0))
+    declare_block_scheme(["a", "b", "c"], globals())
+    print(a, b, c)
+
+
+if __name__ == '__main__':
+    main()

sage/rings/polynomial/pbori/cnf.py

@@ -0,0 +1,241 @@
+# sage_setup: distribution = sagemath-brial
+# sage.doctest: needs sage.rings.polynomial.pbori
+from random import Random
+from sage.rings.polynomial.pbori.pbori import if_then_else as ite
+from .PyPolyBoRi import Polynomial
+from .statistics import used_vars_set
+
+
+class CNFEncoder:
+    def __init__(self, r, random_seed=16):
+        self.random_generator = Random(random_seed)
+        self.one_set = r.one().set()
+        self.empty_set = r.zero().set()
+        self.r = r
+
+    def zero_blocks(self, f):
+        r"""
+        Divide the zero set of f into blocks.
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: r = declare_ring(["x", "y", "z"], dict())
+            sage: from sage.rings.polynomial.pbori.cnf import CNFEncoder
+            sage: e = CNFEncoder(r)
+            sage: e.zero_blocks(r.variable(0)*r.variable(1)*r.variable(2))
+            [{z: 0}, {y: 0}, {x: 0}]
+        """
+        f = Polynomial(f)
+        variables = f.vars_as_monomial()
+
+        space = variables.divisors()
+        variables = list(variables.variables())
+        zeros = f.zeros_in(space)
+        rest = zeros
+        res = []
+
+        # inefficient compared to polynomials lex_lead
+        def choose(s):
+            indices = []
+            assert not s.empty()
+            nav = s.navigation()
+            while not nav.constant():
+                e = nav.else_branch()
+                t = nav.then_branch()
+                if e.constant() and not e.terminal_one():
+                    indices.append(nav.value())
+                    nav = t
+                else:
+                    if self.random_generator.randint(0, 1):
+                        indices.append(nav.value())
+                        nav = t
+
+                    else:
+                        nav = e
+            assert nav.terminal_one()
+            res = self.one_set
+            for i in reversed(indices):
+                res = ite(i, res, self.empty_set)
+            return next(iter(res))
+        while not rest.empty():
+            l = choose(rest)
+            l_variables = set(l.variables())
+
+            def get_val(var):
+                if var in l_variables:
+                    return 1
+                return 0
+            block_dict = {v: get_val(v) for v in variables}
+
+            l = l.set()
+            self.random_generator.shuffle(variables)
+            for v in variables:
+                candidate = l.change(v.index())
+                if candidate.diff(zeros).empty():
+                    l = l.union(candidate)
+                    del block_dict[v]
+            rest = rest.diff(l)
+            res.append(block_dict)
+        return res
+
+    def clauses(self, f):
+        r"""
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: r = declare_ring(["x", "y", "z"], dict())
+            sage: from sage.rings.polynomial.pbori.cnf import CNFEncoder
+            sage: e = CNFEncoder(r)
+            sage: sorted(e.clauses(r.variable(0)*r.variable(1)*r.variable(2)), key=lambda d: sorted(d.items()))
+            [{z: 0, y: 0, x: 0}]
+            sage: sorted(e.clauses(r.variable(Integer(1))+r.variable(Integer(0))), key=lambda d: sorted(d.items()))
+            [{y: 0, x: 1}, {y: 1, x: 0}]
+        """
+        # we form an expression for a var configuration *not* lying in the
+        # block it is evaluated to 0 by f, iff it is not lying in any zero
+        # block of f+1
+        return [{variable: 1 - value for variable, value in b.items()}
+                for b in self.zero_blocks(f + 1)]
+
+    def polynomial_clauses(self, f):
+        r"""
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: r = declare_ring(["x", "y", "z"], dict())
+            sage: from sage.rings.polynomial.pbori.cnf import CNFEncoder
+            sage: e = CNFEncoder(r)
+            sage: e.polynomial_clauses(r.variable(0)*r.variable(1)*r.variable(2))
+            [x*y*z]
+            sage: v = r.variable
+            sage: p = v(1)*v(2)+v(2)*v(0)+1
+            sage: groebner_basis([p], heuristic = False)==groebner_basis(e.polynomial_clauses(p), heuristic = False)
+            True
+        """
+
+        def product(l):
+            res = l[0]
+            for p in l[1:]:
+                res = res * p
+            # please care about the order of these multiplications for
+            # performance
+            return res
+        return [product([variable + value for (variable, value)
+                         in b.items()]) for b in self.clauses(f)]
+
+    def to_dimacs_index(self, v):
+        return v.index() + 1
+
+    def dimacs_encode_clause(self, c):
+        def get_sign(value):
+            if value == 1:
+                return 1
+            return -1
+
+        items = sorted(c.items(), reverse=True)
+        return " ".join([str(v) for v in [
+            get_sign(value) * self.to_dimacs_index(variable)
+            for (variable, value) in items] + [0]])
+
+    def dimacs_encode_polynomial(self, p):
+        r"""
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: d = {}
+            sage: r = declare_ring(["x", "y", "z"], d)
+            sage: from sage.rings.polynomial.pbori.cnf import CNFEncoder
+            sage: e = CNFEncoder(r)
+            sage: sorted(e.dimacs_encode_polynomial(d["x"]+d["y"]+d["z"]))
+            ['-1 -2 -3 0', '-1 2 3 0', '1 -2 3 0', '1 2 -3 0']
+        """
+        return [self.dimacs_encode_clause(c) for c in self.clauses(p)]
+
+    def dimacs_cnf(self, polynomial_system):
+        r"""
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: r = declare_ring(["x", "y", "z"], dict())
+            sage: from sage.rings.polynomial.pbori.cnf import CNFEncoder
+            sage: e = CNFEncoder(r)
+            sage: e.dimacs_cnf([r.variable(0)*r.variable(1)*r.variable(2)])
+            'c cnf generated by PolyBoRi\np cnf 3 1\n-1 -2 -3 0'
+            sage: e.dimacs_cnf([r.variable(1)+r.variable(0)])
+            'c cnf generated by PolyBoRi\np cnf 3 2\n-1 2 0\n1 -2 0'
+            sage: e.dimacs_cnf([r.variable(0)*r.variable(1)*r.variable(2), r.variable(1)+r.variable(0)])
+            'c cnf generated by PolyBoRi\np cnf 3 3\n-1 -2 -3 0\n1 -2 0\n-1 2 0'
+        """
+        clauses_list = [c for p in polynomial_system for c in self.
+                        dimacs_encode_polynomial(p)]
+        res = ["c cnf generated by PolyBoRi"]
+        r = polynomial_system[0].ring()
+        n_variables = r.n_variables()
+        res.append(f"p cnf {n_variables} {len(clauses_list)}")
+        res.extend(clauses_list)
+        return "\n".join(res)
+
+
+class CryptoMiniSatEncoder(CNFEncoder):
+    group_counter = 0
+
+    def dimacs_encode_polynomial(self, p):
+        r"""
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: d=dict()
+            sage: r = declare_ring(["x", "y", "z"], d)
+            sage: from sage.rings.polynomial.pbori.cnf import CryptoMiniSatEncoder
+            sage: e = CryptoMiniSatEncoder(r)
+            sage: p = d["x"]+d["y"]+d["z"]
+            sage: p.deg()
+            1
+            sage: len(p)
+            3
+            sage: e.dimacs_encode_polynomial(p)
+            ['x1 2 3 0\nc g 1 x + y + z']
+            sage: e.dimacs_encode_polynomial(p+1)
+            ['x1 2 -3 0\nc g 2 x + y + z + 1']
+        """
+        if p.deg() != 1 or len(p) <= 1:
+            res = super().dimacs_encode_polynomial(p)
+        else:
+            invert_last = bool(p.has_constant_part())
+            variables = list(p.vars_as_monomial().variables())
+            indices = [self.to_dimacs_index(v) for v in variables]
+            if invert_last:
+                indices[-1] = -indices[-1]
+            indices.append(0)
+            res = ["x" + " ".join(str(v) for v in indices)]
+        self.group_counter = self.group_counter + 1
+        group_comment = f"\nc g {self.group_counter} {str(p)[:30]}"
+        return [c + group_comment for c in res]
+
+    def dimacs_cnf(self, polynomial_system):
+        r"""
+
+        TESTS::
+
+            sage: from sage.rings.polynomial.pbori import *
+            sage: r = declare_ring(["x", "y", "z"], dict())
+            sage: from sage.rings.polynomial.pbori.cnf import CryptoMiniSatEncoder
+            sage: e = CryptoMiniSatEncoder(r)
+            sage: e.dimacs_cnf([r.variable(0)*r.variable(1)*r.variable(2)])
+            'c cnf generated by PolyBoRi\np cnf 3 1\n-1 -2 -3 0\nc g 1 x*y*z\nc v 1 x\nc v 2 y\nc v 3 z'
+            sage: e.dimacs_cnf([r.variable(1)+r.variable(0)])
+            'c cnf generated by PolyBoRi\np cnf 3 1\nx1 2 0\nc g 2 x + y\nc v 1 x\nc v 2 y'
+            sage: e.dimacs_cnf([r.variable(0)*r.variable(1)*r.variable(2), r.variable(1)+r.variable(0)])
+            'c cnf generated by PolyBoRi\np cnf 3 2\n-1 -2 -3 0\nc g 3 x*y*z\nx1 2 0\nc g 4 x + y\nc v 1 x\nc v 2 y\nc v 3 z'
+        """
+        uv = list(used_vars_set(polynomial_system).variables())
+        res = super().dimacs_cnf(polynomial_system)
+        res += "\n" + "\n".join(f"c v {self.to_dimacs_index(v)} {v}"
+                                for v in uv)
+        return res