// File: myudelson[P]hmm-scalable/FitBit.cpp
// Function: FitBit::checkConvergence(FitResult*)
bool FitBit::checkConvergence(FitResult *fr) {
    
    NUMBER criterion = 0;
    NUMBER criterion_oscil = 0; // oscillation between PAR and PARm1, i.e. PAR is close to PARm2
    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 j=0; (this->A != NULL) && j<this->nS; j++) {
                    criterion += pow(this->A[i][j] - this->Am1[i][j],2);
                }
                for(NPAR k=0; (this->B != NULL) && k<this->nO; k++) {
                    criterion += pow(this->B[i][k] - this->Bm1[i][k],2);
                }
            }
            
            if( (!(sqrt(criterion) < this->tol)) && ( this->PIm2 != NULL || this->Am2 != NULL || this->Bm2 != NULL ) ) {
                for(NPAR i=0; i<this->nS; i++)
                {
                    if(this->pi != NULL) criterion_oscil += pow( this->pi[i]-this->PIm2[i], 2 )/*:0*/;
                    for(NPAR j=0; (this->A != NULL) && j<this->nS; j++) {
                        criterion_oscil += pow(this->A[i][j] - this->Am2[i][j],2);
                    }
                    for(NPAR k=0; (this->B != NULL) && k<this->nO; k++) {
                        criterion_oscil += pow(this->B[i][k] - this->Bm2[i][k],2);
                    }
                }
                res = sqrt(criterion_oscil) < this->tol; // double the truth or false
            }
            else
                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 %rbp
