LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O2/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 0xe1180
xorps %xmm0, %xmm0
movq %rsp, %rsi
movaps %xmm0, (%rsi)
andq $0x0, 0x10(%rsi)
leaq 0x20(%rsp), %rdx
movaps %xmm0, (%rdx)
andq $0x0, 0x10(%rdx)
movq %r12, %rdi
callq 0xe1892
movq (%r15), %rax
cmpq %rax, 0x8(%r15)
je 0x59721
movq %rax, 0x8(%r15)
movq (%r14), %rax
cmpq %rax, 0x8(%r14)
je 0x5972e
movq %rax, 0x8(%r14)
movq (%rbx), %rax
cmpq %rax, 0x8(%rbx)
je 0x5973b
movq %rax, 0x8(%rbx)
xorl %ebp, %ebp
pushq $0x2
popq %r12
movl %ebp, %eax
movq 0x28(%rsp), %rcx
subq 0x20(%rsp), %rcx
sarq $0x3, %rcx
cmpq %rax, %rcx
jbe 0x5979d
leal -0x2(%r12), %eax
movl %eax, %esi
shlq $0x2, %rsi
addq (%rsp), %rsi
movq %r15, %rdi
callq 0x9434
leal -0x1(%r12), %esi
shlq $0x2, %rsi
addq (%rsp), %rsi
movq %r14, %rdi
callq 0x9434
movl %r12d, %esi
shlq $0x2, %rsi
addq (%rsp), %rsi
movq %rbx, %rdi
callq 0x9434
incl %ebp
addl $0x3, %r12d
jmp 0x59741
leaq 0x20(%rsp), %rdi
callq 0xb25c
movq %rsp, %rdi
callq 0xb290
leaq 0x68(%rsp), %rdi
callq 0xd6d3c
addq $0x160, %rsp            # imm = 0x160
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
jmp 0x597cb
movq %rax, %rbx
leaq 0x20(%rsp), %rdi
callq 0xb25c
movq %rsp, %rdi
callq 0xb290
leaq 0x68(%rsp), %rdi
callq 0xd6d3c
movq %rbx, %rdi
callq 0x4780
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/mindist.cc
// Function: MIN_PER_DISTANCE::minimum_periodic_distance(double, double, double, double, double, double, double&, double&, double&)
double MIN_PER_DISTANCE::minimum_periodic_distance(double a1, double b1, double c1, double a2, double b2, double c2, double &minDa, double &minDb, double &minDc){
  double va_a = a1, va_b = b1, va_c = c1;
  double vb_a = a2, vb_b = b2, vb_c = c2;
  double vc_a = 0.5, vc_b = 0.5, vc_c = 0.5;
  
  double vd_a = vb_a - va_a + vc_a;
  double vd_b = vb_b - va_b + vc_b;
  double vd_c = vb_c - va_c + vc_c;
  double ve_a = trans_to_origuc(vd_a);
  double ve_b = trans_to_origuc(vd_b);
  double ve_c = trans_to_origuc(vd_c);
  
  double minDist = DBL_MAX;
  
  for(unsigned int index = 0; index < ivals.size(); index++){
    double vi_a = ivals[index], vi_b = jvals[index], vi_c = kvals[index];
    double da = vc_a - (ve_a+vi_a);
    double db = vc_b - (ve_b+vi_b);
    double dc = vc_c - (ve_c+vi_c);
    double dx, dy, dz;
    abc_to_xyz(da, db, dc, dx, dy, dz);
    //cout << "a: " << da << "   b: " << db << "    c: " << dc << endl;
    //cout << "x: " << dx << "   y: " << dy << "    z: " << dz << endl;
    double dist = sqrt(dx*dx + dy*dy + dz*dz);
    if(dist < minDist){
      minDist = dist;
      minDa = -da; minDb = -db; minDc = -dc;
    }
  }
  //cout << "minDist here: " << minDist << endl;
  return minDist;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rcx, %rbx
movq %rdx, %r14
movq %rsi, %r15
movq %rdi, %r12
subsd %xmm0, %xmm3
movsd 0xa2691(%rip), %xmm0    # 0xfc038
addsd %xmm0, %xmm3
subsd %xmm1, %xmm4
addsd %xmm0, %xmm4
movsd %xmm4, (%rsp)
subsd %xmm2, %xmm5
addsd %xmm0, %xmm5
movsd %xmm5, 0x8(%rsp)
movapd %xmm3, %xmm0
callq 0x6377f
movsd %xmm0, 0x10(%rsp)
movsd (%rsp), %xmm0
callq 0x6377f
movsd %xmm0, (%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x6377f
movq (%r12), %rax
movq 0x8(%r12), %rcx
subq %rax, %rcx
sarq $0x2, %rcx
movq 0x18(%r12), %rdx
movq 0x30(%r12), %rsi
movsd 0xa7cd7(%rip), %xmm1    # 0x1016e8
xorl %edi, %edi
movapd 0xa6bb5(%rip), %xmm2    # 0x1005d0
movl %edi, %edi
cmpq %rdi, %rcx
jbe 0x59af9
cvtsi2sdl (%rax,%rdi,4), %xmm4
cvtsi2sdl (%rdx,%rdi,4), %xmm5
cvtsi2sdl (%rsi,%rdi,4), %xmm6
addsd 0x10(%rsp), %xmm4
movsd 0xa25f5(%rip), %xmm7    # 0xfc038
movapd %xmm7, %xmm3
subsd %xmm4, %xmm3
addsd (%rsp), %xmm5
movapd %xmm7, %xmm4
subsd %xmm5, %xmm4
addsd %xmm0, %xmm6
movapd %xmm7, %xmm5
subsd %xmm6, %xmm5
movsd 0x50(%r12), %xmm6
mulsd %xmm4, %xmm6
movsd 0x48(%r12), %xmm7
mulsd %xmm3, %xmm7
addsd %xmm6, %xmm7
movsd 0x60(%r12), %xmm6
mulsd %xmm5, %xmm6
addsd %xmm7, %xmm6
movsd 0x68(%r12), %xmm7
mulsd %xmm5, %xmm7
movsd 0x58(%r12), %xmm8
mulsd %xmm4, %xmm8
addsd %xmm7, %xmm8
movsd 0x70(%r12), %xmm7
mulsd %xmm5, %xmm7
mulsd %xmm8, %xmm8
mulsd %xmm6, %xmm6
addsd %xmm8, %xmm6
mulsd %xmm7, %xmm7
addsd %xmm6, %xmm7
sqrtsd %xmm7, %xmm6
ucomisd %xmm6, %xmm1
jbe 0x59af2
xorpd %xmm2, %xmm3
movlpd %xmm3, (%r15)
xorpd %xmm2, %xmm4
movlpd %xmm4, (%r14)
xorpd %xmm2, %xmm5
movlpd %xmm5, (%rbx)
movapd %xmm6, %xmm1
incl %edi
jmp 0x59a1b
movapd %xmm1, %xmm0
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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 %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
movq %rdi, %rbx
movq 0xf627f(%rip), %r14     # 0x14ffc0
leaq 0xa9e61(%rip), %rsi     # 0x103ba9
movq %r14, %rdi
callq 0x4480
leaq 0xa8827(%rip), %r15     # 0x10257e
movq %rax, %rdi
movq %r15, %rsi
callq 0x4480
xorl %eax, %eax
leaq 0xaa7a7(%rip), %r12     # 0x104512
movl %eax, %r13d
movq (%rbx), %rax
movq 0x8(%rbx), %rcx
subq %rax, %rcx
sarq $0x2, %rcx
cmpq %r13, %rcx
jbe 0x59dd4
movl (%rax,%r13,4), %esi
movq %r14, %rdi
callq 0x4760
movq %rax, %rdi
movq %r12, %rsi
callq 0x4480
movq 0x18(%rbx), %rcx
movl (%rcx,%r13,4), %esi
movq %rax, %rdi
callq 0x4760
movq %rax, %rdi
movq %r12, %rsi
callq 0x4480
movq 0x30(%rbx), %rcx
movl (%rcx,%r13,4), %esi
movq %rax, %rdi
callq 0x4760
movq %rax, %rdi
movq %r15, %rsi
callq 0x4480
leal 0x1(%r13), %eax
jmp 0x59d6b
movq 0xf61e5(%rip), %rdi     # 0x14ffc0
leaq 0xa9dd9(%rip), %rsi     # 0x103bbb
callq 0x4480
leaq 0xa8790(%rip), %rbx     # 0x10257e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
movq %rax, %rdi
movq %rbx, %rsi
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
jmp 0x4480
nopl (%rax)
// File: /workspace/github_cpp_makefile_star2_O2/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 %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x10, %rsp
movq %rdi, %rbx
leaq 0x60(%rdi), %r14
xorpd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
pushq $-0x1
popq %rbp
xorl %r15d, %r15d
pushq $0x48
popq %r12
movq 0x68(%rbx), %rax
subq 0x60(%rbx), %rax
cqto
idivq %r12
cltq
cmpq %rax, %r15
jge 0x5bd7b
movq 0x140(%rbx), %rax
cmpl $-0x1, (%rax,%r15,4)
jne 0x5bd76
movq %r14, %rdi
movq %r15, %rsi
callq 0xabe2
movsd 0x38(%rax), %xmm0
ucomisd 0x8(%rsp), %xmm0
jbe 0x5bd76
movq %r14, %rdi
movq %r15, %rsi
callq 0xabe2
movsd 0x38(%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movl %r15d, %ebp
incq %r15
jmp 0x5bd23
movl %ebp, %eax
addq $0x10, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O2/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
movq %rdi, %r14
xorl %ebp, %ebp
pushq $0x1
popq %rbx
cmpl 0x138(%r14), %ebp
jge 0x5c53b
movq %r14, %rdi
movl %ebp, %esi
callq 0x5c542
incl %ebp
testl %eax, %eax
jne 0x5c520
xorl %ebx, %ebx
movl %ebx, %eax
popq %rbx
popq %r14
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O2/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 %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r14
movl 0x138(%rdi), %eax
movslq %edx, %r15
movq 0x140(%rdi), %rcx
movl %eax, (%rcx,%r15,4)
leaq 0x238(%rdi), %r12
addq $0x60, %rdi
movq %r15, %rsi
callq 0xabe2
leaq 0x38(%rax), %rsi
movq %r12, %rdi
callq 0xac88
movq %r14, %rdi
movq %rbx, %rsi
movl %r15d, %edx
movl %r15d, %ecx
callq 0x5c79c
incl 0x138(%r14)
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/feature.cc
// Function: FEATURE::segment_identify_connections_recurr(std::vector<int, std::allocator<int>>*, std::vector<int, std::allocator<int>>*, int*, int, int, std::vector<SEGCONN, std::allocator<SEGCONN>>*, int)
void FEATURE::segment_identify_connections_recurr(vector <int> *nodeid,vector <int> *edgeinface,int *nodes_togo,int curr_face,int curr_neigh_seg, vector <SEGCONN> *edges,int curr_edge)
{
int i;
// check node at current edge, if not assigned to face, assign, consider its neighbors


for(i=0;i<(int)edges->size();i++)
   {
   if(i!=curr_edge&&edgeinface->at(i)==-1) // only consider if not yet assigned to face
     {
     //only consider pairs of edges that connect to the same segment 
     if(curr_neigh_seg==edges->at(i).to_seg&&edges->at(curr_edge).deltaPos.equals(edges->at(i).deltaPos))
       {
       // for two edges connecting to another segment, two edges are on the same face
       // only if starting nodes of this face overlap
//       int node1id=nodeid->at(edges->at(curr_edge).from);
//       int node2id=nodeid->at(edges->at(i).from);

//       int node1id=(edges->at(curr_edge).from);
//       int node2id=(edges->at(i).from);

//       if(nodes[node1id].max_radius+nodes[node2id].max_radius>calcEuclideanDistance(nodes[node1id].x,nodes[node1id].y,nodes[node1id].z,nodes[node2id].x,nodes[node2id].y,nodes[node2id].z))
          { // starting nodes of edges overlap, assign edge to the current face
          edgeinface->at(i)=curr_face;
          *nodes_togo=(*nodes_togo)--;
          segment_identify_connections_recurr(nodeid,edgeinface,nodes_togo,curr_face,curr_neigh_seg,edges,i);
          };


       }; // end of if - overlapping nodes

     }; // end of if - the same segments
   }; // end of loop over all edges



}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x38, %rsp
movl %r9d, %ebp
movl %r8d, 0xc(%rsp)
movq %rcx, 0x20(%rsp)
movq %rdx, %r15
movq %rsi, 0x18(%rsp)
movq %rdi, 0x10(%rsp)
movl 0x78(%rsp), %r12d
movq 0x70(%rsp), %r14
movslq %r12d, %rax
movq %rax, 0x30(%rsp)
xorl %r13d, %r13d
pushq $0x30
popq %rbx
movq %rdx, 0x28(%rsp)
movq 0x8(%r14), %rax
subq (%r14), %rax
cqto
idivq %rbx
cltq
cmpq %rax, %r13
jge 0x5d15a
cmpq %r13, %r12
je 0x5d152
movq %r15, %rdi
movq %r13, %rsi
callq 0xac00
cmpl $-0x1, (%rax)
jne 0x5d152
movq %r14, %rdi
movq %r13, %rsi
callq 0x5f2b8
cmpl %ebp, 0xc(%rax)
jne 0x5d152
movq %r14, %rdi
movq 0x30(%rsp), %rsi
callq 0x5f2b8
movq %rbx, %r15
movq %r12, %rbx
movl %ebp, %r12d
leaq 0x24(%rax), %rbp
movq %r14, %rdi
movq %r13, %rsi
callq 0x5f2b8
movq 0x24(%rax), %rsi
movl 0x2c(%rax), %edx
movq %rbp, %rdi
movl %r12d, %ebp
movq %rbx, %r12
movq %r15, %rbx
movq 0x28(%rsp), %r15
callq 0x3d332
testb %al, %al
je 0x5d152
movq %r15, %rdi
movq %r13, %rsi
callq 0xac00
movl 0xc(%rsp), %r8d
movl %r8d, (%rax)
movq 0x10(%rsp), %rdi
movq 0x18(%rsp), %rsi
movq %r15, %rdx
movq 0x20(%rsp), %rcx
movl %ebp, %r9d
pushq %r13
pushq %r14
callq 0x5d050
popq %rax
popq %rcx
incq %r13
jmp 0x5d095
addq $0x38, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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
pushq %rax
movq %rdx, %rbx
movq %rsi, %r15
movq %rdi, %r14
movl %ecx, 0x4(%rsp)
testl %ecx, %ecx
movl $0x140, %eax            # imm = 0x140
movl $0x1d8, %ebp            # imm = 0x1D8
cmoveq %rax, %rbp
leaq 0x30(%rdi), %r12
xorl %ecx, %ecx
movl %ecx, (%rsp)
movq 0x68(%r14), %rax
subq 0x60(%r14), %rax
cqto
pushq $0x48
popq %rsi
idivq %rsi
cmpl %eax, %ecx
jge 0x5d1dd
movslq %ecx, %rax
movl (%rbx), %r13d
movq (%r14,%rbp), %rcx
addl (%rcx,%rax,4), %r13d
movq %r12, %rdi
movq %rsp, %rsi
callq 0x5f2d6
movslq (%rax), %rsi
movq %r15, %rdi
callq 0xac00
movl %r13d, (%rax)
movl (%rsp), %ecx
incl %ecx
jmp 0x5d198
xorl %eax, %eax
cmpl $0x0, 0x4(%rsp)
setne %al
shll $0x7, %eax
movl 0x138(%r14,%rax), %eax
addl %eax, (%rbx)
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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)
leaq 0x1c0(%rdi), %r14
leaq 0x250(%rdi), %r15
xorl %ecx, %ecx
leaq 0x4(%rsp), %r12
movl %ecx, 0x4(%rsp)
cmpl %eax, %ecx
jge 0x5d269
movq %r14, %rdi
movq %r12, %rsi
callq 0x9434
movslq 0x4(%rsp), %rsi
shlq $0x3, %rsi
addq 0x238(%rbx), %rsi
movq %r15, %rdi
callq 0xac88
movl 0x4(%rsp), %ecx
incl %ecx
movl 0x138(%rbx), %eax
jmp 0x5d230
leaq 0x1d8(%rbx), %r14
xorl %ecx, %ecx
pushq $0x48
popq %r15
movl %ecx, 0x4(%rsp)
movq 0x68(%rbx), %rax
subq 0x60(%rbx), %rax
cqto
idivq %r15
cmpl %eax, %ecx
jge 0x5d2a9
movslq %ecx, %rsi
shlq $0x2, %rsi
addq 0x140(%rbx), %rsi
movq %r14, %rdi
callq 0x9434
movl 0x4(%rsp), %ecx
incl %ecx
jmp 0x5d276
movq 0xf2d10(%rip), %rdi     # 0x14ffc0
leaq 0xa6af9(%rip), %rsi     # 0x103db0
callq 0x4480
movl 0x1b8(%rbx), %esi
movq %rax, %rdi
callq 0x4760
leaq 0xa52ad(%rip), %rsi     # 0x10257e
movq %rax, %rdi
callq 0x4480
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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
movq 0x140(%rdi), %r8
subq 0x60(%rdi), %rax
pushq $0x48
popq %rcx
cqto
idivq %rcx
xorl %edx, %edx
testl %eax, %eax
cmovlel %edx, %eax
xorl %ecx, %ecx
cmpq %rdx, %rax
je 0x5d31c
xorl %edi, %edi
cmpl %esi, (%r8,%rdx,4)
sete %dil
addl %edi, %ecx
incq %rdx
jmp 0x5d306
movl %ecx, %eax
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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 %esi, %ebx
movq %rdi, %r14
movl 0x1b8(%rdi), %eax
leaq 0x208(%rdi), %rcx
movq %rcx, (%rsp)
movslq %esi, %rcx
movq 0x208(%rdi), %rdx
movl %eax, (%rdx,%rcx,4)
leaq 0x158(%rdi), %r12
xorl %r13d, %r13d
pushq $0x30
popq %r15
movq 0x160(%r14), %rax
subq 0x158(%r14), %rax
cqto
idivq %r15
cltq
cmpq %rax, %r13
jge 0x5d6c1
movq %r12, %rdi
movq %r13, %rsi
callq 0x5f2b8
cmpl %ebx, 0x8(%rax)
je 0x5d66a
movq %r12, %rdi
movq %r13, %rsi
callq 0x5f2b8
cmpl %ebx, 0xc(%rax)
jne 0x5d6b9
movq %r12, %rdi
movq %r13, %rsi
callq 0x5f2b8
cmpl $0x1, 0x20(%rax)
jne 0x5d6b9
movq %r12, %rdi
movq %r13, %rsi
callq 0x5f2b8
xorl %ebp, %ebp
cmpl %ebx, 0x8(%rax)
sete %bpl
movq %r12, %rdi
movq %r13, %rsi
callq 0x5f2b8
movl 0x8(%rax,%rbp,4), %ebp
movslq %ebp, %rsi
movq (%rsp), %rdi
callq 0xac00
cmpl $-0x1, (%rax)
jne 0x5d6b9
movq %r14, %rdi
movl %ebp, %esi
callq 0x5d5f4
incq %r13
jmp 0x5d630
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O2/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 0x5fb44
callq 0x4590
xorl %r14d, %r14d
movslq 0xa4(%rbx), %rax
cmpq %rax, %r14
jge 0x5fb6d
movq 0x98(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x5fb68
callq 0x4590
incq %r14
jmp 0x5fb47
xorl %r14d, %r14d
movslq %eax, %rcx
cmpq %rcx, %r14
jge 0x5fb98
movq 0x90(%rbx), %rcx
movq (%rcx,%r14,8), %rdi
testq %rdi, %rdi
je 0x5fb93
callq 0x4590
movl 0xa4(%rbx), %eax
incq %r14
jmp 0x5fb70
xorl %r14d, %r14d
movslq %eax, %rcx
cmpq %rcx, %r14
jge 0x5fbc3
movq 0x80(%rbx), %rcx
movq (%rcx,%r14,8), %rdi
testq %rdi, %rdi
je 0x5fbbe
callq 0x4590
movl 0xa4(%rbx), %eax
incq %r14
jmp 0x5fb9b
movq 0xc8(%rbx), %rdi
testq %rdi, %rdi
je 0x5fbd4
callq 0x4590
movq 0xc0(%rbx), %rdi
testq %rdi, %rdi
je 0x5fbe5
callq 0x4590
movq 0xb8(%rbx), %rdi
testq %rdi, %rdi
je 0x5fbf6
callq 0x4590
movq 0xb0(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc07
callq 0x4590
movq 0xa8(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc18
callq 0x4590
movq 0x90(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc29
callq 0x4590
movq 0x98(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc3a
callq 0x4590
movq 0x88(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc4b
callq 0x4590
movq 0x80(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc5c
callq 0x4590
xorl %r14d, %r14d
movslq 0x3c(%rbx), %rax
cmpq %rax, %r14
jge 0x5fc91
movq 0x68(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x5fc7a
callq 0x4590
movq 0x60(%rbx), %rax
movq (%rax,%r14,8), %rdi
testq %rdi, %rdi
je 0x5fc8c
callq 0x4590
incq %r14
jmp 0x5fc5f
movq 0x78(%rbx), %rdi
testq %rdi, %rdi
je 0x5fc9f
callq 0x4590
movq 0x70(%rbx), %rdi
testq %rdi, %rdi
je 0x5fcad
callq 0x4590
movq 0x68(%rbx), %rdi
testq %rdi, %rdi
je 0x5fcbb
callq 0x4590
movq 0x60(%rbx), %rdi
addq $0x8, %rsp
testq %rdi, %rdi
je 0x5fcd0
popq %rbx
popq %r14
jmp 0x4590
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O2/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), %r15
movslq %esi, %r12
movl (%r15,%r12,4), %eax
addl %eax, %eax
movl %eax, (%r15,%r12,4)
cmpl $0x10001, %eax          # imm = 0x10001
jl 0x5fd0f
leaq 0xa4221(%rip), %rdi     # 0x103f20
pushq $0x2
popq %rsi
callq 0xe0d0f
movq 0x78(%rbx), %r15
movl (%r15,%r12,4), %eax
leal (,%rax,4), %ecx
movslq %ecx, %rdi
shlq $0x3, %rdi
testl %eax, %eax
pushq $-0x1
popq %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 0x60(%rbx), %rax
movq 0x70(%rbx), %rcx
movslq (%rcx,%r12,4), %rcx
leal (,%rcx,4), %edx
xorl %esi, %esi
testl %edx, %edx
cmovlel %esi, %edx
cmpq %rsi, %rdx
je 0x5fd7d
movq (%rax,%r12,8), %rdi
movsd (%rdi,%rsi,8), %xmm0
movsd %xmm0, (%r14,%rsi,8)
incq %rsi
jmp 0x5fd64
movq 0x68(%rbx), %rdx
xorl %esi, %esi
cmpq %rcx, %rsi
jge 0x5fd98
movq (%rdx,%r12,8), %rdi
movl (%rdi,%rsi,4), %edi
movl %edi, (%r15,%rsi,4)
incq %rsi
jmp 0x5fd83
movq (%rax,%r12,8), %rdi
testq %rdi, %rdi
je 0x5fdaa
callq 0x4590
movq 0x68(%rbx), %rdx
movq (%rdx,%r12,8), %rdi
testq %rdi, %rdi
je 0x5fdb8
callq 0x4590
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
// File: /workspace/github_cpp_makefile_star2_O2/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), %r14
movslq %esi, %r13
movl (%r14,%r13,4), %eax
addl %eax, %eax
movl %eax, (%r14,%r13,4)
cmpl $0x801, %eax            # imm = 0x801
jl 0x6010f
leaq 0xa3e58(%rip), %rdi     # 0x103f54
pushq $0x2
popq %rsi
callq 0xe0d0f
movq 0xb8(%rbx), %r14
movl (%r14,%r13,4), %eax
leal (%rax,%rax), %ecx
movslq %ecx, %rdi
shlq $0x2, %rdi
testl %eax, %eax
pushq $-0x1
popq %rax
cmovsq %rax, %rdi
movq %rax, %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
pushq $-0x1
popq %rax
cmovsq %rax, %rdi
callq 0x4030
movq %rax, %r14
movq 0xa8(%rbx), %rsi
movq 0x80(%rbx), %rax
movq 0x90(%rbx), %rcx
movq 0x98(%rbx), %rdx
movslq (%rsi,%r13,4), %rsi
xorl %edi, %edi
xorl %r8d, %r8d
cmpq %rsi, %r8
jge 0x601be
movq (%rax,%r13,8), %r9
movl (%r9,%rdi), %r9d
movl %r9d, (%r12,%rdi)
movq (%rcx,%r13,8), %r9
movups (%r9,%rdi,4), %xmm0
movups %xmm0, (%r15,%rdi,4)
movq (%rdx,%r13,8), %r9
movl (%r9,%rdi), %r9d
movl %r9d, (%r14,%rdi)
incq %r8
addq $0x4, %rdi
jmp 0x6018a
leaq (%r12,%rbp,4), %rbp
movq 0x88(%rbx), %rcx
movq 0xb0(%rbx), %rdx
movslq (%rdx,%r13,4), %rdx
xorl %esi, %esi
cmpq %rdx, %rsi
jge 0x601eb
movq (%rcx,%r13,8), %rdi
movl (%rdi,%rsi,4), %edi
movl %edi, (%rbp,%rsi,4)
incq %rsi
jmp 0x601d6
movq (%rax,%r13,8), %rdi
testq %rdi, %rdi
je 0x60200
callq 0x4590
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 0x6022b
callq 0x4590
movq 0x90(%rbx), %rax
movq %r15, (%rax,%r13,8)
movq 0x98(%rbx), %rax
movq (%rax,%r13,8), %rdi
testq %rdi, %rdi
je 0x6024b
callq 0x4590
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
// File: /workspace/github_cpp_makefile_star2_O2/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:
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 0xa3c78(%rip), %r15     # 0x103f89
movslq 0xa0(%r14), %rax
cmpq %rax, %r12
jge 0x6046e
movq 0x60(%r14), %rax
movq 0xc0(%r14), %rcx
movslq (%rcx,%r12,4), %rcx
movq (%rax,%rcx,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 0xa8(%r14), %rax
movslq (%rax,%r12,4), %rax
cmpq %rax, %r13
jge 0x60466
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 0x603c6
testl %edx, %edx
jne 0x603c6
testl %ecx, %ecx
jne 0x603c6
testl %eax, %eax
je 0x6045e
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
cvtsi2sd %edx, %xmm0
mulsd (%r14), %xmm0
addsd (%rsi,%rdi), %xmm0
cvtsi2sd %ecx, %xmm1
movsd 0x8(%r14), %xmm2
mulsd %xmm1, %xmm2
addsd %xmm0, %xmm2
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 0x4610
incq %r13
jmp 0x60368
incq %r12
jmp 0x60311
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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 0x6379d
movsd 0x988a7(%rip), %xmm1    # 0xfc038
subsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
negl %eax
jmp 0x637ad
movsd 0x98893(%rip), %xmm1    # 0xfc038
addsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
cvtsi2sd %eax, %xmm1
subsd %xmm1, %xmm0
movsd 0x9889b(%rip), %xmm2    # 0xfc058
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_O2/Luthaf[P]zeopp/geometry.cc
// Function: calcDeterminant(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), %xmm1
movsd 0x20(%rdi), %xmm2
movsd 0x38(%rdi), %xmm3
movsd 0x28(%rdi), %xmm0
movapd %xmm0, %xmm4
mulsd %xmm3, %xmm4
movapd %xmm1, %xmm5
mulsd %xmm2, %xmm5
subsd %xmm4, %xmm5
movsd 0x8(%rdi), %xmm4
movsd 0x10(%rdi), %xmm6
mulsd %xmm6, %xmm3
mulsd %xmm4, %xmm1
subsd %xmm3, %xmm1
mulsd 0x18(%rdi), %xmm1
mulsd (%rdi), %xmm5
subsd %xmm1, %xmm5
mulsd %xmm2, %xmm6
mulsd %xmm4, %xmm0
subsd %xmm6, %xmm0
mulsd 0x30(%rdi), %xmm0
addsd %xmm5, %xmm0
retq
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/geometry.cc
// Function: tryInvertMatrix(double (*) [3], double (*) [3])
bool tryInvertMatrix(double matrix [3][3], double newMatrix[3][3]){
  double determinant = calcDeterminant(matrix);
  double invDet = 0;
  bool ok = false;
  if(determinant!=0) {
    invDet = 1/determinant;
    ok = true;
  }
  newMatrix[0][0] = invDet*   (matrix[2][2]*matrix[1][1]-matrix[2][1]*matrix[1][2]);
  newMatrix[0][1] = invDet*-1*(matrix[2][2]*matrix[0][1]-matrix[2][1]*matrix[0][2]);
  newMatrix[0][2] = invDet*   (matrix[1][2]*matrix[0][1]-matrix[1][1]*matrix[0][2]);
  newMatrix[1][0] = invDet*-1*(matrix[2][2]*matrix[1][0]-matrix[2][0]*matrix[1][2]);
  newMatrix[1][1] = invDet*   (matrix[2][2]*matrix[0][0]-matrix[2][0]*matrix[0][2]);
  newMatrix[1][2] = invDet*-1*(matrix[1][2]*matrix[0][0]-matrix[1][0]*matrix[0][2]);
  newMatrix[2][0] = invDet*   (matrix[2][1]*matrix[1][0]-matrix[2][0]*matrix[1][1]);
  newMatrix[2][1] = invDet*-1*(matrix[2][1]*matrix[0][0]-matrix[2][0]*matrix[0][1]);
  newMatrix[2][2] = invDet*   (matrix[1][1]*matrix[0][0]-matrix[1][0]*matrix[0][1]);
  return ok;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r14
callq 0x63835
xorpd %xmm1, %xmm1
movapd %xmm0, %xmm2
cmpneqsd %xmm1, %xmm2
movq %xmm2, %rax
andl $0x1, %eax
ucomisd %xmm1, %xmm0
movsd 0x9857a(%rip), %xmm1    # 0xfc058
divsd %xmm0, %xmm1
jne 0x63aea
jp 0x63aea
xorpd %xmm1, %xmm1
movsd 0x28(%r14), %xmm0
movsd 0x40(%r14), %xmm2
mulsd 0x38(%r14), %xmm0
mulsd 0x20(%r14), %xmm2
subsd %xmm0, %xmm2
mulsd %xmm1, %xmm2
movsd %xmm2, (%rbx)
movapd 0x9caba(%rip), %xmm0    # 0x1005d0
xorpd %xmm1, %xmm0
movsd 0x10(%r14), %xmm2
mulsd 0x38(%r14), %xmm2
movsd 0x40(%r14), %xmm3
mulsd 0x8(%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x8(%rbx)
movsd 0x10(%r14), %xmm2
movsd 0x28(%r14), %xmm3
mulsd 0x20(%r14), %xmm2
mulsd 0x8(%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x10(%rbx)
movsd 0x28(%r14), %xmm2
movsd 0x40(%r14), %xmm3
mulsd 0x30(%r14), %xmm2
mulsd 0x18(%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x18(%rbx)
movsd 0x10(%r14), %xmm2
mulsd 0x30(%r14), %xmm2
movsd 0x40(%r14), %xmm3
mulsd (%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x20(%rbx)
movsd 0x10(%r14), %xmm2
movsd 0x28(%r14), %xmm3
mulsd 0x18(%r14), %xmm2
mulsd (%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x28(%rbx)
movsd 0x20(%r14), %xmm2
movsd 0x38(%r14), %xmm3
mulsd 0x30(%r14), %xmm2
mulsd 0x18(%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm1, %xmm3
movsd %xmm3, 0x30(%rbx)
movsd 0x8(%r14), %xmm2
mulsd 0x30(%r14), %xmm2
movsd 0x38(%r14), %xmm3
mulsd (%r14), %xmm3
subsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x38(%rbx)
movsd 0x8(%r14), %xmm0
movsd 0x20(%r14), %xmm2
mulsd 0x18(%r14), %xmm0
mulsd (%r14), %xmm2
subsd %xmm0, %xmm2
mulsd %xmm1, %xmm2
movsd %xmm2, 0x40(%rbx)
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O2/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 0x63c84
movl %esi, %eax
leaq (%rdi,%rax,8), %rax
retq
pushq %rbx
movq 0xec364(%rip), %rdi     # 0x14fff0
leaq 0xa03ac(%rip), %rsi     # 0x10403f
callq 0x4480
leaq 0x9e8df(%rip), %rbx     # 0x10257e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
leaq 0x9f971(%rip), %rsi     # 0x103622
movq %rax, %rdi
callq 0x4480
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
pushq $0x1
popq %rdi
callq 0x524c
// File: /workspace/github_cpp_makefile_star2_O2/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 0x98168(%rip), %xmm2    # 0xfc0a0
ucomisd 0x9d2a8(%rip), %xmm2    # 0x1011e8
ja 0x63f94
leaq 0x20(%rsp), %rax
leaq 0x8(%rsp), %rcx
movsd 0x8(%rcx), %xmm0
movsd 0x8(%rax), %xmm1
movapd %xmm0, %xmm2
subsd %xmm1, %xmm2
andpd 0x9813a(%rip), %xmm2    # 0xfc0a0
ucomisd 0x9d27a(%rip), %xmm2    # 0x1011e8
ja 0x63f94
movsd 0x10(%rcx), %xmm0
movsd 0x10(%rax), %xmm1
movapd %xmm0, %xmm2
subsd %xmm1, %xmm2
andpd 0x98116(%rip), %xmm2    # 0xfc0a0
ucomisd 0x9d256(%rip), %xmm2    # 0x1011e8
jbe 0x63f9c
ucomisd %xmm0, %xmm1
seta %al
retq
xorl %eax, %eax
retq
nop
// File: /workspace/github_cpp_makefile_star2_O2/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 0x63fbe
movl %esi, %eax
leaq (%rdi,%rax,4), %rax
retq
pushq %rbx
movq 0xec02a(%rip), %rdi     # 0x14fff0
leaq 0xa00a9(%rip), %rsi     # 0x104076
callq 0x4480
leaq 0x9e5a5(%rip), %rbx     # 0x10257e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
leaq 0x9f637(%rip), %rsi     # 0x103622
movq %rax, %rdi
callq 0x4480
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
pushq $0x1
popq %rdi
callq 0x524c
// File: /workspace/github_cpp_makefile_star2_O2/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 %r14
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r14
callq 0x6418c
movupd (%r14), %xmm1
movupd 0x8(%r14), %xmm2
movhpd 0x8(%rbx), %xmm2        # xmm2 = xmm2[0],mem[0]
movhpd (%rbx), %xmm1           # xmm1 = xmm1[0],mem[0]
mulpd %xmm2, %xmm2
mulpd %xmm1, %xmm1
addpd %xmm2, %xmm1
movsd 0x10(%r14), %xmm2
movhpd 0x10(%rbx), %xmm2       # xmm2 = xmm2[0],mem[0]
mulpd %xmm2, %xmm2
addpd %xmm1, %xmm2
sqrtpd %xmm2, %xmm1
movapd %xmm1, %xmm2
unpckhpd %xmm1, %xmm2            # xmm2 = xmm2[1],xmm1[1]
mulsd %xmm1, %xmm2
divsd %xmm2, %xmm0
movsd 0x97e48(%rip), %xmm1    # 0xfc058
ucomisd %xmm1, %xmm0
ja 0x6422c
movsd 0x9c0ca(%rip), %xmm2    # 0x1002e8
ucomisd %xmm0, %xmm2
movapd %xmm0, %xmm1
jbe 0x6422c
movapd %xmm2, %xmm1
movapd %xmm1, %xmm0
callq 0x4450
movapd %xmm0, %xmm1
cmpordsd %xmm0, %xmm1
andpd %xmm1, %xmm0
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O2/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 0x64262
cmpl $0x2, %esi
je 0x6425e
cmpl $0x1, %esi
jne 0x64266
addq $0x8, %rdi
jmp 0x64262
addq $0x10, %rdi
movq %rdi, %rax
retq
pushq %rbx
movq 0xebd82(%rip), %rdi     # 0x14fff0
leaq 0x9fe3a(%rip), %rsi     # 0x1040af
callq 0x4480
leaq 0x9e2fd(%rip), %rbx     # 0x10257e
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
leaq 0x9f38f(%rip), %rsi     # 0x103622
movq %rax, %rdi
callq 0x4480
movq %rax, %rdi
movq %rbx, %rsi
callq 0x4480
pushq $0x1
popq %rdi
callq 0x524c
// File: /workspace/github_cpp_makefile_star2_O2/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), %xmm0
movsd 0x8(%rdi), %xmm1
subsd (%rsi), %xmm0
subsd 0x8(%rsi), %xmm1
movsd 0x10(%rdi), %xmm2
subsd 0x10(%rsi), %xmm2
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
mulsd %xmm2, %xmm2
addsd %xmm0, %xmm2
sqrtsd %xmm2, %xmm0
retq
// File: /workspace/github_cpp_makefile_star2_O2/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:
movq %rdi, %rax
movsd 0x30(%rsp), %xmm0
movapd %xmm0, %xmm1
subsd 0x18(%rsp), %xmm1
movsd 0x48(%rsp), %xmm2
subsd %xmm0, %xmm2
movupd 0x20(%rsp), %xmm3
movsd 0x8(%rsp), %xmm4
movapd %xmm3, %xmm5
subsd %xmm4, %xmm5
movapd %xmm3, %xmm6
unpckhpd %xmm3, %xmm6            # xmm6 = xmm6[1],xmm3[1]
subsd 0x10(%rsp), %xmm6
movapd 0x38(%rsp), %xmm4
subpd %xmm3, %xmm4
movapd %xmm4, %xmm7
unpckhpd %xmm4, %xmm7            # xmm7 = xmm7[1],xmm4[1]
mulsd %xmm7, %xmm6
mulsd %xmm7, %xmm7
movapd %xmm4, %xmm8
mulsd %xmm4, %xmm8
addsd %xmm7, %xmm8
movapd %xmm2, %xmm7
mulsd %xmm2, %xmm7
addsd %xmm8, %xmm7
mulsd %xmm4, %xmm5
addsd %xmm6, %xmm5
mulsd %xmm2, %xmm1
addsd %xmm5, %xmm1
xorpd 0x9c1a5(%rip), %xmm1    # 0x1005d0
sqrtsd %xmm7, %xmm5
mulsd %xmm5, %xmm5
divsd %xmm5, %xmm1
movapd %xmm1, %xmm5
unpcklpd %xmm1, %xmm5            # xmm5 = xmm5[0],xmm1[0]
mulpd %xmm4, %xmm5
addpd %xmm3, %xmm5
movupd %xmm5, (%rdi)
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rdi)
retq
// File: /workspace/github_cpp_makefile_star2_O2/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
movupd (%rcx), %xmm3
movsd 0x8(%rcx), %xmm4
mulsd %xmm2, %xmm4
unpcklpd %xmm4, %xmm0            # xmm0 = xmm0[0],xmm4[0]
movsd 0x10(%rcx), %xmm4
mulsd %xmm2, %xmm4
unpcklpd %xmm1, %xmm2            # xmm2 = xmm2[0],xmm1[0]
movhpd 0x8(%rdx), %xmm3        # xmm3 = xmm3[0],mem[0]
mulpd %xmm2, %xmm3
addpd %xmm0, %xmm3
movupd %xmm3, (%rdi)
movsd %xmm4, 0x10(%rdi)
retq
// File: /workspace/github_cpp_makefile_star2_O2/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:
subq $0x68, %rsp
movsd %xmm0, 0x60(%rsp)
movq 0x80(%rsp), %rax
movq %rax, 0x10(%rsp)
movaps 0x70(%rsp), %xmm0
movups %xmm0, (%rsp)
callq 0x650f0
movsd %xmm0, 0x20(%rsp)
movsd %xmm1, 0x30(%rsp)
movq 0x98(%rsp), %rax
movq %rax, 0x10(%rsp)
movups 0x88(%rsp), %xmm0
movups %xmm0, (%rsp)
callq 0x650f0
movsd %xmm0, 0x38(%rsp)
movsd 0x30(%rsp), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, 0x30(%rsp)
movsd 0x20(%rsp), %xmm0
callq 0x47e0
movsd %xmm0, 0x28(%rsp)
movsd 0x30(%rsp), %xmm0
callq 0x4110
mulsd 0x28(%rsp), %xmm0
mulsd %xmm0, %xmm0
movsd %xmm0, 0x28(%rsp)
movsd 0x38(%rsp), %xmm0
callq 0x47e0
movsd %xmm0, 0x48(%rsp)
movsd 0x20(%rsp), %xmm0
callq 0x4110
movsd %xmm0, 0x40(%rsp)
movsd 0x38(%rsp), %xmm0
callq 0x4110
movsd %xmm0, 0x58(%rsp)
movsd 0x20(%rsp), %xmm0
callq 0x47e0
movsd %xmm0, 0x50(%rsp)
movsd 0x30(%rsp), %xmm0
callq 0x47e0
movsd 0x50(%rsp), %xmm1
mulsd 0x58(%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd 0x40(%rsp), %xmm0
mulsd 0x48(%rsp), %xmm0
subsd %xmm1, %xmm0
mulsd %xmm0, %xmm0
addsd 0x28(%rsp), %xmm0
sqrtsd %xmm0, %xmm0
movsd %xmm0, 0x28(%rsp)
movsd 0x20(%rsp), %xmm0
callq 0x4110
movsd %xmm0, 0x48(%rsp)
movsd 0x38(%rsp), %xmm0
callq 0x4110
movsd %xmm0, 0x40(%rsp)
movsd 0x20(%rsp), %xmm0
callq 0x47e0
movsd %xmm0, 0x20(%rsp)
movsd 0x38(%rsp), %xmm0
callq 0x47e0
mulsd 0x20(%rsp), %xmm0
movsd %xmm0, 0x20(%rsp)
movsd 0x30(%rsp), %xmm0
callq 0x47e0
mulsd 0x20(%rsp), %xmm0
movsd 0x40(%rsp), %xmm1
mulsd 0x48(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd 0x28(%rsp), %xmm0
divsd %xmm1, %xmm0
callq 0x43c0
mulsd 0x60(%rsp), %xmm0
addq $0x68, %rsp
retq
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/geometry.cc
// Function: findLongAndLat(Point)
pair<double,double> findLongAndLat(Point pt){
  double r = sqrt(pt[0]*pt[0] + pt[1]*pt[1]);
  double phi = atan(pt[2]/r);
  double lambda = atan2(pt[1],pt[0]);
  return pair<double,double> (lambda, phi);
}

// Assembly:
subq $0x18, %rsp
movsd 0x20(%rsp), %xmm1
movsd %xmm1, 0x8(%rsp)
movsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
mulsd %xmm0, %xmm0
mulsd %xmm1, %xmm1
addsd %xmm0, %xmm1
sqrtsd %xmm1, %xmm1
movsd 0x30(%rsp), %xmm0
divsd %xmm1, %xmm0
callq 0x43c0
movsd %xmm0, (%rsp)
movsd 0x10(%rsp), %xmm0
movsd 0x8(%rsp), %xmm1
callq 0x4090
movsd (%rsp), %xmm1
addq $0x18, %rsp
retq
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/geometry.cc
// Function: projectPointOnSphere(Point, double, Point)
Point projectPointOnSphere(Point p, double circRad, Point circCenter){
  double deltaX = p[0] - circCenter[0];
  double deltaY = p[1] - circCenter[1];
  double deltaZ = p[2] - circCenter[2];
  double factor = sqrt((circRad*circRad)/(deltaX*deltaX + deltaY*deltaY + deltaZ*deltaZ));
  return Point(factor*deltaX+circCenter[0], factor*deltaY+circCenter[1], factor*deltaZ+circCenter[2]);
}

// Assembly:
movq %rdi, %rax
movsd 0x18(%rsp), %xmm1
movsd 0x30(%rsp), %xmm2
subsd %xmm2, %xmm1
mulsd %xmm0, %xmm0
movapd 0x8(%rsp), %xmm3
movupd 0x20(%rsp), %xmm4
subpd %xmm4, %xmm3
movapd %xmm3, %xmm5
mulpd %xmm3, %xmm5
unpckhpd %xmm5, %xmm5            # xmm5 = xmm5[1,1]
movapd %xmm3, %xmm6
mulsd %xmm3, %xmm6
addsd %xmm5, %xmm6
movapd %xmm1, %xmm5
mulsd %xmm1, %xmm5
addsd %xmm6, %xmm5
divsd %xmm5, %xmm0
sqrtsd %xmm0, %xmm0
mulsd %xmm0, %xmm1
unpcklpd %xmm0, %xmm0            # xmm0 = xmm0[0,0]
mulpd %xmm3, %xmm0
addpd %xmm4, %xmm0
addsd %xmm2, %xmm1
movupd %xmm0, (%rdi)
movsd %xmm1, 0x10(%rdi)
retq
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/holograms.cc
// Function: get_bin(float, float*)
int get_bin(float measure, float *upper_bounds) {
	int i;
	for(i=0; i<NUM_BINS-1; i++) {
		if(measure<upper_bounds[i]) {
			return i;
		}
	}
	return NUM_BINS-1; //i.e. the final bin will contain everything too large for the others
}

// Assembly:
xorl %ecx, %ecx
pushq $0xf
popq %rax
cmpq $0xf, %rcx
je 0x66109
movss (%rdi,%rcx,4), %xmm1
ucomiss %xmm0, %xmm1
ja 0x66107
incq %rcx
jmp 0x660f2
movl %ecx, %eax
retq
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/grid.cc
// Function: write_distances(_IO_FILE*, double***, int, int, int)
void write_distances(FILE *f, double ***grid, int x_grid_steps, int y_grid_steps, int z_grid_steps) {
	vector <double> visit_friendly_array;
	for(int k=0; k<z_grid_steps; k++) {
		for(int j=0; j<y_grid_steps; j++) {
			for(int i=0; i<x_grid_steps; i++) { //for some reason, the convention is to write out the data in the z,y,x order, i.e., when you read the data, x iterates the fastest, and z the slowest
				visit_friendly_array.push_back(grid[i][j][k]);
			}
		}
	}
	fwrite(&visit_friendly_array[0], sizeof(double), x_grid_steps*y_grid_steps*z_grid_steps, f);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movq %rsi, %r12
movq %rdi, 0x10(%rsp)
xorps %xmm0, %xmm0
leaq 0x30(%rsp), %r13
movaps %xmm0, (%r13)
andq $0x0, 0x10(%r13)
xorl %eax, %eax
testl %edx, %edx
movl $0x0, %ebp
movl %edx, 0x4(%rsp)
cmovgl %edx, %ebp
testl %ecx, %ecx
movl $0x0, %edx
movl %ecx, 0x8(%rsp)
cmovgl %ecx, %edx
movq %rdx, 0x28(%rsp)
testl %r8d, %r8d
movl $0x0, %ecx
movl %r8d, 0xc(%rsp)
cmovgl %r8d, %ecx
movq %rcx, 0x18(%rsp)
cmpq 0x18(%rsp), %rax
je 0x66d22
movq %rax, 0x20(%rsp)
leaq (,%rax,8), %r14
xorl %ebx, %ebx
cmpq 0x28(%rsp), %rbx
je 0x66d18
xorl %r15d, %r15d
cmpq %r15, %rbp
je 0x66d13
movq (%r12,%r15,8), %rax
movq (%rax,%rbx,8), %rsi
addq %r14, %rsi
movq %r13, %rdi
callq 0xac88
incq %r15
jmp 0x66cf6
incq %rbx
jmp 0x66cec
movq 0x20(%rsp), %rax
incq %rax
jmp 0x66cd6
leaq 0x30(%rsp), %r12
movq (%r12), %rdi
movl 0x8(%rsp), %eax
imull 0x4(%rsp), %eax
imull 0xc(%rsp), %eax
movslq %eax, %rdx
pushq $0x8
popq %rsi
movq 0x10(%rsp), %rcx
callq 0x4840
movq %r12, %rdi
callq 0xb25c
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rax, %rbx
leaq 0x30(%rsp), %rdi
callq 0xb25c
movq %rbx, %rdi
callq 0x4780
// File: /workspace/github_cpp_makefile_star2_O2/Luthaf[P]zeopp/grid.cc
// Function: generateGaussianGrid(ATOM_NETWORK*, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char>>, bool, bool)
void generateGaussianGrid(ATOM_NETWORK *atmnet, string cubefilename, bool angstrom_to_bohr, bool useMassFlag){

 GaussianCube cube(atmnet);

 cube.calculateDistanceGrid(atmnet);
 cube.writeGrid(atmnet,cubefilename, angstrom_to_bohr, useMassFlag);
 cube.deinit();

}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xb8, %rsp
movl %ecx, %ebx
movl %edx, %ebp
movq %rsi, %r14
movq %rdi, %r15
leaq 0x28(%rsp), %r12
movq %r12, %rdi
movq %r15, %rsi
callq 0x67e56
movq %r12, %rdi
movq %r15, %rsi
callq 0x66e22
leaq 0x8(%rsp), %r13
movq %r13, %rdi
movq %r14, %rsi
callq 0x45d0
movq %r12, %rdi
movq %r15, %rsi
movq %r13, %rdx
movl %ebp, %ecx
movl %ebx, %r8d
callq 0x66fce
leaq 0x8(%rsp), %rdi
callq 0x4210
leaq 0x28(%rsp), %rdi
callq 0x6811e
addq $0xb8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rax, %rbx
leaq 0x8(%rsp), %rdi
callq 0x4210
movq %rbx, %rdi
callq 0x4780