How can verify a neural network using SMT/ILP/SAT solvers?

Suppose I have a trained network consisting of a fully connected layer and relu activation. How can use the trained weights to verify the accuracy of the model. In my case, I will be classifying MNIST dataset. Is there any package or functions to solve the problem similar to z3 in python or any other SMT/SAT solvers?