Daniel Wennberg

StaticArrays' implementation of Cholesky decomposition assumes that `uplo == 'U'` in `Symmetric` and `Hermitian` wrappers. This gives incorrect results when `uplo == 'L'` and the underlying array is not actually...

Threading kills performance when extracting values from a vector of dual numbers

8

This is strangely specific, but look, a factor ~200 slowdown when enabling multithreading: ```julia julia> using BenchmarkTools, FastBroadcast, ForwardDiff julia> N = 100; x = ForwardDiff.Dual.(randn(N), randn(N)); v = zeros(N);...

Separate RHS and state vectors for `FunctionOperator`

4

Currently, `FunctionOperator` can wrap functions on the form `f(u, p, t)` representing the application of an operator `L(u, p, t)` to the vector `u`, i.e., the operation `L(u, p, t)...

Operator fusion/matrix chain multiplication at constant `(u, p, t)`-slices

1

Say you have `C = ComposedOperator(A::MatrixOperator, B::MatrixOperator)`. You're passing it to an iterative solver. Depending on the shapes of `A` and `B`, precomputing the product `Cfused = A.A * B.A`...

Fix RemoteChannel iterator interface

2

Just a missing import

`multiple call sites` obscures allocation origin

1

Allocations reported by AllocCheck often have stack traces that only contain stack frames from Julia core/Base (array/dict/set constructors, locks, gc code, etc.) before terminating with `multiple call sites`. This is...

Fix fallback for errors in logging

2

This fixes the bug identified in #3012, in which the catch block that's supposed to save you in case an error is thrown during logging ends up throwing an error...

Error when displaying log messages from GPUCompiler

5

GPUCompiler.jl implements slightly tweaked loggers for safe use in generated functions: https://github.com/JuliaGPU/GPUCompiler.jl/blob/27d4b190243f9a4440e8c0bfda876af89dd1140e/src/utils.jl#L55-L73 Pluto.jl fails to display these logs, instead erroring with `ArgumentError: Base.TTY(RawFD(4294967295) invalid status, 0 bytes waiting) is not...

My attempts at forward-over-reverse Hessian-vector products keep crashing Julia. Here's the most minimal MWE I've been able to produce that exhibits this issue. It's essentially a further stripped-down version of...

Sticky log messages

2

I sometimes wish I could use `@info` to update an existing log message. Imagine the following, with each invocation of `@info "..." sticky=true` updating the same log card instead of...