orthopy
A Python library providing orthogonal polynomial classes for lines, triangles, disks, spheres, n-cubes, and nD spaces with stable recurrence schemes.
What is orthopy?
orthopy is a Python library that provides orthogonal polynomial classes for various geometric domains like lines, triangles, disks, spheres, n-cubes, and nD spaces. It solves the problem of accurately computing orthogonal polynomials using numerically stable recurrence schemes, which is essential for numerical integration, approximation, and spectral methods.
Researchers, scientists, and engineers working in numerical analysis, computational physics, finite element methods, or any field requiring orthogonal polynomial bases for approximation and integration.
Developers choose orthopy for its comprehensive coverage of polynomial families across multiple domains, numerically stable implementations, and support for both numerical and symbolic computation, all within a unified Python interface.
Overview
:triangular_ruler: Orthogonal polynomials in all shapes and sizes.
Use Cases
Best For
- Numerical integration using Gaussian quadrature rules
- Spectral methods for partial differential equations
- Approximation theory and function approximation
- Finite element method basis function generation
- Computational physics simulations requiring orthogonal bases
- Symbolic manipulation of orthogonal polynomials
Not Ideal For
- Projects requiring only basic orthogonal polynomials for common domains (e.g., Legendre on [-1,1]) where SciPy's free functions suffice
- Teams with strict open-source or zero-cost licensing requirements, as commercial use requires a paid license
- Applications needing real-time performance with minimal overhead, due to symbolic computation and vectorization latency
- Developers unfamiliar with orthogonal polynomial theory, given the specialized iterator-based interface and recurrence schemes
Pros & Cons
Pros
Uses stable recurrence schemes to ensure accuracy for high-degree polynomials, critical for reliable computations in numerical integration and approximation.
Functions are fully vectorized for efficient evaluation at multiple points, enabling batch processing and performance gains in scientific computing.
Integrates with SymPy for exact arithmetic and symbolic polynomial generation, allowing precise manipulation and analysis for theoretical work.
Supports diverse geometric domains like triangles, disks, and spheres with various weight functions, providing a unified interface for advanced numerical methods.
Cons
Requires a paid license for commercial use via MondayTech, adding cost and administrative complexity compared to fully open-source alternatives like SciPy.
The iterator-based `Eval` pattern and need to manage recurrence coefficients can be less intuitive than direct function calls, posing a steep learning curve.
As a specialized library, it has fewer community resources and examples, with advanced usage often relying on external publications for guidance.
Frequently Asked Questions
Related Projects
quadpy:triangular_ruler: Numerical integration (quadrature, cubature) in Python
HPDDMA framework for high-performance domain decomposition methods.
PyDMDmathLab mirror of Python Dynamic Mode Decomposition
accupy:dart: Accurate sums and dot products for Python.
Found a gem we're missing?
Open-Awesome is built by the community, for the community. Submit a project, suggest an awesome list, or help improve the catalog on GitHub.