Feature Request: Add Import/Export Functions for Boolean Expressions in Z3 and BDD Formats

Question

Feature Request: Add Import/Export Functions for Boolean Expressions in Z3 and BDD Formats

YuriyTigiev opened this issue 22 days ago · comments

Hi SymPy Developers,

Could you please add functions for exporting and importing Boolean expressions in Z3 and BDD formats in the next release? This feature would be incredibly useful.

Thank you for your hard work!

Shekhar Prasad Rajak · Answer 1 · Fri Jun 14 2024 17:30:55 GMT+0800 (China Standard Time)

Can you please share the usecase for the requirement ?

YuriyTigiev · Answer 2 · Fri Jun 14 2024 22:04:21 GMT+0800 (China Standard Time)

They have better performance for optimization and solving boolean expressions. But sympy has better api interface for working with expressions.

…

________________________________ From: Shekhar Prasad Rajak ***@***.***> Sent: Friday, June 14, 2024 1:31:20 PM To: sympy/sympy ***@***.***> Cc: YuriyTigiev ***@***.***>; Author ***@***.***> Subject: Re: [sympy/sympy] Feature Request: Add Import/Export Functions for Boolean Expressions in Z3 and BDD Formats (Issue #26705) Can you please share the usecase for the requirement ? — Reply to this email directly, view it on GitHub<#26705 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AFNSYQSNSY47C7EXKJWTE63ZHKZ6LAVCNFSM6AAAAABJJSG4HKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCNRXGYZTIOBVGM>. You are receiving this because you authored the thread.Message ID: ***@***.***>

Oscar Benjamin · Answer 3 · Sat Jun 15 2024 00:46:37 GMT+0800 (China Standard Time)

The feature request here is reasonable but it is not useful to have an issue open for it if there is no detail about what could be done.

There is no need to have an open issue for someone to work on this if they want. Someone can just open a PR.

I suggest closing this unless some useful information is added.

S.Y. Lee · Answer 4 · Sat Jun 15 2024 13:59:22 GMT+0800 (China Standard Time)

Can you take a look at https://github.com/sympy/sympy/blob/master/sympy/printing/smtlib.py first before that can be still usable for your application?

YuriyTigiev · Answer 5 · Sun Jun 16 2024 00:01:06 GMT+0800 (China Standard Time)

Sorry, but I didn't fully understand the meaning of your advice. Could you please clarify?

…

________________________________ From: S.Y. Lee ***@***.***> Sent: Saturday, June 15, 2024 9:59 AM To: sympy/sympy ***@***.***> Cc: YuriyTigiev ***@***.***>; Author ***@***.***> Subject: Re: [sympy/sympy] Feature Request: Add Import/Export Functions for Boolean Expressions in Z3 and BDD Formats (Issue #26705) Can you take a look at https://github.com/sympy/sympy/blob/master/sympy/printing/smtlib.py first before that can be still usable for your application? — Reply to this email directly, view it on GitHub<#26705 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AFNSYQQTLVBNDXQATZK5YSTZHPJ5DAVCNFSM6AAAAABJJSG4HKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCNRZGE2DQOBUGA>. You are receiving this because you authored the thread.Message ID: ***@***.***>

S.Y. Lee · Answer 6 · Sun Jun 16 2024 00:08:31 GMT+0800 (China Standard Time)

What are Z3 and BDD formats anyway? You have come up with ideas not easily searchable in google, but without any references.

YuriyTigiev · Answer 7 · Sun Jun 16 2024 00:19:39 GMT+0800 (China Standard Time)

Here is some information about the Z3 and dd libraries. I also discovered the pyeda library, which has implemented more powerful algorithms (Espresso and SAT) for optimization. https://github.com/Z3Prover/z3 ### Z3 Z3 is a high-performance theorem prover developed by Microsoft Research. It includes various functionalities for working with Boolean expressions, including simplification and solving. Here are some key functions: #### Simplifying Boolean Expressions - **`simplify(expr)`**: Simplifies a given Boolean expression. ```python from z3 import * x, y, z = Bools('x y z') expr = Or(And(x, y), And(Not(x), z)) simplified_expr = simplify(expr) print(simplified_expr) ``` #### Solving Boolean Expressions - **`Solver()`**: Creates a solver instance. - **`solver.add(expr)`**: Adds a constraint to the solver. - **`solver.check()`**: Checks the satisfiability of the added constraints. - **`solver.model()`**: Returns a model (solution) if the constraints are satisfiable. ```python s = Solver() s.add(And(x, Or(y, z))) if s.check() == sat: m = s.model() print(m) else: print("No solution") ``` https://github.com/tulip-control/dd ### BDD (Binary Decision Diagrams) BDD libraries like `dd` (part of `tulip-control`) are used for efficient representation and manipulation of Boolean functions. Here's a brief overview of some useful functions: #### Creating and Manipulating BDDs - **`bdd.add_expr(expr)`**: Adds a Boolean expression to the BDD. - **`bdd.apply(op, u, v)`**: Applies a binary operation (`op`) to two BDD nodes `u` and `v`. - **`bdd.simplify()`**: Simplifies the BDD representation. ```python from dd import autoref as _bdd bdd = _bdd.BDD() bdd.declare('x', 'y', 'z') u = bdd.add_expr('x & y') v = bdd.add_expr('~x & z') simplified_u = bdd.simplify(u) print(bdd.to_expr(simplified_u)) ``` #### Solving Boolean Expressions BDD libraries typically do not "solve" Boolean expressions like SAT solvers but can be used to check satisfiability and manipulate the expressions efficiently. - **`bdd.let(mapping, u)`**: Substitutes variables in a BDD node `u` according to a given `mapping`. ```python mapping = {'x': True, 'y': False} result = bdd.let(mapping, u) print(bdd.to_expr(result)) ``` These libraries are powerful tools for working with Boolean expressions, each offering unique functionalities and optimizations.

…

________________________________ From: S.Y. Lee ***@***.***> Sent: Saturday, June 15, 2024 8:08 PM To: sympy/sympy ***@***.***> Cc: YuriyTigiev ***@***.***>; Author ***@***.***> Subject: Re: [sympy/sympy] Feature Request: Add Import/Export Functions for Boolean Expressions in Z3 and BDD Formats (Issue #26705) What are Z3 and BDD formats anyway? You have come up with ideas not easily searchable in google, but without any references. — Reply to this email directly, view it on GitHub<#26705 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AFNSYQWQXLWSASOFFNIUAKTZHRRJJAVCNFSM6AAAAABJJSG4HKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCNZQGA4TIMRVHE>. You are receiving this because you authored the thread.

S.Y. Lee · Answer 8 · Sun Jun 16 2024 00:27:45 GMT+0800 (China Standard Time)

The documentation of z3 says that input format is smtlib.
So I'm not sure why my suggestion was not even helpful.

https://github.com/Z3Prover/z3?tab=readme-ov-file#interfaces

Default input format is SMTLIB2

I'm not sure about bdd, however, I'm not even sure about it even can input from some text files
For example, are you suggesting to use their direct Python interface rather than specialized input-output text files?

YuriyTigiev · Answer 9 · Sun Jun 16 2024 00:41:58 GMT+0800 (China Standard Time)

I needed direct conversion from code without using text files.

Below is an example of the code generated by ChatGPT, but I couldn't use it.

I don't quite understand how I can use SMTLIB2.

Below is an example of one of the boolean expressions I am trying to simplify. This is a very small expression. In real tasks, the number of variables will be 100 times greater.


from sympy import symbols
from sympy.logic.boolalg import And, Or, Not, BooleanFunction
from z3 import And as Z3And, Or as Z3Or, Not as Z3Not, Bool, simplify, Solver

# Function to convert a SymPy expression to Z3
def sympy_to_z3(expr):
    if isinstance(expr, BooleanFunction):
        if isinstance(expr, And):
            return Z3And(*[sympy_to_z3(arg) for arg in expr.args])
        elif isinstance(expr, Or):
            return Z3Or(*[sympy_to_z3(arg) for arg in expr.args])
        elif isinstance(expr, Not):
            return Z3Not(sympy_to_z3(expr.args[0]))
    elif expr.is_Symbol:
        return Bool(str(expr))
    else:
        raise ValueError(f"Unsupported SymPy expression: {expr}")

# Function to convert a Z3 expression to SymPy
def z3_to_sympy(expr):
    if expr.decl().kind() == 263:  # Z3_OP_AND
        return And(*[z3_to_sympy(arg) for arg in expr.children()])
    elif expr.decl().kind() == 258:  # Z3_OP_OR
        return Or(*[z3_to_sympy(arg) for arg in expr.children()])
    elif expr.decl().kind() == 260:  # Z3_OP_NOT
        return Not(z3_to_sympy(expr.children()[0]))
    elif expr.decl().kind() == 235:  # Z3_OP_UNINTERPRETED
        return symbols(str(expr))
    else:
        raise ValueError(f"Unsupported Z3 expression: {expr}")

# Function to optimize a Z3 expression using Z3's solver
def optimize_z3_expr(z3_expr):
    s = Solver()
    s.add(z3_expr)
    if s.check() == z3.sat:
        model = s.model()
        return simplify(z3_expr)
    else:
        raise ValueError("The expression is unsatisfiable")

# Example usage
a, b, c = symbols('a b c')
sympy_expr = And(a, Or(b, Not(c)))

# Convert from SymPy to Z3
z3_expr = sympy_to_z3(sympy_expr)
print(f"Original Z3 Expression: {z3_expr}")

# Optimize the Z3 expression
optimized_z3_expr = optimize_z3_expr(z3_expr)
print(f"Optimized Z3 Expression: {optimized_z3_expr}")

# Convert back from Z3 to SymPy
sympy_converted = z3_to_sympy(optimized_z3_expr)
print(f"Optimized SymPy Expression: {sympy_converted}")

((x0 & ~x2) | (x1 & ~x2) | (x2 & ~x0 & ~x1)) & ((y2 & ~y4) | (y3 & ~y4) | (y4 & ~y2 & ~y3) | (y4 & ~y0 & ~y1 & ~y3)) & ((x0 & y4) ^ (x1 & y3) ^ (x2 & y2) ^ ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2)) ^ ((c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3))) & ~(((x0 & x1 & x2 & y2 & y3 & y4 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x0 & x1 & x2 & y2 & y3 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | (x0 & x1 & y3 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x0 & x2 & y2 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x1 & x2 & y2 & y3 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2)))) & ((x0 & x1 & y3 & y4) | (x0 & x2 & y2 & y4) | (x1 & x2 & y2 & y3) | (x0 & y4 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x1 & y3 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x2 & y2 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x0 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | (x1 & y3 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | (x2 & y2 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | ((c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))))) & ~(((x0 & x1 & x2 & y2 & y3 & y4 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x0 & x1 & x2 & y2 & y3 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | (x0 & x1 & y3 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x0 & x2 & y2 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x1 & x2 & y2 & y3 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2)))) ^ ((x0 & x1 & y3 & y4) | (x0 & x2 & y2 & y4) | (x1 & x2 & y2 & y3) | (x0 & y4 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x1 & y3 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x2 & y2 & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2))) | (x0 & y4 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | (x1 & y3 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | (x2 & y2 & (c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3)) | ((c2.0 | c2.1) & (c2.0 | x0) & (c2.0 | x1) & (c2.0 | x2) & (c2.0 | y1) & (c2.0 | y2) & (c2.0 | y3) & (c2.1 | x0) & (c2.1 | x1) & (c2.1 | x2) & (c2.1 | y1) & (c2.1 | y2) & (c2.1 | y3) & (x0 | x1) & (x0 | x2) & (x0 | y1) & (x0 | y2) & (x1 | x2) & (x1 | y1) & (x1 | y3) & (x2 | y2) & (x2 | y3) & (y1 | y2) & (y1 | y3) & (y2 | y3) & ((c2.0 & c2.1) | (c2.0 & x0 & y3) | (c2.0 & x1 & y2) | (c2.0 & x2 & y1) | (c2.1 & x0 & y3) | (c2.1 & x1 & y2) | (c2.1 & x2 & y1) | (x0 & x1 & y2 & y3) | (x0 & x2 & y1 & y3) | (x1 & x2 & y1 & y2)))))

S.Y. Lee · Answer 10 · Sun Jun 16 2024 00:56:36 GMT+0800 (China Standard Time)

It looks like your specification anyway needs extra dependency like z3 python interface.

I'm not even sure if that project have to be inside SymPy,
because it should anyway be dependent on both the versions of SymPy and Z3, for instance.

If you are having some expectation like, to ask us to do something in the next release

Could you please add functions for exporting and importing Boolean expressions in Z3 and BDD formats in the next release?

Unfortunately, the expectation may not be matched for you, sorry.
I don't think that it is high priority project.

And you may better spend time learning about SymPy, and try to fix your issue in your own project,
or try to employ someone else who can do the task.

Oscar Benjamin · Answer 11 · Sun Jun 16 2024 01:08:07 GMT+0800 (China Standard Time)

I'm not even sure if that project have to be inside SymPy,

I think it makes sense to have this in SymPy. There already some things in sympy/logic/algorithms/z3_wrapper.py but just not quite a function for going directly from a SymPy expression.

The ChatGPT code is almost right. It just hallucinated the numbers for e.g. Z3_OP_AND:

In [24]: z3.Z3_OP_AND
Out[24]: 261

It would be better to use z3.Z3_OP_AND in the code rather than hardcoding the numbers.

Note to anyone looking to install z3 since this is the second time I have forgotten this: although the module is called z3 it is name is z3-solver on PyPI.