LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedAB.cpp
// Function: FitBitSlicedAB::init(double*&, double***&, double***&)
void FitBitSlicedAB::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  = init3D<NUMBER>((NDAT)this->nZ, (NDAT)this->nS, (NDAT)this->nO);
        else
            toZero3D<NUMBER>(a_B, (NDAT)this->nZ, (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 0xce04
movq %rsi, %r12
movq (%rsi), %rdi
movsbq 0x1(%r14), %rdx
testq %rdi, %rdi
je 0xcdf3
testb %dl, %dl
jle 0xce04
shlq $0x3, %rdx
xorl %esi, %esi
callq 0x4120
jmp 0xce04
movl $0x8, %esi
movq %rdx, %rdi
callq 0x4130
movq %rax, (%r12)
movq %rbx, 0x20(%rsp)
movq %r14, 0x8(%rsp)
cmpq $0x0, 0x18(%r14)
je 0xce90
movq (%r15), %rcx
movq 0x8(%rsp), %rax
movsbl 0x2(%rax), %edi
movq %rcx, 0x10(%rsp)
testq %rcx, %rcx
je 0xce82
testb %dil, %dil
jle 0xce90
movslq %edi, %rax
movq %rax, 0x18(%rsp)
movq 0x8(%rsp), %rax
movsbq 0x1(%rax), %rbp
movl %ebp, %ebx
leaq (,%rbx,8), %r15
xorl %r14d, %r14d
testb %bpl, %bpl
jle 0xce76
movq 0x10(%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 0xce60
incq %r14
cmpq 0x18(%rsp), %r14
jne 0xce4f
jmp 0xce90
movsbl 0x1(%rax), %edx
movl %edx, %esi
callq 0xcaae
movq %rax, (%r15)
movq 0x8(%rsp), %rax
cmpq $0x0, 0x20(%rax)
je 0xcf18
movq 0x20(%rsp), %rbx
movq (%rbx), %r15
movsbl 0x2(%rax), %edi
testq %r15, %r15
je 0xcf09
testb %dil, %dil
jle 0xcf18
movslq %edi, %rcx
movq %rcx, 0x18(%rsp)
movsbq (%rax), %r13
movsbq 0x1(%rax), %rax
movl %r13d, %ebx
shlq $0x3, %rbx
movq %rax, 0x10(%rsp)
movl %eax, %ebp
xorl %r12d, %r12d
cmpb $0x0, 0x10(%rsp)
jle 0xcefd
xorl %r14d, %r14d
testb %r13b, %r13b
jle 0xcef5
movq (%r15,%r12,8), %rax
movq (%rax,%r14,8), %rdi
xorl %esi, %esi
movq %rbx, %rdx
callq 0x4120
incq %r14
cmpq %r14, %rbp
jne 0xcede
incq %r12
cmpq 0x18(%rsp), %r12
jne 0xced4
jmp 0xcf18
movsbl 0x1(%rax), %esi
movsbl (%rax), %edx
callq 0xcaae
movq %rax, (%rbx)
addq $0x28, %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/FitBitSlicedAB.cpp
// Function: FitBitSlicedAB::negate(double*&, double***&, double***&)
void FitBitSlicedAB::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 z=0; z<this->nZ; z++)
            for(NPAR i=0; i<this->nS; i++)
                for(NPAR m=0; m<this->nO; m++)
                    a_B[z][i][m] = -a_B[z][i][m];
    }
}

// Assembly:
pushq %rbx
cmpq $0x0, 0x10(%rdi)
je 0xcf5d
movsbq 0x1(%rdi), %rax
testq %rax, %rax
jle 0xcf5d
movq (%rsi), %rsi
xorl %r8d, %r8d
movaps 0x269f9(%rip), %xmm0    # 0x33940
movsd (%rsi,%r8,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%rsi,%r8,8)
incq %r8
cmpq %r8, %rax
jne 0xcf47
cmpq $0x0, 0x18(%rdi)
je 0xcfb7
movsbq 0x2(%rdi), %rax
testq %rax, %rax
jle 0xcfb7
movzbl 0x1(%rdi), %esi
xorl %r8d, %r8d
movaps 0x269c4(%rip), %xmm0    # 0x33940
testb %sil, %sil
jle 0xcfaf
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 0xcf91
incq %r10
cmpq %rsi, %r10
jne 0xcf8b
incq %r8
cmpq %rax, %r8
jne 0xcf7c
cmpq $0x0, 0x20(%rdi)
je 0xd018
movsbq 0x2(%rdi), %rax
testq %rax, %rax
jle 0xd018
movzbl 0x1(%rdi), %edx
movzbl (%rdi), %esi
xorl %edi, %edi
movaps 0x26968(%rip), %xmm0    # 0x33940
testb %dl, %dl
jle 0xd010
movq (%rcx), %r8
xorl %r9d, %r9d
testb %sil, %sil
jle 0xd008
movq (%r8,%rdi,8), %r10
movq (%r10,%r9,8), %r10
xorl %r11d, %r11d
movsd (%r10,%r11,8), %xmm1
xorps %xmm0, %xmm1
movlps %xmm1, (%r10,%r11,8)
incq %r11
cmpq %r11, %rsi
jne 0xcff2
incq %r9
cmpq %rdx, %r9
jne 0xcfe2
incq %rdi
cmpq %rax, %rdi
jne 0xcfd8
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedAB.cpp
// Function: FitBitSlicedAB::destroy(double*&, double***&, double***&)
void FitBitSlicedAB::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) free3D<NUMBER>(a_B, (NDAT)this->nZ, (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 0xd364
movq (%r15), %rdi
testq %rdi, %rdi
je 0xd364
callq 0x42a0
cmpq $0x0, 0x18(%r12)
je 0xd385
movq (%r14), %rdi
testq %rdi, %rdi
je 0xd385
movsbl 0x2(%r12), %esi
movsbl 0x1(%r12), %edx
callq 0xca4f
cmpq $0x0, 0x20(%r12)
je 0xd3a6
movq (%rbx), %rdi
testq %rdi, %rdi
je 0xd3a6
movsbl 0x2(%r12), %esi
movsbl 0x1(%r12), %edx
callq 0xca4f
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/FitBitSlicedAB.cpp
// Function: FitBitSlicedAB::copy(FIT_BIT_SLOT, FIT_BIT_SLOT)
void FitBitSlicedAB::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 0x269aa(%rip), %rax     # 0x33eb0
cmpl $0x7, %esi
ja 0xd542
movl %esi, %ecx
movq (%rax,%rcx,8), %rsi
leaq 0x269d8(%rip), %r8      # 0x33ef0
movq (%r8,%rcx,8), %r8
leaq 0x26a0d(%rip), %r9      # 0x33f30
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 0xd591
movl %edx, %ecx
movq (%rax,%rcx,8), %rax
leaq 0x26989(%rip), %rdx     # 0x33ef0
movq (%rdx,%rcx,8), %rdx
leaq 0x269be(%rip), %rsi     # 0x33f30
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 0xd16c
addq $0x38, %rsp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedAB.cpp
// Function: FitBitSlicedAB::add(FIT_BIT_SLOT, FIT_BIT_SLOT)
void FitBitSlicedAB::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 0x268e7(%rip), %r8      # 0x33eb0
cmpl $0x7, %esi
ja 0xd5f8
movl %esi, %ecx
movq (%r8,%rcx,8), %rdx
leaq 0x26915(%rip), %rsi     # 0x33ef0
movq (%rsi,%rcx,8), %r9
leaq 0x2694a(%rip), %rsi     # 0x33f30
movq (%rsi,%rcx,8), %rcx
movq (%rdi,%rdx), %rsi
movq (%rdi,%r9), %rdx
movq (%rdi,%rcx), %rcx
jmp 0xd5fe
xorl %esi, %esi
xorl %edx, %edx
xorl %ecx, %ecx
decl %eax
cmpl $0x7, %eax
ja 0xd62f
movl %eax, %eax
movq (%r8,%rax,8), %r8
leaq 0x268de(%rip), %r9      # 0x33ef0
movq (%r9,%rax,8), %r9
leaq 0x26913(%rip), %r10     # 0x33f30
movq (%r10,%rax,8), %rax
movq (%rdi,%r8), %r8
movq (%rdi,%r9), %r9
movq (%rdi,%rax), %rax
jmp 0xd637
xorl %r8d, %r8d
xorl %r9d, %r9d
xorl %eax, %eax
pushq %rax
movq %rax, (%rsp)
callq 0xd2a0
popq %rax
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/FitBitSlicedAB.cpp
// Function: FitBitSlicedAB::doLog10ScaleGentle(FIT_BIT_SLOT)
void FitBitSlicedAB::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);
//        doLog10Scale3DGentle(a_A,  this->A,  this->nZ, this->nS, this->nS);
	if(this->B  != NULL)
        for(NPAR z=0; z<this->nZ; z++)
            doLog10Scale2DGentle(a_B[z],  this->B[z], this->nS, this->nO);
//        doLog10Scale3DGentle(a_B,  this->B,  this->nZ, this->nS, this->nO);
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rdi, %rbx
decl %esi
cmpl $0x7, %esi
ja 0xd801
movl %esi, %eax
leaq 0x266d7(%rip), %rcx     # 0x33eb0
movq (%rcx,%rax,8), %rcx
leaq 0x2670c(%rip), %rdx     # 0x33ef0
movq (%rdx,%rax,8), %rdx
leaq 0x26741(%rip), %rsi     # 0x33f30
movq (%rsi,%rax,8), %rax
movq (%rbx,%rcx), %rdi
movq (%rbx,%rdx), %r15
movq (%rbx,%rax), %r14
jmp 0xd809
xorl %edi, %edi
xorl %r15d, %r15d
xorl %r14d, %r14d
movq 0x10(%rbx), %rsi
testq %rsi, %rsi
je 0xd81b
movsbl 0x1(%rbx), %edx
callq 0x9418
cmpq $0x0, 0x18(%rbx)
je 0xd84f
cmpb $0x0, 0x2(%rbx)
jle 0xd84f
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 0xd82b
cmpq $0x0, 0x20(%rbx)
je 0xd884
cmpb $0x0, 0x2(%rbx)
jle 0xd884
xorl %r15d, %r15d
movq (%r14,%r15,8), %rdi
movq 0x20(%rbx), %rax
movq (%rax,%r15,8), %rsi
movsbl 0x1(%rbx), %edx
movsbl (%rbx), %ecx
callq 0x96c1
incq %r15
movsbq 0x2(%rbx), %rax
cmpq %rax, %r15
jl 0xd85f
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/./StripedArray.h
// Function: StripedArray<signed char>::add(signed char)
void StripedArray<T>::add(T value) {
	if(this->size == (NDAT_MAX-1)) {
		fprintf(stderr, "Error! Maximum array size reached.\n");
		return;
	}
	
	// add stripe if necessary
	if( this->nstripes==0 || this->size_last_stripe==this->stripe_size)
		this->addStripe();
	// place data
	this->stripes[this->nstripes-1][this->size_last_stripe++] = value;
	this->size++;
}

// Assembly:
pushq %rbp
pushq %rbx
pushq %rax
cmpl $0x7ffffffe, (%rdi)     # imm = 0x7FFFFFFE
je 0xfbd4
movl %esi, %ebp
movq %rdi, %rbx
movl 0x8(%rdi), %eax
testl %eax, %eax
je 0xfba5
movl 0xc(%rbx), %ecx
cmpl 0x4(%rbx), %ecx
jne 0xfbb3
movq %rbx, %rdi
callq 0xff5a
movl 0x8(%rbx), %eax
movl 0xc(%rbx), %ecx
movq 0x10(%rbx), %rdx
cltq
movq -0x8(%rdx,%rax,8), %rax
leal 0x1(%rcx), %edx
movl %edx, 0xc(%rbx)
movslq %ecx, %rcx
movb %bpl, (%rax,%rcx)
incl (%rbx)
addq $0x8, %rsp
popq %rbx
popq %rbp
retq
callq 0x4987
jmp 0xfbcd
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/./StripedArray.h
// Function: StripedArray<int>::toArray()
T* StripedArray<T>::toArray() {
    T *result = (T*)malloc( (size_t)this->size*sizeof(T));
    for(NDAT i=0; i<this->nstripes; i++) {
        T* dest = &result[ i*this->stripe_size ];
        size_t sz = sizeof(T)*( (i<(this->nstripes-1))?this->stripe_size:this->size_last_stripe );
        memcpy( dest, this->stripes[i], sz );  
    }
    return result;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rdi, %r14
movslq (%rdi), %rdi
shlq $0x2, %rdi
callq 0x4300
movq %rax, (%rsp)
movslq 0x8(%r14), %rax
movq %rax, 0x8(%rsp)
testq %rax, %rax
jle 0xfd6c
movslq 0x4(%r14), %r12
movq 0x8(%rsp), %rax
decl %eax
movq %rax, 0x10(%rsp)
movl 0xc(%r14), %ebp
movq 0x10(%r14), %rbx
movq %r12, %r13
shlq $0x2, %r12
xorl %r15d, %r15d
movq (%rsp), %r14
cmpq 0x10(%rsp), %r15
movl %ebp, %eax
cmovbl %r13d, %eax
movslq %eax, %rdx
shlq $0x2, %rdx
movq (%rbx,%r15,8), %rsi
movq %r14, %rdi
callq 0x41c0
incq %r15
addq %r12, %r14
cmpq %r15, 0x8(%rsp)
jne 0xfd41
movq (%rsp), %rax
addq $0x18, %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/./StripedArray.h
// Function: StripedArray<signed char>::~StripedArray()
StripedArray<T>::~StripedArray() {
	for(NDAT i=0; i<this->nstripes;i++) {
//        if(this->complex)
//            for(NDAT j=0; j<( (i<(this->nstripes-1))?this->stripe_size:this->size_last_stripe );j++)
//                free(  (void *)(this->stripes[i][j]) );
		free(this->stripes[i]);
    }
	free(this->stripes);
//	for(NDAT i=0; i<this->nstripes;i++)
//		delete [] this->stripes[i];
//	delete [] this->stripes;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
cmpl $0x0, 0x8(%rdi)
jle 0xfcdb
xorl %r14d, %r14d
movq 0x10(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x42a0
incq %r14
movslq 0x8(%rbx), %rax
cmpq %rax, %r14
jl 0xfcc2
movq 0x10(%rbx), %rdi
addq $0x8, %rsp
popq %rbx
popq %r14
jmp 0x42a0
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/./StripedArray.h
// Function: StripedArray<signed char>::StripedArray(_IO_FILE*, int)
StripedArray<T>::StripedArray(FILE *f, NDAT N) {
	stripe_size = 20000;
//    complex = false;
	size = N;
	nstripes = (NDAT)ceil((double)N/stripe_size);
	size_last_stripe = (NDAT)fmod(N, stripe_size);
	stripes = (T **) calloc((size_t)nstripes, sizeof(T*));
    NDAT num;
    NDAT nread;
    for(NDAT i=0; i<nstripes; i++) {
        num = (i<(nstripes-1)?stripe_size:size_last_stripe);
        stripes[i] = (T*)calloc((size_t)num, sizeof(T)); // alloc data
        nread = (NDAT)fread (stripes[i], sizeof(T), (size_t)num, f);
        if(nread != num) {
            fprintf(stderr,"Error reading data from file\n");
            return;
        }
    }
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r14
movl $0x4e20, 0x4(%rdi)      # imm = 0x4E20
cvtsi2sd %edx, %xmm0
movsd %xmm0, (%rsp)
movl %edx, (%rdi)
divsd 0x2419e(%rip), %xmm0    # 0x33f80
callq 0x4090
cvttsd2si %xmm0, %eax
movl %eax, 0x8(%r14)
movsd (%rsp), %xmm0
movsd 0x24184(%rip), %xmm1    # 0x33f80
callq 0x4310
cvttsd2si %xmm0, %eax
movl %eax, 0xc(%r14)
movslq 0x8(%r14), %r15
movl $0x8, %esi
movq %r15, %rdi
callq 0x4130
movq %rax, 0x10(%r14)
testq %r15, %r15
jle 0xfe77
xorl %r12d, %r12d
decl %r15d
movslq %r15d, %rax
xorl %ecx, %ecx
cmpq %rax, %r12
setge %cl
movslq 0x4(%r14,%rcx,8), %r15
movl $0x1, %esi
movq %r15, %rdi
callq 0x4130
movq 0x10(%r14), %rcx
movq %rax, (%rcx,%r12,8)
movq 0x10(%r14), %rax
movq (%rax,%r12,8), %rdi
movl $0x1, %esi
movq %r15, %rdx
movq %rbx, %rcx
callq 0x4360
cmpl %eax, %r15d
jne 0xfe83
incq %r12
movslq 0x8(%r14), %r15
cmpq %r15, %r12
jl 0xfe26
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
callq 0x49db
jmp 0xfe77
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/./StripedArray.h
// Function: StripedArray<signed char>::addStripe()
void StripedArray<T>::addStripe() {
	if(this->size_last_stripe<this->stripe_size && this->nstripes>0) {
		fprintf(stderr, "Exception! You can only add stripes when none exits or previous is full.\n");
		return;
	}
	
	this->stripes = (T **) realloc(this->stripes,(size_t)(++this->nstripes)*sizeof(T*)); // add pointer
//	T** new_stripes = Calloc(T*, (this->nstripes));
//	memcpy(new_stripes, this->stripes, sizeof(T*)*this->nstripes );
//	free(this->stripes);
//	this->stripes = Malloc(T*,this->nstripes+1);
//	memcpy(this->stripes, new_stripes, sizeof(T*)*this->nstripes );
//	this->nstripes++;
//	free(new_stripes);
	
//	this->stripes[this->nstripes-1] = Calloc(T, (size_t)this->stripe_size); // alloc data
	this->stripes[this->nstripes-1] = (T*)calloc((size_t)this->stripe_size, sizeof(T)); // alloc data
	this->size_last_stripe = 0; // reset counter
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movl 0x8(%rdi), %eax
movl 0xc(%rdi), %ecx
cmpl 0x4(%rdi), %ecx
setl %cl
testl %eax, %eax
setg %dl
andb %cl, %dl
cmpb $0x1, %dl
je 0xffc0
movq %rdi, %rbx
cltq
movq 0x10(%rdi), %rdi
leaq 0x1(%rax), %rcx
movl %ecx, 0x8(%rbx)
leaq 0x8(,%rax,8), %rsi
callq 0x43f0
movq %rax, %r14
movq %rax, 0x10(%rbx)
movslq 0x4(%rbx), %rdi
movl $0x1, %esi
callq 0x4130
movslq 0x8(%rbx), %rcx
movq %rax, -0x8(%r14,%rcx,8)
movl $0x0, 0xc(%rbx)
addq $0x8, %rsp
popq %rbx
popq %r14
retq
callq 0x4a13
jmp 0xffb8
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/PULogistic.cpp
// Function: PULogistic::unite(double*, int, double*, int)
NUMBER PULogistic::unite(NUMBER* standard, NCAT nstandard, NUMBER* pooled, NCAT npooled) {
	NUMBER num = 1;
	NUMBER den = 1;
	NUMBER safe = 0;
	if(standard != NULL && nstandard > 0) {
		for(NCAT i=0; i<=nstandard; i++) {
			safe = safe01num(standard[i]);
			num *= 1 - safe;
			den *=     safe;
		}
	}
	if(pooled != NULL && npooled > 0) {
		for(NCAT i=0; i<=nstandard; i++) {
			safe = safe01num(pooled[i]);
			num *= 1 - safe;
			den *=     safe;
		}
	}
	return 1/( 1 + num/safe0num(den) );
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movl %r8d, %ebp
movq %rcx, %rbx
movl %edx, %r14d
testq %rsi, %rsi
sete %al
testl %edx, %edx
setle %cl
orb %al, %cl
jne 0x104e3
movq %rsi, %r15
leal 0x1(%r14), %r12d
movapd 0x234c3(%rip), %xmm0    # 0x33970
xorl %r13d, %r13d
movapd %xmm0, (%rsp)
movsd (%r15,%r13,8), %xmm0
callq 0x8dd2
movsd 0x23428(%rip), %xmm1    # 0x338f0
subsd %xmm0, %xmm1
unpcklpd %xmm0, %xmm1            # xmm1 = xmm1[0],xmm0[0]
movapd (%rsp), %xmm0
mulpd %xmm1, %xmm0
incq %r13
cmpq %r13, %r12
jne 0x104b0
jmp 0x104eb
movapd 0x23485(%rip), %xmm0    # 0x33970
testq %rbx, %rbx
sete %al
testl %ebp, %ebp
setle %cl
testl %r14d, %r14d
sets %dl
orb %cl, %dl
orb %al, %dl
jne 0x10539
incl %r14d
xorl %r15d, %r15d
movapd %xmm0, (%rsp)
movsd (%rbx,%r15,8), %xmm0
callq 0x8dd2
movsd 0x233d0(%rip), %xmm1    # 0x338f0
subsd %xmm0, %xmm1
unpcklpd %xmm0, %xmm1            # xmm1 = xmm1[0],xmm0[0]
movapd (%rsp), %xmm0
mulpd %xmm1, %xmm0
incq %r15
cmpq %r15, %r14
jne 0x10508
movapd %xmm0, (%rsp)
unpckhpd %xmm0, %xmm0            # xmm0 = xmm0[1,1]
callq 0x8df8
movapd (%rsp), %xmm1
divsd %xmm0, %xmm1
movsd 0x23398(%rip), %xmm0    # 0x338f0
addsd %xmm0, %xmm1
divsd %xmm1, %xmm0
addq $0x18, %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/PULogistic.cpp
// Function: PULogistic::derivativePair(double, double, signed char)
NUMBER PULogistic::derivativePair(NUMBER standard, NUMBER pooled, NPAR which) {
	NUMBER united = pair(standard, pooled);
	// which 0-standard, 1-pooled
	NUMBER deriv_param = 0;
	if(which==0) {
		NUMBER safe_standard = safe01num(standard);
		deriv_param = 1 / ( safe_standard * (1 - safe_standard) );
	} else if(which==1) {
		NUMBER safe_pooled   = safe01num(pooled);
		deriv_param = 1 / (   safe_pooled * (1 -   safe_pooled) );
	}
	return united * (1-united) * deriv_param;
}

// Assembly:
pushq %rbx
subq $0x10, %rsp
movl %esi, %ebx
movq (%rdi), %rax
movsd %xmm0, (%rsp)
movsd %xmm1, 0x8(%rsp)
callq *(%rax)
movapd %xmm0, %xmm2
testl %ebx, %ebx
je 0x105a3
xorpd %xmm1, %xmm1
cmpl $0x1, %ebx
jne 0x105cf
movsd 0x8(%rsp), %xmm0
movsd %xmm0, (%rsp)
movsd (%rsp), %xmm0
movsd %xmm2, (%rsp)
callq 0x8dd2
movsd 0x23336(%rip), %xmm1    # 0x338f0
movapd %xmm1, %xmm2
subsd %xmm0, %xmm2
mulsd %xmm0, %xmm2
divsd %xmm2, %xmm1
movsd (%rsp), %xmm2
movsd 0x23319(%rip), %xmm0    # 0x338f0
subsd %xmm2, %xmm0
mulsd %xmm2, %xmm0
mulsd %xmm1, %xmm0
addq $0x10, %rsp
popq %rbx
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/PULogistic.cpp
// Function: PULogistic::derivativeUnite(double*, int, double*, int, signed char, int)
NUMBER PULogistic::derivativeUnite(NUMBER* standard, NCAT nstandard, NUMBER* pooled, NCAT npooled, NPAR which, NCAT index) {
	NUMBER united = unite(standard, nstandard, pooled, npooled);
	// which 0-standard, 1-pooled
	NUMBER deriv_param = 0;
	if(which==0) {
		NUMBER safe_standard = safe01num(standard[index]);
		deriv_param = 1 / ( safe_standard * (1 - safe_standard) );
	} else if(which==1) {
		NUMBER safe_pooled   = safe01num(pooled[index]);
		deriv_param = 1 / (   safe_pooled * (1 -   safe_pooled) );
	}
	return united * (1-united) * deriv_param;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x10, %rsp
movl %r9d, %ebp
movq %rcx, %rbx
movq %rsi, %r14
movq (%rdi), %rax
callq *0x8(%rax)
movapd %xmm0, %xmm2
testl %ebp, %ebp
je 0x10615
xorpd %xmm1, %xmm1
cmpl $0x1, %ebp
jne 0x10649
movq %rbx, %r14
movslq 0x30(%rsp), %rax
movsd (%r14,%rax,8), %xmm0
movsd %xmm2, 0x8(%rsp)
callq 0x8dd2
movsd 0x232bd(%rip), %xmm1    # 0x338f0
movapd %xmm1, %xmm2
subsd %xmm0, %xmm2
mulsd %xmm0, %xmm2
divsd %xmm2, %xmm1
movsd 0x8(%rsp), %xmm2
movsd 0x2329f(%rip), %xmm0    # 0x338f0
subsd %xmm2, %xmm0
mulsd %xmm2, %xmm0
mulsd %xmm1, %xmm0
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/PULogistic.cpp
// Function: PULogistic::derivativeUnite(double, double, double, signed char)
NUMBER PULogistic::derivativeUnite(NUMBER united, NUMBER standard, NUMBER pooled, NPAR which) {
	// which 0-standard, 1-pooled
	NUMBER deriv_param = 0;
	if(which==0) {
		NUMBER safe_standard = safe01num(standard);
		deriv_param = 1 / ( safe_standard * (1 - safe_standard) );
	} else if(which==1) {
		NUMBER safe_pooled   = safe01num(pooled);
		deriv_param = 1 / (   safe_pooled * (1 -   safe_pooled) );
	}
	return united * (1-united) * deriv_param;
}

// Assembly:
pushq %rax
movapd %xmm0, %xmm3
testl %esi, %esi
je 0x1067c
xorpd %xmm1, %xmm1
cmpl $0x1, %esi
jne 0x106a7
movapd %xmm2, %xmm1
movapd %xmm1, %xmm0
movsd %xmm3, (%rsp)
callq 0x8dd2
movsd (%rsp), %xmm3
movsd 0x23259(%rip), %xmm1    # 0x338f0
movapd %xmm1, %xmm2
subsd %xmm0, %xmm2
mulsd %xmm0, %xmm2
divsd %xmm2, %xmm1
movsd 0x23241(%rip), %xmm0    # 0x338f0
subsd %xmm3, %xmm0
mulsd %xmm3, %xmm0
mulsd %xmm1, %xmm0
popq %rax
retq
nopl (%rax)
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/PUCorbettianAdditive.cpp
// Function: PUCorbettianAdditive::sumOdds(double*, int, double*, int)
NUMBER PUCorbettianAdditive::sumOdds(NUMBER* standard, NCAT nstandard, NUMBER* pooled, NCAT npooled) {
	NUMBER sum_odds = 0;
	if(standard != NULL && nstandard > 0) {
		for(NCAT i=0; i<=nstandard; i++) {
			sum_odds += standard[i] / safe0num(1-standard[i]);
		}
	}
	if(pooled != NULL && npooled > 0) {
		for(NCAT i=0; i<=nstandard; i++) {
			sum_odds += pooled[i];
		}
	}
	return sum_odds;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movl %r8d, %ebp
movq %rcx, %rbx
movl %edx, %r14d
testq %rsi, %rsi
sete %al
testl %edx, %edx
setle %cl
xorpd %xmm0, %xmm0
orb %al, %cl
jne 0x10781
movq %rsi, %r15
leal 0x1(%r14), %r12d
xorl %r13d, %r13d
movsd %xmm0, 0x10(%rsp)
movsd (%r15,%r13,8), %xmm1
movsd %xmm1, 0x8(%rsp)
movsd 0x23194(%rip), %xmm0    # 0x338f0
subsd %xmm1, %xmm0
callq 0x8df8
movsd 0x8(%rsp), %xmm1
divsd %xmm0, %xmm1
movsd 0x10(%rsp), %xmm0
addsd %xmm1, %xmm0
incq %r13
cmpq %r13, %r12
jne 0x10742
testq %rbx, %rbx
sete %al
testl %ebp, %ebp
setle %cl
testl %r14d, %r14d
sets %dl
orb %cl, %dl
orb %al, %dl
jne 0x107aa
incl %r14d
xorl %eax, %eax
addsd (%rbx,%rax,8), %xmm0
incq %rax
cmpq %rax, %r14
jne 0x1079d
addq $0x18, %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/PUCorbettianAdditive.cpp
// Function: PUCorbettianAdditive::derivativeUnite(double*, int, double*, int, signed char, int)
NUMBER PUCorbettianAdditive::derivativeUnite(NUMBER* standard, NCAT nstandard, NUMBER* pooled, NCAT npooled, NPAR which, NCAT index) {
	NUMBER sum_odds = sumOdds(standard, nstandard, pooled, npooled);
	// which 0-standard
	// \frac{1}{ (1+O_{ x_{ i },g_{ j } }^{ add } )^2 } \frac { 1 }{ (1-x_{ i })^2 }
	// which 1-pooled
	// \frac{1}{ (1+O_{ x_{ i },g_{ j } }^{ add } )2 }
	return 1 / ( pow(1+sum_odds,2) ) * ((which==0) ? ( 1/(pow( safe0num(1-standard[index]) ,2)) ) : 1 );
}

// Assembly:
pushq %rbp
pushq %rbx
pushq %rax
movl %r9d, %ebp
movq %rsi, %rbx
callq 0x1070e
movapd %xmm0, %xmm1
movsd 0x23090(%rip), %xmm0    # 0x338f0
movapd %xmm0, %xmm2
testb %bpl, %bpl
jne 0x108a2
movslq 0x20(%rsp), %rax
movsd 0x2307a(%rip), %xmm0    # 0x338f0
subsd (%rbx,%rax,8), %xmm0
movsd %xmm1, (%rsp)
callq 0x8df8
movsd 0x23063(%rip), %xmm2    # 0x338f0
movsd (%rsp), %xmm1
mulsd %xmm0, %xmm0
divsd %xmm0, %xmm2
movsd 0x2304e(%rip), %xmm0    # 0x338f0
addsd %xmm0, %xmm1
mulsd %xmm1, %xmm1
divsd %xmm1, %xmm0
mulsd %xmm2, %xmm0
addq $0x8, %rsp
popq %rbx
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::HMMProblem(param*)
HMMProblem::HMMProblem(struct param *param) {
    NPAR i;
    switch (param->structure) {
        case STRUCTURE_SKILL: // Expectation Maximization (Baum-Welch)
            for(i=0; i<3; i++) this->sizes[i] = param->nK;
            this->n_params = param->nK * 4;
            break;
//        case STRUCTURE_GROUP: // Gradient Descent by group
//            for(i=0; i<3; i++) this->sizes[i] = param->nG;
//            this->n_params = param->nG * 4;
//            break;
        default:
            fprintf(stderr,"Structure specified is not supported and should have been caught earlier\n");
            break;
    }
    init(param);
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r14
leaq 0x2dc25(%rip), %rax     # 0x3e5a8
movq %rax, (%rdi)
cmpb $0x1, 0x3e(%rsi)
jne 0x109be
movl 0x8c0(%rbx), %ecx
movl %ecx, 0xc(%r14)
movl %ecx, 0x10(%r14)
movl %ecx, 0x14(%r14)
movl 0x8c0(%rbx), %ecx
shll $0x2, %ecx
movl %ecx, 0x8(%r14)
movq 0x70(%rax), %rax
movq %r14, %rdi
movq %rbx, %rsi
addq $0x8, %rsp
popq %rbx
popq %r14
jmpq *%rax
movq 0x2e613(%rip), %rax     # 0x3efd8
movq (%rax), %rcx
leaq 0x23a12(%rip), %rdi     # 0x343e1
movl $0x49, %esi
movl $0x1, %edx
callq 0x4480
movq (%r14), %rax
jmp 0x109ab
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::checkPIABConstraints(double*, double**, double**)
bool HMMProblem::checkPIABConstraints(NUMBER* a_PI, NUMBER** a_A, NUMBER** a_B) {
	NPAR i, j;
	// check values
	NUMBER sum_pi = 0.0;
	NUMBER sum_a_row[this->p->nS];
	NUMBER sum_b_row[this->p->nS];
	for(i=0; i<this->p->nS; i++) {
		sum_a_row[i] = 0.0;
		sum_b_row[i] = 0.0;
	}
	
	for(i=0; i<this->p->nS; i++) {
		if( a_PI[i]>1.0 || a_PI[i]<0.0)
			return false;
		sum_pi += a_PI[i];
		for(j=0; j<this->p->nS; j++) {
			if( a_A[i][j]>1.0 || a_A[i][j]<0.0)
				return false;
			sum_a_row[i] += a_A[i][j];
		}// all states 2
		for(int m=0; m<this->p->nO; m++) {
			if( a_B[i][m]>1.0 || a_B[i][m]<0.0)
				return false;
			sum_b_row[i] += a_B[i][m];
		}// all observations
	}// all states
	if(sum_pi!=1.0)
		return false;
	for(i=0; i<this->p->nS; i++)
		if( sum_a_row[i]!=1.0 || sum_b_row[i]!=1.0)
			return false;
	return true;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %rdx, -0x48(%rbp)
movq 0x88(%rdi), %rax
movq %rax, -0x40(%rbp)
movzbl 0x8b8(%rax), %r12d
leal 0xf(,%r12,8), %eax
andl $-0x10, %eax
movq %rsp, %rbx
subq %rax, %rbx
movq %rbx, %rsp
movq %rsp, %r14
subq %rax, %r14
movq %r14, %rsp
testb %r12b, %r12b
jle 0x11170
movq %rsi, %r8
movq %rcx, -0x30(%rbp)
leal (,%r12,8), %r15d
xorl %r13d, %r13d
movq %rbx, %rdi
xorl %esi, %esi
movq %r15, %rdx
movq %r8, -0x38(%rbp)
callq 0x4120
movq %r14, %rdi
xorl %esi, %esi
movq %r15, %rdx
callq 0x4120
movq -0x38(%rbp), %r8
xorpd %xmm2, %xmm2
movsd 0x22872(%rip), %xmm0    # 0x338f0
xorpd %xmm1, %xmm1
movq -0x48(%rbp), %rdi
movsd (%r8,%r13,8), %xmm3
xorl %eax, %eax
ucomisd %xmm0, %xmm3
ja 0x11172
ucomisd %xmm3, %xmm2
ja 0x11172
addsd %xmm3, %xmm1
movq (%rdi,%r13,8), %rcx
xorl %edx, %edx
movsd (%rcx,%rdx,8), %xmm3
xorl %eax, %eax
ucomisd %xmm0, %xmm3
ja 0x11172
ucomisd %xmm3, %xmm2
ja 0x11172
addsd (%rbx,%r13,8), %xmm3
movsd %xmm3, (%rbx,%r13,8)
incq %rdx
cmpq %rdx, %r12
jne 0x110ac
movq -0x40(%rbp), %rax
movsbq 0x8b9(%rax), %rcx
testq %rcx, %rcx
jle 0x1111d
movq -0x30(%rbp), %rax
movq (%rax,%r13,8), %rdx
xorl %esi, %esi
movsd (%rdx,%rsi,8), %xmm3
xorl %eax, %eax
ucomisd %xmm0, %xmm3
ja 0x11172
ucomisd %xmm3, %xmm2
ja 0x11172
addsd (%r14,%r13,8), %xmm3
movsd %xmm3, (%r14,%r13,8)
incq %rsi
cmpq %rsi, %rcx
jne 0x110f6
incq %r13
cmpq %r12, %r13
jne 0x11086
movapd %xmm1, %xmm2
cmpeqsd %xmm0, %xmm2
movq %xmm2, %rax
andl $0x1, %eax
ucomisd %xmm0, %xmm1
jne 0x11172
jp 0x11172
testb %r12b, %r12b
jle 0x11172
xorl %ecx, %ecx
movsd (%rbx,%rcx,8), %xmm1
ucomisd %xmm0, %xmm1
jne 0x11170
jp 0x11170
movsd (%r14,%rcx,8), %xmm1
ucomisd %xmm0, %xmm1
jne 0x11170
jp 0x11170
incq %rcx
movb $0x1, %al
cmpq %rcx, %r12
jne 0x11149
jmp 0x11172
xorl %eax, %eax
leaq -0x28(%rbp), %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/HMMProblem.cpp
// Function: HMMProblem::readModel(char const*, bool)
void HMMProblem::readModel(const char *filename, bool overwrite) {
	FILE *fid = fopen(filename,"r");
	if(fid == NULL)
	{
		fprintf(stderr,"Can't read model file %s\n",filename);
		exit(1);
	}
	int max_line_length = 1024;
	char *line = Malloc(char,(size_t)max_line_length);
	NDAT line_no = 0;
    struct param initparam;
    set_param_defaults(&initparam);
    
    //
    // read solver info
    //
    if(overwrite)
        readSolverInfo(fid, this->p, &line_no);
    else
        readSolverInfo(fid, &initparam, &line_no);
    //
    // read model
    //
    readModelBody(fid, &initparam, &line_no, overwrite);
		
	fclose(fid);
	free(line);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x9b8, %rsp            # imm = 0x9B8
movl %edx, %ebp
movq %rsi, %r15
movq %rdi, %r14
leaq 0x21ea6(%rip), %rsi     # 0x33044
movq %r15, %rdi
callq 0x4290
testq %rax, %rax
je 0x11210
movq %rax, %rbx
movl $0x0, 0xc(%rsp)
leaq 0x10(%rsp), %rdi
callq 0x9c82
testb %bpl, %bpl
je 0x111ce
movq 0x88(%r14), %rsi
jmp 0x111d3
leaq 0x10(%rsp), %rsi
leaq 0xc(%rsp), %rdx
movq %rbx, %rdi
callq 0xa1f3
movq (%r14), %rax
movzbl %bpl, %r8d
leaq 0x10(%rsp), %rdx
leaq 0xc(%rsp), %rcx
movq %r14, %rdi
movq %rbx, %rsi
callq *0x68(%rax)
movq %rbx, %rdi
callq 0x4200
addq $0x9b8, %rsp            # imm = 0x9B8
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
movq %r15, %rdi
callq 0x4a70
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::destroy()
void HMMProblem::destroy() {
	// destroy model data
    if(this->null_obs_ratio != NULL) free(this->null_obs_ratio);
	if(this->pi != NULL) free2D<NUMBER>(this->pi, this->sizes[0]);
	if(this->A  != NULL) free3D<NUMBER>(this->A,  this->sizes[1], this->p->nS);
	if(this->B  != NULL) free3D<NUMBER>(this->B,  this->sizes[2], this->p->nS);
	if(this->lbPI!=NULL) free(this->lbPI);
	if(this->ubPI!=NULL) free(this->ubPI);
	if(this->lbA!=NULL) free2D<NUMBER>(this->lbA, this->p->nS);
	if(this->ubA!=NULL) free2D<NUMBER>(this->ubA, this->p->nS);
	if(this->lbB!=NULL) free2D<NUMBER>(this->lbB, this->p->nS);
	if(this->ubB!=NULL) free2D<NUMBER>(this->ubB, this->p->nS);
}

// Assembly:
pushq %rbx
movq %rdi, %rbx
movq 0x18(%rdi), %rdi
testq %rdi, %rdi
je 0x1125c
callq 0x42a0
movq 0x38(%rbx), %rdi
testq %rdi, %rdi
je 0x1126d
movl 0xc(%rbx), %esi
callq 0xa5dd
movq 0x40(%rbx), %rdi
testq %rdi, %rdi
je 0x1128c
movl 0x10(%rbx), %esi
movq 0x88(%rbx), %rax
movsbl 0x8b8(%rax), %edx
callq 0xca4f
movq 0x48(%rbx), %rdi
testq %rdi, %rdi
je 0x112ab
movl 0x14(%rbx), %esi
movq 0x88(%rbx), %rax
movsbl 0x8b8(%rax), %edx
callq 0xca4f
movq 0x50(%rbx), %rdi
testq %rdi, %rdi
je 0x112b9
callq 0x42a0
movq 0x68(%rbx), %rdi
testq %rdi, %rdi
je 0x112c7
callq 0x42a0
movq 0x58(%rbx), %rdi
testq %rdi, %rdi
je 0x112e3
movq 0x88(%rbx), %rax
movsbl 0x8b8(%rax), %esi
callq 0xa5dd
movq 0x70(%rbx), %rdi
testq %rdi, %rdi
je 0x112ff
movq 0x88(%rbx), %rax
movsbl 0x8b8(%rax), %esi
callq 0xa5dd
movq 0x60(%rbx), %rdi
testq %rdi, %rdi
je 0x1131b
movq 0x88(%rbx), %rax
movsbl 0x8b8(%rax), %esi
callq 0xa5dd
movq 0x78(%rbx), %rdi
testq %rdi, %rdi
je 0x11338
movq 0x88(%rbx), %rax
movsbl 0x8b8(%rax), %esi
popq %rbx
jmp 0xa5dd
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::getPI(int)
NUMBER* HMMProblem::getPI(NCAT x) {
	if( x > (this->sizes[0]-1) ) {
		fprintf(stderr,"While accessing PI, skill index %d exceeded last index of the data %d.\n", x, this->sizes[0]-1);
		exit(1);
	}
	return this->pi[x];
}

// Assembly:
movq %rdi, %rax
movl 0xc(%rdi), %edi
cmpl %esi, %edi
jle 0x1136a
movq 0x38(%rax), %rax
movslq %esi, %rcx
movq (%rax,%rcx,8), %rax
retq
pushq %rax
callq 0x4a94
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::getA(int)
NUMBER** HMMProblem::getA(NCAT x) {
	if( x > (this->sizes[1]-1) ) {
		fprintf(stderr,"While accessing A, skill index %d exceeded last index of the data %d.\n", x, this->sizes[1]-1);
		exit(1);
	}
	return this->A[x];
}

// Assembly:
movq %rdi, %rax
movl 0x10(%rdi), %edi
cmpl %esi, %edi
jle 0x11386
movq 0x40(%rax), %rax
movslq %esi, %rcx
movq (%rax,%rcx,8), %rax
retq
pushq %rax
callq 0x4abb
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::getB(int)
NUMBER** HMMProblem::getB(NCAT x) {
	if( x > (this->sizes[2]-1) ) {
		fprintf(stderr,"While accessing B, skill index %d exceeded last index of the data %d.\n", x, this->sizes[2]-1);
		exit(1);
	}
	return this->B[x];
}

// Assembly:
movq %rdi, %rax
movl 0x14(%rdi), %edi
cmpl %esi, %edi
jle 0x113a2
movq 0x48(%rax), %rax
movslq %esi, %rcx
movq (%rax,%rcx,8), %rax
retq
pushq %rax
callq 0x4ae2
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::getB(data*, signed char, signed char)
NUMBER HMMProblem::getB (struct data* dt, NPAR i, NPAR m) {
    // special attention for "unknonw" observations, i.e. the observation was there but we do not know what it is
    // in this case we simply return 1, effectively resulting in no change in \alpha or \beta vatiables
    if(m<0)
        return 1;
    switch(this->p->structure)
    {
        case STRUCTURE_SKILL:
            return this->B[dt->k][i][m];
            break;
//        case STRUCTURE_GROUP:
//            return this->B[dt->g][i][m];
//            break;
        default:
            fprintf(stderr,"Solver specified is not supported.\n");
            exit(1);
            break;
    }
}

// Assembly:
testb %cl, %cl
js 0x1149e
movq 0x88(%rdi), %rax
cmpb $0x1, 0x3e(%rax)
jne 0x114a7
movq 0x48(%rdi), %rax
movslq 0x58(%rsi), %rsi
movq (%rax,%rsi,8), %rax
movsbq %dl, %rdx
movq (%rax,%rdx,8), %rax
movzbl %cl, %ecx
movsd (%rax,%rcx,8), %xmm0
retq
movsd 0x2244a(%rip), %xmm0    # 0x338f0
retq
pushq %rax
callq 0x4b53
nop
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::initAlpha(int, data**)
void HMMProblem::initAlpha(NCAT xndat, struct data** x_data) {
	NPAR nS = this->p->nS;
//    int parallel_now = this->p->parallel==2; //PAR
//    #pragma omp parallel for schedule(dynamic) if(parallel_now) //PAR
	for(NCAT x=0; x<xndat; x++) {
        NDAT t;
        NPAR i;
		if( x_data[x]->cnt!=0 ) continue; // ... and the thing has not been computed yet (e.g. from group to skill)
		// alpha
		if( x_data[x]->alpha == NULL ) {
			x_data[x]->alpha = Calloc(NUMBER*, (size_t)x_data[x]->n);
			for(t=0; t<x_data[x]->n; t++)
				x_data[x]->alpha[t] = Calloc(NUMBER, (size_t)nS);
			
		} else {
			for(t=0; t<x_data[x]->n; t++)
				for(i=0; i<nS; i++)
					x_data[x]->alpha[t][i] = 0.0;
		}
		// p_O_param
		x_data[x]->p_O_param = 0.0;
		x_data[x]->loglik = 0.0;
        // c - scaling
		if( x_data[x]->c == NULL ) {
            x_data[x]->c = Calloc(NUMBER, (size_t)x_data[x]->n);
        } else {
			for(t=0; t<x_data[x]->n; t++)
                x_data[x]->c[t] = 0.0;
        }
	} // for all groups in skill
}

// Assembly:
testl %esi, %esi
jle 0x1161e
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %rdx, %rbx
movq 0x88(%rdi), %rax
movzbl 0x8b8(%rax), %r12d
movsbq %r12b, %r15
movl %r12d, %eax
shll $0x3, %eax
movq %rax, 0x20(%rsp)
movl %esi, %eax
movq %rax, 0x18(%rsp)
xorl %ebp, %ebp
movq %rdx, 0x10(%rsp)
movq %r15, 0x8(%rsp)
movq %r12, (%rsp)
movq (%rbx,%rbp,8), %r13
cmpl $0x0, 0x4(%r13)
je 0x11546
incq %rbp
cmpq 0x18(%rsp), %rbp
jne 0x1152c
jmp 0x11610
movq 0x28(%r13), %r15
movl (%r13), %r14d
testq %r15, %r15
je 0x1157e
testl %r14d, %r14d
jle 0x115cd
xorl %ebx, %ebx
testb %r12b, %r12b
jle 0x1156f
movq (%r15,%rbx,8), %rdi
xorl %esi, %esi
movq 0x20(%rsp), %rdx
callq 0x4120
incq %rbx
cmpq %rbx, %r14
jne 0x1155a
movq 0x10(%rsp), %rbx
jmp 0x115cd
movslq %r14d, %rdi
movl $0x8, %esi
callq 0x4130
movq %rax, 0x28(%r13)
movq (%rbx,%rbp,8), %r13
movl (%r13), %r14d
testl %r14d, %r14d
jle 0x115cd
xorl %r15d, %r15d
movq 0x8(%rsp), %r12
movl $0x8, %esi
movq %r12, %rdi
callq 0x4130
movq 0x28(%r13), %rcx
movq %rax, (%rcx,%r15,8)
incq %r15
movq (%rbx,%rbp,8), %r13
movslq (%r13), %r14
cmpq %r14, %r15
jl 0x115a4
movq (%rsp), %r12
xorps %xmm0, %xmm0
movups %xmm0, 0x48(%r13)
movq 0x20(%r13), %rdi
testq %rdi, %rdi
je 0x115fa
testl %r14d, %r14d
jle 0x11537
movl %r14d, %edx
shlq $0x3, %rdx
xorl %esi, %esi
callq 0x4120
jmp 0x11537
movslq %r14d, %rdi
movl $0x8, %esi
callq 0x4130
movq %rax, 0x20(%r13)
jmp 0x11537
addq $0x28, %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/HMMProblem.cpp
// Function: HMMProblem::initBeta(int, data**)
void HMMProblem::initBeta(NCAT xndat, struct data** x_data) {
	NPAR nS = this->p->nS;
//    int parallel_now = this->p->parallel==2; //PAR
//    #pragma omp parallel for schedule(dynamic) if(parallel_now) //PAR
	for(NCAT x=0; x<xndat; x++) {
        NDAT t;
        NPAR i;
		if( x_data[x]->cnt!=0 ) continue; // ... and the thing has not been computed yet (e.g. from group to skill)
		// beta
		if( x_data[x]->beta == NULL ) {
			x_data[x]->beta = Calloc(NUMBER*, (size_t)x_data[x]->n);
			for(t=0; t<x_data[x]->n; t++)
				x_data[x]->beta[t] = Calloc(NUMBER, (size_t)nS);
			
		} else {
			for(t=0; t<x_data[x]->n; t++)
				for(i=0; i<nS; i++)
					x_data[x]->beta[t][i] = 0.0;
		}
	} // for all groups in skill
}

// Assembly:
testl %esi, %esi
jle 0x11906
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rdx, %rbx
movq 0x88(%rdi), %rax
movzbl 0x8b8(%rax), %r12d
movsbq %r12b, %rax
movq %rax, 0x8(%rsp)
movl %r12d, %eax
shll $0x3, %eax
movq %rax, 0x10(%rsp)
movl %esi, %eax
movq %rax, (%rsp)
xorl %ebp, %ebp
movq (%rbx,%rbp,8), %r15
cmpl $0x0, 0x4(%r15)
je 0x1187c
incq %rbp
cmpq (%rsp), %rbp
jne 0x11866
jmp 0x118f8
movq 0x30(%r15), %r14
movl (%r15), %r13d
testq %r14, %r14
je 0x118af
testl %r13d, %r13d
jle 0x11871
xorl %r15d, %r15d
testb %r12b, %r12b
jle 0x118a5
movq (%r14,%r15,8), %rdi
xorl %esi, %esi
movq 0x10(%rsp), %rdx
callq 0x4120
incq %r15
cmpq %r15, %r13
jne 0x11890
jmp 0x11871
movslq %r13d, %rdi
movl $0x8, %esi
callq 0x4130
movq %rax, 0x30(%r15)
movq (%rbx,%rbp,8), %r14
cmpl $0x0, (%r14)
jle 0x11871
xorl %r13d, %r13d
movl $0x8, %esi
movq 0x8(%rsp), %rdi
callq 0x4130
movq 0x30(%r14), %rcx
movq %rax, (%rcx,%r13,8)
incq %r13
movq (%rbx,%rbp,8), %r14
movslq (%r14), %rax
cmpq %rax, %r13
jl 0x118cd
jmp 0x11871
addq $0x18, %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/HMMProblem.cpp
// Function: HMMProblem::setGradPI(FitBit*)
void HMMProblem::setGradPI(FitBit *fb){
    if(this->p->block_fitting[0]>0) return;
    NDAT t = 0, ndat = 0;
    NPAR i, o;
    struct data* dt;
    for(NCAT x=0; x<fb->xndat; x++) {
        dt = fb->x_data[x];
        if( dt->cnt!=0 ) continue;
        ndat += dt->n;
        o = this->p->dat_obs[ dt->ix[t] ];//->get( dt->ix[t] );
        for(i=0; i<this->p->nS; i++) {
            fb->gradPI[i] -= dt->beta[t][i] * ((o<0)?1:fb->B[i][o]) / safe0num(dt->p_O_param);
        }
    }
    if( this->p->Cslices>0 ) // penalty
        fb->addL2Penalty(FBV_PI, this->p, (NUMBER)ndat);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq 0x88(%rdi), %rdx
cmpb $0x0, 0x981(%rdx)
jle 0x12133
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rsi, %rbx
movl 0xe8(%rsi), %eax
testl %eax, %eax
jle 0x1221d
movq %rdi, %r14
xorl %r15d, %r15d
xorl %esi, %esi
movq 0xf0(%rbx), %rcx
movq (%rcx,%r15,8), %r12
cmpl $0x0, 0x4(%r12)
jne 0x12205
addl (%r12), %esi
cmpb $0x0, 0x8b8(%rdx)
jle 0x12205
movl %esi, 0xc(%rsp)
movq 0x870(%rdx), %rax
movq 0x10(%r12), %rcx
movslq (%rcx), %rcx
movzbl (%rax,%rcx), %r13d
xorl %ebp, %ebp
movsd 0x2175a(%rip), %xmm0    # 0x338f0
testb %r13b, %r13b
js 0x121a9
movq 0x20(%rbx), %rax
movq (%rax,%rbp,8), %rax
movsd (%rax,%r13,8), %xmm0
movq 0x30(%r12), %rax
movq (%rax), %rax
mulsd (%rax,%rbp,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%r12), %xmm0
callq 0x8df8
movsd 0x10(%rsp), %xmm1
divsd %xmm0, %xmm1
movq 0x58(%rbx), %rax
movsd (%rax,%rbp,8), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax,%rbp,8)
incq %rbp
movq 0x88(%r14), %rdx
movsbq 0x8b8(%rdx), %rax
cmpq %rax, %rbp
jl 0x1218e
movl 0xe8(%rbx), %eax
movl 0xc(%rsp), %esi
incq %r15
movslq %eax, %rcx
cmpq %rcx, %r15
jl 0x1214c
xorps %xmm0, %xmm0
cvtsi2sd %esi, %xmm0
jmp 0x12221
xorpd %xmm0, %xmm0
cmpb $0x0, 0x42(%rdx)
jle 0x12124
movq %rbx, %rdi
movl $0x1, %esi
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xb90a
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::setGradA(FitBit*)
void HMMProblem::setGradA (FitBit *fb){
    if(this->p->block_fitting[1]>0) return;
    NDAT t, ndat = 0;
    NPAR o, i, j;
    struct data* dt;
    for(NCAT x=0; x<fb->xndat; x++) {
        dt = fb->x_data[x];
        if( dt->cnt!=0 ) continue;
        ndat += dt->n;
        for(t=1; t<dt->n; t++) {
            o = this->p->dat_obs[ dt->ix[t] ];//->get( dt->ix[t] );
            for(i=0; i<this->p->nS; i++)
                for(j=0; j<this->p->nS; j++)
                    fb->gradA[i][j] -= dt->beta[t][j] * ((o<0)?1:fb->B[j][o]) * dt->alpha[t-1][i] / safe0num(dt->p_O_param);
        }
    }
    if( this->p->Cslices>0 ) // penalty
        fb->addL2Penalty(FBV_A, this->p, (NUMBER)ndat);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq 0x88(%rdi), %rdx
cmpb $0x0, 0x982(%rdx)
jle 0x12273
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rsi, %r13
movl 0xe8(%rsi), %eax
testl %eax, %eax
jle 0x123e9
movq %rdi, %rbx
xorl %r8d, %r8d
movl $0x0, 0xc(%rsp)
movq 0xf0(%r13), %rcx
movq %r8, 0x18(%rsp)
movq (%rcx,%r8,8), %r12
cmpl $0x0, 0x4(%r12)
jne 0x123ca
movl (%r12), %ecx
addl %ecx, 0xc(%rsp)
cmpl $0x2, %ecx
jl 0x123ca
movl $0x1, %eax
movq %rax, 0x10(%rsp)
movb 0x8b8(%rdx), %al
testb %al, %al
jle 0x123aa
movq 0x870(%rdx), %rcx
movq 0x10(%r12), %r8
movq 0x10(%rsp), %r9
movslq (%r8,%r9,4), %r8
movzbl (%rcx,%r8), %ebp
xorl %r15d, %r15d
testb %al, %al
jle 0x12390
xorl %r14d, %r14d
movsd 0x215e9(%rip), %xmm0    # 0x338f0
testb %bpl, %bpl
js 0x12319
movq 0x20(%r13), %rax
movq (%rax,%r14,8), %rax
movsd (%rax,%rbp,8), %xmm0
movq 0x30(%r12), %rax
movq 0x10(%rsp), %rcx
movq (%rax,%rcx,8), %rax
mulsd (%rax,%r14,8), %xmm0
movq 0x28(%r12), %rax
movq -0x8(%rax,%rcx,8), %rax
mulsd (%rax,%r15,8), %xmm0
movsd %xmm0, 0x20(%rsp)
movsd 0x48(%r12), %xmm0
callq 0x8df8
movsd 0x20(%rsp), %xmm1
divsd %xmm0, %xmm1
movq 0x60(%r13), %rax
movq (%rax,%r15,8), %rax
movsd (%rax,%r14,8), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax,%r14,8)
incq %r14
movq 0x88(%rbx), %rdx
movb 0x8b8(%rdx), %al
movsbq %al, %rcx
cmpq %rcx, %r14
jl 0x122ff
jmp 0x1239a
movb 0x8b8(%rdx), %al
movsbq %al, %rcx
incq %r15
cmpq %rcx, %r15
jl 0x122f4
movl (%r12), %ecx
movq 0x10(%rsp), %r9
incq %r9
movslq %ecx, %rax
movq %r9, 0x10(%rsp)
cmpq %rax, %r9
jl 0x122c9
movl 0xe8(%r13), %eax
movq 0x18(%rsp), %r8
incq %r8
movslq %eax, %rcx
cmpq %rcx, %r8
jl 0x12292
xorps %xmm0, %xmm0
cvtsi2sdl 0xc(%rsp), %xmm0
jmp 0x123ed
xorpd %xmm0, %xmm0
cmpb $0x0, 0x42(%rdx)
jle 0x12264
movq %r13, %rdi
movl $0x2, %esi
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xb90a
// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblem.cpp
// Function: HMMProblem::setGradB(FitBit*)
void HMMProblem::setGradB (FitBit *fb){
    if(this->p->block_fitting[2]>0) return;
    NDAT t, ndat = 0;
	NPAR o, o0, i;//, j;
    struct data* dt;
    for(NCAT x=0; x<fb->xndat; x++) {
        dt = fb->x_data[x];
        if( dt->cnt!=0 ) continue;
        ndat += dt->n;
//        for(t=0; t<dt->n; t++) { // old
        for(t=0; t<dt->n; t++) { // Levinson MMFST
            o  = this->p->dat_obs[ dt->ix[t] ];//->get( dt->ix[t] );
            o0 = this->p->dat_obs[ dt->ix[0] ];//->get( dt->ix[t] );
            if(o<0) // if no observation -- skip
                continue;
            for(i=0; i<this->p->nS; i++)
                fb->gradB[i][o] -= dt->alpha[t][i] * dt->beta[t][i] / safe0num(dt->p_O_param * fb->B[i][o]); // old
//            for(j=0; j<this->p->nS; j++)
//                if(t==0) {
//                    fb->gradB[j][o] -= (o0==o) * fb->pi[j] * dt->beta[0][j];
//                } else {
//                    for(i=0; i<this->p->nS; i++)
//                        fb->gradB[j][o] -= ( dt->alpha[t-1][i] * fb->A[i][j] * dt->beta[t][j] /*+ (o0==o) * getPI(dt,j) * dt->beta[0][j]*/ ) / safe0num(dt->p_O_param); // Levinson MMFST
//                }
        }
    }
    if( this->p->Cslices>0 ) // penalty
        fb->addL2Penalty(FBV_B, this->p, (NUMBER)ndat);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq 0x88(%rdi), %rdx
cmpb $0x0, 0x983(%rdx)
jle 0x1243f
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rsi, %rbx
movl 0xe8(%rsi), %eax
testl %eax, %eax
jle 0x12560
movq %rdi, %r14
xorl %esi, %esi
xorl %edi, %edi
movq 0xf0(%rbx), %rcx
movq (%rcx,%rsi,8), %r12
cmpl $0x0, 0x4(%r12)
jne 0x12548
movl (%r12), %ecx
addl %ecx, %edi
testl %ecx, %ecx
jle 0x12548
movl %edi, 0x4(%rsp)
movq %rsi, 0x8(%rsp)
xorl %r13d, %r13d
movq 0x870(%rdx), %rax
movq 0x10(%r12), %rsi
movslq (%rsi,%r13,4), %rsi
movsbq (%rax,%rsi), %rbp
testq %rbp, %rbp
js 0x1252a
cmpb $0x0, 0x8b8(%rdx)
jle 0x1252a
xorl %r15d, %r15d
movq 0x28(%r12), %rax
movq 0x30(%r12), %rcx
movq (%rax,%r13,8), %rax
movsd (%rax,%r15,8), %xmm0
movq (%rcx,%r13,8), %rax
mulsd (%rax,%r15,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%r12), %xmm0
movq 0x20(%rbx), %rax
movq (%rax,%r15,8), %rax
mulsd (%rax,%rbp,8), %xmm0
callq 0x8df8
movsd 0x10(%rsp), %xmm1
divsd %xmm0, %xmm1
movq 0x68(%rbx), %rax
movq (%rax,%r15,8), %rax
movsd (%rax,%rbp,8), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax,%rbp,8)
incq %r15
movq 0x88(%r14), %rdx
movsbq 0x8b8(%rdx), %rax
cmpq %rax, %r15
jl 0x124b2
movl (%r12), %ecx
incq %r13
movslq %ecx, %rax
cmpq %rax, %r13
jl 0x12488
movl 0xe8(%rbx), %eax
movq 0x8(%rsp), %rsi
movl 0x4(%rsp), %edi
incq %rsi
movslq %eax, %rcx
cmpq %rcx, %rsi
jl 0x12457
xorps %xmm0, %xmm0
cvtsi2sd %edi, %xmm0
jmp 0x12564
xorpd %xmm0, %xmm0
cmpb $0x0, 0x42(%rdx)
jle 0x12430
movq %rbx, %rdi
movl $0x3, %esi
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xb90a
nop