// File: /workspace/github_cpp_makefile_star2_O2/fjamour[P]icentral/src/qube.cc
// Function: muc_t::i_iteration_1(unsigned int, unsigned int, unsigned int, int, int, std::tr1::unordered_map<unsigned int, double, std::tr1::hash<unsigned int>, std::equal_to<unsigned int>, std::allocator<std::pair<unsigned int const, double>>>&)
void muc_t::i_iteration_1(node_id_t s,
                        node_id_t src,
                        node_id_t dst,
                        int d_src,
                        int d_dst,
                        tr1_map_t(double)& bc_map)
{
        vector<vector<node_id_t> >&  P = muc_fast_subgraph.P; 
        vector<int>&                 path_cnt_map = muc_fast_subgraph.path_cnt_vec;
        vector<int>&                 new_path_cnt_map = muc_fast_subgraph.new_path_cnt_vec;
        vector<int>&                 dist_map = muc_fast_subgraph.dist_vec;
        vector<double>&              pair_dep_map = muc_fast_subgraph.pair_dep_vec;
        vector<double>&              new_pair_dep_map = muc_fast_subgraph.new_pair_dep_vec;
        vector<double>&              sigma_t_map = muc_fast_subgraph.sigma_t_map;
        vector<double>&              new_sigma_t_map = muc_fast_subgraph.new_sigma_t_map; 
        queue<node_id_t>&            Q = muc_fast_subgraph.Q;
        stack<node_id_t>             S;
        
        muc_fast_subgraph.init_maps();           

        path_cnt_map[s] = 1;
        new_path_cnt_map[s] = 1;//NEW
        dist_map[s] = 0;
        
        Q.push(s);
        while(!Q.empty()) {
            node_id_t v_i = Q.front(); Q.pop();
            S.push(v_i);
            if(v_i == dst) {//NEW
                new_path_cnt_map[v_i] += path_cnt_map[src];
            }
            for(int i = 0; i < muc_fast_subgraph.nodes_vec[v_i].size(); ++i) {
                node_id_t v_n = muc_fast_subgraph.nodes_vec[v_i][i];
                if(dist_map[v_n] < 0) {
                    Q.push(v_n);
                    dist_map[v_n] = dist_map[v_i] + 1;
                }
                if(dist_map[v_n] == dist_map[v_i] + 1) {
                    path_cnt_map[v_n] = path_cnt_map[v_n] + path_cnt_map[v_i];
                    new_path_cnt_map[v_n] += new_path_cnt_map[v_i];//NEW
                    P[v_n].push_back(v_i);
                }
            }
        }
        
        while(!S.empty()) {
            node_id_t v_n = S.top(); S.pop();
            if(   tmp_conn_vertex_map.find(s) != tmp_conn_vertex_map.end()
               && tmp_conn_vertex_map.find(v_n) != tmp_conn_vertex_map.end()
               && s != v_n) {
                int VG_s, VG_n;
                VG_s = tmp_subgraph_map[s].size();
                VG_n = tmp_subgraph_map[v_n].size();
                int c_t = VG_s*VG_n;
                sigma_t_map[v_n] =  sigma_t_map[v_n] + c_t;
                new_sigma_t_map[v_n] = new_sigma_t_map[v_n] + c_t;//NEW
                /*this guy must not change!*/
                //bc_map[v_n] = bc_map[v_n] - c_t;
                //bc_map[v_n] = bc_map[v_n] + c_t;
            }
            for(int i = 0; i < P[v_n].size(); ++i) {
                node_id_t v_p = P[v_n][i];
                double sp_sn = ((double)path_cnt_map[v_p]/path_cnt_map[v_n]);
                pair_dep_map[v_p] = pair_dep_map[v_p] + sp_sn*(1+pair_dep_map[v_n]);
                double new_sp_sn = ((double)new_path_cnt_map[v_p]/new_path_cnt_map[v_n]);
                new_pair_dep_map[v_p] = new_pair_dep_map[v_p] + new_sp_sn*(1+new_pair_dep_map[v_n]);
                if(tmp_conn_vertex_map.find(s) != tmp_conn_vertex_map.end()) {
                    sigma_t_map[v_p] = sigma_t_map[v_p] + sigma_t_map[v_n]*sp_sn;
                    new_sigma_t_map[v_p] = new_sigma_t_map[v_p] + new_sigma_t_map[v_n]*new_sp_sn;
                    
                    node_id_t new_v_p = muc_fast_subgraph.inout_label_map[v_p];
                    bc_map[new_v_p] = bc_map[new_v_p] - sigma_t_map[v_n]*sp_sn/2.0;
                    bc_map[new_v_p] = bc_map[new_v_p] + new_sigma_t_map[v_n]*new_sp_sn/2.0;
                }
            }
            //IMP: this is the only change that happens to P, @src should be added as parent for dst
            if(v_n == dst) {
                node_id_t v_p = src;
                double new_sp_sn = ((double)new_path_cnt_map[v_p]/new_path_cnt_map[v_n]);
                new_pair_dep_map[v_p] = new_pair_dep_map[v_p] + new_sp_sn*(1+new_pair_dep_map[v_n]);
                if(tmp_conn_vertex_map.find(s) != tmp_conn_vertex_map.end()) {
                    new_sigma_t_map[v_p] = new_sigma_t_map[v_p] + new_sigma_t_map[v_n]*new_sp_sn;
                    
                    node_id_t new_v_p = muc_fast_subgraph.inout_label_map[v_p];
                    bc_map[new_v_p] = bc_map[new_v_p] + new_sigma_t_map[v_n]*new_sp_sn/2.0;
                }
            }
            if(s != v_n) {
                node_id_t new_v_n = muc_fast_subgraph.inout_label_map[v_n];
                bc_map[new_v_n] = bc_map[new_v_n] - pair_dep_map[v_n]/2.0;
                bc_map[new_v_n] = bc_map[new_v_n] + new_pair_dep_map[v_n]/2.0;
            }
            if(tmp_conn_vertex_map.find(s) != tmp_conn_vertex_map.end()) {
                int VG_s = tmp_subgraph_map[s].size();
                
                node_id_t new_v_n = muc_fast_subgraph.inout_label_map[v_n];
                bc_map[new_v_n] = bc_map[new_v_n] - pair_dep_map[v_n]*VG_s;
                bc_map[new_v_n] = bc_map[new_v_n] + new_pair_dep_map[v_n]*VG_s;
            }
        }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x98, %rsp
movl %ecx, 0x2c(%rsp)
movl %edx, %ebx
movl %esi, %ebp
movq %rdi, %r14
movl %esi, 0xc(%rsp)
leaq 0x310(%rdi), %r12
leaq 0x48(%rsp), %rdi
callq 0x766a
movq %r12, %rdi
callq 0x25fea
leaq 0x4b8(%r14), %r13
movl %ebp, %eax
movq 0x3b8(%r14), %rcx
pushq $0x1
popq %rdx
movl %edx, (%rcx,%rax,4)
leaq 0xc(%rsp), %rsi
movl (%rsi), %eax
movq 0x3d0(%r14), %rcx
movl %edx, (%rcx,%rax,4)
movl (%rsi), %eax
movq 0x3e8(%r14), %rcx
andl $0x0, (%rcx,%rax,4)
movq %r13, %rdi
callq 0x62da
movq 0xd0(%rsp), %r15
movl %ebx, %eax
movq %rax, 0x38(%rsp)
movl 0x2c(%rsp), %eax
movq %rax, 0x10(%rsp)
leaq 0x4(%rsp), %rbx
movq 0x4c8(%r14), %rax
cmpq %rax, 0x4e8(%r14)
je 0x21b32
movl (%rax), %eax
movl %eax, 0x4(%rsp)
movq %r13, %rdi
callq 0x650e
leaq 0x48(%rsp), %rdi
movq %rbx, %rsi
callq 0x62da
movl 0x2c(%rsp), %eax
cmpl %eax, 0x4(%rsp)
jne 0x21a81
movq 0x3b8(%r14), %rax
movq 0x3d0(%r14), %rcx
movq 0x38(%rsp), %rdx
movl (%rax,%rdx,4), %eax
movq 0x10(%rsp), %rdx
addl %eax, (%rcx,%rdx,4)
xorl %ebp, %ebp
movl 0x4(%rsp), %eax
movq (%r12), %rdx
imulq $0x18, %rax, %rsi
movq (%rdx,%rsi), %rcx
movq 0x8(%rdx,%rsi), %rdx
subq %rcx, %rdx
sarq $0x2, %rdx
cmpq %rbp, %rdx
jbe 0x21a2a
movl (%rcx,%rbp,4), %ecx
movl %ecx, 0x8(%rsp)
movq 0x3e8(%r14), %rdx
movl (%rdx,%rcx,4), %esi
testl %esi, %esi
jns 0x21ae8
movq %r13, %rdi
leaq 0x8(%rsp), %rsi
callq 0x62da
movl 0x4(%rsp), %eax
movq 0x3e8(%r14), %rdx
movl (%rdx,%rax,4), %eax
incl %eax
movl 0x8(%rsp), %ecx
movl %eax, (%rdx,%rcx,4)
movl 0x8(%rsp), %ecx
movl (%rdx,%rcx,4), %esi
movl 0x4(%rsp), %eax
movl (%rdx,%rax,4), %edx
incl %edx
cmpl %edx, %esi
jne 0x21b2a
movq 0x3b8(%r14), %rdx
movl (%rdx,%rax,4), %eax
addl %eax, (%rdx,%rcx,4)
movl 0x4(%rsp), %eax
movq 0x3d0(%r14), %rcx
movl (%rcx,%rax,4), %eax
movl 0x8(%rsp), %edx
addl %eax, (%rcx,%rdx,4)
movl 0x8(%rsp), %eax
imulq $0x18, %rax, %rdi
addq 0x3a0(%r14), %rdi
movq %rbx, %rsi
callq 0x478c
incq %rbp
jmp 0x21a83
leaq 0x508(%r14), %r12
leaq 0x538(%r14), %rax
movq %rax, 0x30(%rsp)
leaq 0x8(%rsp), %r13
movq 0x78(%rsp), %rax
cmpq 0x58(%rsp), %rax
je 0x220c1
cmpq 0x80(%rsp), %rax
jne 0x21b78
movq 0x90(%rsp), %rax
movq -0x8(%rax), %rax
movl $0x200, %ecx            # imm = 0x200
addq %rcx, %rax
movl -0x4(%rax), %eax
movl %eax, 0x4(%rsp)
leaq 0x48(%rsp), %rdi
callq 0x85f6
movq %r12, %rdi
leaq 0xc(%rsp), %rsi
callq 0x22ed2
movq 0x510(%r14), %rcx
movq 0x518(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0x21c3b
movq %r12, %rdi
leaq 0x4(%rsp), %rsi
callq 0x22ed2
movq 0x510(%r14), %rcx
movq 0x518(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0x21c3b
movl 0xc(%rsp), %eax
cmpl 0x4(%rsp), %eax
je 0x21c3b
movq 0x30(%rsp), %rdi
leaq 0xc(%rsp), %rsi
callq 0x22dc4
movq %rax, %rdi
callq 0x176f6
movq %rax, %rbx
movq 0x30(%rsp), %rdi
leaq 0x4(%rsp), %rsi
callq 0x22dc4
movq %rax, %rdi
callq 0x176f6
imull %ebx, %eax
movl 0x4(%rsp), %ecx
cvtsi2sd %eax, %xmm0
movq 0x448(%r14), %rax
movsd (%rax,%rcx,8), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax,%rcx,8)
movq 0x460(%r14), %rax
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
xorl %ebx, %ebx
movl 0x4(%rsp), %ecx
movq 0x3a0(%r14), %rdx
imulq $0x18, %rcx, %rsi
movq (%rdx,%rsi), %rax
movq 0x8(%rdx,%rsi), %rdx
subq %rax, %rdx
sarq $0x2, %rdx
cmpq %rbx, %rdx
jbe 0x21e09
movl (%rax,%rbx,4), %ebp
movq 0x3b8(%r14), %rax
movq 0x400(%r14), %rdx
cvtsi2sdl (%rax,%rbp,4), %xmm2
cvtsi2sdl (%rax,%rcx,4), %xmm0
divsd %xmm0, %xmm2
movsd (%rdx,%rcx,8), %xmm0
movsd 0x63af(%rip), %xmm1     # 0x28040
addsd %xmm1, %xmm0
movsd %xmm2, 0x10(%rsp)
mulsd %xmm2, %xmm0
addsd (%rdx,%rbp,8), %xmm0
movsd %xmm0, (%rdx,%rbp,8)
movq 0x3d0(%r14), %rax
cvtsi2sdl (%rax,%rbp,4), %xmm2
cvtsi2sdl (%rax,%rcx,4), %xmm0
movq 0x418(%r14), %rax
divsd %xmm0, %xmm2
movsd (%rax,%rcx,8), %xmm0
addsd %xmm1, %xmm0
movsd %xmm2, 0x18(%rsp)
mulsd %xmm2, %xmm0
addsd (%rax,%rbp,8), %xmm0
movsd %xmm0, (%rax,%rbp,8)
movq %r12, %rdi
leaq 0xc(%rsp), %rsi
callq 0x22ed2
movq 0x510(%r14), %rcx
movq 0x518(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0x21e01
movq 0x448(%r14), %rax
movl 0x4(%rsp), %ecx
movsd (%rax,%rcx,8), %xmm0
mulsd 0x10(%rsp), %xmm0
addsd (%rax,%rbp,8), %xmm0
movsd %xmm0, (%rax,%rbp,8)
movq 0x460(%r14), %rax
movsd (%rax,%rcx,8), %xmm0
mulsd 0x18(%rsp), %xmm0
addsd (%rax,%rbp,8), %xmm0
movsd %xmm0, (%rax,%rbp,8)
movq 0x358(%r14), %rax
movl (%rax,%rbp,4), %eax
movl %eax, 0x8(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x20(%rsp)
movl 0x4(%rsp), %eax
movq 0x448(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x40(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd 0x10(%rsp), %xmm0
mulsd 0x40(%rsp), %xmm0
mulsd 0x629d(%rip), %xmm0     # 0x28038
movsd 0x20(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x4(%rsp), %eax
movq 0x460(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x20(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
mulsd 0x623d(%rip), %xmm0     # 0x28030
movsd 0x10(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax)
incq %rbx
jmp 0x21c3d
cmpl 0x2c(%rsp), %ecx
jne 0x21f09
movq 0x3d0(%r14), %rax
movq 0x38(%rsp), %rdx
cvtsi2sdl (%rax,%rdx,4), %xmm1
cvtsi2sdl (%rax,%rcx,4), %xmm0
movq 0x418(%r14), %rax
divsd %xmm0, %xmm1
movsd (%rax,%rcx,8), %xmm0
addsd 0x61ff(%rip), %xmm0     # 0x28040
movsd %xmm1, 0x10(%rsp)
mulsd %xmm1, %xmm0
addsd (%rax,%rdx,8), %xmm0
movsd %xmm0, (%rax,%rdx,8)
movq %r12, %rdi
leaq 0xc(%rsp), %rsi
callq 0x22ed2
movq 0x510(%r14), %rdx
movq 0x518(%r14), %rsi
movl 0x4(%rsp), %ecx
cmpq (%rdx,%rsi,8), %rax
je 0x21f09
movq 0x460(%r14), %rax
movsd (%rax,%rcx,8), %xmm0
mulsd 0x10(%rsp), %xmm0
movq 0x38(%rsp), %rcx
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq 0x358(%r14), %rax
movl (%rax,%rcx,4), %eax
movl %eax, 0x8(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x18(%rsp)
movl 0x4(%rsp), %eax
movq 0x460(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x20(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd 0x10(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
mulsd 0x6139(%rip), %xmm0     # 0x28030
movsd 0x18(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax)
movl 0x4(%rsp), %ecx
cmpl %ecx, 0xc(%rsp)
je 0x21fc7
movl %ecx, %eax
movq 0x358(%r14), %rcx
movl (%rcx,%rax,4), %eax
movl %eax, 0x8(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x4(%rsp), %eax
movq 0x400(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm1
mulsd 0x60d1(%rip), %xmm1     # 0x28038
movsd 0x10(%rsp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x4(%rsp), %eax
movq 0x418(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm1
mulsd 0x6077(%rip), %xmm1     # 0x28030
movsd 0x10(%rsp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movq %r12, %rdi
leaq 0xc(%rsp), %rsi
callq 0x22ed2
movq 0x510(%r14), %rcx
movq 0x518(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0x21b4a
movq 0x30(%rsp), %rdi
leaq 0xc(%rsp), %rsi
callq 0x22dc4
movq %rax, %rdi
callq 0x176f6
movq %rax, %rbx
movl 0x4(%rsp), %eax
movq 0x358(%r14), %rcx
movl (%rcx,%rax,4), %eax
movl %eax, 0x8(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x4(%rsp), %eax
movq 0x400(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
cvtsi2sd %ebx, %xmm0
movsd %xmm0, 0x20(%rsp)
movsd 0x18(%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x4(%rsp), %eax
movq 0x418(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %r13, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
movsd %xmm0, (%rax)
jmp 0x21b4a
leaq 0x48(%rsp), %rdi
callq 0x5280
addq $0x98, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
jmp 0x220f7
movq %rax, %rbx
leaq 0x48(%rsp), %rdi
callq 0x5280
movq %rbx, %rdi
callq 0x23d0