Klaus Crusius

DynamicalODE with immutable state variables and callbacks failing

4

Running this example with immutatble state and callback fails. Obviously `addsteps!` tries to write into the `SVector` state variables. ``` julia> using StaticArrays julia> using OrdinaryDiffEq, LinearAlgebra, ForwardDiff, Plots; gr()...

Overflow behaviour of Double64 does not match Float64

84

Please re-open this issue: The root cause seems to be in `DoubleFloats`: ``` julia> Double64(1e300)^2 NaN julia> Float64(1e300)^2 Inf ``` _Originally posted by @KlausC in https://github.com/JuliaMath/DoubleFloats.jl/issues/149#issuecomment-1133882519_

I think, each `TranscodingStream` should implement `isreadable(io::TranscodingStream) = isreadable(io.stream)` `iswritable(io::TranscodingStream) = iswritable(io.stream)`

New Algorithms

5

Is it valid to consider the results of this paper from 2011: [Inner-Iteration Krylov Subspace Methods for Least Squares Problems](https://www.nii.ac.jp/TechReports/public_html/11-001E.pdf)? Especially the algorithm called there "BAGMRES-NRSOR" seems to show outstanding...

broadcast(/, ::SparseMatrixCSC, ::AbstractArray) stores unnecessary zeros

10

If `A` is a sparse matrix, `B` is a sparse or dense vector of non-zeros with appropriate size, then `C = A ./ B` and `C = A ./ B'`...

Broadcasting over sparse matrix generating wrong type

1

``` julia> VERSION v"1.7.0-DEV.1108" julia> rem.(sparse([1 2; 3 4]), 2) 2×2 SparseMatrixCSC{Int64, Int64} with 2 stored entries: 1 ⋅ 1 ⋅ julia> big.(rem.(sparse([1 2; 3 4]), 2)) 2×2 SparseMatrixCSC{BigInt, Int64}...

Switch documentation from `readthedocs.org` to `Documenter.jl`. * Setup docs folder and make.jl * Remove redundancies form documentation * Make use of Documenter-features as far as appropriate

Code coverage dropped from 99% to 66%. That is not acceptable. New target: 100%

solves #202 proposed insertion of `copy` methods (see https://github.com/JuliaMatrices/BandedMatrices.jl/issues/202#issuecomment-688312771) reveals bug in `lu!(BandedMatrix)`. Further investigation required.

Type of factors in Cholesky decomposition

1

The types of the factors of a Cholesky decomposition of a banded matrix do not reflect the banded layout, but produce a full matrix in half of the cases. Would...