Ioannis Papadopoulos
Ioannis Papadopoulos
Is there a way to work with Fourier and BigFloat? It goes via FFTW which is causing some issues. ``` using ClassicalOrthogonalPolynomials F = Fourier{BigFloat}() F[:, 1:10] \ sin.(axes(F,1)) ```...
@dlfivefifty Is there a reason why the LaguerreWeight axis was originally all of R rather than just the HalfLine?
More of a question than an issue. The following works just fine: `Weighted(ChebyshevT())\Weighted(ChebyshevU())` but `Weighted(ChebyshevU())\Weighted(ChebyshevT())` throws an` ERROR: OutOfMemoryError()` message. Is this expected behaviour? Or do we want something similar...
I managed to implement the partial derivatives with respect to y. The following code now runs: `using MultivariateOrthogonalPolynomials W = Weighted(Zernike(1)) Δ = Laplacian(axes(W,1)) ∂ʸ = PartialDerivative{2}(axes(W,1)) Z⁰ = Zernike(0)...
@dlfivefifty The lowering `Zernike(0,1) \ Weighted(Zernike(0,1))` is a factor of 2 off. This also seems to be confirmed in the tests..? https://github.com/JuliaApproximation/MultivariateOrthogonalPolynomials.jl/blob/cb0fb241eecc49fbdb873a5cd2b7de098897d355/test/test_disk.jl#L303 Is there a reason for this?
Timon @TSGut has done the heavy lifting and has implemented Jacobi matrices for Zernike polynomials on a disc. Alberto @APaganini could you remind us what exactly you required? If I...