Steps to reproduce:

• Run the UK QSP 2026 ipython notebook
• Reach the penultimate code cell, and fail on sol_func = solve(func_profile_prob, alg_int; reoptimize_init=true, verbose=true)

ERROR: TypeError: in setfield!, expected LinearAlgebra.QRCompactWY{Float64, Matrix{Float64}, Matrix{Float64}}, got a value of type LinearAlgebra.QRPivoted{Float64, ComponentMatrix{Float64, Matrix{Float64}, Tuple{Axis{(P = 1, Ap = 2, Ar = 3)}, Axis{(P = 1, Ap = 2, Ar = 3)}}}, Vector{Float64}, Vector{Int64}}
Stacktrace: ...

The error message indicates that a variable had an incorrect type, but it is the ODE solver itself failing. Changing the ODE solver options from alg = AutoTsit5(Rosenbrock23()) to alg = Tsit5() fixes it.

Error found on Julia 1.12.6, running on Ubuntu via WSL.

Julia info

julia> versioninfo()
Julia Version 1.12.6
Commit 15346901f00 (2026-04-09 19:20 UTC)
Build Info:
Official https://julialang.org release
Platform Info:
OS: Linux (x86_64-linux-gnu)
CPU: 16 × AMD Ryzen 7 PRO 6850U with Radeon Graphics
WORD_SIZE: 64
LLVM: libLLVM-18.1.7 (ORCJIT, znver3)
GC: Built with stock GC
Threads: 1 default, 1 interactive, 1 GC (on 16 virtual cores)

Status ~/.julia/environments/v1.12/Project.toml
[b0b7db55] ComponentArrays v0.15.40
[93acb638] LikelihoodProfiler v1.5.3
[22f7324a] OptimizationLBFGSB v1.4.2
[1344f307] OrdinaryDiffEqLowOrderRK v2.1.1
[43230ef6] OrdinaryDiffEqRosenbrock v2.3.2
[b1df2697] OrdinaryDiffEqTsit5 v2.0.2
[91a5bcdd] Plots v1.41.6
[2913bbd2] StatsBase v0.34.12

Operating system info

Edition Windows 11 Enterprise
Version 24H2
Installed on ‎19/‎06/‎2025
OS build 26100.8655
Experience Windows Feature Experience Pack 1000.26100.315.0

$ lsb_release -a
No LSB modules are available.
Distributor ID: Ubuntu
Description: Ubuntu 20.04.6 LTS
Release: 20.04
Codename: focal

$ wsl --version
WSL version: 2.7.10.0
Kernel version: 6.18.33.2-2
WSLg version: 1.0.73.2
MSRDC version: 1.2.6676
Direct3D version: 1.611.1-81528511
DXCore version: 10.0.26100.1-240331-1435.ge-release
Windows version: 10.0.26100.8655