// File: myudelson[P]hmm-scalable/HMMProblemSlicedAB.cpp
// Function: HMMProblemSlicedAB::checkPIABConstraints(double*, double***, double***)
bool HMMProblemSlicedAB::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[nZ][nS];
    for(z=0; z<nZ; z++) {
        for(i=0; i<nS; i++) {
            sum_a_row[z][i] = 0.0;
            sum_b_row[z][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
            for(int m=0; m<nO; m++) {
                if( a_B[z][i][m]>1.0 || a_B[z][i][m]<0.0)
                    return false;
                sum_b_row[z][i] += a_B[z][i][m];
            }// all observations
        }// all states
    } // all slices
	if(sum_pi!=1.0)
		return false;
    for(z=0; z<nZ; z++) {
        for(i=0; i<nS; i++) {
            if( sum_a_row[z][i]!=1.0 || sum_b_row[z][i]!=1.0) {
                return false;
            }
        }
    }
	return true;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x90, %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, -0x88(%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, -0x80(%rbp)
movq %rdx, %rsp
movq %rcx, -0x48(%rbp)
movq %rax, -0x50(%rbp)
movzbl -0x2e(%rbp), %esi
movl %esi, %ecx
movzbl -0x2c(%rbp), %edx
movl %edx, %eax
movq %rax, -0x78(%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, -0x58(%rbp)
movq %rax, -0x60(%rbp)
movb $0x0, -0x2b(%rbp)
movsbl -0x2b(%rbp), %eax
movsbl -0x2e(%rbp), %ecx
cmpl %ecx, %eax
jge 0x394d0
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x394c4
movq -0x70(%rbp), %rax
movq -0x78(%rbp), %rdx
movq -0x80(%rbp), %rcx
movq -0x88(%rbp), %rdi
movsbq -0x2b(%rbp), %rsi
imulq %rdi, %rsi
shlq $0x3, %rsi
addq %rsi, %rcx
movsbq -0x29(%rbp), %rsi
xorps %xmm0, %xmm0
movsd %xmm0, (%rcx,%rsi,8)
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 -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x39461
jmp 0x394c6
movb -0x2b(%rbp), %al
addb $0x1, %al
movb %al, -0x2b(%rbp)
jmp 0x39451
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x39541
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x16512(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x3950f
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x3951f
movb $0x0, -0x1(%rbp)
movl $0x1, -0x64(%rbp)
jmp 0x3980e
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 0x394d4
movb $0x0, -0x2b(%rbp)
movsbl -0x2b(%rbp), %eax
movsbl -0x2e(%rbp), %ecx
cmpl %ecx, %eax
jge 0x39720
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x39711
movb $0x0, -0x2a(%rbp)
movsbl -0x2a(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x39637
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 0x16463(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x395d0
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 0x395e0
movb $0x0, -0x1(%rbp)
movl $0x1, -0x64(%rbp)
jmp 0x3980e
movq -0x80(%rbp), %rax
movq -0x88(%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 0x3956d
movl $0x0, -0x68(%rbp)
movl -0x68(%rbp), %eax
movsbl -0x2d(%rbp), %ecx
cmpl %ecx, %eax
jge 0x39702
movq -0x28(%rbp), %rax
movsbq -0x2b(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x68(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x16394(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x3969e
movq -0x28(%rbp), %rax
movsbq -0x2b(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x68(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x396ae
movb $0x0, -0x1(%rbp)
movl $0x1, -0x64(%rbp)
jmp 0x3980e
movq -0x70(%rbp), %rax
movq -0x78(%rbp), %rdx
movq -0x28(%rbp), %rcx
movsbq -0x2b(%rbp), %rsi
movq (%rcx,%rsi,8), %rcx
movsbq -0x29(%rbp), %rsi
movq (%rcx,%rsi,8), %rcx
movslq -0x68(%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)
movl -0x68(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x68(%rbp)
jmp 0x3963e
jmp 0x39704
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x39559
jmp 0x39713
movb -0x2b(%rbp), %al
addb $0x1, %al
movb %al, -0x2b(%rbp)
jmp 0x39545
movsd -0x38(%rbp), %xmm0
movsd 0x162db(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x39737
jp 0x39737
jmp 0x39747
movb $0x0, -0x1(%rbp)
movl $0x1, -0x64(%rbp)
jmp 0x3980e
movb $0x0, -0x2b(%rbp)
movsbl -0x2b(%rbp), %eax
movsbl -0x2e(%rbp), %ecx
cmpl %ecx, %eax
jge 0x39803
movb $0x0, -0x29(%rbp)
movsbl -0x29(%rbp), %eax
movsbl -0x2c(%rbp), %ecx
cmpl %ecx, %eax
jge 0x397f4
movq -0x80(%rbp), %rax
movq -0x88(%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 0x1626c(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x397d8
jp 0x397d8
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 0x1623a(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x397d8
jp 0x397d8
jmp 0x397e5
movb $0x0, -0x1(%rbp)
movl $0x1, -0x64(%rbp)
jmp 0x3980e
jmp 0x397e7
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x3975f
jmp 0x397f6
movb -0x2b(%rbp), %al
addb $0x1, %al
movb %al, -0x2b(%rbp)
jmp 0x3974b
movb $0x1, -0x1(%rbp)
movl $0x1, -0x64(%rbp)
movq -0x40(%rbp), %rax
movq %rax, %rsp
movb -0x1(%rbp), %al
andb $0x1, %al
movq %rbp, %rsp
popq %rbp
retq
nop