// File: /workspace/github_cpp_makefile_star2_O2/fjamour[P]icentral/src/bcc_delta.cc
// Function: bcc_delta_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 bcc_delta_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 = bcc_fast_subgraph.P; 
        vector<int>&                 path_cnt_map = bcc_fast_subgraph.path_cnt_vec;
        vector<int>&                 new_path_cnt_map = bcc_fast_subgraph.new_path_cnt_vec;
        vector<int>&                 dist_map = bcc_fast_subgraph.dist_vec;
        vector<double>&              pair_dep_map = bcc_fast_subgraph.pair_dep_vec;
        vector<double>&              new_pair_dep_map = bcc_fast_subgraph.new_pair_dep_vec;
        vector<double>&              sigma_t_map = bcc_fast_subgraph.sigma_t_map;
        vector<double>&              new_sigma_t_map = bcc_fast_subgraph.new_sigma_t_map; 
        queue<node_id_t>&            Q = bcc_fast_subgraph.Q;
        stack<node_id_t>             S;
        bcc_fast_subgraph.init_maps();
        
//    int sz = bcc_fast_subgraph.nodes_vec.size();
//        vector<vector<node_id_t> >  P(sz, vector<node_id_t>());// = bcc_fast_subgraph.P; 
//        vector<int>                 path_cnt_map(sz, 0);// = bcc_fast_subgraph.path_cnt_vec;
//        vector<int>                 new_path_cnt_map(sz, 0);// = bcc_fast_subgraph.new_path_cnt_vec;
//        vector<int>                 dist_map(sz, -1);// = bcc_fast_subgraph.dist_vec;
//        vector<double>              pair_dep_map(sz, 0);// = bcc_fast_subgraph.pair_dep_vec;
//        vector<double>              new_pair_dep_map(sz, 0);// = bcc_fast_subgraph.new_pair_dep_vec;
//        vector<double>              sigma_t_map(sz, 0);// = bcc_fast_subgraph.sigma_t_map;
//        vector<double>              new_sigma_t_map(sz, 0);// = bcc_fast_subgraph.new_sigma_t_map; 
//        queue<node_id_t>            Q;// = bcc_fast_subgraph.Q;
//        stack<node_id_t>            S;
        
        
        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 < bcc_fast_subgraph.nodes_vec[v_i].size(); ++i) {
                node_id_t v_n = bcc_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(   art_pt_map.find(s) != art_pt_map.end()
               && art_pt_map.find(v_n) != art_pt_map.end()
               && s != v_n) {
                int VG_s, VG_n;
                VG_s = accumulate(art_pt_map[s].begin(), art_pt_map[s].end(), 0);
                VG_n = accumulate(art_pt_map[v_n].begin(), art_pt_map[v_n].end(), 0);
                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(art_pt_map.find(s) != art_pt_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 = bcc_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(art_pt_map.find(s) != art_pt_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 = bcc_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 = bcc_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(art_pt_map.find(s) != art_pt_map.end()) {
                int VG_s = accumulate(art_pt_map[s].begin(), art_pt_map[s].end(), 0);
                
                node_id_t new_v_n = bcc_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, 0x34(%rsp)
movl %edx, %ebx
movl %esi, %ebp
movq %rdi, %r14
movl %esi, 0x24(%rsp)
leaq 0x2d8(%rdi), %r12
leaq 0x48(%rsp), %rdi
callq 0x766a
movq %r12, %rdi
callq 0x25fea
leaq 0x480(%r14), %r13
movl %ebp, %eax
movq 0x380(%r14), %rcx
pushq $0x1
popq %rdx
movl %edx, (%rcx,%rax,4)
leaq 0x24(%rsp), %rsi
movl (%rsi), %eax
movq 0x398(%r14), %rcx
movl %edx, (%rcx,%rax,4)
movl (%rsi), %eax
movq 0x3b0(%r14), %rcx
andl $0x0, (%rcx,%rax,4)
movq %r13, %rdi
callq 0x62da
movq 0xd0(%rsp), %r15
movl %ebx, %eax
movq %rax, 0x38(%rsp)
movl 0x34(%rsp), %eax
movq %rax, 0x10(%rsp)
leaq 0x8(%rsp), %rbx
movq 0x490(%r14), %rax
cmpq %rax, 0x4b0(%r14)
je 0xd3d8
movl (%rax), %eax
movl %eax, 0x8(%rsp)
movq %r13, %rdi
callq 0x650e
leaq 0x48(%rsp), %rdi
movq %rbx, %rsi
callq 0x62da
movl 0x34(%rsp), %eax
cmpl %eax, 0x8(%rsp)
jne 0xd327
movq 0x380(%r14), %rax
movq 0x398(%r14), %rcx
movq 0x38(%rsp), %rdx
movl (%rax,%rdx,4), %eax
movq 0x10(%rsp), %rdx
addl %eax, (%rcx,%rdx,4)
xorl %ebp, %ebp
movl 0x8(%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 0xd2d0
movl (%rcx,%rbp,4), %ecx
movl %ecx, 0xc(%rsp)
movq 0x3b0(%r14), %rdx
movl (%rdx,%rcx,4), %esi
testl %esi, %esi
jns 0xd38e
movq %r13, %rdi
leaq 0xc(%rsp), %rsi
callq 0x62da
movl 0x8(%rsp), %eax
movq 0x3b0(%r14), %rdx
movl (%rdx,%rax,4), %eax
incl %eax
movl 0xc(%rsp), %ecx
movl %eax, (%rdx,%rcx,4)
movl 0xc(%rsp), %ecx
movl (%rdx,%rcx,4), %esi
movl 0x8(%rsp), %eax
movl (%rdx,%rax,4), %edx
incl %edx
cmpl %edx, %esi
jne 0xd3d0
movq 0x380(%r14), %rdx
movl (%rdx,%rax,4), %eax
addl %eax, (%rdx,%rcx,4)
movl 0x8(%rsp), %eax
movq 0x398(%r14), %rcx
movl (%rcx,%rax,4), %eax
movl 0xc(%rsp), %edx
addl %eax, (%rcx,%rdx,4)
movl 0xc(%rsp), %eax
imulq $0x18, %rax, %rdi
addq 0x368(%r14), %rdi
movq %rbx, %rsi
callq 0x478c
incq %rbp
jmp 0xd329
leaq 0x24(%rsp), %r13
leaq 0xc(%rsp), %rbp
movq 0x78(%rsp), %rax
cmpq 0x58(%rsp), %rax
je 0xd98a
cmpq 0x80(%rsp), %rax
jne 0xd410
movq 0x90(%rsp), %rax
movq -0x8(%rax), %rax
movl $0x200, %ecx            # imm = 0x200
addq %rcx, %rax
movl -0x4(%rax), %eax
movl %eax, 0x8(%rsp)
leaq 0x48(%rsp), %rdi
callq 0x85f6
movq %r14, %rdi
movq %r13, %rsi
callq 0x4a5c
movq 0x8(%r14), %rcx
movq 0x10(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0xd4fc
movq %r14, %rdi
leaq 0x8(%rsp), %rsi
callq 0x4a5c
movq 0x8(%r14), %rcx
movq 0x10(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0xd4fc
movl 0x24(%rsp), %eax
cmpl 0x8(%rsp), %eax
je 0xd4fc
movq %r14, %rdi
movq %r13, %rsi
callq 0x4aaa
movq (%rax), %r12
movq %r14, %rdi
movq %r13, %rsi
callq 0x4aaa
movq 0x8(%rax), %rax
xorl %ebx, %ebx
cmpq %rax, %r12
je 0xd499
addl (%r12), %ebx
addq $0x4, %r12
jmp 0xd48a
movq %r14, %rdi
leaq 0x8(%rsp), %rsi
callq 0x4aaa
movq (%rax), %r12
movq %r14, %rdi
leaq 0x8(%rsp), %rsi
callq 0x4aaa
movq 0x8(%rax), %rcx
xorl %eax, %eax
cmpq %rcx, %r12
je 0xd4cb
addl (%r12), %eax
addq $0x4, %r12
jmp 0xd4bc
imull %ebx, %eax
movl 0x8(%rsp), %ecx
cvtsi2sd %eax, %xmm0
movq 0x410(%r14), %rax
movsd (%rax,%rcx,8), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax,%rcx,8)
movq 0x428(%r14), %rax
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
xorl %ebx, %ebx
movl 0x8(%rsp), %ecx
movq 0x368(%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 0xd6ce
movl (%rax,%rbx,4), %r12d
movq 0x380(%r14), %rax
movq 0x3c8(%r14), %rdx
cvtsi2sdl (%rax,%r12,4), %xmm2
cvtsi2sdl (%rax,%rcx,4), %xmm0
divsd %xmm0, %xmm2
movsd (%rdx,%rcx,8), %xmm0
movsd 0x1aaec(%rip), %xmm1    # 0x28040
addsd %xmm1, %xmm0
movsd %xmm2, 0x10(%rsp)
mulsd %xmm2, %xmm0
addsd (%rdx,%r12,8), %xmm0
movsd %xmm0, (%rdx,%r12,8)
movq 0x398(%r14), %rax
cvtsi2sdl (%rax,%r12,4), %xmm2
cvtsi2sdl (%rax,%rcx,4), %xmm0
movq 0x3e0(%r14), %rax
divsd %xmm0, %xmm2
movsd (%rax,%rcx,8), %xmm0
addsd %xmm1, %xmm0
movsd %xmm2, 0x18(%rsp)
mulsd %xmm2, %xmm0
addsd (%rax,%r12,8), %xmm0
movsd %xmm0, (%rax,%r12,8)
movq %r14, %rdi
movq %r13, %rsi
callq 0x4a5c
movq 0x8(%r14), %rcx
movq 0x10(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0xd6c6
movq 0x410(%r14), %rax
movl 0x8(%rsp), %ecx
movsd (%rax,%rcx,8), %xmm0
mulsd 0x10(%rsp), %xmm0
addsd (%rax,%r12,8), %xmm0
movsd %xmm0, (%rax,%r12,8)
movq 0x428(%r14), %rax
movsd (%rax,%rcx,8), %xmm0
mulsd 0x18(%rsp), %xmm0
addsd (%rax,%r12,8), %xmm0
movsd %xmm0, (%rax,%r12,8)
movq 0x320(%r14), %rax
movl (%rax,%r12,4), %eax
movl %eax, 0xc(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x28(%rsp)
movl 0x8(%rsp), %eax
movq 0x410(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x40(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd 0x10(%rsp), %xmm0
mulsd 0x40(%rsp), %xmm0
mulsd 0x1a9d8(%rip), %xmm0    # 0x28038
movsd 0x28(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x8(%rsp), %eax
movq 0x428(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x28(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm0
mulsd 0x28(%rsp), %xmm0
mulsd 0x1a978(%rip), %xmm0    # 0x28030
movsd 0x10(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax)
incq %rbx
jmp 0xd4fe
cmpl 0x34(%rsp), %ecx
jne 0xd7c6
movq 0x398(%r14), %rax
movq 0x38(%rsp), %rdx
cvtsi2sdl (%rax,%rdx,4), %xmm1
cvtsi2sdl (%rax,%rcx,4), %xmm0
movq 0x3e0(%r14), %rax
divsd %xmm0, %xmm1
movsd (%rax,%rcx,8), %xmm0
addsd 0x1a93a(%rip), %xmm0    # 0x28040
movsd %xmm1, 0x10(%rsp)
mulsd %xmm1, %xmm0
addsd (%rax,%rdx,8), %xmm0
movsd %xmm0, (%rax,%rdx,8)
movq %r14, %rdi
movq %r13, %rsi
callq 0x4a5c
movq 0x8(%r14), %rdx
movq 0x10(%r14), %rsi
movl 0x8(%rsp), %ecx
cmpq (%rdx,%rsi,8), %rax
je 0xd7c6
movq 0x428(%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 0x320(%r14), %rax
movl (%rax,%rcx,4), %eax
movl %eax, 0xc(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x18(%rsp)
movl 0x8(%rsp), %eax
movq 0x428(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x28(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd 0x10(%rsp), %xmm0
mulsd 0x28(%rsp), %xmm0
mulsd 0x1a87c(%rip), %xmm0    # 0x28030
movsd 0x18(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, (%rax)
movl 0x8(%rsp), %ecx
cmpl %ecx, 0x24(%rsp)
je 0xd884
movl %ecx, %eax
movq 0x320(%r14), %rcx
movl (%rcx,%rax,4), %eax
movl %eax, 0xc(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x8(%rsp), %eax
movq 0x3c8(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm1
mulsd 0x1a814(%rip), %xmm1    # 0x28038
movsd 0x10(%rsp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x8(%rsp), %eax
movq 0x3e0(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm1
mulsd 0x1a7ba(%rip), %xmm1    # 0x28030
movsd 0x10(%rsp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movq %r14, %rdi
movq %r13, %rsi
callq 0x4a5c
movq 0x8(%r14), %rcx
movq 0x10(%r14), %rdx
cmpq (%rcx,%rdx,8), %rax
je 0xd3e2
movq %r14, %rdi
movq %r13, %rsi
callq 0x4aaa
movq (%rax), %r12
movq %r14, %rdi
movq %r13, %rsi
callq 0x4aaa
movq 0x8(%rax), %rax
xorl %ebx, %ebx
cmpq %rax, %r12
je 0xd8cf
addl (%r12), %ebx
addq $0x4, %r12
jmp 0xd8c0
movl 0x8(%rsp), %eax
movq 0x320(%r14), %rcx
movl (%rcx,%rax,4), %eax
movl %eax, 0xc(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x8(%rsp), %eax
movq 0x3c8(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
cvtsi2sd %ebx, %xmm0
movsd %xmm0, 0x28(%rsp)
movsd 0x18(%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd (%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movl 0x8(%rsp), %eax
movq 0x3e0(%r14), %rcx
movsd (%rcx,%rax,8), %xmm0
movsd %xmm0, 0x18(%rsp)
movq %r15, %rdi
movq %rbp, %rsi
callq 0xb00c
movsd 0x18(%rsp), %xmm0
mulsd 0x28(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
movsd %xmm0, (%rax)
jmp 0xd3e2
leaq 0x48(%rsp), %rdi
callq 0x5280
addq $0x98, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
jmp 0xd9c0
movq %rax, %rbx
leaq 0x48(%rsp), %rdi
callq 0x5280
movq %rbx, %rdi
callq 0x23d0
nop