Mason Protter

Feature request: Multiple timesteps in Exponentiate

1

Suppose I have some linear operator `A` and I want to know the result of `exp(t * A) * v` for a bunch of different values of `t`. The naïve...

Memory leak in solving ODE many times

3

The following code has monotically increasing memory usage for me and will just keep on growing until the system runs out of RAM and crashes: ```julia julia> using SparseArrays, LinearAlgebra,...

Supporting covariant derivatives

2

## The problem: Currently, when you compute the reverse mode derivative of some function of a `struct MyStruct` in enzyme, i.e. `autodiff(Reverse, f, Duplicated(A::MyStruct, dA::MyStruct))`, the object `dA` must be...

Not yet implemented: dgebal_64_ (Lapack)

1

Trying to differentiate through a matrix exponential took a very very very long time, and then finally errored. MWE: ```julia julia> let M = rand(2, 2) ∂M = similar(M) autodiff(Reverse,...

Feature request: `do` with infix calls

1

I was disapointed today to learn that I can't write the following macro: ``` @forall x in y do a, b ... end ``` because `(x in y) do` is...

https://github.com/JuliaLang/JuliaSyntax.jl/blob/a6f2d1580f7bbad11822033e8c83e607aa31f100/docs/src/design.md?plain=1#L364 Unfortunately, this seems to no longer be true: ```julia julia> JuliaSyntax.parsestmt(SyntaxNode, "@(a + b)") ERROR: ParseError: # Warning @ line 1:2 @(a + b) #└─────┘ ── parenthesizing macro names...

Add precompile statements to reduce startup latency

5

So @timholy and I have done a little preliminary looking into adding `precompile` statements to reduce the latency of IJulia kernels, and it seems to us that a major problem...

Not using "strong zeros" for deriatives of constant values

1

Compare: ```julia julia> using Diffractor, Zygote julia> var"'"(f) = Zygote.var"'"(f); julia> exp'(Inf) Inf julia> exp''(Inf) Inf ``` with ```julia julia> var"'"(f) = Diffractor.PrimeDerivativeFwd(f); julia> exp'(Inf) Inf julia> exp''(Inf) NaN ```...

Feature request: ability to specify method table in `methodinstance`

1

It'd be useful if `methodinstance` could work with any method table, not just the global default one. The version that takes a worldage is easy, it just needs to become...

The example in https://github.com/JuliaLinearAlgebra/InfiniteLinearAlgebra.jl/blob/master/examples/periodicschrodinger.jl is erroring out for me. MWE (same code, just not using `Plots`) ```julia using InfiniteLinearAlgebra A = BlockTridiagonal(Vcat([[0. 1.; 0. 0.]],Fill([0. 1.; 0. 0.], ∞)), Vcat([[-1....