LinearAlgebra/addmul/mul!(::LinearAlgebra.LowerTriangular{Float64, S} where S<:AbstractMatrix{Float64}, ::LinearAlgebra.Diagonal{Float64, V} where V<:AbstractVector{Float64}, ::LinearAlgebra.Diagonal{Int64, V} where V<:AbstractVector{Int64}, α, β)/α = 0.638064744234355, β = 0.638064744234355/β = 0 ignores C .= NaN @test ≈(collect(returned_mat), exp_val, rtol = rtol, atol = atol)