search

sethaxen / ExponentialAction.jl

Compute the action of the matrix exponential
MIT License
24 stars 0 forks source link

Coefficients for theta doesn't agree with scipy _expm_multiply #27

Closed colinfang closed 6 months ago

colinfang commented 6 months ago

See scipy coefficients here for double tolerance

https://github.com/scipy/scipy/blob/v1.13.0/scipy/sparse/linalg/_expm_multiply.py#L309

It doesn't agree with values here https://github.com/sethaxen/ExponentialAction.jl/blob/v0.2.8/src/coefficients.jl#L69

I also checked the original paper (table 3.1). Julia coefficients match in the case of single tolerance, but not double (where scipy version uses)

colinfang commented 6 months ago

Sry, the issue is invalid. They do match against each other.