// File: myudelson[P]hmm-scalable/HMMProblemPiABGK.cpp
// Function: HMMProblemPiABGK::toFile(char const*)
void HMMProblemPiABGK::toFile(const char *filename) {
	FILE *fid = fopen(filename,"w");
	if(fid == NULL) {
		fprintf(stderr,"Can't write output model file %s\n",filename);
		exit(1);
	}
    
    // write solved id
    writeSolverInfo(fid, this->p);
    
	fprintf(fid,"Null skill ratios\t");
	for(NPAR m=0; m<this->p->nO; m++)
		fprintf(fid," %12.10f%s",this->null_obs_ratio[m],(m==(this->p->nO-1))?"\n":"\t");
	NCAT k, g;
    NPAR i,j,m;
	std::map<NCAT,std::string>::iterator it;
	for(g=0;g<this->p->nG; g++) {
		it = this->p->map_group_bwd->find(g);
		fprintf(fid,"%d\t%s\n",g,it->second.c_str());
		fprintf(fid,"PIg\t");
		for(i=0; i<this->p->nS; i++)
			fprintf(fid,"%12.10f%s",this->PIg[g][i],(i==(this->p->nS-1))?"\n":"\t");
		fprintf(fid,"Ag\t");
		for(i=0; i<this->p->nS; i++)
			for(j=0; j<this->p->nS; j++)
				fprintf(fid,"%12.10f%s",this->Ag[g][i][j],(i==(this->p->nS-1) && j==(this->p->nS-1))?"\n":"\t");
		fprintf(fid,"Bg\t");
		for(i=0; i<this->p->nS; i++)
			for(m=0; m<this->p->nO; m++)
				fprintf(fid,"%12.10f%s",this->Bg[g][i][m],(i==(this->p->nS-1) && m==(this->p->nO-1))?"\n":"\t");
    }
	for(k=0;k<this->p->nK;k++) {
		it = this->p->map_skill_bwd->find(k);
		fprintf(fid,"%d\t%s\n",k,it->second.c_str());
		fprintf(fid,"PIk\t");
		for(i=0; i<this->p->nS; i++)
			fprintf(fid,"%12.10f%s",this->pi[k][i],(i==(this->p->nS-1))?"\n":"\t");
		fprintf(fid,"Ak\t");
		for(i=0; i<this->p->nS; i++)
			for(j=0; j<this->p->nS; j++)
				fprintf(fid,"%12.10f%s",this->A[k][i][j],(i==(this->p->nS-1) && j==(this->p->nS-1))?"\n":"\t");
		fprintf(fid,"Bk\t");
		for(i=0; i<this->p->nS; i++)
			for(m=0; m<this->p->nO; m++)
				fprintf(fid,"%12.10f%s",this->B[k][i][m],(i==(this->p->nS-1) && m==(this->p->nO-1))?"\n":"\t");
	}
	fclose(fid);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x110, %rsp            # imm = 0x110
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x10(%rbp), %rdi
leaq 0x1a5f4(%rip), %rsi     # 0x5058a
callq 0x5370
movq %rax, -0x18(%rbp)
cmpq $0x0, -0x18(%rbp)
jne 0x35fcc
movq 0x2502b(%rip), %rax     # 0x5afd8
movq (%rax), %rdi
movq -0x10(%rbp), %rdx
leaq 0x1a5d1(%rip), %rsi     # 0x5058c
movb $0x0, %al
callq 0x5440
movl $0x1, %edi
callq 0x5390
movq -0x48(%rbp), %rax
movq -0x18(%rbp), %rdi
movq 0x88(%rax), %rsi
callq 0xbc70
movq -0x18(%rbp), %rdi
leaq 0x1a5c3(%rip), %rsi     # 0x505ae
movb $0x0, %al
callq 0x5440
movb $0x0, -0x19(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x19(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
cmpl %ecx, %eax
jge 0x3608a
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0x58(%rbp)
movq 0x18(%rcx), %rax
movsbq -0x19(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0x50(%rbp)
movsbl -0x19(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x36055
leaq 0x19956(%rip), %rax     # 0x4f9a5
movq %rax, -0x60(%rbp)
jmp 0x36062
leaq 0x1ac1e(%rip), %rax     # 0x50c7a
movq %rax, -0x60(%rbp)
jmp 0x36062
movsd -0x50(%rbp), %xmm0
movq -0x58(%rbp), %rdi
movq -0x60(%rbp), %rdx
leaq 0x1a8ec(%rip), %rsi     # 0x50962
movb $0x1, %al
callq 0x5440
movb -0x19(%rbp), %al
addb $0x1, %al
movb %al, -0x19(%rbp)
jmp 0x35ff6
leaq -0x30(%rbp), %rdi
callq 0x170b0
movl $0x0, -0x24(%rbp)
movq -0x48(%rbp), %rcx
movl -0x24(%rbp), %eax
movq 0x88(%rcx), %rcx
cmpl 0x8bc(%rcx), %eax
jge 0x36432
movq -0x48(%rbp), %rax
movq 0x88(%rax), %rax
movq 0x928(%rax), %rdi
leaq -0x24(%rbp), %rsi
callq 0x27d20
movq %rax, -0x38(%rbp)
movq -0x38(%rbp), %rax
movq %rax, -0x30(%rbp)
movq -0x18(%rbp), %rax
movq %rax, -0x70(%rbp)
movl -0x24(%rbp), %eax
movl %eax, -0x64(%rbp)
leaq -0x30(%rbp), %rdi
callq 0x17160
movq %rax, %rdi
addq $0x8, %rdi
callq 0x5110
movq -0x70(%rbp), %rdi
movl -0x64(%rbp), %edx
movq %rax, %rcx
leaq 0x1a4bc(%rip), %rsi     # 0x505cb
movb $0x0, %al
callq 0x5440
movq -0x18(%rbp), %rdi
leaq 0x1a84c(%rip), %rsi     # 0x5096d
movb $0x0, %al
callq 0x5440
movb $0x0, -0x25(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x361d8
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0x80(%rbp)
movq 0x90(%rcx), %rax
movslq -0x24(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x25(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0x78(%rbp)
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x3619d
leaq 0x19811(%rip), %rax     # 0x4f9a5
movq %rax, -0x88(%rbp)
jmp 0x361ad
leaq 0x1aad6(%rip), %rax     # 0x50c7a
movq %rax, -0x88(%rbp)
jmp 0x361ad
movsd -0x78(%rbp), %xmm0
movq -0x80(%rbp), %rdi
movq -0x88(%rbp), %rdx
leaq 0x1a79f(%rip), %rsi     # 0x50963
movb $0x1, %al
callq 0x5440
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0x3612c
movq -0x18(%rbp), %rdi
leaq 0x1a8f9(%rip), %rsi     # 0x50adc
movb $0x0, %al
callq 0x5440
movb $0x0, -0x25(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x362fd
movb $0x0, -0x26(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x26(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x362ee
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0x98(%rbp)
movq 0x98(%rcx), %rax
movslq -0x24(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x25(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x26(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0x90(%rbp)
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x362ad
movq -0x48(%rbp), %rcx
movsbl -0x26(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x362ad
leaq 0x19701(%rip), %rax     # 0x4f9a5
movq %rax, -0xa0(%rbp)
jmp 0x362bd
leaq 0x1a9c6(%rip), %rax     # 0x50c7a
movq %rax, -0xa0(%rbp)
jmp 0x362bd
movsd -0x90(%rbp), %xmm0
movq -0x98(%rbp), %rdi
movq -0xa0(%rbp), %rdx
leaq 0x1a689(%rip), %rsi     # 0x50963
movb $0x1, %al
callq 0x5440
movb -0x26(%rbp), %al
addb $0x1, %al
movb %al, -0x26(%rbp)
jmp 0x36210
jmp 0x362f0
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0x361ee
movq -0x18(%rbp), %rdi
leaq 0x1a922(%rip), %rsi     # 0x50c2a
movb $0x0, %al
callq 0x5440
movb $0x0, -0x25(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x36422
movb $0x0, -0x27(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x27(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
cmpl %ecx, %eax
jge 0x36413
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0xb0(%rbp)
movq 0xa0(%rcx), %rax
movslq -0x24(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x25(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x27(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0xa8(%rbp)
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x363d2
movq -0x48(%rbp), %rcx
movsbl -0x27(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x363d2
leaq 0x195dc(%rip), %rax     # 0x4f9a5
movq %rax, -0xb8(%rbp)
jmp 0x363e2
leaq 0x1a8a1(%rip), %rax     # 0x50c7a
movq %rax, -0xb8(%rbp)
jmp 0x363e2
movsd -0xa8(%rbp), %xmm0
movq -0xb0(%rbp), %rdi
movq -0xb8(%rbp), %rdx
leaq 0x1a564(%rip), %rsi     # 0x50963
movb $0x1, %al
callq 0x5440
movb -0x27(%rbp), %al
addb $0x1, %al
movb %al, -0x27(%rbp)
jmp 0x36335
jmp 0x36415
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0x36313
jmp 0x36424
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
jmp 0x3609a
movl $0x0, -0x20(%rbp)
movq -0x48(%rbp), %rcx
movl -0x20(%rbp), %eax
movq 0x88(%rcx), %rcx
cmpl 0x8c0(%rcx), %eax
jge 0x367e0
movq -0x48(%rbp), %rax
movq 0x88(%rax), %rax
movq 0x948(%rax), %rdi
leaq -0x20(%rbp), %rsi
callq 0x27d20
movq %rax, -0x40(%rbp)
movq -0x40(%rbp), %rax
movq %rax, -0x30(%rbp)
movq -0x18(%rbp), %rax
movq %rax, -0xc8(%rbp)
movl -0x20(%rbp), %eax
movl %eax, -0xbc(%rbp)
leaq -0x30(%rbp), %rdi
callq 0x17160
movq %rax, %rdi
addq $0x8, %rdi
callq 0x5110
movq -0xc8(%rbp), %rdi
movl -0xbc(%rbp), %edx
movq %rax, %rcx
leaq 0x1a111(%rip), %rsi     # 0x505cb
movb $0x0, %al
callq 0x5440
movq -0x18(%rbp), %rdi
leaq 0x1a4a6(%rip), %rsi     # 0x50972
movb $0x0, %al
callq 0x5440
movb $0x0, -0x25(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x3658c
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0xd8(%rbp)
movq 0x38(%rcx), %rax
movslq -0x20(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x25(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0xd0(%rbp)
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x3654b
leaq 0x19463(%rip), %rax     # 0x4f9a5
movq %rax, -0xe0(%rbp)
jmp 0x3655b
leaq 0x1a728(%rip), %rax     # 0x50c7a
movq %rax, -0xe0(%rbp)
jmp 0x3655b
movsd -0xd0(%rbp), %xmm0
movq -0xd8(%rbp), %rdi
movq -0xe0(%rbp), %rdx
leaq 0x1a3eb(%rip), %rsi     # 0x50963
movb $0x1, %al
callq 0x5440
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0x364d7
movq -0x18(%rbp), %rdi
leaq 0x1a549(%rip), %rsi     # 0x50ae0
movb $0x0, %al
callq 0x5440
movb $0x0, -0x25(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x366ae
movb $0x0, -0x26(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x26(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x3669f
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0xf0(%rbp)
movq 0x40(%rcx), %rax
movslq -0x20(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x25(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x26(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0xe8(%rbp)
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x3665e
movq -0x48(%rbp), %rcx
movsbl -0x26(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x3665e
leaq 0x19350(%rip), %rax     # 0x4f9a5
movq %rax, -0xf8(%rbp)
jmp 0x3666e
leaq 0x1a615(%rip), %rax     # 0x50c7a
movq %rax, -0xf8(%rbp)
jmp 0x3666e
movsd -0xe8(%rbp), %xmm0
movq -0xf0(%rbp), %rdi
movq -0xf8(%rbp), %rdx
leaq 0x1a2d8(%rip), %rsi     # 0x50963
movb $0x1, %al
callq 0x5440
movb -0x26(%rbp), %al
addb $0x1, %al
movb %al, -0x26(%rbp)
jmp 0x365c4
jmp 0x366a1
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0x365a2
movq -0x18(%rbp), %rdi
leaq 0x1a575(%rip), %rsi     # 0x50c2e
movb $0x0, %al
callq 0x5440
movb $0x0, -0x25(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
cmpl %ecx, %eax
jge 0x367d0
movb $0x0, -0x27(%rbp)
movq -0x48(%rbp), %rcx
movsbl -0x27(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
cmpl %ecx, %eax
jge 0x367c1
movq -0x48(%rbp), %rcx
movq -0x18(%rbp), %rax
movq %rax, -0x108(%rbp)
movq 0x48(%rcx), %rax
movslq -0x20(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x25(%rbp), %rdx
movq (%rax,%rdx,8), %rax
movsbq -0x27(%rbp), %rdx
movsd (%rax,%rdx,8), %xmm0
movsd %xmm0, -0x100(%rbp)
movsbl -0x25(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b8(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x36780
movq -0x48(%rbp), %rcx
movsbl -0x27(%rbp), %eax
movq 0x88(%rcx), %rcx
movsbl 0x8b9(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jne 0x36780
leaq 0x1922e(%rip), %rax     # 0x4f9a5
movq %rax, -0x110(%rbp)
jmp 0x36790
leaq 0x1a4f3(%rip), %rax     # 0x50c7a
movq %rax, -0x110(%rbp)
jmp 0x36790
movsd -0x100(%rbp), %xmm0
movq -0x108(%rbp), %rdi
movq -0x110(%rbp), %rdx
leaq 0x1a1b6(%rip), %rsi     # 0x50963
movb $0x1, %al
callq 0x5440
movb -0x27(%rbp), %al
addb $0x1, %al
movb %al, -0x27(%rbp)
jmp 0x366e6
jmp 0x367c3
movb -0x25(%rbp), %al
addb $0x1, %al
movb %al, -0x25(%rbp)
jmp 0x366c4
jmp 0x367d2
movl -0x20(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x20(%rbp)
jmp 0x36439
movq -0x18(%rbp), %rdi
callq 0x52c0
addq $0x110, %rsp            # imm = 0x110
popq %rbp
retq
nopw %cs:(%rax,%rax)