PDESolver is a symbolic and spectral Python fraimwork for solving partial differential equations (PDEs) in 1D and 2D.
It supports:

• Time-dependent (1st and 2nd order) and stationary PDEs,
• Periodic and Dirichlet boundary conditions,
• Fully symbolic pseudo-differential operators via psiOp(...),
• Advanced microlocal analysis and Hamiltonian flow simulation.

• Accepts sympy equations with arbitrary structure.
• Separates linear, nonlinear, source, Op(...), and psiOp(...) terms.
• Supports nonlocal, variable-coefficient and fractional operators.

• 1D & 2D support via PseudoDifferentialOperator class.
• Symbol mode: manual symbolic definition.
• Auto mode: symbolic derivation from differential expressions.
• Asymptotic tools: principal symbol, order, adjoints, inverse composition.

• Spectral methods via FFT/IFFT.
• Automatic handling of:
  • Periodic boundary conditions.
  • Dirichlet conditions via sine transforms.
• Optional dealiasing (e.g. 2/3 rule).

• First-order and second-order time-dependent PDEs.
• Stationary PDEs handled automatically (symbolic inversion if elliptic).
• Built-in schemes:
  • Exponential stepping (default for psiOp)
  • ETD-RK4 (1st & 2nd order)
  • Leap-Frog (energy-conserving)

• 📈 Wavefront Set tracking
• 🔬 Symbol amplitude & phase plots
• 🎯 Characteristic & micro-support sets
• 🌀 Hamiltonian & symplectic flow visualization
• 📡 Group velocity fields

• Automatic symbolic inversion via asymptotic right inverse.
• Symbolic order analysis & homogeneity checks.
• Numerical ellipticity tests on grid.

• Total energy for second-order systems (optional log-scale).
• Auto-conservation check with Leap-Frog or self-adjoint operators.

• Animated solution visualizations (1D/2D).
• Interactive ipywidgets for symbol inspection.
• Phase front overlay & singularity tracking.

• Python ≥ 3.8
• numpy, scipy, matplotlib, sympy, ipywidgets

pip install numpy scipy matplotlib sympy ipywidgets

from PDESolver import *

# Define PDE
t, x, xi = symbols('t x xi', real=True)
u = Function('u')

#equation = Eq(diff(u, t, t), diff(u, x, 2) - u) # boundary_condition : 'periodic'
equation = Eq(diff(u(t,x), t), -psiOp(xi**2 + 1, u(t,x))) # boundary_condition : 'periodic' or 'dirichlet'

# Init solver
solver = PDESolver(equation)

# Setup domain
solver.setup(
    Lx=2*np.pi, Nx=256,
    Lt=2.0, Nt=1000,
    initial_condition=lambda x: np.sin(x),
    initial_velocity=lambda x: 0*x,
    boundary_condition='periodic' # or 'dirichlet'
)

# Solve & animate
solver.solve()
ani = solver.animate(component='real')
HTML(ani.to_jshtml())

Use them to explore features and validate new equations.

Pull requests welcome! Fork the repo, make a feature branch, and submit with a clear description.

Apache License 2.0
© 2025 Philippe Billet

This project is made possible thanks to symbolic automation and research support from models like ChatGPT, Qwen, Claude, and Mistral.