movq %rsp, %rbp
subq $0x40, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x38(%rbp)
movq $0x0, -0x18(%rbp)
movq $0x0, -0x20(%rbp)
movb $0x0, -0x21(%rbp)
movsbl 0x108(%rax), %eax
movl %eax, -0x2c(%rbp)
subl $0x6c, %eax
je 0xe9e3
jmp 0xe671
movl -0x2c(%rbp), %eax
subl $0x70, %eax
jne 0xea24
jmp 0xe67f
movb $0x0, -0x22(%rbp)
movq -0x38(%rbp), %rcx
movsbl -0x22(%rbp), %eax
movsbl 0x1(%rcx), %ecx
cmpl %ecx, %eax
jge 0xe7e6
movq -0x38(%rbp), %rax
cmpq $0x0, 0x10(%rax)
je 0xe6d6
movq -0x38(%rbp), %rax
movq 0x10(%rax), %rcx
movsbq -0x22(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movq 0x28(%rax), %rax
movsbq -0x22(%rbp), %rcx
subsd (%rax,%rcx,8), %xmm0
movl $0x2, %edi
callq 0xf6b0
addsd -0x18(%rbp), %xmm0
movsd %xmm0, -0x18(%rbp)
movb $0x0, -0x23(%rbp)
movq -0x38(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, 0x18(%rcx)
movb %al, -0x39(%rbp)
je 0xe6fe
movq -0x38(%rbp), %rcx
movsbl -0x23(%rbp), %eax
movsbl 0x1(%rcx), %ecx
cmpl %ecx, %eax
setl %al
movb %al, -0x39(%rbp)
movb -0x39(%rbp), %al
testb $0x1, %al
jne 0xe707
jmp 0xe757
movq -0x38(%rbp), %rax
movq 0x18(%rax), %rcx
movsbq -0x22(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movsbq -0x23(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movq 0x30(%rax), %rax
movsbq -0x22(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x23(%rbp), %rcx
subsd (%rax,%rcx,8), %xmm0
movl $0x2, %edi
callq 0xf6b0
addsd -0x18(%rbp), %xmm0
movsd %xmm0, -0x18(%rbp)
movb -0x23(%rbp), %al
addb $0x1, %al
movb %al, -0x23(%rbp)
jmp 0xe6da
movb $0x0, -0x24(%rbp)
movq -0x38(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, 0x20(%rcx)
movb %al, -0x3a(%rbp)
je 0xe77e
movq -0x38(%rbp), %rcx
movsbl -0x24(%rbp), %eax
movsbl (%rcx), %ecx
cmpl %ecx, %eax
setl %al
movb %al, -0x3a(%rbp)
movb -0x3a(%rbp), %al
testb $0x1, %al
jne 0xe787
jmp 0xe7d7
movq -0x38(%rbp), %rax
movq 0x20(%rax), %rcx
movsbq -0x22(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movsbq -0x24(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movq 0x38(%rax), %rax
movsbq -0x22(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x24(%rbp), %rcx
subsd (%rax,%rcx,8), %xmm0
movl $0x2, %edi
callq 0xf6b0
addsd -0x18(%rbp), %xmm0
movsd %xmm0, -0x18(%rbp)
movb -0x24(%rbp), %al
addb $0x1, %al
movb %al, -0x24(%rbp)
jmp 0xe75b
jmp 0xe7d9
movb -0x22(%rbp), %al
addb $0x1, %al
movb %al, -0x22(%rbp)
jmp 0xe683
movsd -0x18(%rbp), %xmm0
callq 0x5220
movq -0x38(%rbp), %rax
movaps %xmm0, %xmm1
movsd 0x100(%rax), %xmm0
ucomisd %xmm1, %xmm0
ja 0xe9bc
movq -0x38(%rbp), %rax
cmpq $0x0, 0x40(%rax)
jne 0xe82e
movq -0x38(%rbp), %rax
cmpq $0x0, 0x48(%rax)
jne 0xe82e
movq -0x38(%rbp), %rax
cmpq $0x0, 0x50(%rax)
je 0xe9bc
movb $0x0, -0x25(%rbp)
movq -0x38(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movsbl 0x1(%rcx), %ecx
cmpl %ecx, %eax
jge 0xe995
movq -0x38(%rbp), %rax
cmpq $0x0, 0x10(%rax)
je 0xe885
movq -0x38(%rbp), %rax
movq 0x10(%rax), %rcx
movsbq -0x25(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movq 0x40(%rax), %rax
movsbq -0x25(%rbp), %rcx
subsd (%rax,%rcx,8), %xmm0
movl $0x2, %edi
callq 0xf6b0
addsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x20(%rbp)
movb $0x0, -0x26(%rbp)
movq -0x38(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, 0x18(%rcx)
movb %al, -0x3b(%rbp)
je 0xe8ad
movq -0x38(%rbp), %rcx
movsbl -0x26(%rbp), %eax
movsbl 0x1(%rcx), %ecx
cmpl %ecx, %eax
setl %al
movb %al, -0x3b(%rbp)
movb -0x3b(%rbp), %al
testb $0x1, %al
jne 0xe8b6
jmp 0xe906
movq -0x38(%rbp), %rax
movq 0x18(%rax), %rcx
movsbq -0x25(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movsbq -0x26(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movq 0x48(%rax), %rax
movsbq -0x25(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x26(%rbp), %rcx
subsd (%rax,%rcx,8), %xmm0
movl $0x2, %edi
callq 0xf6b0
addsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x20(%rbp)
movb -0x26(%rbp), %al
addb $0x1, %al
movb %al, -0x26(%rbp)
jmp 0xe889
movb $0x0, -0x27(%rbp)
movq -0x38(%rbp), %rcx
xorl %eax, %eax
cmpq $0x0, 0x20(%rcx)
movb %al, -0x3c(%rbp)
je 0xe92d
movq -0x38(%rbp), %rcx
movsbl -0x27(%rbp), %eax
movsbl (%rcx), %ecx
cmpl %ecx, %eax
setl %al
movb %al, -0x3c(%rbp)
movb -0x3c(%rbp), %al
testb $0x1, %al
jne 0xe936
jmp 0xe986
movq -0x38(%rbp), %rax
movq 0x20(%rax), %rcx
movsbq -0x25(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movsbq -0x27(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movq 0x50(%rax), %rax
movsbq -0x25(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x27(%rbp), %rcx
subsd (%rax,%rcx,8), %xmm0
movl $0x2, %edi
callq 0xf6b0
addsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x20(%rbp)
movb -0x27(%rbp), %al
addb $0x1, %al
movb %al, -0x27(%rbp)
jmp 0xe90a
jmp 0xe988
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0xe832
movsd -0x20(%rbp), %xmm0
callq 0x5220
movq -0x38(%rbp), %rax
movaps %xmm0, %xmm1
movsd 0x100(%rax), %xmm0
ucomisd %xmm1, %xmm0
seta %al
andb $0x1, %al
movb %al, -0x21(%rbp)
jmp 0xe9e1
movsd -0x18(%rbp), %xmm0
callq 0x5220
movq -0x38(%rbp), %rax
movaps %xmm0, %xmm1
movsd 0x100(%rax), %xmm0
ucomisd %xmm1, %xmm0
seta %al
andb $0x1, %al
movb %al, -0x21(%rbp)
jmp 0xea24
movq -0x38(%rbp), %rax
movq -0x10(%rbp), %rcx
movsd 0x10(%rcx), %xmm0
movq -0x10(%rbp), %rcx
subsd 0x18(%rcx), %xmm0
movq -0x10(%rbp), %rcx
cvtsi2sdl 0x24(%rcx), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0x18(%rbp)
movsd -0x18(%rbp), %xmm1
movsd 0x100(%rax), %xmm0
ucomisd %xmm1, %xmm0
seta %al
andb $0x1, %al
movb %al, -0x21(%rbp)
movb -0x21(%rbp), %al
andb $0x1, %al
addq $0x40, %rsp
popq %rbp
retq
nop