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

function opti3(fM, qn, nFlows)
    fsource = zeros(Int, nFlows)
    ftarget = zeros(Int, nFlows)
    flows = zeros(eltype(fM), nFlows)
    u = 1
    for p in 1:qn
        for q in 1:qn
            if fM[p, q] > 0
                fsource[u] = p
                ftarget[u] = q
                flows[u] = fM[p, q]
                u += 1
            end
        end
    end
    return fsource, ftarget, flows
end