Оптимизация СМО

function opti1(nM, qn, nEdges)
    target = zeros(Int, nEdges)
    source = zeros(Int, nEdges)  
    Mu = zeros(Float64, nEdges)
    u = 1
    
    for p in 1:qn
        for q in 1:qn
            if nM[p, q] > 0
                source[u] = p
                target[u] = q
                Mu[u] = nM[p, q]
                u += 1
            end
        end
    end
    
    return source, target, Mu
end