// File: 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
subq $0x70, %rsp
movq %rdi, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq %rcx, -0x28(%rbp)
movq -0x10(%rbp), %rax
movq %rax, -0x70(%rbp)
movq $0x0, -0x38(%rbp)
movq 0x88(%rax), %rcx
movzbl 0x8b8(%rcx), %esi
movl %esi, %ecx
movq %rsp, %rdx
movq %rdx, -0x40(%rbp)
movl %esi, %edx
leal 0xf(,%rdx,8), %edx
andl $-0x10, %edx
movl %edx, %esi
movq %rsp, %rdx
subq %rsi, %rdx
movq %rdx, -0x68(%rbp)
movq %rdx, %rsp
movq %rcx, -0x48(%rbp)
movq 0x88(%rax), %rax
movzbl 0x8b8(%rax), %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, -0x60(%rbp)
movq %rcx, %rsp
movq %rax, -0x50(%rbp)
movb $0x0, -0x29(%rbp)
movq -0x70(%rbp), %rcx
movsbl -0x29(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x19f57
movq -0x60(%rbp), %rax
movq -0x68(%rbp), %rcx
movsbq -0x29(%rbp), %rdx
xorps %xmm0, %xmm0
movsd %xmm0, (%rcx,%rdx,8)
movsbq -0x29(%rbp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x19f11
movb $0x0, -0x29(%rbp)
movq -0x70(%rbp), %rcx
movsbl -0x29(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x1a133
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x35a79(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x19fa8
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x19fb8
movb $0x0, -0x1(%rbp)
movl $0x1, -0x54(%rbp)
jmp 0x1a1da
movq -0x18(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
addsd -0x38(%rbp), %xmm0
movsd %xmm0, -0x38(%rbp)
movb $0x0, -0x2a(%rbp)
movq -0x70(%rbp), %rcx
movsbl -0x2a(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x1a07a
movq -0x20(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x2a(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x359f7(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x1a033
movq -0x20(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movsbq -0x2a(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x1a043
movb $0x0, -0x1(%rbp)
movl $0x1, -0x54(%rbp)
jmp 0x1a1da
movq -0x68(%rbp), %rax
movq -0x20(%rbp), %rcx
movsbq -0x29(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movsbq -0x2a(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
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 0x19fd4
movl $0x0, -0x58(%rbp)
movq -0x70(%rbp), %rcx
movl -0x58(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
cmpl %ecx, %eax
jge 0x1a124
movq -0x28(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x58(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x3594c(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
ja 0x1a0dd
movq -0x28(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movslq -0x58(%rbp), %rcx
xorps %xmm0, %xmm0
ucomisd (%rax,%rcx,8), %xmm0
jbe 0x1a0ed
movb $0x0, -0x1(%rbp)
movl $0x1, -0x54(%rbp)
jmp 0x1a1da
movq -0x60(%rbp), %rax
movq -0x28(%rbp), %rcx
movsbq -0x29(%rbp), %rdx
movq (%rcx,%rdx,8), %rcx
movslq -0x58(%rbp), %rdx
movsd (%rcx,%rdx,8), %xmm0
movsbq -0x29(%rbp), %rcx
addsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rax,%rcx,8)
movl -0x58(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x58(%rbp)
jmp 0x1a081
jmp 0x1a126
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x19f5b
movsd -0x38(%rbp), %xmm0
movsd 0x358c8(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x1a14a
jp 0x1a14a
jmp 0x1a15a
movb $0x0, -0x1(%rbp)
movl $0x1, -0x54(%rbp)
jmp 0x1a1da
movb $0x0, -0x29(%rbp)
movq -0x70(%rbp), %rcx
movsbl -0x29(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x1a1cf
movq -0x68(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x3587a(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x1a1b6
jp 0x1a1b6
movq -0x60(%rbp), %rax
movsbq -0x29(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd 0x3585c(%rip), %xmm1    # 0x4fa08
ucomisd %xmm1, %xmm0
jne 0x1a1b6
jp 0x1a1b6
jmp 0x1a1c3
movb $0x0, -0x1(%rbp)
movl $0x1, -0x54(%rbp)
jmp 0x1a1da
jmp 0x1a1c5
movb -0x29(%rbp), %al
addb $0x1, %al
movb %al, -0x29(%rbp)
jmp 0x1a15e
movb $0x1, -0x1(%rbp)
movl $0x1, -0x54(%rbp)
movq -0x40(%rbp), %rax
movq %rax, %rsp
movb -0x1(%rbp), %al
andb $0x1, %al
movq %rbp, %rsp
popq %rbp
retq
nopl (%rax,%rax)