FastLSQ 0.2.5__tar.gz → 0.3.0__tar.gz

{fastlsq-0.2.5 → fastlsq-0.3.0}/CHANGELOG.md

@@ -2,6 +2,45 @@
 
 All notable changes to FastLSQ will be documented in this file.
 
+## [0.3.0] - 2026-06-21
+
+### Added
+
+- **Exact closed-form integral & integro-differential operators.** The cyclic
+  derivative identity `D^α sin(W·x+b) = (∏ W^{α_k}) Φ_{|α| mod 4}` runs *backwards* --
+  integration is differentiation of *negative* order, `∫ sin(wx+b) dx = −(1/w)cos(wx+b)`.
+  `SinusoidalBasis.derivative` now accepts **signed** multi-indices (negative entries =
+  indefinite integration), with a DC guard (`dc_eps`) that zeros features whose frequency
+  along an integrated axis is ~0 (their primitive is a ramp that leaves the sinusoidal
+  family). `DiffOperator.antiderivative(dim, order, d)` exposes this as a factory, so
+  `Op.partial(0,1,1) - k*Op.antiderivative(0,1,1)` composes integro-differential operators.
+- **`IntegralOperator`** -- definite and running (Volterra) integrals with limits,
+  via the new `SinusoidalBasis.definite_integral`. Evaluated with a numerically stable
+  `sinc` identity (no `1/w` division), so the running integral of a near-DC feature is
+  exact rather than singular. Factories `IntegralOperator.volterra(...)` /
+  `IntegralOperator.definite(...)`.
+- **`IntegroDifferentialOperator`** -- the common roof under which differential and
+  integral terms compose (`+`, `−`, scalar/`nn.Parameter` `*`) into one `(M, N)`
+  linear-least-squares design matrix; coefficients stay differentiable so learnable
+  integral-term coefficients train through the solve. All three exported from `fastlsq`.
+- **Examples** `examples/integro_differential_demo.py` (one-shot forward solve of
+  `u'(x) + ∫_0^x u ds = f`, rel-L2 ~3e-11) and
+  `examples/inverse/inverse_memory_kernel.py` (recover an unknown memory strength λ from
+  noisy data via AdamW through the differentiable solve). New `tests/test_integral.py`
+  asserts ∫-then-∂ round-trips to identity, Volterra/definite match quadrature, the DC
+  guard stays finite, and gradients flow to learnable coefficients.
+
+## [0.2.6] - 2026-06-09
+
+### Changed
+
+- **Faster `_auto_solve` fallback for ill-conditioned systems.** When the
+  Cholesky path fails on CPU with no ridge (`mu = 0`), `_auto_solve` now goes
+  straight to the rank-deficient-safe SVD solve (LAPACK `gelsd`) instead of
+  first attempting Householder QR with the blow-up guard -- on CPU `gelsd` is
+  faster than QR, so the detour only added a full extra factorization. The
+  QR-then-SVD path is unchanged for ridge solves and non-CPU devices.
+
 ## [0.2.5] - 2026-06-04
 
 ### Fixed

{fastlsq-0.2.5 → fastlsq-0.3.0}/FastLSQ.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: FastLSQ
-Version: 0.2.5
+Version: 0.3.0
 Summary: One-shot PDE solving via Fourier features with exact analytical derivatives; rank-revealing solvers, learnable anisotropic bandwidth, and CPU/CUDA/MPS support
 Author: Antonin Sulc
 License-Expression: MIT

{fastlsq-0.2.5 → fastlsq-0.3.0}/FastLSQ.egg-info/SOURCES.txt

@@ -18,6 +18,7 @@ examples/grid_rl_control.py
 examples/grid_swing.py
 examples/gs_inverse.py
 examples/gs_rl_control.py
+examples/integro_differential_demo.py
 examples/inverse_heat_source.py
 examples/inverse_magnetostatics.py
 examples/inverse_source_position.py
@@ -32,6 +33,17 @@ examples/run_nonlinear.py
 examples/tutorial_basic.py
 examples/tutorial_nonlinear.py
 examples/vector_basis_stream_vorticity.py
+examples/digital_twins/darcy_heat.py
+examples/digital_twins/pendulum.py
+examples/digital_twins/pendulum_benchmark.py
+examples/digital_twins/plasma_wakefield.py
+examples/digital_twins/plasma_wakefield_2D_1.py
+examples/digital_twins/plasma_wakefield_2D_2.py
+examples/digital_twins/plasma_wakefield_2d_3.py
+examples/digital_twins/plasma_wakefield_parameteric.py
+examples/digital_twins/plot_utils.py
+examples/digital_twins/structural_health_simple.py
+examples/digital_twins/turbulence_gravity_cooling.py
 examples/extras/fred_sde.py
 examples/extras/fred_sde_fastlsq.py
 examples/extras/gaia_potential.py
@@ -73,6 +85,19 @@ examples/extras/scenarios/s12_mems_resonator.py
 examples/extras/scenarios/s13_variable_stars_kepler.py
 examples/extras/scenarios/s14_eeg.py
 examples/extras/scenarios/s15_circadian.py
+examples/inverse/aero_.py
+examples/inverse/denoising_parameter_estimation.py
+examples/inverse/elastic_wave_animation.py
+examples/inverse/heat_from_video.py
+examples/inverse/inverse_memory_kernel.py
+examples/inverse/inverse_turbulence.py
+examples/inverse/shape_ns.py
+examples/inverse/subsurface_imaging.py
+examples/inverse/wing_optimize_simple.py
+examples/sindy/compare_sindy_methods.py
+examples/sindy/sindy_benchmarks.py
+examples/sindy/sindy_differentiable.py
+examples/sindy/sindy_minimal_diff.py
 fastlsq/__init__.py
 fastlsq/api.py
 fastlsq/basis.py
@@ -95,6 +120,16 @@ fastlsq/problems/__init__.py
 fastlsq/problems/linear.py
 fastlsq/problems/nonlinear.py
 fastlsq/problems/regression.py
+misc/fastlsq_teaser.png
+misc/ideal_quadrupole.png
+misc/inverse_heat_source.gif
+misc/inverse_heat_source.png
+misc/inverse_magnetostatics.png
+misc/inverse_magnetostatics_convergence.png
+misc/quadrupole_convergence.png
+misc/quadrupole_optimization.png
+misc/tutorial_nlpoisson_convergence.png
+misc/tutorial_nlpoisson_solution.png
 tests/test_basic.py
 tests/test_benchmarks_inverse.py
 tests/test_block.py
@@ -102,6 +137,7 @@ tests/test_derivatives.py
 tests/test_device.py
 tests/test_grad_shafranov.py
 tests/test_grid_swing.py
+tests/test_integral.py
 tests/test_learnable.py
 tests/test_orbit_hill.py
 tests/test_vector_basis.py

{fastlsq-0.2.5 → fastlsq-0.3.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: FastLSQ
-Version: 0.2.5
+Version: 0.3.0
 Summary: One-shot PDE solving via Fourier features with exact analytical derivatives; rank-revealing solvers, learnable anisotropic bandwidth, and CPU/CUDA/MPS support
 Author: Antonin Sulc
 License-Expression: MIT