LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O3/hariguchi[P]rtacl/./rtacl.hpp
// Function: rtacl::db<long>::makeMin(sockaddr_in const&, sockaddr_in const&, unsigned char, unsigned char, boost::geometry::model::point<long, 6ul, boost::geometry::cs::cartesian>&)
inline void
db<ADDR>::makeTuple (const sockaddr_in& src,
                     const sockaddr_in& dst,
                     const u8 proto,
                     const u8 dscp,
                     const s32 offset,
                     tuple<rtacl::ipv4a>& result)
{
    assert(af == AF_INET);

    s64 val;

    val = static_cast<s64>(ntohl(src.sin_addr.s_addr)) + offset;
    bg::set<0>(result, val);
    val = static_cast<s64>(ntohl(dst.sin_addr.s_addr)) + offset;
    bg::set<1>(result, val);
    val = static_cast<s64>(ntohs(src.sin_port)) + offset;
    bg::set<2>(result, val);
    val = static_cast<s64>(ntohs(dst.sin_port)) + offset;
    bg::set<3>(result, val);
    val = static_cast<s64>(proto) + offset;
    bg::set<4>(result, val);
    val = static_cast<s64>(dscp) + offset;
    bg::set<5>(result, val);
}

// Assembly:
cmpw $0x2, 0x18(%rdi)
jne 0x14d06
movl 0x4(%rsi), %eax
bswapl %eax
decq %rax
movq %rax, (%r9)
movl 0x4(%rdx), %eax
bswapl %eax
decq %rax
movq %rax, 0x8(%r9)
movzwl 0x2(%rsi), %eax
rolw $0x8, %ax
movzwl %ax, %eax
decq %rax
movq %rax, 0x10(%r9)
movzwl 0x2(%rdx), %eax
rolw $0x8, %ax
movzwl %ax, %eax
decq %rax
movq %rax, 0x18(%r9)
movzbl %cl, %eax
decq %rax
movq %rax, 0x20(%r9)
movzbl %r8b, %eax
decq %rax
movq %rax, 0x28(%r9)
retq
pushq %rax
leaq 0x27544(%rip), %rdi     # 0x3c252
leaq 0x20841(%rip), %rsi     # 0x35556
leaq 0x27544(%rip), %rcx     # 0x3c260
movl $0x25b, %edx            # imm = 0x25B
callq 0x4140
// File: /workspace/github_cpp_makefile_star2_O3/hariguchi[P]rtacl/./rtacl.hpp
// Function: rtacl::db<boost::multiprecision::number<boost::multiprecision::backends::cpp_int_backend<256u, 256u, (boost::multiprecision::cpp_integer_type)1, (boost::multiprecision::cpp_int_check_type)0, void>, (boost::multiprecision::expression_template_option)0>>::find(boost::geometry::model::point<boost::multiprecision::number<boost::multiprecision::backends::cpp_int_backend<256u, 256u, (boost::multiprecision::cpp_integer_type)1, (boost::multiprecision::cpp_int_check_type)0, void>, (boost::multiprecision::expression_template_option)0>, 6ul, boost::geometry::cs::cartesian> const&)
inline result<ADDR>
db<ADDR>::find (const tuple<ADDR>& key)
{
    result<ADDR> r;
    size_t n = rtree.query(bgi::contains(key), std::back_inserter(r));
    if (n != r.size()) {
        fprintf(stderr, "n(%ld) != r.size()(%ld)\n", n, r.size());
    }
    assert(n == r.size());

    return r;
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
subq $0x130, %rsp            # imm = 0x130
movq %rdi, %rbx
xorps %xmm0, %xmm0
movups %xmm0, (%rdi)
movq $0x0, 0x10(%rdi)
movl $0x22, %eax
movaps -0x22(%rdx,%rax), %xmm0
movaps -0x12(%rdx,%rax), %xmm1
movaps %xmm1, -0x2(%rsp,%rax)
movaps %xmm0, -0x12(%rsp,%rax)
movzwl -0x2(%rdx,%rax), %ecx
movw %cx, 0xe(%rsp,%rax)
movb (%rdx,%rax), %cl
movb %cl, 0x10(%rsp,%rax)
addq $0x30, %rax
cmpq $0x142, %rax            # imm = 0x142
jne 0x157f2
cmpq $0x0, 0x10(%rsi)
je 0x1584e
leaq 0x10(%rsp), %rax
leaq 0xf(%rsp), %rcx
movq %rsi, %rdi
movq %rax, %rsi
movq %rbx, %rdx
callq 0x21a1c
movq %rax, %r14
movq (%rbx), %rax
movq 0x8(%rbx), %rcx
jmp 0x15855
xorl %eax, %eax
xorl %ecx, %ecx
xorl %r14d, %r14d
subq %rax, %rcx
sarq $0x4, %rcx
movabsq $0x14c1bacf914c1bad, %r15 # imm = 0x14C1BACF914C1BAD
imulq %r15, %rcx
cmpq %rcx, %r14
jne 0x15893
movq 0x8(%rbx), %rcx
subq %rax, %rcx
sarq $0x4, %rcx
imulq %r15, %rcx
cmpq %rcx, %r14
jne 0x158b3
movq %rbx, %rax
addq $0x130, %rsp            # imm = 0x130
popq %rbx
popq %r14
popq %r15
retq
movq 0x5872e(%rip), %rax     # 0x6dfc8
movq (%rax), %rdi
leaq 0x1fc8b(%rip), %rsi     # 0x3552f
movq %r14, %rdx
xorl %eax, %eax
callq 0x42d0
movq (%rbx), %rax
jmp 0x1586f
leaq 0x1fc8e(%rip), %rdi     # 0x35548
leaq 0x1fc95(%rip), %rsi     # 0x35556
leaq 0x428d3(%rip), %rcx     # 0x5819b
movl $0x23d, %edx            # imm = 0x23D
callq 0x4140
movq %rax, %r14
movq (%rbx), %rdi
testq %rdi, %rdi
je 0x158e9
movq 0x10(%rbx), %rsi
subq %rdi, %rsi
callq 0x4210
movq %r14, %rdi
callq 0x4370
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: projsimplex(double*, signed char)
void projsimplex(NUMBER* y, NPAR size) {
    
	bool bget = false;
	NUMBER *s = init1D<NUMBER>((NDAT)size);
	cpy1D<NUMBER>(y, s, (NDAT)size);
	
	qsortNumberRev(s, size);
	NUMBER tmpsum = 0, tmax = 0;
    
	for(NPAR i=0; i<(size-1); i++) {
		tmpsum = tmpsum + s[i];
		tmax = (tmpsum - 1)/(i+1);
		if(tmax >= s[i+1]) {
			bget = true;
			break;
		}
	}
	if(!bget) tmax = (tmpsum + s[size-1] -1)/size;
	free(s);
    
	for(NPAR i=0; i<size; i++)
		y[i] = ((y[i]-tmax)<0)?0:(y[i]-tmax);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
pushq %rax
movl %esi, %ebp
movq %rdi, %rbx
movslq %esi, %r15
movl $0x8, %esi
movq %r15, %rdi
callq 0x4130
movq %rax, %r14
leaq (,%r15,8), %rdx
movq %rax, %rdi
movq %rbx, %rsi
callq 0x41c0
leaq -0xa8(%rip), %rcx       # 0x63c1
movl $0x8, %edx
movq %r14, %rdi
movq %r15, %rsi
callq 0x4470
decl %r15d
xorpd %xmm0, %xmm0
cmpb $0x2, %bpl
jl 0x64c8
xorl %ecx, %ecx
movsd 0x2d458(%rip), %xmm1    # 0x338e8
xorl %eax, %eax
movsbq %al, %rdx
addsd (%r14,%rdx,8), %xmm0
movapd %xmm0, %xmm3
addsd %xmm1, %xmm3
leal 0x1(%rcx), %edx
xorps %xmm2, %xmm2
cvtsi2sd %edx, %xmm2
divsd %xmm2, %xmm3
movslq %ecx, %rcx
ucomisd 0x8(%r14,%rcx,8), %xmm3
jae 0x64ec
incb %al
movsbl %al, %ecx
cmpl %ecx, %r15d
jg 0x6492
movslq %r15d, %rax
addsd (%r14,%rax,8), %xmm0
addsd 0x2d40f(%rip), %xmm0    # 0x338e8
movsbl %bpl, %eax
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
divsd %xmm1, %xmm0
movapd %xmm0, %xmm3
movsd %xmm3, (%rsp)
movq %r14, %rdi
callq 0x42a0
movsd (%rsp), %xmm2
testb %bpl, %bpl
jle 0x6527
movzbl %bpl, %eax
xorl %ecx, %ecx
movsd (%rbx,%rcx,8), %xmm0
subsd %xmm2, %xmm0
xorpd %xmm1, %xmm1
maxsd %xmm0, %xmm1
movsd %xmm1, (%rbx,%rcx,8)
incq %rcx
cmpq %rcx, %rax
jne 0x6509
addq $0x8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: issimplex(double*, signed char)
bool issimplex(NUMBER* ar, NPAR size) {
	NUMBER sum = 1;
	for(NPAR i=0; i<size; i++) {
		sum -= ar[i];
		if( ar[i] < 0 || ar[i] > 1)
			return false;
	}
//    for(NPAR i=0; i<size && fabs(sum)<SAFETY && fabs(sum)>0; i++) {
//        if( (ar[i]+sum) >=0 ) {
//            ar[i]+=sum;
//            break;
//        }
//        if( (ar[i]+sum) <=1 ) {
//            ar[i]+=sum;
//            break;
//        }
//    }
	return fabs(sum)<SAFETY;
}

// Assembly:
testb %sil, %sil
jle 0x6597
movsd (%rdi), %xmm2
movsd 0x2d3ad(%rip), %xmm1    # 0x338f0
movapd %xmm1, %xmm0
subsd %xmm2, %xmm0
xorl %ecx, %ecx
xorpd %xmm3, %xmm3
ucomisd %xmm2, %xmm3
ja 0x65a1
ucomisd 0x2d391(%rip), %xmm2    # 0x338f0
ja 0x65a1
movzbl %sil, %eax
movl $0x1, %edx
xorpd %xmm2, %xmm2
movq %rdx, %rcx
cmpq %rdx, %rax
je 0x658f
movsd (%rdi,%rcx,8), %xmm3
subsd %xmm3, %xmm0
ucomisd %xmm3, %xmm2
ja 0x658f
leaq 0x1(%rcx), %rdx
ucomisd %xmm1, %xmm3
jbe 0x656e
cmpq %rax, %rcx
setae %cl
jmp 0x65a1
movsd 0x2d351(%rip), %xmm0    # 0x338f0
movb $0x1, %cl
andpd 0x2d387(%rip), %xmm0    # 0x33930
movsd 0x2d347(%rip), %xmm1    # 0x338f8
ucomisd %xmm0, %xmm1
seta %al
andb %cl, %al
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: issimplexbounded(double*, double*, double*, signed char)
bool issimplexbounded(NUMBER* ar, NUMBER *lb, NUMBER *ub, NPAR size) {
	NUMBER sum = 1;
	for(NPAR i=0; i<size; i++) {
		sum -= ar[i];
		if( ar[i] < lb[i] || ar[i] > ub[i])
			return false;
	}
//    for(NPAR i=0; i<size && fabs(sum)<SAFETY && fabs(sum)>0; i++) {
//        if( (ar[i]+sum) >=lb[i] ) {
//            ar[i]+=sum;
//            break;
//        }
//        if( (ar[i]+sum) <=ub[i] ) {
//            ar[i]+=sum;
//            break;
//        }
//    }
	return fabs(sum)<SAFETY;
}

// Assembly:
testb %cl, %cl
jle 0x65dd
movsd (%rdi), %xmm1
movsd 0x2d325(%rip), %xmm0    # 0x338f0
subsd %xmm1, %xmm0
movsd (%rsi), %xmm2
ucomisd %xmm1, %xmm2
jbe 0x65e9
xorl %ecx, %ecx
jmp 0x661e
movsd 0x2d30b(%rip), %xmm0    # 0x338f0
movb $0x1, %cl
jmp 0x661e
movzbl %cl, %eax
movl $0x1, %r8d
xorl %ecx, %ecx
ucomisd -0x8(%rdx,%r8,8), %xmm1
ja 0x661e
cmpq %rax, %r8
setae %cl
je 0x661e
movsd (%rdi,%r8,8), %xmm1
subsd %xmm1, %xmm0
movsd (%rsi,%r8,8), %xmm2
incq %r8
ucomisd %xmm1, %xmm2
jbe 0x65f4
andpd 0x2d30a(%rip), %xmm0    # 0x33930
movsd 0x2d2ca(%rip), %xmm1    # 0x338f8
ucomisd %xmm0, %xmm1
seta %al
andb %cl, %al
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: maxn(double*, int)
NUMBER maxn(NUMBER *ar, NDAT n) {
	NUMBER mx = ar[0];
	for(NDAT i=1; i<n; i++)
		if(ar[i]>mx) mx = ar[i];
	return mx;
}

// Assembly:
movsd (%rdi), %xmm1
cmpl $0x2, %esi
jl 0x8fd8
movl %esi, %eax
movl $0x1, %ecx
movsd (%rdi,%rcx,8), %xmm0
maxsd %xmm1, %xmm0
incq %rcx
movapd %xmm0, %xmm1
cmpq %rcx, %rax
jne 0x8fc2
retq
movapd %xmm1, %xmm0
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: add3DNumbersWeighted(double***, double***, signed char, signed char, signed char, double)
void add3DNumbersWeighted(NUMBER*** sourse, NUMBER*** target, NPAR size1, NPAR size2, NPAR size3, NUMBER weight) {
    for(NPAR i=0; i<size1; i++)
        for(NPAR j=0; j<size2; j++)
            for(NPAR l=0; l<size3; l++)
                target[i][j][l] = target[i][j][l] + sourse[i][j][l]*weight;
}

// Assembly:
testb %dl, %dl
jle 0x90cd
pushq %r15
pushq %r14
pushq %rbx
movzbl %dl, %eax
movzbl %cl, %edx
movzbl %r8b, %r9d
xorl %r10d, %r10d
testb %cl, %cl
jle 0x90c0
xorl %r11d, %r11d
testb %r8b, %r8b
jle 0x90b8
movq (%rsi,%r10,8), %rbx
movq (%rbx,%r11,8), %rbx
movq (%rdi,%r10,8), %r14
movq (%r14,%r11,8), %r14
xorl %r15d, %r15d
movsd (%r14,%r15,8), %xmm1
mulsd %xmm0, %xmm1
addsd (%rbx,%r15,8), %xmm1
movsd %xmm1, (%rbx,%r15,8)
incq %r15
cmpq %r15, %r9
jne 0x909a
incq %r11
cmpq %rdx, %r11
jne 0x9082
incq %r10
cmpq %rax, %r10
jne 0x907b
popq %rbx
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: isPasses(double*, signed char)
bool isPasses(NUMBER* ar, NPAR size) {
	NUMBER sum = 0;
	for(NPAR i=0; i<size; i++) {
		if( ar[i]<0 || ar[i]>1)
			return false;
		sum += ar[i];
	}
	return sum==1;
}

// Assembly:
testb %sil, %sil
jle 0x9144
movsd (%rdi), %xmm0
xorl %eax, %eax
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
ja 0x9143
ucomisd 0x2a805(%rip), %xmm0    # 0x338f0
ja 0x9143
movzbl %sil, %eax
movl $0x1, %edx
movsd 0x2a7f2(%rip), %xmm2    # 0x338f0
xorpd %xmm3, %xmm3
movapd %xmm3, %xmm4
movq %rdx, %rcx
movapd %xmm0, %xmm3
addsd %xmm4, %xmm3
cmpq %rdx, %rax
je 0x912b
movsd (%rdi,%rcx,8), %xmm0
ucomisd %xmm0, %xmm1
ja 0x912b
leaq 0x1(%rcx), %rdx
ucomisd %xmm2, %xmm0
jbe 0x9102
cmpq %rax, %rcx
setae %cl
ucomisd 0x2a7b7(%rip), %xmm3    # 0x338f0
setnp %dl
sete %al
andb %dl, %al
andb %cl, %al
retq
xorl %eax, %eax
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: isPassesLim(double*, signed char, double*, double*)
bool isPassesLim(NUMBER* ar, NPAR size, NUMBER *lb, NUMBER* ub) {
	NUMBER sum = 0;
	for(NPAR i=0; i<size; i++) {
		if( (lb[i]-ar[i])>SAFETY || (ar[i]-ub[i])>SAFETY )
			return false;
		sum += ar[i];
	}
	return fabs(sum-1)<SAFETY;
}

// Assembly:
testb %sil, %sil
jle 0x919a
movzbl %sil, %eax
xorpd %xmm0, %xmm0
xorl %esi, %esi
movsd 0x2a79a(%rip), %xmm1    # 0x338f8
movsd (%rdx,%rsi,8), %xmm3
movsd (%rdi,%rsi,8), %xmm2
subsd %xmm2, %xmm3
ucomisd %xmm1, %xmm3
ja 0x9197
movapd %xmm2, %xmm3
subsd (%rcx,%rsi,8), %xmm3
ucomisd %xmm1, %xmm3
ja 0x9197
addsd %xmm2, %xmm0
incq %rsi
cmpq %rsi, %rax
jne 0x915e
addsd 0x2a753(%rip), %xmm0    # 0x338e8
jmp 0x91a2
xorl %eax, %eax
retq
movsd 0x2a746(%rip), %xmm0    # 0x338e8
andpd 0x2a786(%rip), %xmm0    # 0x33930
movsd 0x2a746(%rip), %xmm1    # 0x338f8
ucomisd %xmm0, %xmm1
seta %al
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: doLog10Scale1D(double*, signed char)
NUMBER doLog10Scale1D(NUMBER *ar, NPAR size) {
	NPAR i;
	NUMBER min_10_scale = 1000, max_10_scale = 0, candidate;
	for(i=0; i<size; i++) {
		if( fabs(ar[i]) < SAFETY ) // 0 gradient
			continue;
		candidate = floor( log10( fabs(ar[i]) ) );
		if(candidate < min_10_scale)
			min_10_scale = candidate;
		candidate = ceil( log10( fabs(ar[i]) ) );
		if(candidate > max_10_scale)
			max_10_scale = candidate;
	}
	min_10_scale++;
	max_10_scale++;
	if(max_10_scale > 0)
		for(i=0; i<size; i++)
			ar[i] = ar[i] / pow(10, max_10_scale);
    //	if(min_10_scale<1000)
    //		for(i=0; i<size; i++)
    //			ar[i] = ar[i] / pow(10, min_10_scale);
	return pow(10, max_10_scale);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x18, %rsp
movl %esi, %ebp
testb %bpl, %bpl
jle 0x92ad
movq %rdi, %rbx
movzbl %bpl, %r14d
xorpd %xmm1, %xmm1
xorl %r15d, %r15d
movapd 0x2a74b(%rip), %xmm2    # 0x33930
movsd 0x2a70b(%rip), %xmm3    # 0x338f8
movsd (%rbx,%r15,8), %xmm0
andpd %xmm2, %xmm0
ucomisd %xmm0, %xmm3
ja 0x9240
movsd %xmm1, 0x8(%rsp)
callq 0x4110
movsd (%rbx,%r15,8), %xmm0
andpd 0x2a71a(%rip), %xmm0    # 0x33930
callq 0x4110
callq 0x4090
movsd 0x2a6d0(%rip), %xmm3    # 0x338f8
movapd 0x2a700(%rip), %xmm2    # 0x33930
movsd 0x8(%rsp), %xmm1
ucomisd %xmm1, %xmm0
jbe 0x9240
movapd %xmm0, %xmm1
incq %r15
cmpq %r15, %r14
jne 0x91ed
testb %bpl, %bpl
addsd 0x2a69d(%rip), %xmm1    # 0x338f0
setle %al
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
setbe %cl
orb %al, %cl
jne 0x92b5
xorl %r15d, %r15d
movsd %xmm1, 0x8(%rsp)
movsd (%rbx,%r15,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x2a68e(%rip), %xmm0    # 0x33910
movsd 0x8(%rsp), %xmm1
callq 0x4420
movsd 0x8(%rsp), %xmm1
movsd 0x10(%rsp), %xmm2
divsd %xmm0, %xmm2
movsd %xmm2, (%rbx,%r15,8)
incq %r15
cmpq %r15, %r14
jne 0x926e
jmp 0x92b5
movsd 0x2a63b(%rip), %xmm1    # 0x338f0
movsd 0x2a653(%rip), %xmm0    # 0x33910
addq $0x18, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
jmp 0x4420
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: doLog10Scale2D(double**, signed char, signed char)
NUMBER doLog10Scale2D(NUMBER **ar, NPAR size1, NPAR size2) {
	NPAR i,j;
	NUMBER min_10_scale = 1000, max_10_scale = 0, candidate;
	for(i=0; i<size1; i++)
		for(j=0; j<size2; j++) {
			if( fabs(ar[i][j]) < SAFETY ) // 0 gradient
				continue;
			candidate = floor( log10( fabs(ar[i][j]) ) );
			if(candidate < min_10_scale)
				min_10_scale = candidate;
			candidate = ceil( log10( fabs(ar[i][j]) ) );
			if(candidate > max_10_scale)
				max_10_scale = candidate;
		}
	min_10_scale++;
	max_10_scale++;
	if(max_10_scale >0 )
		for(i=0; i<size1; i++)
			for(j=0; j<size2; j++)
				ar[i][j] = ar[i][j] / pow(10, max_10_scale);
    //	if(min_10_scale<1000)
    //		for(i=0; i<size1; i++)
    //			for(j=0; j<size2; j++)
    //				ar[i][j] = ar[i][j] / pow(10, min_10_scale);
	return pow(10, max_10_scale);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
testb %sil, %sil
jle 0x93f5
movl %edx, %ebx
movq %rdi, %r14
movl %esi, 0x10(%rsp)
movzbl %sil, %r15d
movzbl %bl, %r12d
xorpd %xmm1, %xmm1
xorl %r13d, %r13d
movapd 0x2a62d(%rip), %xmm2    # 0x33930
movsd 0x2a5ed(%rip), %xmm3    # 0x338f8
testb %bl, %bl
jle 0x9372
xorl %ebp, %ebp
movq (%r14,%r13,8), %rax
movsd (%rax,%rbp,8), %xmm0
andpd %xmm2, %xmm0
ucomisd %xmm0, %xmm3
ja 0x936a
movsd %xmm1, 0x8(%rsp)
callq 0x4110
movq (%r14,%r13,8), %rax
movsd (%rax,%rbp,8), %xmm0
andpd 0x2a5f0(%rip), %xmm0    # 0x33930
callq 0x4110
callq 0x4090
movsd 0x2a5a6(%rip), %xmm3    # 0x338f8
movapd 0x2a5d6(%rip), %xmm2    # 0x33930
movsd 0x8(%rsp), %xmm1
ucomisd %xmm1, %xmm0
jbe 0x936a
movapd %xmm0, %xmm1
incq %rbp
cmpq %rbp, %r12
jne 0x9311
incq %r13
cmpq %r15, %r13
jne 0x930b
cmpb $0x0, 0x10(%rsp)
addsd 0x2a569(%rip), %xmm1    # 0x338f0
setle %al
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
setbe %cl
orb %al, %cl
jne 0x93fd
xorl %r13d, %r13d
movsd %xmm1, 0x8(%rsp)
testb %bl, %bl
jle 0x93eb
movq (%r14,%r13,8), %rax
xorl %ebp, %ebp
movsd (%rax,%rbp,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x2a551(%rip), %xmm0    # 0x33910
movsd 0x8(%rsp), %xmm1
callq 0x4420
movsd 0x8(%rsp), %xmm1
movsd 0x10(%rsp), %xmm2
divsd %xmm0, %xmm2
movq (%r14,%r13,8), %rax
movsd %xmm2, (%rax,%rbp,8)
incq %rbp
cmpq %rbp, %r12
jne 0x93ac
incq %r13
cmpq %r15, %r13
jne 0x93a2
jmp 0x93fd
movsd 0x2a4f3(%rip), %xmm1    # 0x338f0
movsd 0x2a50b(%rip), %xmm0    # 0x33910
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x4420
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: writeSolverInfo(_IO_FILE*, param*)
void writeSolverInfo(FILE *fid, struct param *param) {
    // solver id
    if( param->solver_setting>0 ) {
        fprintf(fid,"SolverId\t%d.%d.%d\n",param->structure,param->solver,param->solver_setting);
    } else {
        fprintf(fid,"SolverId\t%d.%d\n",param->structure,param->solver);
    }
	// nK
    fprintf(fid,"nK\t%d\n",param->nK);
    // nG
    fprintf(fid,"nG\t%d\n",param->nG);
    // nS
    fprintf(fid,"nS\t%d\n",param->nS);
    // nO
    fprintf(fid,"nO\t%d\n",param->nO);
    // nZ
    fprintf(fid,"nZ\t%d\n",param->nZ);
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rsi, %r14
movq %rdi, %rbx
movsbl 0x40(%rsi), %r8d
testl %r8d, %r8d
jle 0xa156
movsbl 0x3e(%r14), %edx
movsbl 0x3f(%r14), %ecx
leaq 0x299b3(%rip), %rsi     # 0x33afd
movq %rbx, %rdi
xorl %eax, %eax
callq 0x4330
jmp 0xa171
movsbl 0x3e(%r14), %edx
movsbl 0x3f(%r14), %ecx
leaq 0x299a9(%rip), %rsi     # 0x33b10
movq %rbx, %rdi
xorl %eax, %eax
callq 0x4330
movl 0x8c0(%r14), %edx
leaq 0x299a1(%rip), %rsi     # 0x33b20
movq %rbx, %rdi
xorl %eax, %eax
callq 0x4330
movl 0x8bc(%r14), %edx
leaq 0x29990(%rip), %rsi     # 0x33b27
movq %rbx, %rdi
xorl %eax, %eax
callq 0x4330
movsbl 0x8b8(%r14), %edx
leaq 0x2997e(%rip), %rsi     # 0x33b2e
movq %rbx, %rdi
xorl %eax, %eax
callq 0x4330
movsbl 0x8b9(%r14), %edx
leaq 0x2996c(%rip), %rsi     # 0x33b35
movq %rbx, %rdi
xorl %eax, %eax
callq 0x4330
movsbl 0x8c8(%r14), %edx
leaq 0x2995a(%rip), %rsi     # 0x33b3c
movq %rbx, %rdi
xorl %eax, %eax
addq $0x8, %rsp
popq %rbx
popq %r14
jmp 0x4330
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: readSolverInfo(_IO_FILE*, param*, int*)
void readSolverInfo(FILE *fid, struct param *param, NDAT *line_no) {
    string s;
    int c, i1, i2;
    // SolverId
    fscanf(fid,"SolverId\t%i.%i", &i1, &i2);
    param->structure = (NPAR) i1;
    param->solver    = (NPAR) i2;
    fscanf(fid,"SolverId\t%hhu.%hhu", &param->structure, &param->solver);
    c = fscanf(fid,".%hhu\n", &param->solver_setting);
    if( c<1 ) {
        fscanf(fid,"\n");
        param->solver_setting = -1;
    }
    (*line_no)++;
	// nK
    fscanf(fid,"nK\t%i\n",&param->nK);
    (*line_no)++;
    // nG
    fscanf(fid,"nG\t%i\n",&param->nG);
    (*line_no)++;
    // nS
    fscanf(fid,"nS\t%hhu\n",&param->nS);
    (*line_no)++;
    // nO
    fscanf(fid,"nO\t%hhu\n",&param->nO);
    (*line_no)++;
    // nZ
    fscanf(fid,"nZ\t%hhu\n",&param->nZ);
    (*line_no)++;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x30, %rsp
movq %rdx, %rbx
movq %rsi, %r14
movq %rdi, %r15
leaq 0x20(%rsp), %r13
movq %r13, -0x10(%r13)
movq $0x0, -0x8(%r13)
movb $0x0, (%r13)
leaq 0x2991d(%rip), %rsi     # 0x33b43
leaq 0xc(%rsp), %rdx
leaq 0x8(%rsp), %rcx
xorl %eax, %eax
callq 0x4460
movb 0xc(%rsp), %al
leaq 0x3e(%r14), %rdx
movb %al, 0x3e(%r14)
movb 0x8(%rsp), %al
leaq 0x3f(%r14), %rcx
movb %al, 0x3f(%r14)
leaq 0x298fc(%rip), %rsi     # 0x33b52
movq %r15, %rdi
xorl %eax, %eax
callq 0x4460
leaq 0x40(%r14), %r12
leaq 0x298fa(%rip), %rsi     # 0x33b65
movq %r15, %rdi
movq %r12, %rdx
xorl %eax, %eax
callq 0x4460
testl %eax, %eax
jg 0xa292
leaq 0x28fda(%rip), %rsi     # 0x3325d
movq %r15, %rdi
xorl %eax, %eax
callq 0x4460
movb $-0x1, (%r12)
incl (%rbx)
leaq 0x8c0(%r14), %rdx
leaq 0x298ca(%rip), %rsi     # 0x33b6c
movq %r15, %rdi
xorl %eax, %eax
callq 0x4460
incl (%rbx)
leaq 0x8bc(%r14), %rdx
leaq 0x298b7(%rip), %rsi     # 0x33b73
movq %r15, %rdi
xorl %eax, %eax
callq 0x4460
incl (%rbx)
leaq 0x8b8(%r14), %rdx
leaq 0x298a4(%rip), %rsi     # 0x33b7a
movq %r15, %rdi
xorl %eax, %eax
callq 0x4460
incl (%rbx)
leaq 0x8b9(%r14), %rdx
leaq 0x29893(%rip), %rsi     # 0x33b83
movq %r15, %rdi
xorl %eax, %eax
callq 0x4460
incl (%rbx)
addq $0x8c8, %r14            # imm = 0x8C8
leaq 0x29882(%rip), %rsi     # 0x33b8c
movq %r15, %rdi
movq %r14, %rdx
xorl %eax, %eax
callq 0x4460
incl (%rbx)
movq 0x10(%rsp), %rdi
cmpq %r13, %rdi
je 0xa330
movq 0x20(%rsp), %rsi
incq %rsi
callq 0x4220
addq $0x30, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
retq
movq %rax, %rbx
movq 0x10(%rsp), %rdi
cmpq %r13, %rdi
je 0xa358
movq 0x20(%rsp), %rsi
incq %rsi
callq 0x4220
movq %rbx, %rdi
callq 0x4430
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: zeroLabels(param*)
void zeroLabels(struct param* param) { // set counts in data sequences to zero
	NCAT k,g;
	for(k=0; k<param->nK; k++)
		for(g=0; g<param->k_numg[k]; g++)
			param->k_g_data[k][g]->cnt = 0;
}

// Assembly:
movslq 0x8c0(%rdi), %rax
testq %rax, %rax
jle 0xa3a9
movq 0x8e0(%rdi), %rcx
xorl %edx, %edx
cmpl $0x0, (%rcx,%rdx,4)
jle 0xa3a1
movq 0x8f8(%rdi), %rsi
movq (%rsi,%rdx,8), %rsi
xorl %r8d, %r8d
movq (%rsi,%r8,8), %r9
movl $0x0, 0x4(%r9)
incq %r8
movslq (%rcx,%rdx,4), %r9
cmpq %r9, %r8
jl 0xa389
incq %rdx
cmpq %rax, %rdx
jne 0xa375
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: RecycleFitData(int, data**, param*)
void RecycleFitData(NCAT xndat, struct data** x_data, struct param *param) {
	NCAT x;
	NDAT t;
	for(x=0; x<xndat; x++) {
        //        if( x_data[x][0].cnt != 0)
        //            continue;
		if( x_data[x][0].alpha != NULL ) {
			free2D<NUMBER>(x_data[x][0].alpha, x_data[x][0].n);  // only free data here
			x_data[x][0].alpha = NULL;
		}
		if( x_data[x][0].c != NULL ) {
			free(x_data[x][0].c);  // only free data here
			x_data[x][0].c = NULL;
		}
		if( x_data[x][0].beta != NULL ) {
			free2D<NUMBER>(x_data[x][0].beta,  x_data[x][0].n); // only free data here
			x_data[x][0].beta = NULL;
		}
		if( x_data[x][0].gamma != NULL ) {
			free2D<NUMBER>(x_data[x][0].gamma, x_data[x][0].n); // only free data here
			x_data[x][0].gamma = NULL;
		}
		if( x_data[x][0].xi != NULL ) {
			for(t=0;t<x_data[x][0].n; t++)
				free2D<NUMBER>(x_data[x][0].xi[t],  (NDAT)param->nS); // only free data here
			x_data[x][0].xi = NULL;
		}
	}
}

// Assembly:
testl %edi, %edi
jle 0xa498
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
movq %rdx, %rbx
movq %rsi, %r14
movl %edi, %r15d
xorl %r12d, %r12d
movq (%r14,%r12,8), %rax
movq 0x28(%rax), %rdi
testq %rdi, %rdi
je 0xa3eb
movl (%rax), %esi
callq 0xa5dd
movq (%r14,%r12,8), %rax
movq $0x0, 0x28(%rax)
movq (%r14,%r12,8), %rax
movq 0x20(%rax), %rdi
testq %rdi, %rdi
je 0xa409
callq 0x42a0
movq (%r14,%r12,8), %rax
movq $0x0, 0x20(%rax)
movq (%r14,%r12,8), %rax
movq 0x30(%rax), %rdi
testq %rdi, %rdi
je 0xa429
movl (%rax), %esi
callq 0xa5dd
movq (%r14,%r12,8), %rax
movq $0x0, 0x30(%rax)
movq (%r14,%r12,8), %rax
movq 0x38(%rax), %rdi
testq %rdi, %rdi
je 0xa449
movl (%rax), %esi
callq 0xa5dd
movq (%r14,%r12,8), %rax
movq $0x0, 0x38(%rax)
movq (%r14,%r12,8), %rax
cmpq $0x0, 0x40(%rax)
je 0xa483
cmpl $0x0, (%rax)
jle 0xa47b
xorl %r13d, %r13d
movq 0x40(%rax), %rax
movq (%rax,%r13,8), %rdi
movsbl 0x8b8(%rbx), %esi
callq 0xa5dd
incq %r13
movq (%r14,%r12,8), %rax
movslq (%rax), %rcx
cmpq %rcx, %r13
jl 0xa458
movq $0x0, 0x40(%rax)
incq %r12
cmpq %r15, %r12
jne 0xa3c7
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/utils.cpp
// Function: sec_to_9cat(int, int, int*, signed char)
NPAR sec_to_9cat(int time1, int time2, int *limits, NPAR nlimits) {
    int diff = time2 - time1;
    
    if(diff <0 ) {
        fprintf(stderr,"ERROR! time 1 should be smaller than time 2\n");
        return 0;
    } else if ( diff < (2*60) ) { // 20min
        return 0;
    } else if ( diff < (20*60) ) { // 20min
        return 1;
    } else if (diff < (60*60) ) { // 1h
        return 2;
    } else if(diff < (14*24*60*60)) { // detect date structures
        time_t t1 = (time_t)time1;
        time_t t2 = (time_t)time2;
        struct tm * ttm1 = localtime (&t1);
        struct tm * ttm2 = localtime (&t2);
        if( diff < (24*60*60) && ttm1->tm_mday==ttm2->tm_mday ) { // same day
            return 3;
        } else if(    diff < (2*24*60*60) &&                      // next day
                  ( (ttm2->tm_wday == (ttm1->tm_wday+1)) ||
                   (ttm1->tm_wday==6 && ttm2->tm_wday==0)
                   )
                  ) {
            return 4;
        } else if( diff < (7*24*60*60) && (ttm1->tm_wday<ttm2->tm_wday) ) { // this week
            return 5;
        } else { // next week
            return 6;
        }
    } else if(diff < (30*24*60*60)) { // less than 30 days
        return 7;
    } else if(diff >= (30*24*60*60)) { // 30 days or more
        return 8;
    }
    return 0;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x10, %rsp
movl %esi, %ebp
subl %edi, %ebp
js 0xa53b
cmpl $0x78, %ebp
jae 0xa4d1
xorl %eax, %eax
jmp 0xa532
movb $0x1, %al
cmpl $0x4b0, %ebp            # imm = 0x4B0
jb 0xa532
movb $0x2, %al
cmpl $0xe10, %ebp            # imm = 0xE10
jb 0xa532
cmpl $0x1274ff, %ebp         # imm = 0x1274FF
ja 0xa528
movslq %edi, %rax
leaq 0x8(%rsp), %rdi
movq %rax, (%rdi)
movslq %esi, %rax
movq %rsp, %r14
movq %rax, (%r14)
callq 0x4190
movq %rax, %rbx
movq %r14, %rdi
callq 0x4190
movq %rax, %rcx
cmpl $0x1517f, %ebp          # imm = 0x1517F
ja 0xa542
movl 0xc(%rbx), %edx
movb $0x3, %al
cmpl 0xc(%rcx), %edx
je 0xa532
jmp 0xa54a
cmpl $0x278d00, %ebp         # imm = 0x278D00
movb $0x8, %al
sbbb $0x0, %al
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
retq
callq 0x4803
jmp 0xa4cd
cmpl $0x2a2ff, %ebp          # imm = 0x2A2FF
ja 0xa564
movl 0x18(%rcx), %ecx
movl 0x18(%rbx), %edx
leal 0x1(%rdx), %esi
movb $0x4, %al
cmpl %esi, %ecx
je 0xa532
movl %edx, %esi
xorl $0x6, %esi
orl %ecx, %esi
je 0xa532
jmp 0xa572
cmpl $0x93a7f, %ebp          # imm = 0x93A7F
ja 0xa578
movl 0x18(%rbx), %edx
movl 0x18(%rcx), %ecx
movb $0x5, %al
cmpl %ecx, %edx
jl 0xa532
movb $0x6, %al
jmp 0xa532
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/./utils.h
// Function: void free2D<double>(double**, int)
void free2D(T** ar, NDAT size1) {
	for(NDAT i=0; i<size1; i++)
		free(ar[i]);
	free(ar);
    //    &ar = NULL;
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
movq %rdi, %rbx
testl %esi, %esi
jle 0xa600
movl %esi, %r14d
xorl %r15d, %r15d
movq (%rbx,%r15,8), %rdi
callq 0x42a0
incq %r15
cmpq %r15, %r14
jne 0xa5ef
movq %rbx, %rdi
popq %rbx
popq %r14
popq %r15
jmp 0x42a0
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBit.cpp
// Function: FitBit::init(double*&, double**&, double**&)
void FitBit::init(NUMBER* &a_PI, NUMBER** &a_A, NUMBER** &a_B) {
    if(this->pi != NULL) {
        if(a_PI == NULL)
            a_PI = init1D<NUMBER>((NDAT)this->nS); // init1DNumber(this->nS);
        else
            toZero1D<NUMBER>(a_PI, (NDAT)this->nS);
    }
    if(this->A  != NULL) {
        if(a_A == NULL)
            a_A  = init2D<NUMBER>((NDAT)this->nS, (NDAT)this->nS);
        else
            toZero2D<NUMBER>(a_A,  (NDAT)this->nS, (NDAT)this->nS);
    }
    if(this->B  != NULL) {
        if(a_B == NULL)
            a_B  = init2D<NUMBER>((NDAT)this->nS, (NDAT)this->nO);
        else
            toZero2D<NUMBER>(a_B, (NDAT)this->nS, (NDAT)this->nO);
    }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rcx, %rbx
movq %rdx, %r15
movq %rdi, %r14
cmpq $0x0, 0x10(%rdi)
je 0xab05
movq %rsi, %r12
movq (%rsi), %rdi
movsbq 0x1(%r14), %rdx
testq %rdi, %rdi
je 0xaaf4
testb %dl, %dl
jle 0xab05
shlq $0x3, %rdx
xorl %esi, %esi
callq 0x4120
jmp 0xab05
movl $0x8, %esi
movq %rdx, %rdi
callq 0x4130
movq %rax, (%r12)
cmpq $0x0, 0x18(%r14)
je 0xab4d
movq (%r15), %r12
movsbl 0x1(%r14), %edi
testq %r12, %r12
je 0xab43
testb %dil, %dil
jle 0xab4d
movslq %edi, %r13
leaq (,%r13,8), %r15
xorl %ebp, %ebp
movq (%r12,%rbp,8), %rdi
xorl %esi, %esi
movq %r15, %rdx
callq 0x4120
incq %rbp
cmpq %rbp, %r13
jne 0xab2b
jmp 0xab4d
movl %edi, %esi
callq 0xbb36
movq %rax, (%r15)
cmpq $0x0, 0x20(%r14)
je 0xaba0
movq (%rbx), %r15
movsbl 0x1(%r14), %edi
testq %r15, %r15
je 0xab94
testb %dil, %dil
jle 0xaba0
movslq %edi, %r12
movsbq (%r14), %r14
movl %r14d, %ebx
shlq $0x3, %rbx
xorl %r13d, %r13d
testb %r14b, %r14b
jle 0xab8a
movq (%r15,%r13,8), %rdi
xorl %esi, %esi
movq %rbx, %rdx
callq 0x4120
incq %r13
cmpq %r13, %r12
jne 0xab77
jmp 0xaba0
movsbl (%r14), %esi
callq 0xbb36
movq %rax, (%rbx)
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBit.cpp
// Function: FitBit::negate(double*&, double**&, double**&)
void FitBit::negate(NUMBER* &a_PI, NUMBER** &a_A, NUMBER** &a_B) {
    if(this->pi != NULL) {
        for(NPAR i=0; i<this->nS; i++) a_PI[i] = -a_PI[i];
    }
    if(this->A  != NULL) {
        for(NPAR i=0; i<this->nS; i++)
            for(NPAR j=0; j<this->nS; j++)
                a_A[i][j] = -a_A[i][j];
    }
    if(this->B  != NULL) {
        for(NPAR i=0; i<this->nS; i++)
            for(NPAR m=0; m<this->nO; m++)
                a_B[i][m] = -a_B[i][m];
    }
}

// Assembly:
cmpq $0x0, 0x10(%rdi)
je 0xabe4
movsbq 0x1(%rdi), %rax
testq %rax, %rax
jle 0xabe4
movq (%rsi), %rsi
xorl %r8d, %r8d
movaps 0x28d72(%rip), %xmm0    # 0x33940
movsd (%rsi,%r8,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%rsi,%r8,8)
incq %r8
cmpq %r8, %rax
jne 0xabce
cmpq $0x0, 0x18(%rdi)
je 0xac26
movsbq 0x1(%rdi), %rax
testq %rax, %rax
jle 0xac26
movq (%rdx), %rdx
xorl %esi, %esi
movaps 0x28d3f(%rip), %xmm0    # 0x33940
movq (%rdx,%rsi,8), %r8
xorl %r9d, %r9d
movsd (%r8,%r9,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%r8,%r9,8)
incq %r9
cmpq %r9, %rax
jne 0xac08
incq %rsi
cmpq %rax, %rsi
jne 0xac01
cmpq $0x0, 0x20(%rdi)
je 0xac6f
movsbq 0x1(%rdi), %rax
testq %rax, %rax
jle 0xac6f
movzbl (%rdi), %edx
movq (%rcx), %rcx
xorl %esi, %esi
movaps 0x28cfa(%rip), %xmm0    # 0x33940
testb %dl, %dl
jle 0xac67
movq (%rcx,%rsi,8), %rdi
xorl %r8d, %r8d
movsd (%rdi,%r8,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%rdi,%r8,8)
incq %r8
cmpq %r8, %rdx
jne 0xac51
incq %rsi
cmpq %rax, %rsi
jne 0xac46
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBit.cpp
// Function: FitBit::destroy(double*&, double**&, double**&)
void FitBit::destroy(NUMBER* &a_PI, NUMBER** &a_A, NUMBER** &a_B) {
    if(this->pi != NULL && a_PI != NULL) free(a_PI);
    if(this->A  != NULL && a_A  != NULL) free2D<NUMBER>(a_A, (NDAT)this->nS);
    if(this->B  != NULL && a_B  != NULL) free2D<NUMBER>(a_B, (NDAT)this->nS);
    a_PI = NULL;
    a_A  = NULL;
    a_B  = NULL;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rcx, %rbx
movq %rdx, %r14
movq %rsi, %r15
movq %rdi, %r12
cmpq $0x0, 0x10(%rdi)
je 0xaec2
movq (%r15), %rdi
testq %rdi, %rdi
je 0xaec2
callq 0x42a0
cmpq $0x0, 0x18(%r12)
je 0xaedd
movq (%r14), %rdi
testq %rdi, %rdi
je 0xaedd
movsbl 0x1(%r12), %esi
callq 0xa5dd
cmpq $0x0, 0x20(%r12)
je 0xaef8
movq (%rbx), %rdi
testq %rdi, %rdi
je 0xaef8
movsbl 0x1(%r12), %esi
callq 0xa5dd
xorl %eax, %eax
movq %rax, (%r15)
movq %rax, (%r14)
movq %rax, (%rbx)
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBit.cpp
// Function: FitBit::copy(FIT_BIT_SLOT, FIT_BIT_SLOT)
void FitBit::copy(enum FIT_BIT_SLOT sourse_fbs, enum FIT_BIT_SLOT target_fbs) {
    NUMBER *soursePI = NULL;
    NUMBER **sourseA = NULL;
    NUMBER **sourseB = NULL;
    get(sourse_fbs, soursePI, sourseA, sourseB);
    NUMBER *targetPI = NULL;
    NUMBER **targetA = NULL;
    NUMBER **targetB = NULL;
    get(target_fbs, targetPI, targetA, targetB);
    
    copy(soursePI, sourseA, sourseB, targetPI, targetA, targetB);
}

// Assembly:
subq $0x38, %rsp
xorl %eax, %eax
movq %rax, 0x30(%rsp)
movq %rax, 0x28(%rsp)
movq %rax, 0x20(%rsp)
decl %esi
leaq 0x28c90(%rip), %rax     # 0x33ce8
cmpl $0x8, %esi
ja 0xb094
movl %esi, %ecx
movq (%rax,%rcx,8), %rsi
leaq 0x28cc6(%rip), %r8      # 0x33d30
movq (%r8,%rcx,8), %r8
leaq 0x28d03(%rip), %r9      # 0x33d78
movq (%r9,%rcx,8), %rcx
movq (%rdi,%rsi), %rsi
movq %rsi, 0x30(%rsp)
movq (%rdi,%r8), %rsi
movq %rsi, 0x28(%rsp)
movq (%rdi,%rcx), %rcx
movq %rcx, 0x20(%rsp)
xorl %ecx, %ecx
movq %rcx, 0x18(%rsp)
movq %rcx, 0x10(%rsp)
movq %rcx, 0x8(%rsp)
decl %edx
cmpl $0x8, %edx
ja 0xb0e3
movl %edx, %ecx
movq (%rax,%rcx,8), %rax
leaq 0x28c77(%rip), %rdx     # 0x33d30
movq (%rdx,%rcx,8), %rdx
leaq 0x28cb4(%rip), %rsi     # 0x33d78
movq (%rsi,%rcx,8), %rcx
movq (%rdi,%rax), %rax
movq %rax, 0x18(%rsp)
movq (%rdi,%rdx), %rax
movq %rax, 0x10(%rsp)
movq (%rdi,%rcx), %rax
movq %rax, 0x8(%rsp)
leaq 0x8(%rsp), %rax
movq %rax, (%rsp)
leaq 0x30(%rsp), %rsi
leaq 0x28(%rsp), %rdx
leaq 0x20(%rsp), %rcx
leaq 0x18(%rsp), %r8
leaq 0x10(%rsp), %r9
callq 0xad48
addq $0x38, %rsp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBit.cpp
// Function: FitBit::add(FIT_BIT_SLOT, FIT_BIT_SLOT)
void FitBit::add(enum FIT_BIT_SLOT sourse_fbs, enum FIT_BIT_SLOT target_fbs) {
    NUMBER *soursePI = NULL;
    NUMBER **sourseA = NULL;
    NUMBER **sourseB = NULL;
    get(sourse_fbs, soursePI, sourseA, sourseB);
    NUMBER *targetPI = NULL;
    NUMBER **targetA = NULL;
    NUMBER **targetB = NULL;
    get(target_fbs, targetPI, targetA, targetB);
    
    add(soursePI, sourseA, sourseB, targetPI, targetA, targetB);
}

// Assembly:
movl %edx, %eax
decl %esi
leaq 0x28bcd(%rip), %r8      # 0x33ce8
cmpl $0x8, %esi
ja 0xb14a
movl %esi, %ecx
movq (%r8,%rcx,8), %rdx
leaq 0x28c03(%rip), %rsi     # 0x33d30
movq (%rsi,%rcx,8), %r9
leaq 0x28c40(%rip), %rsi     # 0x33d78
movq (%rsi,%rcx,8), %rcx
movq (%rdi,%rdx), %rsi
movq (%rdi,%r9), %rdx
movq (%rdi,%rcx), %rcx
jmp 0xb150
xorl %esi, %esi
xorl %edx, %edx
xorl %ecx, %ecx
decl %eax
cmpl $0x8, %eax
ja 0xb181
movl %eax, %eax
movq (%r8,%rax,8), %r8
leaq 0x28bcc(%rip), %r9      # 0x33d30
movq (%r9,%rax,8), %r9
leaq 0x28c09(%rip), %r10     # 0x33d78
movq (%r10,%rax,8), %rax
movq (%rdi,%r8), %r8
movq (%rdi,%r9), %r9
movq (%rdi,%rax), %rax
jmp 0xb189
xorl %r8d, %r8d
xorl %r9d, %r9d
xorl %eax, %eax
pushq %rax
movq %rax, (%rsp)
callq 0xae0a
popq %rax
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/./utils.h
// Function: double** init2D<double>(int, int)
T** init2D(NDAT size1, NDAT size2) {
	T** ar = (T **)Calloc(T *, (size_t)size1);
	for(NDAT i=0; i<size1; i++)
		ar[i] = (T *)Calloc(T, (size_t)size2);
	return ar;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movl %esi, %r14d
movl %edi, %ebp
movslq %edi, %r15
movl $0x8, %esi
movq %r15, %rdi
callq 0x4130
movq %rax, %rbx
testl %r15d, %r15d
jle 0xbb7d
movslq %r14d, %r14
movl %ebp, %r15d
xorl %r12d, %r12d
movl $0x8, %esi
movq %r14, %rdi
callq 0x4130
movq %rax, (%rbx,%r12,8)
incq %r12
cmpq %r12, %r15
jne 0xbb64
movq %rbx, %rax
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
nopl (%rax)
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::init(double*&, double***&, double**&)
void FitBitSlicedA::init(NUMBER* &a_PI, NUMBER*** &a_A, NUMBER** &a_B) {
    if(this->pi != NULL) {
        if(a_PI == NULL)
            a_PI = init1D<NUMBER>((NDAT)this->nS); // init1DNumber(this->nS);
        else
            toZero1D<NUMBER>(a_PI, (NDAT)this->nS);
    }
    if(this->A  != NULL) {
        if(a_A == NULL)
            a_A  = init3D<NUMBER>((NDAT)this->nZ, (NDAT)this->nS, (NDAT)this->nS);
        else
            toZero3D<NUMBER>(a_A, (NDAT)this->nZ, (NDAT)this->nS, (NDAT)this->nS);
    }
    if(this->B  != NULL) {
        if(a_B == NULL)
            a_B  = init2D<NUMBER>( (NDAT)this->nS, (NDAT)this->nO);
        else
            toZero2D<NUMBER>(a_B, (NDAT)this->nS, (NDAT)this->nO);
    }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %rcx, %rbx
movq %rdx, %r15
movq %rdi, %r14
cmpq $0x0, 0x10(%rdi)
je 0xbe28
movq %rsi, %r12
movq (%rsi), %rdi
movsbq 0x1(%r14), %rdx
testq %rdi, %rdi
je 0xbe17
testb %dl, %dl
jle 0xbe28
shlq $0x3, %rdx
xorl %esi, %esi
callq 0x4120
jmp 0xbe28
movl $0x8, %esi
movq %rdx, %rdi
callq 0x4130
movq %rax, (%r12)
movq %rbx, 0x10(%rsp)
movq %r14, 0x8(%rsp)
cmpq $0x0, 0x18(%r14)
je 0xbeb4
movq (%r15), %rcx
movq 0x8(%rsp), %rax
movsbl 0x2(%rax), %edi
movq %rcx, 0x18(%rsp)
testq %rcx, %rcx
je 0xbea6
testb %dil, %dil
jle 0xbeb4
movslq %edi, %rax
movq %rax, 0x20(%rsp)
movq 0x8(%rsp), %rax
movsbq 0x1(%rax), %rbp
movl %ebp, %ebx
leaq (,%rbx,8), %r15
xorl %r14d, %r14d
testb %bpl, %bpl
jle 0xbe9a
movq 0x18(%rsp), %rax
movq (%rax,%r14,8), %r12
xorl %r13d, %r13d
movq (%r12,%r13,8), %rdi
xorl %esi, %esi
movq %r15, %rdx
callq 0x4120
incq %r13
cmpq %r13, %rbx
jne 0xbe84
incq %r14
cmpq 0x20(%rsp), %r14
jne 0xbe73
jmp 0xbeb4
movsbl 0x1(%rax), %edx
movl %edx, %esi
callq 0xcaae
movq %rax, (%r15)
movq 0x8(%rsp), %rax
cmpq $0x0, 0x20(%rax)
je 0xbf0f
movq 0x10(%rsp), %rbx
movq (%rbx), %r15
movsbl 0x1(%rax), %edi
testq %r15, %r15
je 0xbf04
testb %dil, %dil
jle 0xbf0f
movslq %edi, %r14
movsbq (%rax), %r12
movl %r12d, %ebx
shlq $0x3, %rbx
xorl %r13d, %r13d
testb %r12b, %r12b
jle 0xbefa
movq (%r15,%r13,8), %rdi
xorl %esi, %esi
movq %rbx, %rdx
callq 0x4120
incq %r13
cmpq %r13, %r14
jne 0xbee7
jmp 0xbf0f
movsbl (%rax), %esi
callq 0xbb36
movq %rax, (%rbx)
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::negate(double*&, double***&, double**&)
void FitBitSlicedA::negate(NUMBER* &a_PI, NUMBER*** &a_A, NUMBER** &a_B) {
    if(this->pi != NULL) {
        for(NPAR i=0; i<this->nS; i++) a_PI[i] = -a_PI[i];
    }
    if(this->A  != NULL) {
        for(NPAR z=0; z<this->nZ; z++)
            for(NPAR i=0; i<this->nS; i++)
                for(NPAR j=0; j<this->nS; j++)
                    a_A[z][i][j] = -a_A[z][i][j];
    }
    if(this->B  != NULL) {
        for(NPAR i=0; i<this->nS; i++)
            for(NPAR m=0; m<this->nO; m++)
                a_B[i][m] = -a_B[i][m];
    }
}

// Assembly:
pushq %rbx
cmpq $0x0, 0x10(%rdi)
je 0xbf53
movsbq 0x1(%rdi), %rax
testq %rax, %rax
jle 0xbf53
movq (%rsi), %rsi
xorl %r8d, %r8d
movaps 0x27a03(%rip), %xmm0    # 0x33940
movsd (%rsi,%r8,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%rsi,%r8,8)
incq %r8
cmpq %r8, %rax
jne 0xbf3d
cmpq $0x0, 0x18(%rdi)
je 0xbfad
movsbq 0x2(%rdi), %rax
testq %rax, %rax
jle 0xbfad
movzbl 0x1(%rdi), %esi
xorl %r8d, %r8d
movaps 0x279ce(%rip), %xmm0    # 0x33940
testb %sil, %sil
jle 0xbfa5
movq (%rdx), %r9
movq (%r9,%r8,8), %r9
xorl %r10d, %r10d
movq (%r9,%r10,8), %r11
xorl %ebx, %ebx
movsd (%r11,%rbx,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%r11,%rbx,8)
incq %rbx
cmpq %rbx, %rsi
jne 0xbf87
incq %r10
cmpq %rsi, %r10
jne 0xbf81
incq %r8
cmpq %rax, %r8
jne 0xbf72
cmpq $0x0, 0x20(%rdi)
je 0xbff6
movsbq 0x1(%rdi), %rax
testq %rax, %rax
jle 0xbff6
movzbl (%rdi), %edx
movq (%rcx), %rcx
xorl %esi, %esi
movaps 0x27973(%rip), %xmm0    # 0x33940
testb %dl, %dl
jle 0xbfee
movq (%rcx,%rsi,8), %rdi
xorl %r8d, %r8d
movsd (%rdi,%r8,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%rdi,%r8,8)
incq %r8
cmpq %r8, %rdx
jne 0xbfd8
incq %rsi
cmpq %rax, %rsi
jne 0xbfcd
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::destroy(double*&, double***&, double**&)
void FitBitSlicedA::destroy(NUMBER* &a_PI, NUMBER*** &a_A, NUMBER** &a_B) {
    if(this->pi != NULL && a_PI != NULL) free(a_PI);
    if(this->A  != NULL && a_A  != NULL) free3D<NUMBER>(a_A, (NDAT)this->nZ, (NDAT)this->nS);
    if(this->B  != NULL && a_B  != NULL) free2D<NUMBER>(a_B, (NDAT)this->nS);
    a_PI = NULL;
    a_A  = NULL;
    a_B  = NULL;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rcx, %rbx
movq %rdx, %r14
movq %rsi, %r15
movq %rdi, %r12
cmpq $0x0, 0x10(%rdi)
je 0xc2e6
movq (%r15), %rdi
testq %rdi, %rdi
je 0xc2e6
callq 0x42a0
cmpq $0x0, 0x18(%r12)
je 0xc307
movq (%r14), %rdi
testq %rdi, %rdi
je 0xc307
movsbl 0x2(%r12), %esi
movsbl 0x1(%r12), %edx
callq 0xca4f
cmpq $0x0, 0x20(%r12)
je 0xc322
movq (%rbx), %rdi
testq %rdi, %rdi
je 0xc322
movsbl 0x1(%r12), %esi
callq 0xa5dd
xorl %eax, %eax
movq %rax, (%r15)
movq %rax, (%r14)
movq %rax, (%rbx)
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::copy(FIT_BIT_SLOT, FIT_BIT_SLOT)
void FitBitSlicedA::copy(enum FIT_BIT_SLOT sourse_fbs, enum FIT_BIT_SLOT target_fbs) {
    NUMBER *soursePI = NULL;
    NUMBER ***sourseA = NULL;
    NUMBER **sourseB = NULL;
    get(sourse_fbs, soursePI, sourseA, sourseB);
    NUMBER *targetPI = NULL;
    NUMBER ***targetA = NULL;
    NUMBER **targetB = NULL;
    get(target_fbs, targetPI, targetA, targetB);
    
    copy(soursePI, sourseA, sourseB, targetPI, targetA, targetB);
}

// Assembly:
subq $0x38, %rsp
xorl %eax, %eax
movq %rax, 0x30(%rsp)
movq %rax, 0x28(%rsp)
movq %rax, 0x20(%rsp)
decl %esi
leaq 0x2796e(%rip), %rax     # 0x33df0
cmpl $0x7, %esi
ja 0xc4be
movl %esi, %ecx
movq (%rax,%rcx,8), %rsi
leaq 0x2799c(%rip), %r8      # 0x33e30
movq (%r8,%rcx,8), %r8
leaq 0x279d1(%rip), %r9      # 0x33e70
movq (%r9,%rcx,8), %rcx
movq (%rdi,%rsi), %rsi
movq %rsi, 0x30(%rsp)
movq (%rdi,%r8), %rsi
movq %rsi, 0x28(%rsp)
movq (%rdi,%rcx), %rcx
movq %rcx, 0x20(%rsp)
xorl %ecx, %ecx
movq %rcx, 0x18(%rsp)
movq %rcx, 0x10(%rsp)
movq %rcx, 0x8(%rsp)
decl %edx
cmpl $0x7, %edx
ja 0xc50d
movl %edx, %ecx
movq (%rax,%rcx,8), %rax
leaq 0x2794d(%rip), %rdx     # 0x33e30
movq (%rdx,%rcx,8), %rdx
leaq 0x27982(%rip), %rsi     # 0x33e70
movq (%rsi,%rcx,8), %rcx
movq (%rdi,%rax), %rax
movq %rax, 0x18(%rsp)
movq (%rdi,%rdx), %rax
movq %rax, 0x10(%rsp)
movq (%rdi,%rcx), %rax
movq %rax, 0x8(%rsp)
leaq 0x8(%rsp), %rax
movq %rax, (%rsp)
leaq 0x30(%rsp), %rsi
leaq 0x28(%rsp), %rdx
leaq 0x20(%rsp), %rcx
leaq 0x18(%rsp), %r8
leaq 0x10(%rsp), %r9
callq 0xc120
addq $0x38, %rsp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::add(FIT_BIT_SLOT, FIT_BIT_SLOT)
void FitBitSlicedA::add(enum FIT_BIT_SLOT sourse_fbs, enum FIT_BIT_SLOT target_fbs) {
    NUMBER *soursePI = NULL;
    NUMBER ***sourseA = NULL;
    NUMBER **sourseB = NULL;
    get(sourse_fbs, soursePI, sourseA, sourseB);
    NUMBER *targetPI = NULL;
    NUMBER ***targetA = NULL;
    NUMBER **targetB = NULL;
    get(target_fbs, targetPI, targetA, targetB);
    
    add(soursePI, sourseA, sourseB, targetPI, targetA, targetB);
}

// Assembly:
movl %edx, %eax
decl %esi
leaq 0x278ab(%rip), %r8      # 0x33df0
cmpl $0x7, %esi
ja 0xc574
movl %esi, %ecx
movq (%r8,%rcx,8), %rdx
leaq 0x278d9(%rip), %rsi     # 0x33e30
movq (%rsi,%rcx,8), %r9
leaq 0x2790e(%rip), %rsi     # 0x33e70
movq (%rsi,%rcx,8), %rcx
movq (%rdi,%rdx), %rsi
movq (%rdi,%r9), %rdx
movq (%rdi,%rcx), %rcx
jmp 0xc57a
xorl %esi, %esi
xorl %edx, %edx
xorl %ecx, %ecx
decl %eax
cmpl $0x7, %eax
ja 0xc5ab
movl %eax, %eax
movq (%r8,%rax,8), %r8
leaq 0x278a2(%rip), %r9      # 0x33e30
movq (%r9,%rax,8), %r9
leaq 0x278d7(%rip), %r10     # 0x33e70
movq (%r10,%rax,8), %rax
movq (%rdi,%r8), %r8
movq (%rdi,%r9), %r9
movq (%rdi,%rax), %rax
jmp 0xc5b3
xorl %r8d, %r8d
xorl %r9d, %r9d
xorl %eax, %eax
pushq %rax
movq %rax, (%rsp)
callq 0xc228
popq %rax
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::checkConvergence(FitResult*)
bool FitBitSlicedA::checkConvergence(FitResult *fr) {

	NUMBER criterion = 0;
    bool res = false;
    switch (this->tol_mode) {
        case 'p':
            for(NPAR i=0; i<this->nS; i++)
                if(this->pi != NULL) criterion += pow( this->pi[i]-this->PIm1[i], 2 )/*:0*/;
            
            for(NPAR z=0; z<this->nZ; z++)
                for(NPAR i=0; i<this->nS; i++)
                {
                    for(NPAR j=0; (this->A != NULL) && j<this->nS; j++) {
                        criterion += pow(this->A[z][i][j] - this->Am1[z][i][j],2);
                    }
                    
                }
            for(NPAR i=0; i<this->nS; i++)
            {
                for(NPAR k=0; (this->B != NULL) && k<this->nO; k++) {
                    criterion += pow(this->B[i][k] - this->Bm1[i][k],2);
                }
            }
            res = sqrt(criterion) < this->tol; // double the truth or false
            break;
        case 'l':
            criterion = (fr->pOmid-fr->pO)/fr->ndat;
            res = criterion < this->tol;
            break;
    }
    return res;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
movzbl 0xf0(%rdi), %eax
cmpl $0x6c, %eax
je 0xc611
cmpl $0x70, %eax
jne 0xc629
movsbq 0x1(%rbx), %rax
testq %rax, %rax
jle 0xc630
movq 0x10(%rbx), %rcx
movq 0x28(%rbx), %rdx
xorpd %xmm0, %xmm0
xorl %esi, %esi
testq %rcx, %rcx
je 0xc607
movsd (%rcx,%rsi,8), %xmm1
subsd (%rdx,%rsi,8), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm1, %xmm0
incq %rsi
cmpq %rsi, %rax
jne 0xc5f0
jmp 0xc634
movsd 0x10(%rsi), %xmm0
subsd 0x18(%rsi), %xmm0
cvtsi2sdl 0x24(%rsi), %xmm1
divsd %xmm1, %xmm0
jmp 0xc711
xorl %eax, %eax
jmp 0xc720
xorpd %xmm0, %xmm0
movsbq 0x2(%rbx), %rcx
testq %rcx, %rcx
jle 0xc6a4
testb %al, %al
setg %r8b
movq 0x18(%rbx), %rdx
movq 0x30(%rbx), %rsi
testq %rdx, %rdx
setne %dil
andb %r8b, %dil
xorl %r8d, %r8d
testb %al, %al
jle 0xc69c
xorl %r9d, %r9d
testb %dil, %dil
je 0xc694
movq (%rdx,%r8,8), %r10
movq (%r10,%r9,8), %r10
movq (%rsi,%r8,8), %r11
movq (%r11,%r9,8), %r11
xorl %r14d, %r14d
movsd (%r10,%r14,8), %xmm1
subsd (%r11,%r14,8), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm1, %xmm0
incq %r14
cmpq %rax, %r14
jl 0xc678
incq %r9
cmpq %rax, %r9
jne 0xc660
incq %r8
cmpq %rcx, %r8
jne 0xc659
testb %al, %al
jle 0xc6fc
movq 0x20(%rbx), %rcx
movq 0x38(%rbx), %rdx
testq %rcx, %rcx
setne %r8b
movsbq (%rbx), %rsi
testq %rsi, %rsi
setg %dil
andb %r8b, %dil
xorl %r8d, %r8d
testb %dil, %dil
je 0xc6f4
movq (%rcx,%r8,8), %r9
movq (%rdx,%r8,8), %r10
xorl %r11d, %r11d
movsd (%r9,%r11,8), %xmm1
subsd (%r10,%r11,8), %xmm1
mulsd %xmm1, %xmm1
addsd %xmm1, %xmm0
incq %r11
cmpq %rsi, %r11
jl 0xc6d8
incq %r8
cmpq %rax, %r8
jne 0xc6c8
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0xc70c
sqrtsd %xmm0, %xmm0
jmp 0xc711
callq 0x41a0
movsd 0xe8(%rbx), %xmm1
ucomisd %xmm0, %xmm1
seta %al
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedA.cpp
// Function: FitBitSlicedA::doLog10ScaleGentle(FIT_BIT_SLOT)
void FitBitSlicedA::doLog10ScaleGentle(enum FIT_BIT_SLOT fbs) {
    // fbs - gradient or direction
    NUMBER *a_PI = NULL;
    NUMBER ***a_A = NULL;
    NUMBER **a_B = NULL;
    get(fbs, a_PI, a_A, a_B);
	if(this->pi != NULL) doLog10Scale1DGentle(a_PI, this->pi, this->nS);
    if(this->A  != NULL) {
        for(NPAR z=0; z<this->nZ; z++) {
            doLog10Scale2DGentle(a_A[z],  this->A[z], this->nS, this->nS);
        }
    }
    if(this->B  != NULL) {
        doLog10Scale2DGentle(a_B,  this->B, this->nS, this->nO);
    }
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
decl %esi
cmpl $0x7, %esi
ja 0xc76b
movl %esi, %eax
leaq 0x276ad(%rip), %rcx     # 0x33df0
movq (%rcx,%rax,8), %rcx
leaq 0x276e2(%rip), %rdx     # 0x33e30
movq (%rdx,%rax,8), %rdx
leaq 0x27717(%rip), %rsi     # 0x33e70
movq (%rsi,%rax,8), %rax
movq (%rbx,%rcx), %rdi
movq (%rbx,%rdx), %r15
movq (%rbx,%rax), %r14
jmp 0xc773
xorl %edi, %edi
xorl %r15d, %r15d
xorl %r14d, %r14d
movq 0x10(%rbx), %rsi
testq %rsi, %rsi
je 0xc785
movsbl 0x1(%rbx), %edx
callq 0x9418
cmpq $0x0, 0x18(%rbx)
je 0xc7b9
cmpb $0x0, 0x2(%rbx)
jle 0xc7b9
xorl %r12d, %r12d
movq (%r15,%r12,8), %rdi
movq 0x18(%rbx), %rax
movq (%rax,%r12,8), %rsi
movsbl 0x1(%rbx), %ecx
movl %ecx, %edx
callq 0x96c1
incq %r12
movsbq 0x2(%rbx), %rax
cmpq %rax, %r12
jl 0xc795
movq 0x20(%rbx), %rsi
testq %rsi, %rsi
je 0xc7dc
movsbl 0x1(%rbx), %edx
movsbl (%rbx), %ecx
movq %r14, %rdi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
jmp 0x96c1
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq