// File: /workspace/github_cpp_makefile_star2_O3/myudelson[P]hmm-scalable/HMMProblemAGK.cpp
// Function: HMMProblemAGK::GradientDescent()
NUMBER HMMProblemAGK::GradientDescent() {
	NCAT k=0, g=0, /*ki, gi, nX, */x;
    NCAT nK = this->p->nK, nG = this->p->nG;
    
	//
	// fit all as 1 skill first, set group gradients to 0, and do not fit them
	//
	if(this->p->single_skill>0) {
        FitBit *fb = new FitBit(this->p->nS, this->p->nO, this->p->nK, this->p->nG, this->p->tol, this->p->tol_mode);
        fb->link( HMMProblem::getPI(0), HMMProblem::getA(0), HMMProblem::getB(0), this->p->nSeq, this->p->k_data);// link skill 0 (we'll copy fit parameters to others
        if(this->p->block_fitting[0]!=0) fb->pi = NULL;
        if(this->p->block_fitting[1]!=0) fb->A  = NULL;
        if(this->p->block_fitting[2]!=0) fb->B  = NULL;
        
        fb->init(FBS_PARm1);
        fb->init(FBS_PARm2);
        fb->init(FBS_GRAD);
        if(this->p->solver==METHOD_CGD) {
            fb->init(FBS_GRADm1);
            fb->init(FBS_DIRm1);
        }
        
        NCAT* original_ks = Calloc(NCAT, (size_t)this->p->nSeq);
        for(x=0; x<this->p->nSeq; x++) { original_ks[x] = this->p->all_data[x].k; this->p->all_data[x].k = 0; } // save progonal k's
        FitResult fr = GradientDescentBit(fb);
        for(x=0; x<this->p->nSeq; x++) { this->p->all_data[x].k = original_ks[x]; } // restore original k's
        free(original_ks);
        if( !this->p->quiet )
            printf("single skill iter#%3d p(O|param)= %15.7f >> %15.7f, conv=%d\n", fr.iter,fr.pO0,fr.pO,fr.conv);
        delete fb;
    }
	//
	// Main fit
	//
    if( this->p->single_skill!=2 ) { // if not "force single skill"
        NCAT first_iteration_qualify = this->p->first_iteration_qualify; // at what iteration, qualification for skill/group convergence should start
        NCAT iterations_to_qualify = this->p->iterations_to_qualify; // how many concecutive iterations necessary for skill/group to qualify as converged
        NCAT* iter_qual_skill = Calloc(NCAT, (size_t)nK);
        NCAT* iter_qual_group = Calloc(NCAT, (size_t)nG);
        NCAT* iter_fail_skill = Calloc(NCAT, (size_t)nK); // counting concecutive number of failures to converge for skills
        NCAT* iter_fail_group = Calloc(NCAT, (size_t)nG); // counting concecutive number of failures to converge for groups
        int skip_k = 0, skip_g = 0;
        // utilize fitting larger data first
        
        int i = 0; // count runs
//        int parallel_now = this->p->parallel==1; //PAR
//        #pragma omp parallel if(parallel_now) shared(iter_qual_group,iter_qual_skill,iter_fail_skill,iter_fail_group)//PAR
//        {//PAR
            while(skip_k<nK || skip_g<nG) {
                //
                // Skills first
                //
//                printf("... while in,  run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
                if(skip_k<nK) {
//                    printf("... still k,   run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
//                    #pragma omp for schedule(dynamic) //PAR
                    for(k=0; k<nK; k++) { // for all A,B-by-skill
                        if(iter_qual_skill[k]==iterations_to_qualify || iter_fail_skill[k]==iterations_to_qualify)
                            continue;
//                        printf("... doing k,   run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
                        FitBit *fb = new FitBit(this->p->nS, this->p->nO, this->p->nK, this->p->nG, this->p->tol, this->p->tol_mode);
                        // link
                        fb->link( HMMProblem::getPI(k), HMMProblem::getA(k), HMMProblem::getB(k), this->p->k_numg[k], this->p->k_g_data[k]);// link skill 0 (we'll copy fit parameters to others
                        if(this->p->block_fitting[0]!=0) fb->pi = NULL;
                        if(this->p->block_fitting[1]!=0) fb->A  = NULL;
                        if(this->p->block_fitting[2]!=0) fb->B  = NULL;
                        
                        fb->Cslice = 0;
                        fb->init(FBS_PARm1);
                        fb->init(FBS_PARm2);
                        fb->init(FBS_GRAD);
                        if(this->p->solver==METHOD_CGD) {
                            fb->init(FBS_GRADm1);
                            fb->init(FBS_DIRm1);
                        }
                        
                        FitResult fr = GradientDescentBit(fb);
						
						// count number of concecutive failures
						if(fr.iter==this->p->maxiter) {
							iter_fail_skill[k]++;
						} else {
							iter_fail_skill[k]=0;
						}
						
                        // decide on convergence
                        if(i>=first_iteration_qualify || fb->xndat==0) {
                            if(fr.iter==1 /*e<=this->p->tol*/ || skip_g==nG || fb->xndat==0) { // converged quick, or don't care (others all converged
                                iter_qual_skill[k]++;
                                if(fb->xndat==0) {
                                    iter_qual_skill[k]=iterations_to_qualify;
                                    fr.conv = 1;
                                }
                                if(iter_qual_skill[k]==iterations_to_qualify || skip_g==nG) {// criterion met, or don't care (others all converged)
                                    if(skip_g==nG) iter_qual_skill[k]=iterations_to_qualify; // G not changing anymore
//                                    #pragma omp critical(update_skip_k)//PAR
//                                    {//PAR
                                        skip_k++;
//                                    }//PAR
                                    if( !this->p->quiet && ( /*(!conv && iter<this->p->maxiter) ||*/ (fr.conv || fr.iter==this->p->maxiter) )) {
                                        printf("run %2d skipK %4d skill %4d iter#%3d p(O|param)= %15.7f >> %15.7f, conv=%d\n",i,skip_k,k,fr.iter,fr.pO0,fr.pO,fr.conv);
                                    }
                                }
                            }
                            else
                                iter_qual_skill[k]=0;
                        } // decide on convergence
                        delete fb;
                    } // for all skills
                }
                //
                // PIg second
                //
                if(skip_g<nG){
//                    printf("... still g,   run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
//                    #pragma omp for schedule(dynamic)//PAR
                    for(g=0; g<nG; g++) { // for all PI-by-user
                        if(iter_qual_group[g]==iterations_to_qualify || iter_fail_group[g]==iterations_to_qualify)
                            continue;
//                        printf("... doing g,   run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
                        FitBit *fb = new FitBit(this->p->nS, this->p->nO, this->p->nK, this->p->nG, this->p->tol, this->p->tol_mode);
                        // link
                        // vvvvvvvvvvvvvvvvvvvvv ONLY PART THAT IS DIFFERENT FROM HMMProblemPiGK
                        fb->link(NULL, this->getAg(g), NULL, this->p->g_numk[g], this->p->g_k_data[g]);
                        fb->observe_non01_boundaries = 0;
                        if(this->p->block_fitting[0]!=0) fb->pi = NULL;
                        if(this->p->block_fitting[1]!=0) fb->A  = NULL;
                        if(this->p->block_fitting[2]!=0) fb->B  = NULL;
                        // ^^^^^^^^^^^^^^^^^^^^^
                        fb->Cslice = 1;
                        fb->init(FBS_PARm1);
                        fb->init(FBS_PARm2);
                        fb->init(FBS_GRAD);
                        if(this->p->solver==METHOD_CGD) {
                            fb->init(FBS_GRADm1);
                            fb->init(FBS_DIRm1);
                        }
                        FitResult fr = GradientDescentBit(fb);
						
						// count number of concecutive failures
						if(fr.iter==this->p->maxiter) {
							iter_fail_group[g]++;
						} else {
							iter_fail_group[g]=0;
						}
						
                        // decide on convergence
                        if(i>=first_iteration_qualify || fb->xndat==0) { //can qualify or  student had no skill labelled rows
                            if(fr.iter==1 /*e<=this->p->tol*/ || skip_k==nK || fb->xndat==0) { // converged quick, or don't care (others all converged), or
                                iter_qual_group[g]++;
                                if(fb->xndat==0) {
                                    iter_qual_group[g]=iterations_to_qualify;
                                    fr.conv = 1;
                                }
                                if(iter_qual_group[g]==iterations_to_qualify || skip_k==nK) {// criterion met, or don't care (others all converged)
                                    if(skip_k==nK) iter_qual_group[g]=iterations_to_qualify; // K not changing anymore
//                                    #pragma omp critical(update_skip_g)//PAR
//                                    {//PAR
                                        skip_g++;
//                                    }//PAR
                                    if( !this->p->quiet && ( /*(!conv && iter<this->p->maxiter) ||*/ (fr.conv || fr.iter==this->p->maxiter) )) {
                                        printf("run %2d skipG %4d group %4d iter#%3d p(O|param)= %15.7f >> %15.7f, conv=%d\n",i,skip_g,g,fr.iter,fr.pO0,fr.pO,fr.conv);
                                    }
                                }
                            }
                            else
                                iter_qual_group[g]=0;
                        } // decide on convergence
                        delete fb;
                    } // for all groups
                }
                //            fr.conv = checkConvergenceBig(fbs, nK+nG, this->p->tol, &crit);
                //
                //            computeAlphaAndPOParam(this->p->nSeq, this->p->k_data);
                //            ll = HMMProblem::getSumLogPOPara(this->p->nSeq, this->p->k_data);
                //            printf("*%i ll=%15.7f, crit=%12.10f\n",i,ll,crit);
//                #pragma omp single//PAR
//                {//PAR
                    i++;
//                    printf("... done run,  run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
//                }//PAR
//                printf("... while out, run=%3d, k=%6d, skippedK=%6d, g=%6d, skippedG=%6d, thread id=%d\n",i,k,skip_k,g,skip_g, omp_get_thread_num());//PAR
            }
//        }//PAR
        // recycle qualifications
        if( iter_qual_skill != NULL ) free(iter_qual_skill);
        if( iter_qual_group != NULL) free(iter_qual_group);
        if( iter_fail_skill != NULL ) free(iter_fail_skill);
        if( iter_fail_group != NULL) free(iter_fail_group);
        
    } // if not "force single skill
    
    // compute loglik
    return getSumLogPOPara(this->p->nSeq, this->p->k_data);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x98, %rsp
movq %rdi, %r13
movq 0x88(%rdi), %rbx
movslq 0x8c0(%rbx), %r15
movslq 0x8bc(%rbx), %r14
cmpb $0x0, 0x3d(%rbx)
movq %r14, 0x88(%rsp)
jle 0x1ec0a
movl $0x110, %edi            # imm = 0x110
callq 0x4210
movq %r14, %r8
movq %rax, %r14
movsbl 0x8b8(%rbx), %esi
movsbl 0x8b9(%rbx), %edx
movsd 0x28(%rbx), %xmm0
movsbl 0x30(%rbx), %r9d
movq %rax, %rdi
movq %r15, %rbp
movl %r15d, %ecx
callq 0xa858
movq %r13, %rdi
xorl %esi, %esi
callq 0x11354
movq %rax, %r15
movq %r13, %rdi
xorl %esi, %esi
callq 0x11370
movq %rax, %r12
movq %r13, %rdi
xorl %esi, %esi
callq 0x1138c
movq 0x88(%r13), %rcx
movl 0x8dc(%rcx), %r8d
movq 0x8f0(%rcx), %r9
movq %r14, %rdi
movq %r15, %rsi
movq %r12, %rdx
movq %rax, %rcx
callq 0xac70
movq 0x88(%r13), %rax
cmpb $0x0, 0x981(%rax)
je 0x1eaaf
movq $0x0, 0x10(%r14)
cmpb $0x0, 0x982(%rax)
je 0x1eac0
movq $0x0, 0x18(%r14)
cmpb $0x0, 0x983(%rax)
je 0x1ead1
movq $0x0, 0x20(%r14)
movq %r14, %rdi
movl $0x2, %esi
callq 0xaf10
movq %r14, %rdi
movl $0x7, %esi
callq 0xaf10
movq %r14, %rdi
movl $0x3, %esi
callq 0xaf10
movq 0x88(%r13), %rbx
cmpb $0x3, 0x3f(%rbx)
jne 0x1eb26
movq %r14, %rdi
movl $0x4, %esi
callq 0xaf10
movq %r14, %rdi
movl $0x6, %esi
callq 0xaf10
movq 0x88(%r13), %rbx
movslq 0x8dc(%rbx), %r12
movl $0x4, %esi
movq %r12, %rdi
callq 0x4130
movq %rax, %r15
testq %r12, %r12
jle 0x1eb67
movq 0x8e8(%rbx), %rax
addq $0x58, %rax
xorl %ecx, %ecx
movl (%rax), %edx
movl %edx, (%r15,%rcx,4)
movl $0x0, (%rax)
incq %rcx
addq $0x60, %rax
cmpq %rcx, %r12
jne 0x1eb4f
leaq 0x58(%rsp), %rdi
movq %r13, %rsi
movq %r14, %rdx
callq 0x14c44
movq 0x88(%r13), %rbx
movslq 0x8dc(%rbx), %rax
testq %rax, %rax
jle 0x1eba9
movq 0x8e8(%rbx), %rcx
addq $0x58, %rcx
xorl %edx, %edx
movl (%r15,%rdx,4), %esi
movl %esi, (%rcx)
incq %rdx
addq $0x60, %rcx
cmpq %rdx, %rax
jne 0x1eb97
movq %r15, %rdi
callq 0x42a0
cmpb $0x0, 0x3c(%rbx)
jne 0x1ebd9
movl 0x58(%rsp), %esi
movsd 0x60(%rsp), %xmm0
movsd 0x70(%rsp), %xmm1
movl 0x78(%rsp), %edx
leaq 0x15d97(%rip), %rdi     # 0x34969
movb $0x2, %al
callq 0x4050
movq %r14, %rdi
callq 0xa8fa
movl $0x110, %esi            # imm = 0x110
movq %r14, %rdi
callq 0x4220
movq 0x88(%r13), %rbx
cmpb $0x2, 0x3d(%rbx)
movq %rbp, %r15
movq 0x88(%rsp), %r14
je 0x1f260
movsbl 0x978(%rbx), %eax
movl %eax, 0x2c(%rsp)
movsbl 0x979(%rbx), %ebx
movl $0x4, %esi
movq %r15, %rdi
callq 0x4130
movq %rax, %r12
movl $0x4, %esi
movq %r14, %rdi
callq 0x4130
movq %rax, 0x20(%rsp)
movl $0x4, %esi
movq %r15, %rdi
callq 0x4130
movq %rax, 0x38(%rsp)
movl $0x4, %esi
movq %r14, %rdi
callq 0x4130
movq %rax, 0x30(%rsp)
movq %r15, 0x80(%rsp)
testl %r15d, %r15d
setg %al
testl %r14d, %r14d
setg %r9b
movl %eax, %ecx
orb %r9b, %cl
testb $0x1, %cl
je 0x1f226
movl %ebx, %edx
movl 0x80(%rsp), %esi
movl %r14d, %r8d
movl $0x0, 0x1c(%rsp)
xorl %ecx, %ecx
movl $0x0, 0x10(%rsp)
movl %ebx, 0x18(%rsp)
movq %r12, 0x40(%rsp)
movq %rsi, 0x48(%rsp)
movq %r8, 0x50(%rsp)
movb %r9b, 0xf(%rsp)
movl %ecx, 0x14(%rsp)
testl %esi, %esi
setle %cl
notb %al
orb %cl, %al
testb $0x1, %al
jne 0x1ef7f
xorl %r15d, %r15d
movq 0x38(%rsp), %rax
movq 0x20(%rsp), %r14
movl (%r12,%r15,4), %ebx
cmpl %edx, %ebx
je 0x1ef71
movl (%rax,%r15,4), %r14d
cmpl %edx, %r14d
jne 0x1ecfc
movq 0x20(%rsp), %r14
jmp 0x1ef71
movl $0x110, %edi            # imm = 0x110
callq 0x4210
movq %rax, %r12
movq 0x88(%r13), %rax
movsbl 0x8b8(%rax), %esi
movsbl 0x8b9(%rax), %edx
movl 0x8c0(%rax), %ecx
movl 0x8bc(%rax), %r8d
movsd 0x28(%rax), %xmm0
movsbl 0x30(%rax), %r9d
movq %r12, %rdi
callq 0xa858
movq %r13, %rdi
movl %r15d, %esi
callq 0x11354
movq %rax, 0x90(%rsp)
movq %r13, %rdi
movl %r15d, %esi
callq 0x11370
movq %rax, %rbp
movq %r13, %rdi
movl %r15d, %esi
callq 0x1138c
movq 0x88(%r13), %rcx
movq 0x8e0(%rcx), %rdx
movq 0x8f8(%rcx), %rcx
movl (%rdx,%r15,4), %r8d
movq (%rcx,%r15,8), %r9
movq %r12, %rdi
movq 0x90(%rsp), %rsi
movq %rbp, %rdx
movq %rax, %rcx
callq 0xac70
movq 0x88(%r13), %rax
cmpb $0x0, 0x981(%rax)
je 0x1edb5
movq $0x0, 0x10(%r12)
cmpb $0x0, 0x982(%rax)
je 0x1edc7
movq $0x0, 0x18(%r12)
cmpb $0x0, 0x983(%rax)
je 0x1edd9
movq $0x0, 0x20(%r12)
movb $0x0, 0xfa(%r12)
movq %r12, %rdi
movl $0x2, %esi
callq 0xaf10
movq %r12, %rdi
movl $0x7, %esi
callq 0xaf10
movq %r12, %rdi
movl $0x3, %esi
callq 0xaf10
movq 0x88(%r13), %rax
cmpb $0x3, 0x3f(%rax)
jne 0x1ee30
movq %r12, %rdi
movl $0x4, %esi
callq 0xaf10
movq %r12, %rdi
movl $0x6, %esi
callq 0xaf10
leaq 0x58(%rsp), %rdi
movq %r13, %rsi
movq %r12, %rdx
callq 0x14c44
movl 0x58(%rsp), %r8d
movq 0x88(%r13), %rcx
movl 0x38(%rcx), %eax
incl %r14d
cmpl %eax, %r8d
movl $0x0, %edx
cmovnel %edx, %r14d
movq 0x38(%rsp), %rdx
movl %r14d, (%rdx,%r15,4)
movl 0x1c(%rsp), %edx
cmpl 0x2c(%rsp), %edx
jge 0x1ee86
cmpl $0x0, 0xe8(%r12)
je 0x1ee86
movq 0x20(%rsp), %r14
jmp 0x1ef44
movl 0xe8(%r12), %edx
cmpl $0x1, %r8d
movq 0x88(%rsp), %rsi
movq 0x20(%rsp), %r14
movq 0x50(%rsp), %r9
je 0x1eec3
cmpl %r9d, 0x14(%rsp)
je 0x1eec3
testl %edx, %edx
movq 0x40(%rsp), %rdi
je 0x1eee2
movl $0x0, (%rdi,%r15,4)
jmp 0x1ef44
incl %ebx
movq 0x40(%rsp), %rdi
movl %ebx, (%rdi,%r15,4)
testl %edx, %edx
je 0x1eee2
movl 0x18(%rsp), %edx
cmpl %edx, %ebx
je 0x1eef2
cmpl %esi, 0x14(%rsp)
jne 0x1ef44
jmp 0x1eef2
movl 0x18(%rsp), %edx
movl %edx, (%rdi,%r15,4)
movl $0x1, 0x78(%rsp)
cmpl %r9d, 0x14(%rsp)
jne 0x1eefd
movl %edx, (%rdi,%r15,4)
incl 0x10(%rsp)
cmpb $0x0, 0x3c(%rcx)
jne 0x1ef44
cmpl %eax, %r8d
sete %al
movl 0x78(%rsp), %r9d
testl %r9d, %r9d
setne %cl
orb %al, %cl
cmpb $0x1, %cl
jne 0x1ef44
movsd 0x60(%rsp), %xmm0
movsd 0x70(%rsp), %xmm1
leaq 0x15a74(%rip), %rdi     # 0x349a6
movl 0x1c(%rsp), %esi
movl 0x10(%rsp), %edx
movl %r15d, %ecx
movb $0x2, %al
callq 0x4050
movq %r12, %rdi
callq 0xa8fa
movl $0x110, %esi            # imm = 0x110
movq %r12, %rdi
callq 0x4220
movq 0x38(%rsp), %rax
movl 0x18(%rsp), %edx
movq 0x40(%rsp), %r12
movq 0x48(%rsp), %rsi
movq 0x50(%rsp), %r8
incq %r15
cmpq %r15, %rsi
jne 0x1ecdd
jmp 0x1ef84
movq 0x20(%rsp), %r14
testl %r8d, %r8d
setg %al
movb 0xf(%rsp), %cl
andb %al, %cl
cmpb $0x1, %cl
jne 0x1f201
xorl %r15d, %r15d
movq 0x30(%rsp), %rax
movl (%r14,%r15,4), %ebx
cmpl %edx, %ebx
je 0x1f1f5
movl (%rax,%r15,4), %ebp
cmpl %edx, %ebp
je 0x1f1f5
movl $0x110, %edi            # imm = 0x110
callq 0x4210
movq %rax, %r12
movq 0x88(%r13), %rax
movsbl 0x8b8(%rax), %esi
movsbl 0x8b9(%rax), %edx
movl 0x8c0(%rax), %ecx
movl 0x8bc(%rax), %r8d
movsd 0x28(%rax), %xmm0
movsbl 0x30(%rax), %r9d
movq %r12, %rdi
callq 0xa858
movq %r13, %rdi
movl %r15d, %esi
callq 0x1e1d8
movq 0x88(%r13), %rcx
movq 0x900(%rcx), %rdx
movq 0x910(%rcx), %rcx
movl (%rdx,%r15,4), %r8d
movq (%rcx,%r15,8), %r9
movq %r12, %rdi
xorl %esi, %esi
movq %rax, %rdx
xorl %ecx, %ecx
callq 0xac70
movb $0x0, 0xf9(%r12)
movq 0x88(%r13), %rax
cmpb $0x0, 0x981(%rax)
je 0x1f053
movq $0x0, 0x10(%r12)
cmpb $0x0, 0x982(%rax)
je 0x1f065
movq $0x0, 0x18(%r12)
cmpb $0x0, 0x983(%rax)
je 0x1f077
movq $0x0, 0x20(%r12)
movb $0x1, 0xfa(%r12)
movq %r12, %rdi
movl $0x2, %esi
callq 0xaf10
movq %r12, %rdi
movl $0x7, %esi
callq 0xaf10
movq %r12, %rdi
movl $0x3, %esi
callq 0xaf10
movq 0x88(%r13), %rax
cmpb $0x3, 0x3f(%rax)
jne 0x1f0ce
movq %r12, %rdi
movl $0x4, %esi
callq 0xaf10
movq %r12, %rdi
movl $0x6, %esi
callq 0xaf10
leaq 0x58(%rsp), %rdi
movq %r13, %rsi
movq %r12, %rdx
callq 0x14c44
movl 0x58(%rsp), %r8d
movq 0x88(%r13), %rcx
movl 0x38(%rcx), %eax
incl %ebp
cmpl %eax, %r8d
movl $0x0, %edx
cmovnel %edx, %ebp
movq 0x30(%rsp), %rdx
movl %ebp, (%rdx,%r15,4)
movl 0x1c(%rsp), %edx
cmpl 0x2c(%rsp), %edx
jge 0x1f11c
cmpl $0x0, 0xe8(%r12)
jne 0x1f1c8
movl 0xe8(%r12), %edx
cmpl $0x1, %r8d
movq 0x48(%rsp), %rsi
je 0x1f146
cmpl %esi, 0x10(%rsp)
je 0x1f146
testl %edx, %edx
je 0x1f167
movl $0x0, (%r14,%r15,4)
jmp 0x1f1c8
incl %ebx
movl %ebx, (%r14,%r15,4)
testl %edx, %edx
je 0x1f167
movl 0x18(%rsp), %edx
cmpl %edx, %ebx
je 0x1f177
movl 0x10(%rsp), %edi
cmpl 0x80(%rsp), %edi
jne 0x1f1c8
jmp 0x1f177
movl 0x18(%rsp), %edx
movl %edx, (%r14,%r15,4)
movl $0x1, 0x78(%rsp)
cmpl %esi, 0x10(%rsp)
jne 0x1f181
movl %edx, (%r14,%r15,4)
incl 0x14(%rsp)
cmpb $0x0, 0x3c(%rcx)
jne 0x1f1c8
cmpl %eax, %r8d
sete %al
movl 0x78(%rsp), %r9d
testl %r9d, %r9d
setne %cl
orb %al, %cl
cmpb $0x1, %cl
jne 0x1f1c8
movsd 0x60(%rsp), %xmm0
movsd 0x70(%rsp), %xmm1
leaq 0x1583c(%rip), %rdi     # 0x349f2
movl 0x1c(%rsp), %esi
movl 0x14(%rsp), %edx
movl %r15d, %ecx
movb $0x2, %al
callq 0x4050
movq %r12, %rdi
callq 0xa8fa
movl $0x110, %esi            # imm = 0x110
movq %r12, %rdi
callq 0x4220
movq 0x30(%rsp), %rax
movl 0x18(%rsp), %edx
movq 0x40(%rsp), %r12
movq 0x48(%rsp), %rsi
movq 0x50(%rsp), %r8
incq %r15
cmpq %r15, %r8
jne 0x1efa1
incl 0x1c(%rsp)
movl 0x14(%rsp), %ecx
cmpl %r8d, %ecx
setl %r9b
cmpl %esi, 0x10(%rsp)
setl %al
jl 0x1ecb6
cmpl %r8d, %ecx
jl 0x1ecb6
testq %r12, %r12
je 0x1f233
movq %r12, %rdi
callq 0x42a0
movq 0x20(%rsp), %rdi
testq %rdi, %rdi
je 0x1f242
callq 0x42a0
movq 0x38(%rsp), %rdi
testq %rdi, %rdi
je 0x1f251
callq 0x42a0
movq 0x30(%rsp), %rdi
testq %rdi, %rdi
je 0x1f260
callq 0x42a0
movq 0x88(%r13), %rax
movl 0x8dc(%rax), %esi
movq 0x8f0(%rax), %rdx
movq %r13, %rdi
callq 0x114ae
addq $0x98, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rax, %rbx
movl $0x110, %esi            # imm = 0x110
movq %r14, %rdi
jmp 0x1f2a8
jmp 0x1f29d
movq %rax, %rbx
movl $0x110, %esi            # imm = 0x110
movq %r12, %rdi
callq 0x4220
movq %rbx, %rdi
callq 0x4430
nop