ExponentialAction.jl
ExponentialAction.jl copied to clipboard
sethaxen
Compute the action of the matrix exponential
In https://github.com/sethaxen/ExponentialAction.jl/blob/v0.2.8/src/parameters.jl#L54 ``` Aᵖ⁺¹ = A * A d = t_norm * sqrt(opnormest1(A)) ``` Should `A` here be `Aᵖ⁺¹` instead? Because we are calculating the 1-norm for `A^2` as in...
I just checked, and on Julia v1.10.0, Enzyme v0.11.15 can differentiate `expv` with no troubles. It fails on `expv_sequence` though. Here's the code: ```julia using Enzyme, ExponentialAction f(t, A, v)...
This pull request changes the compat entry for the `AbstractDifferentiation` package from `0.4, 0.5` to `0.4, 0.5, 0.6`. This keeps the compat entries for earlier versions. Note: I have not...
![github-actions[bot] avatar](https://avatars.githubusercontent.com/in/15368?v=4 "github-actions[bot] avatar"){kind=avatar}
This PR adds enzyme compatibility by marking as inactivate using EnzymeRules the numerical methods used only for control flow. Since Enzyme supports mutation, a future PR could add mutating variants...
This PR adds support for computing the action of the exponential at half-precision tolerance. It also ensures that the Taylor parameters are computed using double precision even if the eltype...