Skip to main content

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.

Project details


Download files

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

Source Distribution

Diophantine-0.2.0.tar.gz (7.2 kB view hashes)

Uploaded Source

Built Distribution

Diophantine-0.2.0-py2.py3-none-any.whl (9.8 kB view hashes)

Uploaded Python 2 Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page