LinearAlgebra/addmul/mul!(::LinearAlgebra.UpperTriangular{Float64, S} where S<:AbstractMatrix{Float64}, ::LinearAlgebra.UpperTriangular{Float32, S} where S<:AbstractMatrix{Float32}, ::LinearAlgebra.UnitUpperTriangular{Float64, S} where S<:AbstractMatrix{Float64}, α, β)/α = true, β = 0.946568898140176/β = 0 ignores C .= NaN @test ≈(collect(returned_mat), exp_val, rtol = rtol, atol = atol)