SuiteSparseGraphBLAS.jl

SuiteSparseGraphBLAS.jl copied to clipboard

function zero_out(mat, bus_index)
    mat[bus_index, :] .= 0
    mat[:, bus_index] .= 0
    mat[bus_index, bus_index] = -1
    return mat
end

For the above code, GBMatrix is > 13.5 times slower than the default sparse matrix for a 20k * 20k sparse matrix with 0.018% sparsity.