LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/channel.cc
// Function: CHANNEL::findBoundingAtoms(ATOM_NETWORK*, std::vector<BASIC_VCELL, std::allocator<BASIC_VCELL>>&, std::vector<int, std::allocator<int>>&)
void CHANNEL::findBoundingAtoms(ATOM_NETWORK *atmnet, vector<BASIC_VCELL> &vcells,
                                vector<int> &atomIDs)
{
    atomIDs.clear();
    
    for(unsigned int i = 0; i < vcells.size(); i++){
        BASIC_VCELL cell = vcells[i];
        for(int j = 0; j < cell.getNumNodes(); j++){
            if(idMappings.find(cell.getNodeID(j)) != idMappings.end()){
                atomIDs.push_back(i);
                break;
            }
        }
    }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movq %rcx, %r14
movq %rdx, 0x8(%rsp)
movq %rdi, %r15
movq (%rcx), %rax
cmpq %rax, 0x8(%rcx)
je 0x7498a
movq %rax, 0x8(%r14)
movq 0x8(%rsp), %rax
movq (%rax), %rbp
cmpq %rbp, 0x8(%rax)
je 0x74abe
leaq 0x8(%r15), %r12
xorl %ebx, %ebx
leaq 0x18(%rsp), %r13
movq %r15, 0x10(%rsp)
movq %rbx, %rax
movq %r14, %rbx
movq %rax, %r15
leaq (%rax,%rax,2), %r14
shlq $0x4, %r14
leaq (%r14,%rbp), %rsi
movq %r13, %rdi
callq 0xe7f6
leaq (%r14,%rbp), %rsi
addq $0x18, %rsi
leaq 0x30(%rsp), %rdi
callq 0xe886
xorl %ebp, %ebp
movq %rbx, %r14
movq %r15, %rbx
movq 0x10(%rsp), %r15
movq %r13, %rdi
callq 0x6c80a
cmpl %eax, %ebp
jge 0x74a64
movq %r13, %rdi
movl %ebp, %esi
callq 0x6c83a
movq 0x10(%r15), %rdx
movq %r12, %rcx
testq %rdx, %rdx
je 0x74a24
movq %r12, %rcx
xorl %esi, %esi
cmpl %eax, 0x20(%rdx)
setl %sil
cmovgeq %rdx, %rcx
movq 0x10(%rdx,%rsi,8), %rdx
testq %rdx, %rdx
jne 0x74a0d
movq %r12, %rdx
cmpq %r12, %rcx
je 0x74a36
cmpl 0x20(%rcx), %eax
cmovlq %r12, %rcx
movq %rcx, %rdx
incl %ebp
cmpq %r12, %rdx
je 0x749e8
movl %ebx, 0x4(%rsp)
movq 0x8(%r14), %rsi
cmpq 0x10(%r14), %rsi
je 0x74a57
movl %ebx, (%rsi)
addq $0x4, %rsi
movq %rsi, 0x8(%r14)
jmp 0x74a64
movq %r14, %rdi
leaq 0x4(%rsp), %rdx
callq 0x1422c
movq 0x30(%rsp), %rdi
testq %rdi, %rdi
je 0x74a7b
movq 0x40(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
movq 0x18(%rsp), %rdi
testq %rdi, %rdi
je 0x74a92
movq 0x28(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
incl %ebx
movq 0x8(%rsp), %rax
movq (%rax), %rbp
movq 0x8(%rax), %rax
subq %rbp, %rax
sarq $0x4, %rax
movabsq $-0x5555555555555555, %rcx # imm = 0xAAAAAAAAAAAAAAAB
imulq %rcx, %rax
cmpq %rbx, %rax
ja 0x749ac
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0x74aed
movq %rax, %rbx
movq 0x18(%rsp), %rdi
testq %rdi, %rdi
je 0x74afa
movq 0x28(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
jmp 0x74afa
jmp 0x74aed
movq %rax, %rbx
leaq 0x18(%rsp), %rdi
callq 0xf420
movq %rbx, %rdi
callq 0x4780
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/channel.cc
// Function: CHANNEL::findFreeIncludedSphereDiameter()
pair<double, pair<double,double> > CHANNEL::findFreeIncludedSphereDiameter(){
    
    pair<double, pair<double,double> > maxdidfdif;
    maxdidfdif.first=0;
    maxdidfdif.second.first=0;
    maxdidfdif.second.second=0;
    
    
    for(int i=0;i<nodes.size();i++)
    {
        // loop over all the nodes in the current channel
        if(i==0) {
            maxdidfdif = findFreeIncludedSphereDiameterforNode(i, maxdidfdif);
        }
        else
        {
            pair<double, pair<double,double> > didfdif = findFreeIncludedSphereDiameterforNode(i, maxdidfdif);
            if(didfdif.second.first > maxdidfdif.second.first) maxdidfdif.second = didfdif.second;
            if(didfdif.first > maxdidfdif.first) maxdidfdif.first = didfdif.first;
        };
    };
    
    return maxdidfdif;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x40, %rsp
movq %rdi, %rbx
xorpd %xmm0, %xmm0
movupd %xmm0, (%rdi)
movq $0x0, 0x10(%rdi)
movq 0x68(%rsi), %rax
cmpq 0x60(%rsi), %rax
je 0x74c86
movq %rsi, %r14
movsd (%rbx), %xmm1
movsd 0x8(%rbx), %xmm2
xorl %r15d, %r15d
leaq 0x18(%rsp), %r12
movabsq $-0x71c71c71c71c71c7, %r13 # imm = 0x8E38E38E38E38E39
testq %r15, %r15
je 0x74c2d
movsd %xmm1, 0x38(%rsp)
movq 0x10(%rbx), %rax
movq %rax, 0x10(%rsp)
movups (%rbx), %xmm0
movups %xmm0, (%rsp)
movq %r12, %rdi
movq %r14, %rsi
movl %r15d, %edx
movsd %xmm2, 0x30(%rsp)
callq 0x74c98
movsd 0x30(%rsp), %xmm2
movsd 0x20(%rsp), %xmm0
ucomisd %xmm2, %xmm0
jbe 0x74c11
movsd %xmm0, 0x8(%rbx)
movsd 0x28(%rsp), %xmm1
movsd %xmm1, 0x10(%rbx)
movapd %xmm0, %xmm2
movsd 0x18(%rsp), %xmm0
movsd 0x38(%rsp), %xmm1
ucomisd %xmm1, %xmm0
jbe 0x74c6a
movsd %xmm0, (%rbx)
movapd %xmm0, %xmm1
jmp 0x74c6a
movq 0x10(%rbx), %rax
movq %rax, 0x10(%rsp)
movups (%rbx), %xmm0
movups %xmm0, (%rsp)
movq %r12, %rdi
movq %r14, %rsi
xorl %edx, %edx
callq 0x74c98
movsd 0x18(%rsp), %xmm1
movsd %xmm1, (%rbx)
movsd 0x20(%rsp), %xmm2
movsd %xmm2, 0x8(%rbx)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x10(%rbx)
movq 0x68(%r14), %rax
subq 0x60(%r14), %rax
incq %r15
sarq $0x3, %rax
imulq %r13, %rax
cmpq %r15, %rax
ja 0x74bbc
movq %rbx, %rax
addq $0x40, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/channel.cc
// Function: PORE::getIncludedSphereDiameter()
double PORE::getIncludedSphereDiameter(){
    
    double maxdi;
    
    for(unsigned int i=0; i<nodes.size(); i++)
    {
        if(i==0)
        {
            maxdi = nodes.at(i).max_radius;
        }
        else{
            if(nodes.at(i).max_radius > maxdi) maxdi = nodes.at(i).max_radius;
        };
    };
    
    maxdi=maxdi*2.0;
    
    return maxdi;
    
}

// Assembly:
movq 0x60(%rdi), %rax
movq 0x68(%rdi), %rdx
subq %rax, %rdx
je 0x75507
sarq $0x3, %rdx
movabsq $-0x71c71c71c71c71c7, %rcx # imm = 0x8E38E38E38E38E39
imulq %rdx, %rcx
addq $0x38, %rax
movl $0x1, %edx
movsd (%rax), %xmm1
cmpl $0x1, %edx
je 0x754f1
ucomisd %xmm0, %xmm1
jbe 0x754f5
movapd %xmm1, %xmm0
movl %edx, %esi
addq $0x48, %rax
incl %edx
cmpq %rsi, %rcx
ja 0x754e2
addsd %xmm0, %xmm0
retq
movsd 0xe1011(%rip), %xmm0    # 0x156520
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/mindist.cc
// Function: compute_images(double, double, double, double, double, double, std::vector<int, std::allocator<int>>&, std::vector<int, std::allocator<int>>&, std::vector<int, std::allocator<int>>&)
void compute_images(double in_bx, double in_bxy, double in_by, double in_bxz, double in_byz, double in_bz,  
		    vector<int> &i_vals, vector<int> &j_vals, vector<int> &k_vals){

  unitcell uc(in_bx, in_bxy, in_by, in_bxz, in_byz, in_bz);
  vector<int> ijk_vals; vector<double> vd;
  uc.images(ijk_vals, vd);
  unsigned int i, j;
  i_vals.clear(); j_vals.clear(); k_vals.clear();
  for(i = j = 0; i < vd.size(); i++, j+=3){
    i_vals.push_back(ijk_vals[j]);
    j_vals.push_back(ijk_vals[j+1]);
    k_vals.push_back(ijk_vals[j+2]);
  }
  
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x160, %rsp            # imm = 0x160
movq %rdx, %rbx
movq %rsi, %r14
movq %rdi, %r15
leaq 0x38(%rsp), %r12
movq %r12, %rdi
callq 0x12f184
xorl %eax, %eax
leaq 0x20(%rsp), %rsi
movq %rax, 0x10(%rsi)
xorps %xmm0, %xmm0
movaps %xmm0, (%rsi)
movq %rsp, %rdx
movq %rax, 0x10(%rdx)
movaps %xmm0, (%rdx)
movq %r12, %rdi
callq 0x1304fe
movq (%r15), %rax
cmpq %rax, 0x8(%r15)
je 0x787d5
movq %rax, 0x8(%r15)
movq (%r14), %rax
cmpq %rax, 0x8(%r14)
je 0x787e2
movq %rax, 0x8(%r14)
movq (%rbx), %rax
cmpq %rax, 0x8(%rbx)
je 0x787ef
movq %rax, 0x8(%rbx)
movq (%rsp), %rdi
cmpq %rdi, 0x8(%rsp)
je 0x788af
xorl %r12d, %r12d
movl $0x1, %ebp
movl %r12d, %edx
shlq $0x2, %rdx
addq 0x20(%rsp), %rdx
movq 0x8(%r15), %rsi
cmpq 0x10(%r15), %rsi
je 0x7882a
movl (%rdx), %eax
movl %eax, (%rsi)
addq $0x4, %rsi
movq %rsi, 0x8(%r15)
jmp 0x78832
movq %r15, %rdi
callq 0xf7b4
leal 0x1(%r12), %edx
shlq $0x2, %rdx
addq 0x20(%rsp), %rdx
movq 0x8(%r14), %rsi
cmpq 0x10(%r14), %rsi
je 0x78858
movl (%rdx), %eax
movl %eax, (%rsi)
addq $0x4, %rsi
movq %rsi, 0x8(%r14)
jmp 0x78860
movq %r14, %rdi
callq 0xf7b4
leal 0x2(%r12), %edx
shlq $0x2, %rdx
addq 0x20(%rsp), %rdx
movq 0x8(%rbx), %rsi
cmpq 0x10(%rbx), %rsi
je 0x78886
movl (%rdx), %eax
movl %eax, (%rsi)
addq $0x4, %rsi
movq %rsi, 0x8(%rbx)
jmp 0x7888e
movq %rbx, %rdi
callq 0xf7b4
addl $0x3, %r12d
movl %ebp, %eax
movq (%rsp), %rdi
movq 0x8(%rsp), %rcx
subq %rdi, %rcx
sarq $0x3, %rcx
incl %ebp
cmpq %rax, %rcx
ja 0x78806
testq %rdi, %rdi
je 0x788c1
movq 0x10(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
movq 0x20(%rsp), %rdi
testq %rdi, %rdi
je 0x788d8
movq 0x30(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
leaq 0x68(%rsp), %rdi
callq 0x122f4a
addq $0x160, %rsp            # imm = 0x160
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
jmp 0x788f4
movq %rax, %rbx
movq (%rsp), %rdi
testq %rdi, %rdi
je 0x7890d
movq 0x10(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
movq 0x20(%rsp), %rdi
testq %rdi, %rdi
je 0x78924
movq 0x30(%rsp), %rsi
subq %rdi, %rsi
callq 0x4470
leaq 0x68(%rsp), %rdi
callq 0x122f4a
movq %rbx, %rdi
callq 0x4780
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/mindist.cc
// Function: abc_to_xyz(double, double, double, double&, double&, double&)
void abc_to_xyz(double a, double b, double c, double &x, double &y, double &z){
    //cout << "bx: " << bx << "     bxy: " << bxy << "    bxz: " << bxz << endl;
    //cout << "by: " << by << "     byz: " << byz << endl;
    //cout << "bz: " << bz << endl;
  x = a*bx+b*bxy+c*bxz;
  y = b*by+c*byz;
  z = c*bz;
}

// Assembly:
xorpd %xmm3, %xmm3
mulsd %xmm3, %xmm1
mulsd %xmm3, %xmm0
addsd %xmm1, %xmm0
mulsd %xmm3, %xmm2
addsd %xmm2, %xmm0
movsd %xmm0, (%rdi)
addsd %xmm2, %xmm1
movsd %xmm1, (%rsi)
movsd %xmm2, (%rdx)
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/mindist.cc
// Function: MIN_PER_DISTANCE::print_images()
void MIN_PER_DISTANCE::print_images(){
  cout << "Printing images: " << "\n";
  for(unsigned int i = 0; i < ivals.size(); i++){
    cout << ivals[i] << " " << jvals[i] << " " << kvals[i] << "\n";
  }
  cout << "Images printed" << "\n" << "\n";
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
movq 0x11e02b(%rip), %r14    # 0x196fc0
leaq 0xddb9d(%rip), %rsi     # 0x156b39
movl $0x11, %edx
movq %r14, %rdi
callq 0x44e0
leaq 0xdc57e(%rip), %rsi     # 0x15552e
movl $0x1, %edx
movq %r14, %rdi
callq 0x44e0
movq (%rbx), %rax
cmpq %rax, 0x8(%rbx)
je 0x7905e
movl $0x1, %ebp
xorl %r14d, %r14d
leaq 0xde4ad(%rip), %r15     # 0x157486
leaq 0xdc54e(%rip), %r12     # 0x15552e
movl (%rax,%r14,4), %esi
movq 0x11dfd5(%rip), %rdi    # 0x196fc0
callq 0x4750
movq %rax, %r13
movl $0x1, %edx
movq %rax, %rdi
movq %r15, %rsi
callq 0x44e0
movq 0x18(%rbx), %rax
movl (%rax,%r14,4), %esi
movq %r13, %rdi
callq 0x4750
movq %rax, %r13
movl $0x1, %edx
movq %rax, %rdi
movq %r15, %rsi
callq 0x44e0
movq 0x30(%rbx), %rax
movl (%rax,%r14,4), %esi
movq %r13, %rdi
callq 0x4750
movl $0x1, %edx
movq %rax, %rdi
movq %r12, %rsi
callq 0x44e0
movl %ebp, %r14d
movq (%rbx), %rax
movq 0x8(%rbx), %rcx
subq %rax, %rcx
sarq $0x2, %rcx
incl %ebp
cmpq %r14, %rcx
ja 0x78fe0
movq 0x11df5b(%rip), %rbx    # 0x196fc0
leaq 0xddadf(%rip), %rsi     # 0x156b4b
movl $0xe, %edx
movq %rbx, %rdi
callq 0x44e0
leaq 0xdc4ae(%rip), %r14     # 0x15552e
movl $0x1, %edx
movq %rbx, %rdi
movq %r14, %rsi
callq 0x44e0
movl $0x1, %edx
movq %rbx, %rdi
movq %r14, %rsi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x44e0
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::segment_findMaxNode()
int FEATURE::segment_findMaxNode(){
int i;
int nd=-1;
double min=0;
for(i=0;i<(int)nodes.size();i++)
   {
   if(segments[i]==-1&&nodes.at(i).max_radius>min) 
       {
       nd=i;
       min=nodes.at(i).max_radius;
       };
   };
return nd;
}

// Assembly:
pushq %rax
movq 0x60(%rdi), %rcx
movq 0x68(%rdi), %rax
subq %rcx, %rax
sarq $0x3, %rax
movabsq $-0x71c71c71c71c71c7, %rdx # imm = 0x8E38E38E38E38E39
imulq %rax, %rdx
testl %edx, %edx
jle 0x7ba7b
movq 0x140(%rdi), %rdi
movl %edx, %r8d
andl $0x7fffffff, %r8d       # imm = 0x7FFFFFFF
addq $0x38, %rcx
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
xorpd %xmm0, %xmm0
xorl %esi, %esi
cmpl $-0x1, (%rdi,%rsi,4)
jne 0x7ba6d
cmpq %rsi, %rdx
jbe 0x7ba82
movsd (%rcx), %xmm1
ucomisd %xmm0, %xmm1
jbe 0x7ba6d
movl %esi, %eax
movapd %xmm1, %xmm0
incq %rsi
addq $0x48, %rcx
cmpq %rsi, %r8
jne 0x7ba52
jmp 0x7ba80
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
popq %rcx
retq
leaq 0xd3c1c(%rip), %rdi     # 0x14f6a5
xorl %eax, %eax
callq 0x4630
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::validateSegmentation()
int FEATURE::validateSegmentation(){
int i;
int valid=1;
for(i=0;i<nsegments;i++)
   {
   if(validateSegment(i)==0)
     {
     valid=0;
     break;
     };

   }

return valid;

}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
movl $0x1, %ebx
cmpl $0x0, 0x138(%rdi)
jle 0x7c39e
movq %rdi, %r14
xorl %ebp, %ebp
movq %r14, %rdi
movl %ebp, %esi
callq 0x7c3a6
testl %eax, %eax
je 0x7c39c
incl %ebp
cmpl 0x138(%r14), %ebp
jl 0x7c381
jmp 0x7c39e
xorl %ebx, %ebx
movl %ebx, %eax
popq %rbx
popq %r14
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::segment_growSegment_start(ATOM_NETWORK*, int)
void FEATURE::segment_growSegment_start(ATOM_NETWORK *uc, int nodeid){

segments[nodeid]=nsegments; // assign new segment number to a nw segment
segment_max_r.push_back(nodes.at(nodeid).max_radius); //store maximum r in a segment
segment_growSegment_cont(uc,nodeid,nodeid);

nsegments++;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
movl %edx, %ebp
movq %rsi, %r14
movl 0x138(%rdi), %eax
movslq %edx, %rsi
movq 0x140(%rdi), %rcx
movl %eax, (%rcx,%rsi,4)
movq 0x60(%rdi), %rax
movq 0x68(%rdi), %rcx
subq %rax, %rcx
sarq $0x3, %rcx
movabsq $-0x71c71c71c71c71c7, %rdx # imm = 0x8E38E38E38E38E39
imulq %rcx, %rdx
cmpq %rsi, %rdx
jbe 0x7c786
movq %rdi, %rbx
leaq (%rsi,%rsi,8), %rcx
leaq (%rax,%rcx,8), %rdx
addq $0x38, %rdx
movq 0x240(%rdi), %rsi
cmpq 0x248(%rdi), %rsi
je 0x7c760
movsd (%rdx), %xmm0
movsd %xmm0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x240(%rbx)
jmp 0x7c76c
leaq 0x238(%rbx), %rdi
callq 0x148c4
movq %rbx, %rdi
movq %r14, %rsi
movl %ebp, %edx
movl %ebp, %ecx
callq 0x7c794
incl 0x138(%rbx)
popq %rbx
popq %r14
popq %rbp
retq
leaq 0xd2f18(%rip), %rdi     # 0x14f6a5
xorl %eax, %eax
callq 0x4630
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::update_node_segmentinfo(std::vector<int, std::allocator<int>>*, int*, int)
void FEATURE::update_node_segmentinfo(vector <int> *segid,int *curr_id,int flag)
{
	int i;
	//loop over all nodes in channel system and update voronoi network nodes

	for(i=0;i<(int)nodes.size();i++)
	{
		if(flag==0)
		{
			segid->at(reverseIDMappings[i])=segments[i]+(*curr_id);
		}else{
			segid->at(reverseIDMappings[i])=feature_node_mapping[i]+(*curr_id);
		};
	};

	if(flag==0) {
		*curr_id+=nsegments;
	}else{
		*curr_id+=nfeatures;
	};

}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movl %ecx, %ebp
movq %rdx, 0x10(%rsp)
movq %rdi, %r14
movl $0x0, 0xc(%rsp)
movq 0x68(%rdi), %rax
subq 0x60(%rdi), %rax
shrq $0x3, %rax
imull $0x38e38e39, %eax, %eax # imm = 0x38E38E39
testl %eax, %eax
jle 0x7d494
movq %rsi, %r15
leaq 0x30(%r14), %r12
xorl %eax, %eax
cltq
movq 0x10(%rsp), %rcx
movl (%rcx), %r13d
testl %ebp, %ebp
je 0x7d440
movq 0x1d8(%r14), %rcx
jmp 0x7d447
movq 0x140(%r14), %rcx
movl (%rcx,%rax,4), %ebx
movq %r12, %rdi
leaq 0xc(%rsp), %rsi
callq 0x8000c
movslq (%rax), %rsi
movq (%r15), %rax
movq 0x8(%r15), %rdx
subq %rax, %rdx
sarq $0x2, %rdx
cmpq %rsi, %rdx
jbe 0x7d4bc
addl %ebx, %r13d
movl %r13d, (%rax,%rsi,4)
movl 0xc(%rsp), %eax
incl %eax
movl %eax, 0xc(%rsp)
movq 0x68(%r14), %rcx
subq 0x60(%r14), %rcx
shrq $0x3, %rcx
imull $0x38e38e39, %ecx, %ecx # imm = 0x38E38E39
cmpl %ecx, %eax
jl 0x7d429
xorl %eax, %eax
testl %ebp, %ebp
setne %al
shll $0x7, %eax
movl 0x138(%r14,%rax), %eax
movq 0x10(%rsp), %rcx
addl %eax, (%rcx)
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0xd21e2(%rip), %rdi     # 0x14f6a5
xorl %eax, %eax
callq 0x4630
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::merge_nothing()
void FEATURE::merge_nothing()
{
	int i;

	nfeatures=nsegments;
	for(i=0;i<nsegments;i++)
	{
		features.push_back(i);
		feature_max_r.push_back(segment_max_r[i]);
	};

	for(i=0;i<(int)nodes.size();i++)
	{
		feature_node_mapping.push_back(segments[i]);
	};

	cout << "After merging nothing, nfeatures = " << nfeatures << "\n";

}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
movl 0x138(%rdi), %eax
movl %eax, 0x1b8(%rdi)
movl $0x0, 0x4(%rsp)
testl %eax, %eax
jle 0x7d57f
leaq 0x1c0(%rbx), %r14
leaq 0x250(%rbx), %r15
xorl %eax, %eax
leaq 0x4(%rsp), %r12
movq 0x1c8(%rbx), %rsi
cmpq 0x1d0(%rbx), %rsi
je 0x7d525
movl %eax, (%rsi)
addq $0x4, %rsi
movq %rsi, 0x1c8(%rbx)
jmp 0x7d530
movq %r14, %rdi
movq %r12, %rdx
callq 0xf7b4
movslq 0x4(%rsp), %rdx
shlq $0x3, %rdx
addq 0x238(%rbx), %rdx
movq 0x258(%rbx), %rsi
cmpq 0x260(%rbx), %rsi
je 0x7d565
movsd (%rdx), %xmm0
movsd %xmm0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x258(%rbx)
jmp 0x7d56d
movq %r15, %rdi
callq 0x148c4
movl 0x4(%rsp), %eax
incl %eax
movl %eax, 0x4(%rsp)
cmpl 0x138(%rbx), %eax
jl 0x7d506
movl $0x0, 0x4(%rsp)
movq 0x68(%rbx), %rax
subq 0x60(%rbx), %rax
shrq $0x3, %rax
imull $0x38e38e39, %eax, %eax # imm = 0x38E38E39
testl %eax, %eax
jle 0x7d5fd
leaq 0x1d8(%rbx), %r14
xorl %eax, %eax
movslq %eax, %rdx
shlq $0x2, %rdx
addq 0x140(%rbx), %rdx
movq 0x1e0(%rbx), %rsi
cmpq 0x1e8(%rbx), %rsi
je 0x7d5d5
movl (%rdx), %eax
movl %eax, (%rsi)
addq $0x4, %rsi
movq %rsi, 0x1e0(%rbx)
jmp 0x7d5dd
movq %r14, %rdi
callq 0xf7b4
movl 0x4(%rsp), %eax
incl %eax
movl %eax, 0x4(%rsp)
movq 0x68(%rbx), %rcx
subq 0x60(%rbx), %rcx
shrq $0x3, %rcx
imull $0x38e38e39, %ecx, %ecx # imm = 0x38E38E39
cmpl %ecx, %eax
jl 0x7d5a6
movq 0x1199bc(%rip), %r14    # 0x196fc0
leaq 0xd9735(%rip), %rsi     # 0x156d40
movl $0x23, %edx
movq %r14, %rdi
callq 0x44e0
movl 0x1b8(%rbx), %esi
movq %r14, %rdi
callq 0x4750
leaq 0xd7f01(%rip), %rsi     # 0x15552e
movl $0x1, %edx
movq %rax, %rdi
callq 0x44e0
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::number_nodes_in_segment(int)
int FEATURE::number_nodes_in_segment(int segid)
{
	int i;
	int count=0;

	for(i=0;i<(int)nodes.size();i++)
	{
		if(segments[i]==segid) count++;
	};

	return count;
}

// Assembly:
movq 0x68(%rdi), %rax
subq 0x60(%rdi), %rax
sarq $0x3, %rax
movabsq $-0x71c71c71c71c71c7, %rcx # imm = 0x8E38E38E38E38E39
imulq %rax, %rcx
testl %ecx, %ecx
jle 0x7d68b
movq 0x140(%rdi), %rdx
andl $0x7fffffff, %ecx       # imm = 0x7FFFFFFF
xorl %edi, %edi
xorl %eax, %eax
xorl %r8d, %r8d
cmpl %esi, (%rdx,%rdi,4)
sete %r8b
addl %r8d, %eax
incq %rdi
cmpq %rdi, %rcx
jne 0x7d675
retq
xorl %eax, %eax
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::merge_findNotAssignedFragment()
int FEATURE::merge_findNotAssignedFragment()
{
	int i;
	int r=-1;
	for(i=segment_merge_list.size()-1;i>-1;i--)
	{
		if(segment_merge_list.at(i)==-1) r=i;
	};
	return r;
}

// Assembly:
pushq %rax
movq 0x208(%rdi), %rcx
movq 0x210(%rdi), %rdx
subq %rcx, %rdx
sarq $0x2, %rdx
movl %edx, %eax
decl %eax
js 0x7d9dd
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
movl %edx, %esi
decl %esi
cmpq %rsi, %rdx
jbe 0x7d9e4
cmpl $-0x1, (%rcx,%rsi,4)
cmovel %esi, %eax
testl %esi, %esi
jg 0x7d9c9
jmp 0x7d9e2
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
popq %rcx
retq
leaq 0xd1cba(%rip), %rdi     # 0x14f6a5
xorl %eax, %eax
callq 0x4630
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::merge_newSegmentGrow(int)
void FEATURE::merge_newSegmentGrow(int segid)
{
	int k;
	int conn_to;

	segment_merge_list[segid]=nfeatures;


	for(k=0;k<(int)segment_conn.size();k++)
	{
		if((segment_conn.at(k).from_seg==segid||segment_conn.at(k).to_seg==segid)&&segment_conn.at(k).merged==1)
		{
			if(segment_conn.at(k).from_seg==segid) 
				{conn_to=segment_conn.at(k).to_seg;}
			else {conn_to=segment_conn.at(k).from_seg;};

			if(segment_merge_list.at(conn_to)==-1) merge_newSegmentGrow(conn_to);


		}; // end if
	}; //end for loop


}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movl 0x1b8(%rdi), %eax
movslq %esi, %rcx
movq 0x208(%rdi), %rdx
movl %eax, (%rdx,%rcx,4)
movq 0x158(%rdi), %rax
movq 0x160(%rdi), %rdx
subq %rax, %rdx
sarq $0x4, %rdx
movabsq $-0x5555555555555555, %r12 # imm = 0xAAAAAAAAAAAAAAAB
imulq %r12, %rdx
testl %edx, %edx
jle 0x7dad5
movl %esi, %ebx
movq %rdi, %r14
movl $0x8, %r13d
movl $0xc, %ebp
xorl %r15d, %r15d
cmpq %r15, %rdx
jbe 0x7dae4
movl (%rax,%r13), %ecx
cmpl %ebx, %ecx
je 0x7da66
cmpl %ebx, 0x4(%rax,%r13)
jne 0x7daa5
cmpl $0x1, 0x18(%rax,%r13)
jne 0x7daa5
cmpl %ebx, %ecx
je 0x7da77
addq %r13, %rax
jmp 0x7da7a
addq %rbp, %rax
movslq (%rax), %rsi
movq 0x208(%r14), %rax
movq 0x210(%r14), %rdx
subq %rax, %rdx
sarq $0x2, %rdx
cmpq %rsi, %rdx
jbe 0x7daf0
cmpl $-0x1, (%rax,%rsi,4)
jne 0x7daa5
movq %r14, %rdi
callq 0x7d9f2
incq %r15
movq 0x158(%r14), %rax
movq 0x160(%r14), %rdx
subq %rax, %rdx
sarq $0x4, %rdx
imulq %r12, %rdx
movslq %edx, %rcx
addq $0x30, %r13
addq $0x30, %rbp
cmpq %rcx, %r15
jl 0x7da4e
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0xd1bba(%rip), %rdi     # 0x14f6a5
movq %r15, %rsi
jmp 0x7daf7
leaq 0xd1bae(%rip), %rdi     # 0x14f6a5
xorl %eax, %eax
callq 0x4630
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/v_network.cc
// Function: voronoi_network::~voronoi_network()
voronoi_network::~voronoi_network() {
	int l;

	// Remove Voronoi mapping array
	delete [] vmap;

	// Remove individual edge arrays
	for(l=0;l<edmem;l++) delete [] pered[l];
	for(l=0;l<edmem;l++) delete [] raded[l];
	for(l=0;l<edmem;l++) delete [] ed[l];

	// Remove back pointers
	delete [] regp;delete [] reg;

	// Remove memory for edges and related statistics
	delete [] numem;delete [] nec;delete [] nu;
	delete [] raded;delete [] pered;
	delete [] ne;delete [] ed;

	// Remove vertex structure arrays
	for(l=0;l<nxyz;l++) {
		delete [] idmem[l];
		delete [] pts[l];
	}
	delete [] ptsmem;delete [] ptsc;
	delete [] idmem;delete [] pts;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
movq 0xd0(%rdi), %rdi
testq %rdi, %rdi
je 0x80934
callq 0x4570
cmpl $0x0, 0xa4(%rbx)
jle 0x80964
xorl %r14d, %r14d
movq 0x98(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x80955
callq 0x4570
incq %r14
movslq 0xa4(%rbx), %rax
cmpq %rax, %r14
jl 0x80940
cmpl $0x0, 0xa4(%rbx)
jle 0x80994
xorl %r14d, %r14d
movq 0x90(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x80985
callq 0x4570
incq %r14
movslq 0xa4(%rbx), %rax
cmpq %rax, %r14
jl 0x80970
cmpl $0x0, 0xa4(%rbx)
jle 0x809c4
xorl %r14d, %r14d
movq 0x80(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x809b5
callq 0x4570
incq %r14
movslq 0xa4(%rbx), %rax
cmpq %rax, %r14
jl 0x809a0
movq 0xc8(%rbx), %rdi
testq %rdi, %rdi
je 0x809d5
callq 0x4570
movq 0xc0(%rbx), %rdi
testq %rdi, %rdi
je 0x809e6
callq 0x4570
movq 0xb8(%rbx), %rdi
testq %rdi, %rdi
je 0x809f7
callq 0x4570
movq 0xb0(%rbx), %rdi
testq %rdi, %rdi
je 0x80a08
callq 0x4570
movq 0xa8(%rbx), %rdi
testq %rdi, %rdi
je 0x80a19
callq 0x4570
movq 0x90(%rbx), %rdi
testq %rdi, %rdi
je 0x80a2a
callq 0x4570
movq 0x98(%rbx), %rdi
testq %rdi, %rdi
je 0x80a3b
callq 0x4570
movq 0x88(%rbx), %rdi
testq %rdi, %rdi
je 0x80a4c
callq 0x4570
movq 0x80(%rbx), %rdi
testq %rdi, %rdi
je 0x80a5d
callq 0x4570
cmpl $0x0, 0x3c(%rbx)
jle 0x80a96
xorl %r14d, %r14d
movq 0x68(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x80a78
callq 0x4570
movq 0x60(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x80a8a
callq 0x4570
incq %r14
movslq 0x3c(%rbx), %rax
cmpq %rax, %r14
jl 0x80a66
movq 0x78(%rbx), %rdi
testq %rdi, %rdi
je 0x80aa4
callq 0x4570
movq 0x70(%rbx), %rdi
testq %rdi, %rdi
je 0x80ab2
callq 0x4570
movq 0x68(%rbx), %rdi
testq %rdi, %rdi
je 0x80ac0
callq 0x4570
movq 0x60(%rbx), %rdi
addq $0x8, %rsp
testq %rdi, %rdi
je 0x80ad5
popq %rbx
popq %r14
jmp 0x4570
popq %rbx
popq %r14
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/v_network.cc
// Function: voronoi_network::add_network_memory(int)
void voronoi_network::add_network_memory(int l) {
	ptsmem[l]<<=1;
	
	// Check to see that an absolute maximum in memory allocation
	// has not been reached, to prevent runaway allocation
	if(ptsmem[l]>max_network_vertex_memory)
		voro_fatal_error("Container vertex maximum memory allocation exceeded",VOROPP_MEMORY_ERROR);
	
	// Allocate new arrays
	double *npts(new double[4*ptsmem[l]]);
	int *nidmem(new int[ptsmem[l]]);

	// Copy the contents of the old arrays into the new ones
	for(int i=0;i<4*ptsc[l];i++) npts[i]=pts[l][i];
	for(int i=0;i<ptsc[l];i++) nidmem[i]=idmem[l][i];

	// Delete old arrays and update pointers to the new ones
	delete [] pts[l];delete [] idmem[l];
	pts[l]=npts;idmem[l]=nidmem;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
movq %rdi, %rbx
movq 0x78(%rdi), %rax
movslq %esi, %r12
movl (%rax,%r12,4), %ecx
addl %ecx, %ecx
movl %ecx, (%rax,%r12,4)
cmpl $0x10001, %ecx          # imm = 0x10001
jl 0x80b10
leaq 0xd6395(%rip), %rdi     # 0x156e9b
movl $0x2, %esi
callq 0x12ecbd
movq 0x78(%rbx), %r15
movslq (%r15,%r12,4), %rax
movq %rax, %rdi
shlq $0x5, %rdi
testq %rax, %rax
movq $-0x1, %r13
cmovsq %r13, %rdi
callq 0x4030
movq %rax, %r14
movslq (%r15,%r12,4), %rax
leaq (,%rax,4), %rdi
testq %rax, %rax
cmovsq %r13, %rdi
callq 0x4030
movq %rax, %r15
movq 0x70(%rbx), %rax
movl (%rax,%r12,4), %esi
testl %esi, %esi
jle 0x80b87
shll $0x2, %esi
movq 0x60(%rbx), %rcx
movq (%rcx,%r12,8), %rcx
cmpl $0x2, %esi
movl $0x1, %edx
cmovgel %esi, %edx
xorl %esi, %esi
movsd (%rcx,%rsi,8), %xmm0
movsd %xmm0, (%r14,%rsi,8)
incq %rsi
cmpq %rsi, %rdx
jne 0x80b74
cmpl $0x0, (%rax,%r12,4)
jle 0x80bab
movq 0x68(%rbx), %rcx
movq (%rcx,%r12,8), %rcx
xorl %edx, %edx
movl (%rcx,%rdx,4), %esi
movl %esi, (%r15,%rdx,4)
incq %rdx
movslq (%rax,%r12,4), %rsi
cmpq %rsi, %rdx
jl 0x80b98
movq 0x60(%rbx), %rax
movq (%rax,%r12,8), %rdi
testq %rdi, %rdi
je 0x80bbd
callq 0x4570
movq 0x68(%rbx), %rax
movq (%rax,%r12,8), %rdi
testq %rdi, %rdi
je 0x80bcf
callq 0x4570
movq 0x60(%rbx), %rax
movq %r14, (%rax,%r12,8)
movq 0x68(%rbx), %rax
movq %r15, (%rax,%r12,8)
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/v_network.cc
// Function: voronoi_network::add_particular_vertex_memory(int)
void voronoi_network::add_particular_vertex_memory(int l) {
	numem[l]<<=1;

	// Check that the vertex allocation does not exceed a maximum safe
	// limit
	if(numem[l]>max_vertex_order)
		voro_fatal_error("Particular vertex maximum memory allocation exceeded",VOROPP_MEMORY_ERROR);
	
	// Allocate new arrays
	int *ned(new int[2*numem[l]]);
	int *nne(ned+numem[l]);
	block *nraded(new block[numem[l]]);
	unsigned int *npered(new unsigned int[numem[l]]);

	// Copy the contents of the old arrays into the new ones 
	for(int i=0;i<nu[l];i++) {
		ned[i]=ed[l][i];
		nraded[i]=raded[l][i];
		npered[i]=pered[l][i];
	}
	for(int i=0;i<nec[l];i++) nne[i]=ne[l][i];
	
	// Delete old arrays and update pointers to the new ones
	delete [] ed[l];ed[l]=ned;ne[l]=nne;
	delete [] raded[l];raded[l]=nraded;
	delete [] pered[l];pered[l]=npered;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
movq 0xb8(%rdi), %rax
movslq %esi, %r13
movl (%rax,%r13,4), %ecx
addl %ecx, %ecx
movl %ecx, (%rax,%r13,4)
cmpl $0x801, %ecx            # imm = 0x801
jl 0x80fe1
leaq 0xd5ef8(%rip), %rdi     # 0x156ecf
movl $0x2, %esi
callq 0x12ecbd
movq 0xb8(%rbx), %r14
movslq (%r14,%r13,4), %rax
leaq (,%rax,8), %rdi
testq %rax, %rax
movq $-0x1, %rax
cmovsq %rax, %rdi
movq $-0x1, %r15
callq 0x4030
movq %rax, %r12
movslq (%r14,%r13,4), %rbp
movq %rbp, %rdi
shlq $0x4, %rdi
testq %rbp, %rbp
cmovsq %r15, %rdi
callq 0x4030
movq %rax, %r15
movslq (%r14,%r13,4), %rax
leaq (,%rax,4), %rdi
testq %rax, %rax
movq $-0x1, %rax
cmovsq %rax, %rdi
callq 0x4030
movq %rax, %r14
movq 0xa8(%rbx), %rax
cmpl $0x0, (%rax,%r13,4)
jle 0x810ab
movq 0x80(%rbx), %rcx
movq 0x90(%rbx), %rdx
movq 0x98(%rbx), %rsi
xorl %edi, %edi
xorl %r8d, %r8d
movq (%rcx,%r13,8), %r9
movl (%r9,%rdi), %r9d
movl %r9d, (%r12,%rdi)
movq (%rdx,%r13,8), %r9
movups (%r9,%rdi,4), %xmm0
movups %xmm0, (%r15,%rdi,4)
movq (%rsi,%r13,8), %r9
movl (%r9,%rdi), %r9d
movl %r9d, (%r14,%rdi)
incq %r8
movslq (%rax,%r13,4), %r9
addq $0x4, %rdi
cmpq %r9, %r8
jl 0x81075
leaq (%r12,%rbp,4), %rbp
movq 0xb0(%rbx), %rax
cmpl $0x0, (%rax,%r13,4)
jle 0x810dd
movq 0x88(%rbx), %rcx
movq (%rcx,%r13,8), %rcx
xorl %edx, %edx
movl (%rcx,%rdx,4), %esi
movl %esi, (%rbp,%rdx,4)
incq %rdx
movslq (%rax,%r13,4), %rsi
cmpq %rsi, %rdx
jl 0x810ca
movq 0x80(%rbx), %rax
movq (%rax,%r13,8), %rdi
testq %rdi, %rdi
je 0x810f2
callq 0x4570
movq 0x80(%rbx), %rax
movq %r12, (%rax,%r13,8)
movq 0x88(%rbx), %rax
movq %rbp, (%rax,%r13,8)
movq 0x90(%rbx), %rax
movq (%rax,%r13,8), %rdi
testq %rdi, %rdi
je 0x8111d
callq 0x4570
movq 0x90(%rbx), %rax
movq %r15, (%rax,%r13,8)
movq 0x98(%rbx), %rax
movq (%rax,%r13,8), %rdi
testq %rdi, %rdi
je 0x8113d
callq 0x4570
movq 0x98(%rbx), %rax
movq %r14, (%rax,%r13,8)
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/v_network.cc
// Function: voronoi_network::clear_network()
void voronoi_network::clear_network() {
	int l;
	edc=0;
	for(l=0;l<nxyz;l++) ptsc[l]=0;
	for(l=0;l<edmem;l++) nu[l]=0;
}

// Assembly:
movl $0x0, 0xa0(%rdi)
cmpl $0x0, 0x3c(%rdi)
jle 0x811d3
movq 0x70(%rdi), %rax
xorl %ecx, %ecx
movl $0x0, (%rax,%rcx,4)
incq %rcx
movslq 0x3c(%rdi), %rdx
cmpq %rdx, %rcx
jl 0x811c0
cmpl $0x0, 0xa4(%rdi)
jle 0x811fb
movq 0xa8(%rdi), %rax
xorl %ecx, %ecx
movl $0x0, (%rax,%rcx,4)
incq %rcx
movslq 0xa4(%rdi), %rdx
cmpq %rdx, %rcx
jl 0x811e5
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/v_network.cc
// Function: voronoi_network::draw_network(_IO_FILE*)
void voronoi_network::draw_network(FILE *fp) {
	int l,q,ai,aj,ak;
	double x,y,z,*ptsp;
	for(l=0;l<edc;l++) {
		ptsp=pts[reg[l]]+4*regp[l];
		x=*(ptsp++);y=*(ptsp++);z=*ptsp;
		for(q=0;q<nu[l];q++) {
			unpack_periodicity(pered[l][q],ai,aj,ak);
			if(ed[l][q]<l&&ai==0&&aj==0&&ak==0) continue;
			ptsp=pts[reg[ed[l][q]]]+4*regp[ed[l][q]];
			fprintf(fp,"%g %g %g\n%g %g %g\n\n\n",x,y,z,
				*ptsp+bx*ai+bxy*aj+bxz*ak,
				ptsp[1]+by*aj+byz*ak,ptsp[2]+bz*ak);
		}
	}
}

// Assembly:
cmpl $0x0, 0xa0(%rdi)
jle 0x813a8
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rsi, %rbx
movq %rdi, %r14
xorl %r12d, %r12d
movl $0xff, %ebp
leaq 0xd5cd8(%rip), %r15     # 0x156f04
movq 0xa8(%r14), %rax
cmpl $0x0, (%rax,%r12,4)
jle 0x81387
movq 0xc0(%r14), %rax
movslq (%rax,%r12,4), %rax
movq 0x60(%r14), %rcx
movq (%rcx,%rax,8), %rax
movq 0xc8(%r14), %rcx
movslq (%rcx,%r12,4), %rcx
shlq $0x5, %rcx
movsd (%rax,%rcx), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rax,%rcx), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rax,%rcx), %xmm0
movsd %xmm0, (%rsp)
xorl %r13d, %r13d
movq 0x80(%r14), %rsi
movq 0x98(%r14), %rax
movq (%rax,%r12,8), %rax
movl (%rax,%r13,4), %eax
movl %eax, %edx
shrl $0x10, %edx
addl $-0x7f, %edx
movl %eax, %ecx
shrl $0x8, %ecx
andl %ebp, %ecx
addl $-0x7f, %ecx
andl %ebp, %eax
addl $-0x7f, %eax
movq (%rsi,%r12,8), %rsi
movslq (%rsi,%r13,4), %rdi
cmpq %rdi, %r12
jle 0x812cf
testl %edx, %edx
jne 0x812cf
testl %ecx, %ecx
jne 0x812cf
testl %eax, %eax
je 0x81370
movq 0x60(%r14), %rsi
movq 0xc0(%r14), %r8
movslq (%r8,%rdi,4), %r8
movq (%rsi,%r8,8), %rsi
movq 0xc8(%r14), %r8
movslq (%r8,%rdi,4), %rdi
shlq $0x5, %rdi
xorps %xmm0, %xmm0
cvtsi2sd %edx, %xmm0
mulsd (%r14), %xmm0
addsd (%rsi,%rdi), %xmm0
xorps %xmm1, %xmm1
cvtsi2sd %ecx, %xmm1
movsd 0x8(%r14), %xmm2
mulsd %xmm1, %xmm2
addsd %xmm0, %xmm2
xorps %xmm5, %xmm5
cvtsi2sd %eax, %xmm5
movsd 0x18(%r14), %xmm3
mulsd %xmm5, %xmm3
mulsd 0x10(%r14), %xmm1
addsd 0x8(%rsi,%rdi), %xmm1
addsd %xmm2, %xmm3
movsd 0x20(%r14), %xmm4
mulsd %xmm5, %xmm4
addsd %xmm1, %xmm4
mulsd 0x28(%r14), %xmm5
addsd 0x10(%rsi,%rdi), %xmm5
movq %rbx, %rdi
movq %r15, %rsi
movsd 0x10(%rsp), %xmm0
movsd 0x8(%rsp), %xmm1
movsd (%rsp), %xmm2
movb $0x6, %al
callq 0x45f0
incq %r13
movq 0xa8(%r14), %rax
movslq (%rax,%r12,4), %rax
cmpq %rax, %r13
jl 0x81285
incq %r12
movslq 0xa0(%r14), %rax
cmpq %rax, %r12
jl 0x8122c
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: trans_to_origuc(double)
double trans_to_origuc(double x){
  double newX =  x - (x < 0.0 ?(-(int)(0.5-x)):((int)(0.5+x)));
  if(newX < 0)
    return newX+1;
  else
    return newX;
}

// Assembly:
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jbe 0x84e89
movsd 0xca1bb(%rip), %xmm1    # 0x14f038
subsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
negl %eax
jmp 0x84e99
movsd 0xca1a7(%rip), %xmm1    # 0x14f038
addsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
subsd %xmm1, %xmm0
movsd 0xca1ac(%rip), %xmm2    # 0x14f058
addsd %xmm0, %xmm2
xorpd %xmm3, %xmm3
movapd %xmm0, %xmm1
cmpltsd %xmm3, %xmm1
andpd %xmm1, %xmm2
andnpd %xmm0, %xmm1
orpd %xmm2, %xmm1
movapd %xmm1, %xmm0
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: invertMatrix(double (*) [3], double (*) [3])
double calcDeterminant(double matrix[3][3]){
    return
    matrix[0][0]*(matrix[2][2]*matrix[1][1] - matrix[2][1]*matrix[1][2])
    - matrix[1][0]*(matrix[2][2]*matrix[0][1] - matrix[2][1]*matrix[0][2])
    + matrix[2][0]*(matrix[1][2]*matrix[0][1] - matrix[1][1]*matrix[0][2]);
}

// Assembly:
movsd 0x40(%rdi), %xmm0
movsd 0x20(%rdi), %xmm3
movsd 0x38(%rdi), %xmm4
movsd 0x28(%rdi), %xmm2
movapd %xmm2, %xmm5
mulsd %xmm4, %xmm5
movapd %xmm0, %xmm1
mulsd %xmm3, %xmm1
subsd %xmm5, %xmm1
movsd 0x8(%rdi), %xmm5
movsd 0x10(%rdi), %xmm6
mulsd %xmm6, %xmm4
mulsd %xmm5, %xmm0
subsd %xmm4, %xmm0
mulsd 0x18(%rdi), %xmm0
movsd (%rdi), %xmm4
mulsd %xmm1, %xmm4
subsd %xmm0, %xmm4
mulsd %xmm3, %xmm6
mulsd %xmm5, %xmm2
subsd %xmm6, %xmm2
mulsd 0x30(%rdi), %xmm2
addsd %xmm4, %xmm2
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm2
jne 0x8511c
jnp 0x85241
movsd 0xc9f34(%rip), %xmm0    # 0x14f058
divsd %xmm2, %xmm0
mulsd %xmm0, %xmm1
movsd %xmm1, (%rsi)
movapd 0xce478(%rip), %xmm1    # 0x1535b0
xorpd %xmm0, %xmm1
movsd 0x10(%rdi), %xmm2
mulsd 0x38(%rdi), %xmm2
movsd 0x40(%rdi), %xmm3
mulsd 0x8(%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x8(%rsi)
movsd 0x10(%rdi), %xmm2
movsd 0x28(%rdi), %xmm3
mulsd 0x20(%rdi), %xmm2
mulsd 0x8(%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x10(%rsi)
movsd 0x28(%rdi), %xmm2
movsd 0x40(%rdi), %xmm3
mulsd 0x30(%rdi), %xmm2
mulsd 0x18(%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x18(%rsi)
movsd 0x10(%rdi), %xmm2
mulsd 0x30(%rdi), %xmm2
movsd 0x40(%rdi), %xmm3
mulsd (%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x20(%rsi)
movsd 0x10(%rdi), %xmm2
movsd 0x28(%rdi), %xmm3
mulsd 0x18(%rdi), %xmm2
mulsd (%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x28(%rsi)
movsd 0x20(%rdi), %xmm2
movsd 0x38(%rdi), %xmm3
mulsd 0x30(%rdi), %xmm2
mulsd 0x18(%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x30(%rsi)
movsd 0x8(%rdi), %xmm2
mulsd 0x30(%rdi), %xmm2
movsd 0x38(%rdi), %xmm3
mulsd (%rdi), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x38(%rsi)
movsd 0x8(%rdi), %xmm1
movsd 0x20(%rdi), %xmm2
mulsd 0x18(%rdi), %xmm1
mulsd (%rdi), %xmm2
subsd %xmm1, %xmm2
mulsd %xmm0, %xmm2
movsd %xmm2, 0x40(%rsi)
retq
pushq %rbx
movq 0x111da7(%rip), %rdi    # 0x196ff0
leaq 0xd1d28(%rip), %rsi     # 0x156f78
callq 0x4440
leaq 0xd02d2(%rip), %rbx     # 0x15552e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
leaq 0xca8d5(%rip), %rsi     # 0x14fb43
movq %rax, %rdi
callq 0x4440
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
movl $0x1, %edi
callq 0x710c
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: Point::operator[](int)
double& Point::operator[](const int index) {
     if(index < 0 || index > 2){
      cerr << "Error: Invalid index to [] operator for Point instance" << "\n"
	   << "Exiting..." << "\n";
      exit(1);
    }
    return vals[index];
}

// Assembly:
cmpl $0x3, %esi
jae 0x85484
movl %esi, %eax
leaq (%rdi,%rax,8), %rax
retq
pushq %rbx
movq 0x111b64(%rip), %rdi    # 0x196ff0
leaq 0xd1b24(%rip), %rsi     # 0x156fb7
callq 0x4440
leaq 0xd008f(%rip), %rbx     # 0x15552e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
leaq 0xd1101(%rip), %rsi     # 0x1565b2
movq %rax, %rdi
callq 0x4440
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
movl $0x1, %edi
callq 0x710c
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: pointIsLess(Point, Point)
bool pointIsLess(Point p1, Point p2){
  double tol = 0.0000001;
  if (abs(p1[0] - p2[0]) > tol)
    return p1[0] < p2[0];
  else if (abs(p1[1] - p2[1]) > tol)
    return p1[1] < p2[1];
  else if (abs(p1[2] - p2[2]) > tol)
    return p1[2] < p2[2];
  else
    return false;
}

// Assembly:
movsd 0x8(%rsp), %xmm0
movsd 0x20(%rsp), %xmm1
movapd %xmm0, %xmm2
subsd %xmm1, %xmm2
andpd 0xc9809(%rip), %xmm2    # 0x14f090
ucomisd 0xce939(%rip), %xmm2    # 0x1541c8
ja 0x858e3
leaq 0x20(%rsp), %rax
leaq 0x8(%rsp), %rcx
movsd 0x8(%rcx), %xmm0
movsd 0x8(%rax), %xmm1
movapd %xmm0, %xmm2
subsd %xmm1, %xmm2
andpd 0xc97db(%rip), %xmm2    # 0x14f090
ucomisd 0xce90b(%rip), %xmm2    # 0x1541c8
ja 0x858e3
movsd 0x10(%rcx), %xmm0
movsd 0x10(%rax), %xmm1
movapd %xmm0, %xmm2
subsd %xmm1, %xmm2
andpd 0xc97b7(%rip), %xmm2    # 0x14f090
ucomisd 0xce8e7(%rip), %xmm2    # 0x1541c8
jbe 0x858eb
ucomisd %xmm0, %xmm1
seta %al
retq
xorl %eax, %eax
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: TRIPLET::operator[](int)
int& TRIPLET::operator[](int index){
  if(index < 0 || index > 2){
    cerr << "Error: Invalid index to [] operator for TRIPLET instance" << "\n"
	 << "Exiting..." << "\n";
    exit(1);
  }
  return vals[index];
}

// Assembly:
cmpl $0x3, %esi
jae 0x8590c
movl %esi, %eax
leaq (%rdi,%rax,4), %rax
retq
pushq %rbx
movq 0x1116dc(%rip), %rdi    # 0x196ff0
leaq 0xd16d3(%rip), %rsi     # 0x156fee
callq 0x4440
leaq 0xcfc07(%rip), %rbx     # 0x15552e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
leaq 0xd0c79(%rip), %rsi     # 0x1565b2
movq %rax, %rdi
callq 0x4440
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
movl $0x1, %edi
callq 0x710c
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: XYZ::angle_between(XYZ const&) const
double XYZ::angle_between(const XYZ& other)  const {
  double cos_angle = dot_product(other)/(magnitude()*other.magnitude());
  // These steps are necessary for rounding issues when cos_angle should be 
  // exactly 1 or -1, but falls very slightly on the wrong side
  if(cos_angle>1) {
    cos_angle = 1;
  } else if(cos_angle<-1) {
    cos_angle = -1; 
  }
  double output = acos(cos_angle);
  if(isnan(output)) {
    return 0;
  } else return output; //this step is necessary in the case of identity, which produces NAN
}

// Assembly:
pushq %rbx
subq $0x10, %rsp
movq %rsi, %rbx
movsd (%rdi), %xmm1
movsd 0x8(%rdi), %xmm2
movsd 0x8(%rsi), %xmm0
mulsd %xmm2, %xmm0
movsd (%rsi), %xmm4
mulsd %xmm1, %xmm4
addsd %xmm0, %xmm4
movsd 0x10(%rdi), %xmm0
movsd 0x10(%rsi), %xmm3
mulsd %xmm0, %xmm3
mulsd %xmm2, %xmm2
mulsd %xmm1, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
xorpd %xmm5, %xmm5
ucomisd %xmm5, %xmm0
jb 0x85c8b
xorps %xmm2, %xmm2
sqrtsd %xmm0, %xmm2
jmp 0x85cae
movsd %xmm3, 0x8(%rsp)
movsd %xmm4, (%rsp)
callq 0x4350
xorpd %xmm5, %xmm5
movsd (%rsp), %xmm4
movsd 0x8(%rsp), %xmm3
movapd %xmm0, %xmm2
addsd %xmm4, %xmm3
movsd (%rbx), %xmm1
movsd 0x8(%rbx), %xmm0
mulsd %xmm0, %xmm0
mulsd %xmm1, %xmm1
addsd %xmm0, %xmm1
movsd 0x10(%rbx), %xmm0
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
ucomisd %xmm5, %xmm0
jb 0x85ce0
sqrtsd %xmm0, %xmm0
jmp 0x85cfb
movsd %xmm3, 0x8(%rsp)
movsd %xmm2, (%rsp)
callq 0x4350
movsd (%rsp), %xmm2
movsd 0x8(%rsp), %xmm3
mulsd %xmm0, %xmm2
divsd %xmm2, %xmm3
movsd 0xc934d(%rip), %xmm0    # 0x14f058
ucomisd %xmm0, %xmm3
ja 0x85d27
movsd 0xcd5b7(%rip), %xmm1    # 0x1532d0
ucomisd %xmm3, %xmm1
movapd %xmm3, %xmm0
jbe 0x85d27
movapd %xmm1, %xmm0
callq 0x4400
movapd %xmm0, %xmm1
cmpordsd %xmm0, %xmm1
andpd %xmm1, %xmm0
addq $0x10, %rsp
popq %rbx
retq
nop
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: XYZ::operator[](int)
double& XYZ::operator[](const int index){
     if (index==0)
       return x;
     if (index==1)
       return y;
     if (index==2)
       return z;
     cerr << "Error: Invalid index to [] operator for XYZ instance" << "\n"
	  << "Exiting..." << "\n";
     exit(1);
}

// Assembly:
testl %esi, %esi
je 0x85d58
cmpl $0x2, %esi
je 0x85d54
cmpl $0x1, %esi
jne 0x85d5c
addq $0x8, %rdi
jmp 0x85d58
addq $0x10, %rdi
movq %rdi, %rax
retq
pushq %rbx
movq 0x11128c(%rip), %rdi    # 0x196ff0
leaq 0xd12bc(%rip), %rsi     # 0x157027
callq 0x4440
leaq 0xcf7b7(%rip), %rbx     # 0x15552e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
leaq 0xd0829(%rip), %rsi     # 0x1565b2
movq %rax, %rdi
callq 0x4440
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4440
movl $0x1, %edi
callq 0x710c
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: XYZ::euclid_dist(XYZ const&) const
double XYZ::euclid_dist(const XYZ& other) const {  //Distance between two XYZ's
    double dx = x-other.x;  
    double dy = y-other.y; 
    double dz = z-other.z;
    return sqrt(dx*dx + dy*dy + dz*dz);
}

// Assembly:
movsd (%rdi), %xmm1
movsd 0x8(%rdi), %xmm2
subsd (%rsi), %xmm1
subsd 0x8(%rsi), %xmm2
movsd 0x10(%rdi), %xmm0
subsd 0x10(%rsi), %xmm0
mulsd %xmm2, %xmm2
mulsd %xmm1, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x4350
sqrtsd %xmm0, %xmm0
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: project_onto_line(XYZ, XYZ, XYZ)
XYZ project_onto_line(XYZ initial, XYZ line_start, XYZ line_end) {
  XYZ projection;
  XYZ initial_to_line_start_vector = get_vector_from_to(initial, line_start);
  XYZ line_vector = get_vector_from_to(line_start, line_end);
  double line_length = line_vector.magnitude();
  double movement_along_line = (-1*initial_to_line_start_vector.dot_product(line_vector))/(line_length*line_length);
  projection.x = line_start.x + (movement_along_line*line_vector.x);
  projection.y = line_start.y + (movement_along_line*line_vector.y);
  projection.z = line_start.z + (movement_along_line*line_vector.z);
  return projection;
}

// Assembly:
pushq %rbx
subq $0x50, %rsp
movq %rdi, %rbx
xorps %xmm0, %xmm0
movups %xmm0, (%rdi)
movq $0x0, 0x10(%rdi)
movsd 0x78(%rsp), %xmm5
movsd 0x80(%rsp), %xmm3
movapd %xmm5, %xmm9
subsd 0x60(%rsp), %xmm9
movsd 0x88(%rsp), %xmm2
movapd %xmm3, %xmm10
subsd 0x68(%rsp), %xmm10
movapd %xmm2, %xmm4
subsd 0x70(%rsp), %xmm4
movsd 0x90(%rsp), %xmm8
movsd 0x98(%rsp), %xmm7
movsd 0xa0(%rsp), %xmm6
subsd %xmm5, %xmm8
subsd %xmm3, %xmm7
subsd %xmm2, %xmm6
movapd %xmm7, %xmm0
mulsd %xmm7, %xmm0
movapd %xmm8, %xmm1
mulsd %xmm8, %xmm1
addsd %xmm0, %xmm1
movapd %xmm6, %xmm0
mulsd %xmm6, %xmm0
addsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x85f96
sqrtsd %xmm0, %xmm0
jmp 0x8600b
movsd %xmm2, 0x38(%rsp)
movsd %xmm3, 0x30(%rsp)
movapd %xmm4, 0x40(%rsp)
movsd %xmm5, 0x28(%rsp)
movsd %xmm6, 0x20(%rsp)
movsd %xmm7, 0x18(%rsp)
movsd %xmm8, 0x10(%rsp)
movsd %xmm9, 0x8(%rsp)
movsd %xmm10, (%rsp)
callq 0x4350
movsd (%rsp), %xmm10
movsd 0x8(%rsp), %xmm9
movsd 0x10(%rsp), %xmm8
movsd 0x18(%rsp), %xmm7
movsd 0x20(%rsp), %xmm6
movsd 0x28(%rsp), %xmm5
movapd 0x40(%rsp), %xmm4
movsd 0x30(%rsp), %xmm3
movsd 0x38(%rsp), %xmm2
mulsd %xmm7, %xmm10
mulsd %xmm8, %xmm9
addsd %xmm10, %xmm9
mulsd %xmm6, %xmm4
addsd %xmm9, %xmm4
xorpd 0xcd585(%rip), %xmm4    # 0x1535b0
mulsd %xmm0, %xmm0
divsd %xmm0, %xmm4
mulsd %xmm4, %xmm8
addsd %xmm5, %xmm8
movsd %xmm8, (%rbx)
mulsd %xmm4, %xmm7
addsd %xmm3, %xmm7
movsd %xmm7, 0x8(%rbx)
mulsd %xmm6, %xmm4
addsd %xmm2, %xmm4
movsd %xmm4, 0x10(%rbx)
movq %rbx, %rax
addq $0x50, %rsp
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: abc_to_xyz(double, double, double, XYZ const&, XYZ const&, XYZ const&)
Point abc_to_xyz(double a, double b, double c, const XYZ& v_a, const XYZ& v_b, const XYZ& v_c){

  double x = v_a.x+b*v_b.x+c*v_c.x;
  double y = b*v_b.y+c*v_c.y;
  double z = c*v_c.z;
  return Point (x, y, z);
}

// Assembly:
movq %rdi, %rax
movsd (%rdx), %xmm0
mulsd %xmm1, %xmm0
addsd (%rsi), %xmm0
movsd (%rcx), %xmm3
mulsd %xmm2, %xmm3
addsd %xmm0, %xmm3
movsd 0x8(%rcx), %xmm0
mulsd 0x8(%rdx), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm0, %xmm1
mulsd 0x10(%rcx), %xmm2
movsd %xmm3, (%rdi)
movsd %xmm1, 0x8(%rdi)
movsd %xmm2, 0x10(%rdi)
retq
// File: /workspace/github_cpp_makefile_star2_O1/Luthaf[P]zeopp/geometry.cc
// Function: calcSphereDistance(Point, Point, double)
double calcSphereDistance(Point p1, Point p2, double circRad){
  pair<double,double> coords1 = findLongAndLat(p1);
  pair<double,double> coords2 = findLongAndLat(p2);
  double phi1 = coords1.first; double phi2 = coords2.first;
  double lambda1 = coords1.second; double lambda2 = coords2.second;
  double deltaLambda = lambda1-lambda2;

  double deltaSigma = atan( (sqrt (pow(cos(phi1)*sin(deltaLambda),2.0) +  pow(cos(phi2)*sin(phi1)-sin(phi2)*cos(phi1)*cos(deltaLambda),2))) 
			    / (sin(phi1)*sin(phi2) + cos(phi1)*cos(phi2)*cos(deltaLambda)));
  return circRad*deltaSigma;
}

// Assembly:
pushq %rbx
subq $0x50, %rsp
movsd %xmm0, 0x48(%rsp)
leaq 0x78(%rsp), %rbx
movsd 0x60(%rsp), %xmm0
movsd 0x68(%rsp), %xmm1
movsd 0x70(%rsp), %xmm2
movsd %xmm1, 0x18(%rsp)
mulsd %xmm1, %xmm1
movsd %xmm0, 0x8(%rsp)
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x871bf
sqrtsd %xmm0, %xmm0
jmp 0x871d0
movsd %xmm2, 0x28(%rsp)
callq 0x4350
movsd 0x28(%rsp), %xmm2
divsd %xmm0, %xmm2
movapd %xmm2, %xmm0
callq 0x4380
movsd %xmm0, 0x28(%rsp)
movsd 0x18(%rsp), %xmm0
movsd 0x8(%rsp), %xmm1
callq 0x40a0
movsd %xmm0, 0x8(%rsp)
movsd (%rbx), %xmm0
movsd 0x8(%rbx), %xmm1
movsd 0x10(%rbx), %xmm2
movsd %xmm1, 0x10(%rsp)
mulsd %xmm1, %xmm1
movsd %xmm0, 0x18(%rsp)
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
ucomisd 0xcc098(%rip), %xmm0    # 0x1532c0
jb 0x87230
sqrtsd %xmm0, %xmm0
jmp 0x87241
movsd %xmm2, 0x20(%rsp)
callq 0x4350
movsd 0x20(%rsp), %xmm2
divsd %xmm0, %xmm2
movapd %xmm2, %xmm0
callq 0x4380
movsd %xmm0, 0x20(%rsp)
movsd 0x10(%rsp), %xmm0
movsd 0x18(%rsp), %xmm1
callq 0x40a0
movsd %xmm0, 0x18(%rsp)
movsd 0x28(%rsp), %xmm0
subsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x28(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x4800
movsd %xmm0, 0x10(%rsp)
movsd 0x28(%rsp), %xmm0
callq 0x4100
mulsd 0x10(%rsp), %xmm0
mulsd %xmm0, %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x4800
movsd %xmm0, 0x20(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x4100
movsd %xmm0, 0x30(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x4100
movsd %xmm0, 0x40(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x4800
movsd %xmm0, 0x38(%rsp)
movsd 0x28(%rsp), %xmm0
callq 0x4800
movsd 0x38(%rsp), %xmm1
mulsd 0x40(%rsp), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, %xmm0
movsd 0x30(%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
subsd %xmm0, %xmm1
mulsd %xmm1, %xmm1
addsd 0x10(%rsp), %xmm1
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
movapd %xmm1, %xmm0
jb 0x8733a
sqrtsd %xmm0, %xmm0
jmp 0x8733f
callq 0x4350
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x4100
movsd %xmm0, 0x20(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x4100
movsd %xmm0, 0x30(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x4800
movsd %xmm0, 0x8(%rsp)
movsd 0x18(%rsp), %xmm0
callq 0x4800
mulsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x28(%rsp), %xmm0
callq 0x4800
mulsd 0x8(%rsp), %xmm0
movsd 0x30(%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
divsd %xmm1, %xmm0
callq 0x4380
mulsd 0x48(%rsp), %xmm0
addq $0x50, %rsp
popq %rbx
retq