// File: /workspace/github_cpp_makefile_star2_O1/bobverity[P]MavericK/writeOut.cpp
// Function: printQmatrix_gene(globals&, int)
void printQmatrix_gene(globals &globals, int Kindex) {
    int K = globals.Kmin+Kindex;
    
    // open outputQmatrix_gene filestream for this K
    string s = globals.outputQmatrix_gene_filePath;
    s.erase(end(s)-4,end(s));
    s = s + "_K" + to_string((long long)K) + ".csv";
    globals.outputQmatrix_gene_fileStream = safe_ofstream(s,globals.outputLog_on, globals.outputLog_fileStream);
    
    // write headers
    globals.outputQmatrix_gene_fileStream << "index,label";
    if (globals.popCol_on)
        globals.outputQmatrix_gene_fileStream << ",given_population";
    globals.outputQmatrix_gene_fileStream << ",locus,gene_copy";
    for (int k=0; k<K; k++)
        globals.outputQmatrix_gene_fileStream << ",deme" << k+1;
    globals.outputQmatrix_gene_fileStream << "\n";
    
    // write results
    int groupIndex=-1;
    for (int ind=0; ind<globals.n; ind++) {
        for (int l=0; l<globals.loci; l++) {
            for (int p=0; p<globals.ploidy_vec[ind]; p++) {
                groupIndex++;
                
                // individual index
                globals.outputQmatrix_gene_fileStream << ind+1 << ",";
                
                // individual label
                globals.outputQmatrix_gene_fileStream << globals.indLabels_vec[ind] << ",";
                
                // population data
                if (globals.popCol_on)
                    globals.outputQmatrix_gene_fileStream << globals.pop_vec[ind] << ",";
                
                // locus
                globals.outputQmatrix_gene_fileStream << l+1 << ",";
                
                // gene copy
                globals.outputQmatrix_gene_fileStream << p+1 << ",";
                
                // probabilities
                for (int k=0; k<K; k++) {
                    char * buffer = new char[256];
                    sprintf(buffer, "%.3f", my_round(globals.Qmatrix_gene[Kindex][groupIndex][k]*1000)/1000.0);
                    string s = buffer;
                    delete [] buffer;
                    globals.outputQmatrix_gene_fileStream << s;
                    if (k<(K-1))
                        globals.outputQmatrix_gene_fileStream << ",";
                }
                
                globals.outputQmatrix_gene_fileStream << "\n";
                
            }
        }
    }
    globals.outputQmatrix_gene_fileStream.flush();
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x2e8, %rsp            # imm = 0x2E8
movl 0x2bc8(%rdi), %eax
movl %esi, 0x30(%rsp)
addl %esi, %eax
movq %rax, 0x18(%rsp)
leaq 0xb8(%rsp), %rax
movq %rax, -0x10(%rax)
movq 0x3a0(%rdi), %rsi
movq %rdi, 0x8(%rsp)
movq 0x3a8(%rdi), %rdx
addq %rsi, %rdx
leaq 0xa8(%rsp), %r14
movq %r14, %rdi
callq 0xbe10
movq (%r14), %rax
movq 0x8(%r14), %rcx
leaq -0x4(%rcx), %rdx
movq %rdx, 0x8(%r14)
movb $0x0, -0x4(%rax,%rcx)
movq (%r14), %rsi
leaq 0x98(%rsp), %r13
movq %r13, -0x10(%r13)
movq 0x8(%r14), %rdx
addq %rsi, %rdx
leaq 0x88(%rsp), %rdi
callq 0xbe10
leaq 0x6e9e(%rip), %rsi      # 0x3987e
leaq 0x88(%rsp), %rdi
callq 0x35e0
movslq 0x18(%rsp), %r15
movq %r15, %r14
negq %r14
cmovsq %r15, %r14
movl $0x1, %r12d
cmpq $0xa, %r14
jb 0x32a5c
movl $0x4, %r12d
movabsq $0x346dc5d63886594b, %rsi # imm = 0x346DC5D63886594B
movq %r14, %rcx
cmpq $0x63, %rcx
jbe 0x32a53
cmpq $0x3e7, %rcx            # imm = 0x3E7
jbe 0x32a59
cmpq $0x2710, %rcx           # imm = 0x2710
jb 0x32a5c
movq %rcx, %rax
mulq %rsi
shrq $0xb, %rdx
addl $0x4, %r12d
cmpq $0x1869f, %rcx          # imm = 0x1869F
movq %rdx, %rcx
ja 0x32a1b
addl $-0x3, %r12d
jmp 0x32a5c
addl $-0x2, %r12d
jmp 0x32a5c
decl %r12d
shrq $0x3f, %r15
leal (%r15,%r12), %esi
leaq 0x78(%rsp), %rbp
movq %rbp, -0x10(%rbp)
leaq 0x68(%rsp), %rdi
movl $0x2d, %edx
callq 0x33a0
addq 0x68(%rsp), %r15
movq %r15, %rdi
movl %r12d, %esi
movq %r14, %rdx
callq 0x6dda
movq 0x88(%rsp), %rcx
movq 0x90(%rsp), %r8
movq 0x70(%rsp), %rdx
leaq (%rdx,%r8), %rax
movl $0xf, %esi
cmpq %r13, %rcx
je 0x32aba
movq 0x98(%rsp), %rsi
cmpq %rsi, %rax
jbe 0x32ad5
movl $0xf, %esi
cmpq %rbp, 0x68(%rsp)
je 0x32ad0
movq 0x78(%rsp), %rsi
cmpq %rsi, %rax
jbe 0x32ae9
movq 0x68(%rsp), %rsi
leaq 0x88(%rsp), %rdi
callq 0x30e0
jmp 0x32af7
leaq 0x68(%rsp), %rdi
xorl %esi, %esi
xorl %edx, %edx
callq 0x3370
leaq 0x58(%rsp), %rdx
movq %rdx, -0x10(%rdx)
movq (%rax), %rsi
leaq 0x10(%rax), %rcx
cmpq %rcx, %rsi
je 0x32b1b
movq %rsi, 0x48(%rsp)
movq (%rcx), %rdx
movq %rdx, 0x58(%rsp)
jmp 0x32b23
movupd (%rcx), %xmm0
movupd %xmm0, (%rdx)
movq %rax, %rdx
addq $0x8, %rdx
movq 0x8(%rax), %rsi
movq %rsi, 0x50(%rsp)
movq %rcx, (%rax)
movq $0x0, (%rdx)
movb $0x0, (%rcx)
leaq 0x489e(%rip), %rsi      # 0x373e5
leaq 0x48(%rsp), %rdi
callq 0x35e0
leaq 0xf8(%rsp), %rbx
movq %rbx, -0x10(%rbx)
movq (%rax), %rdx
movq %rax, %rcx
addq $0x10, %rcx
cmpq %rcx, %rdx
je 0x32b81
movq %rdx, 0xe8(%rsp)
movq (%rcx), %rdx
movq %rdx, 0xf8(%rsp)
jmp 0x32b89
movupd (%rcx), %xmm0
movupd %xmm0, (%rbx)
movq 0x8(%rax), %rdx
leaq 0xe8(%rsp), %r14
movq %rdx, 0x8(%r14)
movq %rcx, (%rax)
movq $0x0, 0x8(%rax)
movb $0x0, 0x10(%rax)
leaq 0xa8(%rsp), %rdi
movq %r14, %rsi
callq 0x33e0
movq (%r14), %rdi
cmpq %rbx, %rdi
je 0x32bc5
callq 0x3310
leaq 0x58(%rsp), %rax
movq -0x10(%rax), %rdi
cmpq %rax, %rdi
je 0x32bd8
callq 0x3310
movq 0x68(%rsp), %rdi
cmpq %rbp, %rdi
je 0x32be7
callq 0x3310
movq 0x88(%rsp), %rdi
cmpq %r13, %rdi
je 0x32bf9
callq 0x3310
leaq 0xd8(%rsp), %rbx
movq %rbx, -0x10(%rbx)
movq 0xa8(%rsp), %rsi
movq 0xb0(%rsp), %rdx
addq %rsi, %rdx
leaq 0xc8(%rsp), %rdi
callq 0xbe10
movq 0x8(%rsp), %r13
movzbl 0x2c3c(%r13), %edx
leaq 0x970(%r13), %rcx
leaq 0xe8(%rsp), %r12
leaq 0xc8(%rsp), %rsi
movq %r12, %rdi
callq 0x2f394
leaq 0x1170(%r13), %r14
movq 0x1170(%r13), %rax
movq -0x18(%rax), %rdi
addq %r14, %rdi
leaq 0xf0(%rsp), %r15
movq -0x8(%r15), %rax
addq -0x18(%rax), %r12
movq %r12, %rsi
callq 0x35d0
leaq 0x1178(%r13), %rdi
movq %r15, %rsi
callq 0x3340
movq 0xf30c(%rip), %rax      # 0x41fa0
movq (%rax), %rcx
movq 0x18(%rax), %rax
leaq 0x1e0(%rsp), %r12
movq %rcx, -0xf8(%r12)
movq -0x18(%rcx), %rcx
movq %rax, 0xe8(%rsp,%rcx)
movq %r15, %rdi
callq 0x32b0
movq %r12, %rdi
callq 0x30f0
movq 0xc8(%rsp), %rdi
cmpq %rbx, %rdi
je 0x32cd9
callq 0x3310
leaq 0x6d56(%rip), %rsi      # 0x39a36
movl $0xb, %edx
movq %r14, %rdi
callq 0x3390
cmpb $0x1, 0x2ba1(%r13)
jne 0x32d0b
leaq 0x6d44(%rip), %rsi      # 0x39a42
movl $0x11, %edx
movq %r14, %rdi
callq 0x3390
leaq 0x6d42(%rip), %rsi      # 0x39a54
movl $0x10, %edx
movq %r14, %rdi
callq 0x3390
xorl %ebx, %ebx
movq 0x18(%rsp), %rax
testl %eax, %eax
cmovgl %eax, %ebx
negl %ebx
movl $0x1, %r12d
leaq 0x6d2b(%rip), %r15      # 0x39a65
leal (%rbx,%r12), %eax
cmpl $0x1, %eax
je 0x32d68
movl $0x5, %edx
movq %r14, %rdi
movq %r15, %rsi
callq 0x3390
leal 0x1(%r12), %ebp
movq %r14, %rdi
movl %r12d, %esi
callq 0x3550
movl %ebp, %r12d
jmp 0x32d3a
leaq 0x4bb3(%rip), %rsi      # 0x37922
movl $0x1, %edx
movq %r14, %rdi
callq 0x3390
cmpl $0x0, 0x2d18(%r13)
jle 0x33052
movq %r13, %rdx
movslq 0x30(%rsp), %rax
movq 0x18(%rsp), %rsi
leal -0x1(%rsi), %ecx
movslq %ecx, %rcx
movq %rcx, 0x30(%rsp)
movl %esi, %r12d
movl $0xffffffff, %ecx       # imm = 0xFFFFFFFF
movq %rcx, 0x28(%rsp)
movq $0x0, 0x20(%rsp)
leaq (%rax,%rax,2), %r13
cmpl $0x0, 0x2d1c(%rdx)
jle 0x33030
movq 0x20(%rsp), %rax
leal 0x1(%rax), %ecx
movl %ecx, 0x3c(%rsp)
movq %rax, %rcx
xorl %eax, %eax
shlq $0x5, %rcx
movq %rcx, 0x40(%rsp)
incl %eax
movl %eax, 0x14(%rsp)
movq 0x28(%rsp), %rax
cltq
movq %rax, 0x28(%rsp)
movl $0x0, 0x10(%rsp)
movq 0x8(%rsp), %rax
movq 0x2cd0(%rax), %rax
movq 0x20(%rsp), %rcx
movl 0x10(%rsp), %edx
cmpl (%rax,%rcx,4), %edx
jge 0x33019
movq %r14, %rdi
movl 0x3c(%rsp), %esi
callq 0x3550
movl $0x1, %edx
movq %rax, %rdi
leaq 0x445a(%rip), %rsi      # 0x37290
callq 0x3390
movq 0x8(%rsp), %rax
movq 0x2c58(%rax), %rax
movq 0x40(%rsp), %rcx
movq (%rax,%rcx), %rsi
movq 0x8(%rax,%rcx), %rdx
movq %r14, %rdi
callq 0x3390
movl $0x1, %edx
movq %rax, %rdi
leaq 0x4424(%rip), %rsi      # 0x37290
callq 0x3390
movq 0x8(%rsp), %rax
cmpb $0x1, 0x2ba1(%rax)
jne 0x32eb5
movq 0x8(%rsp), %rax
movq 0x2c70(%rax), %rax
movq 0x40(%rsp), %rcx
movq (%rax,%rcx), %rsi
movq 0x8(%rax,%rcx), %rdx
movq %r14, %rdi
callq 0x3390
movl $0x1, %edx
movq %rax, %rdi
leaq 0x43e0(%rip), %rsi      # 0x37290
callq 0x3390
movq %r14, %rdi
movl 0x14(%rsp), %esi
callq 0x3550
movl $0x1, %edx
movq %rax, %rdi
leaq 0x43c0(%rip), %rsi      # 0x37290
callq 0x3390
movl 0x10(%rsp), %esi
incl %esi
movq %r14, %rdi
movl %esi, 0x10(%rsp)
callq 0x3550
movl $0x1, %edx
movq %rax, %rdi
leaq 0x439a(%rip), %rsi      # 0x37290
callq 0x3390
incq 0x28(%rsp)
cmpl $0x0, 0x18(%rsp)
jle 0x33000
movq 0x28(%rsp), %rax
leaq (%rax,%rax,2), %rbx
xorl %ebp, %ebp
movl $0x100, %edi            # imm = 0x100
callq 0x3030
movsd 0x67e0(%rip), %xmm1     # 0x39708
movq %rax, %r15
movq 0x8(%rsp), %rax
movq 0x2d70(%rax), %rax
movq (%rax,%r13,8), %rax
movq (%rax,%rbx,8), %rax
movsd (%rax,%rbp,8), %xmm0
mulsd %xmm1, %xmm0
callq 0x1aa2c
movsd 0x67b3(%rip), %xmm1     # 0x39708
xorps %xmm0, %xmm0
cvtsi2sd %eax, %xmm0
divsd %xmm1, %xmm0
movq %r15, %rdi
leaq 0x6b01(%rip), %rsi      # 0x39a6b
movb $0x1, %al
callq 0x3080
leaq 0xf8(%rsp), %rax
movq %rax, 0xe8(%rsp)
movq %r15, %rdi
callq 0x3140
movq %r15, %rdx
addq %rax, %rdx
leaq 0xe8(%rsp), %rdi
movq %r15, %rsi
callq 0x6e56
movq %r15, %rdi
callq 0x33f0
movq 0xe8(%rsp), %rsi
movq 0xf0(%rsp), %rdx
movq %r14, %rdi
callq 0x3390
cmpq 0x30(%rsp), %rbp
jge 0x32fda
movl $0x1, %edx
movq %r14, %rdi
leaq 0x42bb(%rip), %rsi      # 0x37290
callq 0x3390
movq 0xe8(%rsp), %rdi
leaq 0xf8(%rsp), %rax
cmpq %rax, %rdi
je 0x32ff4
callq 0x3310
incq %rbp
cmpq %rbp, %r12
jne 0x32f16
movl $0x1, %edx
movq %r14, %rdi
leaq 0x4913(%rip), %rsi      # 0x37922
callq 0x3390
jmp 0x32dfd
movq 0x8(%rsp), %rax
movl 0x14(%rsp), %ecx
cmpl 0x2d1c(%rax), %ecx
movl %ecx, %eax
jl 0x32de3
movq 0x20(%rsp), %rsi
incq %rsi
movq 0x8(%rsp), %rdx
movslq 0x2d18(%rdx), %rax
movq %rsi, 0x20(%rsp)
cmpq %rax, %rsi
jl 0x32dbc
movq %r14, %rdi
callq 0x31f0
movq 0xa8(%rsp), %rdi
leaq 0xb8(%rsp), %rax
cmpq %rax, %rdi
je 0x33074
callq 0x3310
addq $0x2e8, %rsp            # imm = 0x2E8
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0x3311a
movq %rax, %r14
leaq 0xe8(%rsp), %rdi
callq 0x3440
jmp 0x330a0
movq %rax, %r14
movq 0xc8(%rsp), %rdi
cmpq %rbx, %rdi
jne 0x33113
jmp 0x3311d
jmp 0x3311a
movq %rax, %r14
leaq 0x58(%rsp), %rax
movq -0x10(%rax), %rdi
cmpq %rax, %rdi
je 0x330cc
callq 0x3310
jmp 0x330cc
movq %rax, %r14
movq 0x68(%rsp), %rdi
cmpq %rbp, %rdi
je 0x330e2
callq 0x3310
jmp 0x330e2
jmp 0x330df
movq %rax, %r14
movq 0x88(%rsp), %rdi
cmpq %r13, %rdi
jne 0x33113
jmp 0x3311d
jmp 0x3311a
jmp 0x3311a
jmp 0x3311a
jmp 0x3311a
jmp 0x3311a
movq %rax, %r14
movq 0xe8(%rsp), %rdi
leaq 0xf8(%rsp), %rax
cmpq %rax, %rdi
je 0x3311d
callq 0x3310
jmp 0x3311d
movq %rax, %r14
movq 0xa8(%rsp), %rdi
leaq 0xb8(%rsp), %rax
cmpq %rax, %rdi
je 0x33137
callq 0x3310
movq %r14, %rdi
callq 0x3580