// File: myudelson[P]hmm-scalable/HMMProblemSlicedA.cpp
// Function: HMMProblemSlicedA::checkPIABConstraints(double*, double***, double**)
bool HMMProblemSlicedA::checkPIABConstraints(NUMBER* a_PI, NUMBER*** a_A, NUMBER** a_B) {
	NPAR i, j, z;
    NPAR nS = this->p->nS, nO = this->p->nO, nZ = this->p->nZ;
	// check values
	NUMBER sum_pi = 0.0;
	NUMBER sum_a_row[nZ][nS];
	NUMBER sum_b_row[nS];
    for(i=0; i<nS; i++) {
        for(z=0; z<nZ; z++) {
            sum_a_row[z][i] = 0.0;
        }
        sum_b_row[i] = 0.0;
    }
	for(i=0; i<nS; i++) {
		if( a_PI[i]>1.0 || a_PI[i]<0.0)
			return false;
		sum_pi += a_PI[i];
    }
    
    for(z=0; z<nZ; z++) {
        for(i=0; i<nS; i++) {
            for(j=0; j<nS; j++) {
                if( a_A[z][i][j]>1.0 || a_A[z][i][j]<0.0)
                    return false;
                sum_a_row[z][i] += a_A[z][i][j];
            }// all states 2
        }// all states
    } // all slices

    for(i=0; i<nS; i++) {
        for(int m=0; m<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<nS; i++) {
        if( sum_b_row[i]!=1.0) return false;
        for(z=0; z<nZ; z++)
            if( sum_a_row[z][i]!=1.0 ) return false;
    }
	return true;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x80, %rsp
movq %rdi, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq %rcx, -0x28(%rbp)
movq -0x10(%rbp), %rax
movq 0x88(%rax), %rcx
movb 0x8b8(%rcx), %cl
movb %cl, -0x2c(%rbp)
movq 0x88(%rax), %rcx
movb 0x8b9(%rcx), %cl
movb %cl, -0x2d(%rbp)
movq 0x88(%rax), %rax
movb 0x8c8(%rax), %al
movb %al, -0x2e(%rbp)
movq $0x0, -0x38(%rbp)
movzbl -0x2e(%rbp), %esi
movl %esi, %ecx
movzbl -0x2c(%rbp), %edx
movl %edx, %eax
movq %rax, -0x78(%rbp)
movq %rsp, %rdi
movq %rdi, -0x40(%rbp)
imull %edx, %esi
movl %esi, %edx
leal 0xf(,%rdx,8), %edx
andl $-0x10, %edx
movl %edx, %esi
movq %rsp, %rdx
subq %rsi, %rdx
movq %rdx, -0x70(%rbp)
movq %rdx, %rsp
movq %rcx, -0x48(%rbp)
movq %rax, -0x50(%rbp)
movzbl -0x2c(%rbp), %edx
movl %edx, %eax
movl %edx, %ecx
leal 0xf(,%rcx,8), %ecx
andl $-0x10, %ecx
movl %ecx, %edx
movq %rsp, %rcx
subq %rdx, %rcx
movq %rcx, -0x68(%rbp)
movq %rcx, %rsp
movq %rax, -0x58(%rbp)
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x43043
movb $0x0, -0x2b(%rbp)
movsbl -0x2b(%rbp), %eax
movsbl -0x2e(%rbp), %ecx
cmpl %ecx, %eax
jge 0x43028
movq -0x70(%rbp), %rax
movq -0x78(%rbp), %rdx
movsbq -0x2b(%rbp), %rcx
imulq %rdx, %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movsbq -0x29(%rbp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movb -0x2b(%rbp), %al
addb $0x1, %al
movb %al, -0x2b(%rbp)
jmp 0x42fed
movq -0x68(%rbp), %rax
movsbq -0x29(%rbp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x42fdd
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x430b4
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0xc99f(%rip), %xmm1     # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x43082
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x43092
movb $0x0, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
jmp 0x43363
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
addsd -0x38(%rbp), %xmm0
movsd %xmm0, -0x38(%rbp)
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x43047
movb $0x0, -0x2b(%rbp)
movsbl -0x2b(%rbp), %eax
movsbl -0x2e(%rbp), %ecx
cmpl %ecx, %eax
jge 0x431c5
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x431b6
movb $0x0, -0x2a(%rbp)
movsbl -0x2a(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x431a7
movq -0x20(%rbp), %rax
movsbq -0x2b(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x2a(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0xc8f0(%rip), %xmm1     # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x43143
movq -0x20(%rbp), %rax
movsbq -0x2b(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x2a(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x43153
movb $0x0, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
jmp 0x43363
movq -0x70(%rbp), %rax
movq -0x78(%rbp), %rdx
movq -0x20(%rbp), %rcx
movsbq -0x2b(%rbp), %rsi
movq (%rcx,%rsi,8), %rcx
movsbq -0x29(%rbp), %rsi
movq (%rcx,%rsi,8), %rcx
movsbq -0x2a(%rbp), %rsi
movsd (%rcx,%rsi,8), %xmm0
movsbq -0x2b(%rbp), %rcx
imulq %rdx, %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movsbq -0x29(%rbp), %rcx
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
movb -0x2a(%rbp), %al
addb $0x1, %al
movb %al, -0x2a(%rbp)
jmp 0x430e0
jmp 0x431a9
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x430cc
jmp 0x431b8
movb -0x2b(%rbp), %al
addb $0x1, %al
movb %al, -0x2b(%rbp)
jmp 0x430b8
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x43284
movl $0x0, -0x60(%rbp)
movl -0x60(%rbp), %eax
movsbl -0x2d(%rbp), %ecx
cmpl %ecx, %eax
jge 0x43275
movq -0x28(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x60(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0xc7fb(%rip), %xmm1     # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x4322e
movq -0x28(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x60(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x4323e
movb $0x0, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
jmp 0x43363
movq -0x68(%rbp), %rax
movq -0x28(%rbp), %rcx
movsbq -0x29(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movslq -0x60(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movsbq -0x29(%rbp), %rcx
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
movl -0x60(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x60(%rbp)
jmp 0x431e0
jmp 0x43277
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x431c9
movsd -0x38(%rbp), %xmm0
movsd 0xc777(%rip), %xmm1     # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x4329b
jp 0x4329b
jmp 0x432ab
movb $0x0, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
jmp 0x43363
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x43358
movq -0x68(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0xc733(%rip), %xmm1     # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x432df
jp 0x432df
jmp 0x432ec
movb $0x0, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
jmp 0x43363
movb $0x0, -0x2b(%rbp)
movsbl -0x2b(%rbp), %eax
movsbl -0x2e(%rbp), %ecx
cmpl %ecx, %eax
jge 0x43349
movq -0x70(%rbp), %rax
movq -0x78(%rbp), %rdx
movsbq -0x2b(%rbp), %rcx
imulq %rdx, %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0xc6e2(%rip), %xmm1     # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x43330
jp 0x43330
jmp 0x4333d
movb $0x0, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
jmp 0x43363
jmp 0x4333f
movb -0x2b(%rbp), %al
addb $0x1, %al
movb %al, -0x2b(%rbp)
jmp 0x432f0
jmp 0x4334b
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x432af
movb $0x1, -0x1(%rbp)
movl $0x1, -0x5c(%rbp)
movq -0x40(%rbp), %rax
movq %rax, %rsp
movb -0x1(%rbp), %al
andb $0x1, %al
movq %rbp, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)