Skip to main content
Help improve PyPI by participating in a 5-minute user interface survey!

A python package for finding small solutions of systems of diophantine (integer algebra) equations

Project Description

Unit Test Status Unit Test Coverage

Author: Thomas G. Close (tom.g.close@gmail.com)

Diophantine (http://github.com/tclose/Diophantine) is a Python package for finding small (integer) solutions of systems of diophantine equations (see http://en.wikipedia.org/wiki/Diophantine_equation). It is based on PHP code by Keith Matthews (see www.number-theory.org) that implements the algorithm described in https://github.com/tclose/Diophantine/blob/master/algorithm.pdf (see http://www.numbertheory.org/lll.html for a list of associated publications), which uses the LLL algorithm to calculate the Hermite-normal-form described in the paper:

Extended gcd and Hermite normal form algorithms via lattice basis reduction, G. Havas, B.S. Majewski, K.R. Matthews, Experimental Mathematics, Vol 7 (1998) 125-136

(please cite this paper if you use this code in a scientific publication)

There are two branches of this code in the GitHub repository (see https://github.com/tclose/Diophantine.git), ‘master’, which uses the sympy library and therefore uses arbitrarily long integer representations, and ‘numpy’, which uses the numpy library, which is faster but can suffer from integer overflow errors despite using int64 representations

To find small solutions to a system of diophantine equations, A x = b, where A is a M x N matrix of coefficents, b is a M x 1 vector and x is the N x 1 vector, use the ‘solve’ method in the module, e.g.

>>> from sympy import Matrix
>>> from diophantine import solve
>>> A = Matrix([[1, 0, 0, 2], [0, 2, 3, 5], [2, 0, 3, 1], [-6, -1, 0, 2],
                [0, 1, 1, 1], [-1, 2, 0,1], [-1, -2, 1, 0]]).T
>>> b = Matrix([1, 1, 1, 1])
>>> solve(A, b)
[Matrix([
[-1],
[ 1],
[ 0],
[ 0],
[-1],
[-1],
[-1]])]

The returned solution vector will tend to be one with the smallest norms. If multiple solutions with the same norm are found they will all be returned. If there are no solutions the empty list will be returned.

Diophantine is released under the MIT Licence (see Licence for details)

Installation

Diophantine is available from the Python Package Index (http://pypi.python.org) and can be installed with the command

pip install diophantine

Alternatively, the master branch can be installed from using the setuptools install command

python setup.py install

from a clone of the GitHub repository (https://github.com/tclose/Diophantine) or simply adding the cloned directory to your PYTHONPATH.

Release history Release notifications

This version
History Node

0.2.0

History Node

0.1.1

History Node

0.1

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Filename, size & hash SHA256 hash help File type Python version Upload date
Diophantine-0.2.0-py2.py3-none-any.whl (9.8 kB) Copy SHA256 hash SHA256 Wheel py2.py3 Nov 20, 2017
Diophantine-0.2.0.tar.gz (7.2 kB) Copy SHA256 hash SHA256 Source None Nov 20, 2017

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging CloudAMQP CloudAMQP RabbitMQ AWS AWS Cloud computing Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page