LinearAlgebra/addmul/mul!(::Matrix{BigFloat}, ::LinearAlgebra.Bidiagonal{Float32, V} where V<:AbstractVector{Float32}, ::LinearAlgebra.LowerTriangular{Float32, S} where S<:AbstractMatrix{Float32}, α, β)/α = 1.0, β = 1.04796756784327982359172892756760120391845703125/β = 0 ignores C .= NaN @test ≈(collect(returned_mat), exp_val, rtol = rtol, atol = atol)