Chris Elrod

A few comments: > I don't have any background in integer optimization, so odds are a lot of things could be done better / I am open to suggestions. I...

```julia julia> format_text(""" begin a = rand(100); sum(a) end """) |> println begin a = rand(100) sum(a) end ``` When debugging, I often run the contents of blocks in my...

Could we also get have an option we can set in our `.JuliaFormatter.toml` to end all lines with a semicolon?

> Do we have an efficient Julia implementation of `dot` somewhere to easily compare with BLAS versions? ```julia function dot_julia(x,y) s = zero(promote_eltype(x, y)) @fastmath for i in eachindex(x, y)...

Updated script with some fixes: ```julia using BenchmarkTools, LinearAlgebra, LoopVectorization # `@noinline`, otherwise the compiler optimizes it away @noinline function dot_julia(x,y) s = zero(Base.promote_eltype(x, y)) # I thought Julia's bound...

> Maybe part of it is that OpenBLAS is just exceedingly awful on Ryzen If you were on Julia 1.6, this would be the case. Julia 1.6.5 used OpenBLAS 0.3.10:...

I'll take a look at this.

I can't reproduce this. 1.5: ```julia julia> @time using OrdinaryDiffEq, Base.Threads, Random 0.000234 seconds (999 allocations: 65.375 KiB) julia> N = 512; julia> Random.seed!(42); u = rand(2N-1, N); du =...

> And I'm the complete opposite! For me, 1.5 was 3.8% slower than 1.6 and not just 2% slower ;) . I also have a very similar CPU to OP...

I feel like "FastBroadcast.jl" is the new "thanks, Obama". My docs aren't building? Thanks, Fastbroadcast.jl. I'll try an official binary.