// File: /workspace/github_c_makefile_star2_O0/sc932[P]ALE/src/ALElike.c
// Function: computeDepthStats
int computeDepthStats(assemblyT *theAssembly, libraryParametersT *libParams){
    // 1. Find the GC content of each read
    int i, j, k;
    int GCpct;
    long place;
    double depthNormalizer[102];
    long depthNormalizerCount[102];

    double negBinomParam_p[102];
    double negBinomParam_r[102];
    double negBinomParamZnorm_r[102];
    double tempLogLike;
    long tooLowCoverageBases = 0;
    long noGCInformation = 0;
    long unmappableRegions = 0;
    long unmappableBases = 0;
    long tmpunmappable = 0;
    unsigned char *GCcont = NULL;

    // initialize
    for(j = 0; j < 101; j++){
        depthNormalizer[j] = 0.0;
        negBinomParam_p[j] = 0.0;
        negBinomParam_r[j] = 0.0;
        negBinomParamZnorm_r[j] = 0.0;
        depthNormalizerCount[j] = 0;
    }
    if (!isMetagenome()) {
        for(i = 0; i < theAssembly->numContigs; i++){ // for each contig
            contig_t *contig = theAssembly->contigs[i];
            if (GCcont != NULL) free(GCcont);
            GCcont = calculateContigGCcont(contig, libParams->avgReadSize);
            tmpunmappable = 0;
            for(j = 0; j < contig->seqLen; j++){
                float depth = contig->depth[j];
                if (depth < 0.1) {
                    tooLowCoverageBases++;
        			if (contig->seq[j] == 'N')
        				tmpunmappable++;
                } else
                	tmpunmappable = 0;

                // Do not record this depth and normalization if this region is not possible to map to.
                if (tmpunmappable > libParams->avgReadSize / 3) {
                	if (tmpunmappable == (libParams->avgReadSize / 3) + 1)
                		unmappableRegions++;
                	unmappableBases++;
                	continue;
                }

                theAssembly->depthAvgSum += depth;
                theAssembly->depthAvgNorm += 1.0;
                GCpct = GCcont[j];
                if (GCpct > 100) {
                    noGCInformation++;
                    continue;
                }
                depthNormalizer[GCpct] += depth;
                depthNormalizerCount[GCpct] += 1;
            }
        }

        // 2. Find the parameters for the distributions

        computeNormaliziedDepthGCParameters(depthNormalizer, depthNormalizerCount, negBinomParam_r, negBinomParam_p, negBinomParamZnorm_r, theAssembly->depthAvgSum / theAssembly->depthAvgNorm);
    }

    for(i = 0; i < theAssembly->numContigs; i++){ // for each contig
        contig_t *contig = theAssembly->contigs[i];
        if (GCcont != NULL) free(GCcont);
        GCcont = calculateContigGCcont(contig, libParams->avgReadSize);
        if (isMetagenome()) {
        	// re-initialize depth for this contig
        	for(j = 0; j < 101; j++){
        		depthNormalizer[j] = 0.0;
        		negBinomParam_p[j] = 0.0;
        		negBinomParam_r[j] = 0.0;
        		negBinomParamZnorm_r[j] = 0.0;
        		depthNormalizerCount[j] = 0;
        	}

            tmpunmappable = 0;
            double avgDepth = 0.0;
            double countDepth = 0.0;
        	for(j = 0; j < contig->seqLen; j++){
        		float depth = contig->depth[j];
        		if (depth < 0.1) {
        			tooLowCoverageBases++;
        			if (contig->seq[j] == 'N')
        				tmpunmappable++;
        		} else
        			tmpunmappable = 0;

                // Do not record this depth and normalization if this region is not possible to map to.
                if (tmpunmappable > libParams->avgReadSize / 3) {
                	if (tmpunmappable == libParams->avgReadSize / 3 + 1)
                		unmappableRegions++;
                	unmappableBases++;
                	continue;
                }

        		theAssembly->depthAvgSum += depth;
        		theAssembly->depthAvgNorm += 1.0;
        		avgDepth += depth;
        		countDepth += 1.0;
        		GCpct = GCcont[j];
        		if (GCpct > 100) {
        			noGCInformation++;
        			continue;
        		}
        		depthNormalizer[GCpct] += depth;
        		depthNormalizerCount[GCpct] += 1;
        	}

        	// 2. Find the parameters for the distributions

            computeNormaliziedDepthGCParameters(depthNormalizer, depthNormalizerCount, negBinomParam_r, negBinomParam_p, negBinomParamZnorm_r, avgDepth / countDepth);
        }

        //printf("Calculating likelihoods for %d positions\n", contig->seqLen);
        // 3. Find the depth likelihood
        for(j = 0; j < contig->seqLen; j++){
            GCpct = GCcont[j];
            if (GCpct > 100){
                //printf("GC correction failed for contig %d at %d '%c'\n", i, j, contig->seq[j]);
            	contig->depthLogLikelihood[j]  = getMinLogLike();
            	theAssembly->totalScore += getMinLogLike();
            	theAssembly->depthScoreAvgSum += getMinLogLike();
            	theAssembly->depthScoreAvgNorm += 1.0;
            } else {

            	// compute the depth likelihood using poisson or negBinomial
            	// contig->depth[j] is the depth at position j
            	// depthNormalizer[GCpct] is avg depth for that GC content

            	// tempLike = poissonPMF(contig->depth[j], depthNormalizer[GCpct]); // log poisson pmf
            	tempLogLike = negBinomPMF((int)floor(contig->depth[j]), negBinomParam_r[GCpct], 0.5); // log neg binom pmf
            	assert(tempLogLike <= 0.0);

            	// z normalization
            	tempLogLike -= negBinomParamZnorm_r[GCpct];
            	//if((int)floor(negBinomParam_r[GCpct]) < 2047){
            	//    tempLike -= negBinomZ[(int)floor(negBinomParam_r[GCpct])];
            	//}else{
            	//    // not in lookup table, compute
            	//    tempLike -= getNegBinomZnorm(negBinomParam_r[GCpct]);
            	//}
            	tempLogLike = validateLogLikelihood( tempLogLike );

            	// apply to assembly and totalScore
            	contig->depthLogLikelihood[j] = tempLogLike;
            	theAssembly->totalScore += tempLogLike; // depth contribution to totalScore at this position
            	theAssembly->depthScoreAvgSum += tempLogLike;
            	theAssembly->depthScoreAvgNorm += 1.0;
            }
            // at this point contig->matchLikelihood[j] contains the sum of the logs of the TOTAL likelihood of all the reads that map over position j
            // then we take the gemetric average by dividing by the depth and change it (if it is a valid likelihood)
            
            if (contig->depth[j] == 0) { // avoid div by zero...
            	contig->matchLogLikelihood[j] = getMinLogLike();
            	contig->insertLogLikelihood[j] = getMinLogLike();
            } else {
                // match
            	contig->matchLogLikelihood[j] =  validateLogLikelihood( contig->matchLogLikelihood[j]/contig->depth[j] );;
                // insert
                contig->insertLogLikelihood[j] = validateLogLikelihood( contig->insertLogLikelihood[j]/contig->depth[j] );;
            }
            
        }
    }
    if (GCcont != NULL) free(GCcont);
    printf("bases with no GC depth correction metric (small contigs): %ld\n", noGCInformation);
    printf("unmappable bases: %ld across %ld different regions\n", unmappableBases, unmappableRegions);
    printf("bases with 0 coverage: %ld\n", tooLowCoverageBases);
    return 1;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x10b0, %rsp           # imm = 0x10B0
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq $0x0, -0x1030(%rbp)
movq $0x0, -0x1038(%rbp)
movq $0x0, -0x1040(%rbp)
movq $0x0, -0x1048(%rbp)
movq $0x0, -0x1050(%rbp)
movq $0x0, -0x1058(%rbp)
movl $0x0, -0x18(%rbp)
cmpl $0x65, -0x18(%rbp)
jge 0x81fd
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0x360(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0x9c0(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0xcf0(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0x1020(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
movq $0x0, -0x690(%rbp,%rax,8)
movl -0x18(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x18(%rbp)
jmp 0x819c
movb $0x0, %al
callq 0xbec0
cmpl $0x0, %eax
jne 0x8496
movl $0x0, -0x14(%rbp)
movl -0x14(%rbp), %eax
movq -0x8(%rbp), %rcx
cmpl (%rcx), %eax
jge 0x845c
movq -0x8(%rbp), %rax
movq 0xa0(%rax), %rax
movslq -0x14(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0x1060(%rbp)
cmpq $0x0, -0x1058(%rbp)
je 0x8253
movq -0x1058(%rbp), %rdi
callq 0x3050
movq -0x1060(%rbp), %rdi
movq -0x10(%rbp), %rax
movq 0x380(%rax), %rax
movl %eax, %esi
callq 0xd810
movq %rax, -0x1058(%rbp)
movq $0x0, -0x1050(%rbp)
movl $0x0, -0x18(%rbp)
movl -0x18(%rbp), %eax
movq -0x1060(%rbp), %rcx
cmpl 0x8(%rcx), %eax
jge 0x844c
movq -0x1060(%rbp), %rax
movq 0x30(%rax), %rax
movslq -0x18(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
movss %xmm0, -0x1064(%rbp)
movss -0x1064(%rbp), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x23808(%rip), %xmm0    # 0x2bad0
ucomisd %xmm1, %xmm0
jbe 0x830c
movq -0x1030(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1030(%rbp)
movq -0x1060(%rbp), %rax
movq 0x10(%rax), %rax
movslq -0x18(%rbp), %rcx
movsbl (%rax,%rcx), %eax
cmpl $0x4e, %eax
jne 0x830a
movq -0x1050(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1050(%rbp)
jmp 0x8317
movq $0x0, -0x1050(%rbp)
movq -0x1050(%rbp), %rax
movq %rax, -0x1090(%rbp)
movq -0x10(%rbp), %rax
movq 0x380(%rax), %rax
movl $0x3, %ecx
cqto
idivq %rcx
movq %rax, %rcx
movq -0x1090(%rbp), %rax
cmpq %rcx, %rax
jle 0x83a8
movq -0x1050(%rbp), %rax
movq %rax, -0x1098(%rbp)
movq -0x10(%rbp), %rax
movq 0x380(%rax), %rax
movl $0x3, %ecx
cqto
idivq %rcx
movq %rax, %rcx
movq -0x1098(%rbp), %rax
addq $0x1, %rcx
cmpq %rcx, %rax
jne 0x8391
movq -0x1040(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1040(%rbp)
movq -0x1048(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1048(%rbp)
jmp 0x843e
movss -0x1064(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
movq -0x8(%rbp), %rax
addsd 0x58(%rax), %xmm0
movsd %xmm0, 0x58(%rax)
movq -0x8(%rbp), %rax
movsd 0x23672(%rip), %xmm0    # 0x2ba40
addsd 0x60(%rax), %xmm0
movsd %xmm0, 0x60(%rax)
movq -0x1058(%rbp), %rax
movslq -0x18(%rbp), %rcx
movzbl (%rax,%rcx), %eax
movl %eax, -0x20(%rbp)
cmpl $0x64, -0x20(%rbp)
jle 0x8404
movq -0x1038(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1038(%rbp)
jmp 0x843e
movss -0x1064(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
movslq -0x20(%rbp), %rax
addsd -0x360(%rbp,%rax,8), %xmm0
movsd %xmm0, -0x360(%rbp,%rax,8)
movslq -0x20(%rbp), %rax
movq -0x690(%rbp,%rax,8), %rcx
addq $0x1, %rcx
movq %rcx, -0x690(%rbp,%rax,8)
movl -0x18(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x18(%rbp)
jmp 0x8285
jmp 0x844e
movl -0x14(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x14(%rbp)
jmp 0x8214
leaq -0x360(%rbp), %rdi
leaq -0x690(%rbp), %rsi
leaq -0xcf0(%rbp), %rdx
leaq -0x9c0(%rbp), %rcx
leaq -0x1020(%rbp), %r8
movq -0x8(%rbp), %rax
movsd 0x58(%rax), %xmm0
movq -0x8(%rbp), %rax
divsd 0x60(%rax), %xmm0
callq 0x7fc0
movl $0x0, -0x14(%rbp)
movl -0x14(%rbp), %eax
movq -0x8(%rbp), %rcx
cmpl (%rcx), %eax
jge 0x8a81
movq -0x8(%rbp), %rax
movq 0xa0(%rax), %rax
movslq -0x14(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0x1070(%rbp)
cmpq $0x0, -0x1058(%rbp)
je 0x84dc
movq -0x1058(%rbp), %rdi
callq 0x3050
movq -0x1070(%rbp), %rdi
movq -0x10(%rbp), %rax
movq 0x380(%rax), %rax
movl %eax, %esi
callq 0xd810
movq %rax, -0x1058(%rbp)
movb $0x0, %al
callq 0xbec0
cmpl $0x0, %eax
je 0x87cf
movl $0x0, -0x18(%rbp)
cmpl $0x65, -0x18(%rbp)
jge 0x8574
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0x360(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0x9c0(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0xcf0(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
xorps %xmm0, %xmm0
movsd %xmm0, -0x1020(%rbp,%rax,8)
movslq -0x18(%rbp), %rax
movq $0x0, -0x690(%rbp,%rax,8)
movl -0x18(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x18(%rbp)
jmp 0x8513
movq $0x0, -0x1050(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x1078(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x1080(%rbp)
movl $0x0, -0x18(%rbp)
movl -0x18(%rbp), %eax
movq -0x1070(%rbp), %rcx
cmpl 0x8(%rcx), %eax
jge 0x8797
movq -0x1070(%rbp), %rax
movq 0x30(%rax), %rax
movslq -0x18(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
movss %xmm0, -0x1084(%rbp)
movss -0x1084(%rbp), %xmm0
cvtss2sd %xmm0, %xmm1
movsd 0x234f1(%rip), %xmm0    # 0x2bad0
ucomisd %xmm1, %xmm0
jbe 0x8623
movq -0x1030(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1030(%rbp)
movq -0x1070(%rbp), %rax
movq 0x10(%rax), %rax
movslq -0x18(%rbp), %rcx
movsbl (%rax,%rcx), %eax
cmpl $0x4e, %eax
jne 0x8621
movq -0x1050(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1050(%rbp)
jmp 0x862e
movq $0x0, -0x1050(%rbp)
movq -0x1050(%rbp), %rax
movq %rax, -0x10a0(%rbp)
movq -0x10(%rbp), %rax
movq 0x380(%rax), %rax
movl $0x3, %ecx
cqto
idivq %rcx
movq %rax, %rcx
movq -0x10a0(%rbp), %rax
cmpq %rcx, %rax
jle 0x86bf
movq -0x1050(%rbp), %rax
movq %rax, -0x10a8(%rbp)
movq -0x10(%rbp), %rax
movq 0x380(%rax), %rax
movl $0x3, %ecx
cqto
idivq %rcx
movq %rax, %rcx
movq -0x10a8(%rbp), %rax
addq $0x1, %rcx
cmpq %rcx, %rax
jne 0x86a8
movq -0x1040(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1040(%rbp)
movq -0x1048(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1048(%rbp)
jmp 0x8789
movss -0x1084(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
movq -0x8(%rbp), %rax
addsd 0x58(%rax), %xmm0
movsd %xmm0, 0x58(%rax)
movq -0x8(%rbp), %rax
movsd 0x2335b(%rip), %xmm0    # 0x2ba40
addsd 0x60(%rax), %xmm0
movsd %xmm0, 0x60(%rax)
movss -0x1084(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
addsd -0x1078(%rbp), %xmm0
movsd %xmm0, -0x1078(%rbp)
movsd 0x2332d(%rip), %xmm0    # 0x2ba40
addsd -0x1080(%rbp), %xmm0
movsd %xmm0, -0x1080(%rbp)
movq -0x1058(%rbp), %rax
movslq -0x18(%rbp), %rcx
movzbl (%rax,%rcx), %eax
movl %eax, -0x20(%rbp)
cmpl $0x64, -0x20(%rbp)
jle 0x874f
movq -0x1038(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x1038(%rbp)
jmp 0x8789
movss -0x1084(%rbp), %xmm0
cvtss2sd %xmm0, %xmm0
movslq -0x20(%rbp), %rax
addsd -0x360(%rbp,%rax,8), %xmm0
movsd %xmm0, -0x360(%rbp,%rax,8)
movslq -0x20(%rbp), %rax
movq -0x690(%rbp,%rax,8), %rcx
addq $0x1, %rcx
movq %rcx, -0x690(%rbp,%rax,8)
movl -0x18(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x18(%rbp)
jmp 0x859c
leaq -0x360(%rbp), %rdi
leaq -0x690(%rbp), %rsi
leaq -0xcf0(%rbp), %rdx
leaq -0x9c0(%rbp), %rcx
leaq -0x1020(%rbp), %r8
movsd -0x1078(%rbp), %xmm0
divsd -0x1080(%rbp), %xmm0
callq 0x7fc0
movl $0x0, -0x18(%rbp)
movl -0x18(%rbp), %eax
movq -0x1070(%rbp), %rcx
cmpl 0x8(%rcx), %eax
jge 0x8a71
movq -0x1058(%rbp), %rax
movslq -0x18(%rbp), %rcx
movzbl (%rax,%rcx), %eax
movl %eax, -0x20(%rbp)
cmpl $0x64, -0x20(%rbp)
jle 0x8865
movb $0x0, %al
callq 0xbdf0
cvtsd2ss %xmm0, %xmm0
movq -0x1070(%rbp), %rax
movq 0x48(%rax), %rax
movslq -0x18(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movb $0x0, %al
callq 0xbdf0
movq -0x8(%rbp), %rax
addsd 0x10(%rax), %xmm0
movsd %xmm0, 0x10(%rax)
movb $0x0, %al
callq 0xbdf0
movq -0x8(%rbp), %rax
addsd 0x48(%rax), %xmm0
movsd %xmm0, 0x48(%rax)
movq -0x8(%rbp), %rax
movsd 0x231ea(%rip), %xmm0    # 0x2ba40
addsd 0x50(%rax), %xmm0
movsd %xmm0, 0x50(%rax)
jmp 0x896e
movq -0x1070(%rbp), %rax
movq 0x30(%rax), %rax
movslq -0x18(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0x3560
cvttsd2si %xmm0, %edi
movslq -0x20(%rbp), %rax
movsd -0xcf0(%rbp,%rax,8), %xmm0
movsd 0x231d5(%rip), %xmm1    # 0x2ba70
callq 0x7dd0
movsd %xmm0, -0x1028(%rbp)
xorps %xmm0, %xmm0
ucomisd -0x1028(%rbp), %xmm0
jb 0x88b7
jmp 0x88d6
leaq 0x27cdc(%rip), %rdi     # 0x3059a
leaq 0x27500(%rip), %rsi     # 0x2fdc5
movl $0x68e, %edx            # imm = 0x68E
leaq 0x27cdc(%rip), %rcx     # 0x305ad
callq 0x3220
movslq -0x20(%rbp), %rax
movsd -0x1020(%rbp,%rax,8), %xmm1
movsd -0x1028(%rbp), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, -0x1028(%rbp)
movsd -0x1028(%rbp), %xmm0
callq 0xbe00
movsd %xmm0, -0x1028(%rbp)
movsd -0x1028(%rbp), %xmm0
cvtsd2ss %xmm0, %xmm0
movq -0x1070(%rbp), %rax
movq 0x48(%rax), %rax
movslq -0x18(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movsd -0x1028(%rbp), %xmm0
movq -0x8(%rbp), %rax
addsd 0x10(%rax), %xmm0
movsd %xmm0, 0x10(%rax)
movsd -0x1028(%rbp), %xmm0
movq -0x8(%rbp), %rax
addsd 0x48(%rax), %xmm0
movsd %xmm0, 0x48(%rax)
movq -0x8(%rbp), %rax
movsd 0x230dc(%rip), %xmm0    # 0x2ba40
addsd 0x50(%rax), %xmm0
movsd %xmm0, 0x50(%rax)
movq -0x1070(%rbp), %rax
movq 0x30(%rax), %rax
movslq -0x18(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
xorps %xmm1, %xmm1
ucomiss %xmm1, %xmm0
jne 0x89cf
jp 0x89cf
movb $0x0, %al
callq 0xbdf0
cvtsd2ss %xmm0, %xmm0
movq -0x1070(%rbp), %rax
movq 0x38(%rax), %rax
movslq -0x18(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movb $0x0, %al
callq 0xbdf0
cvtsd2ss %xmm0, %xmm0
movq -0x1070(%rbp), %rax
movq 0x40(%rax), %rax
movslq -0x18(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
jmp 0x8a61
movq -0x1070(%rbp), %rax
movq 0x38(%rax), %rax
movslq -0x18(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
movq -0x1070(%rbp), %rax
movq 0x30(%rax), %rax
movslq -0x18(%rbp), %rcx
divss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0xbe00
cvtsd2ss %xmm0, %xmm0
movq -0x1070(%rbp), %rax
movq 0x38(%rax), %rax
movslq -0x18(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
movq -0x1070(%rbp), %rax
movq 0x40(%rax), %rax
movslq -0x18(%rbp), %rcx
movss (%rax,%rcx,4), %xmm0
movq -0x1070(%rbp), %rax
movq 0x30(%rax), %rax
movslq -0x18(%rbp), %rcx
divss (%rax,%rcx,4), %xmm0
cvtss2sd %xmm0, %xmm0
callq 0xbe00
cvtsd2ss %xmm0, %xmm0
movq -0x1070(%rbp), %rax
movq 0x40(%rax), %rax
movslq -0x18(%rbp), %rcx
movss %xmm0, (%rax,%rcx,4)
jmp 0x8a63
movl -0x18(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x18(%rbp)
jmp 0x87d6
jmp 0x8a73
movl -0x14(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x14(%rbp)
jmp 0x849d
cmpq $0x0, -0x1058(%rbp)
je 0x8a97
movq -0x1058(%rbp), %rdi
callq 0x3050
movq -0x1038(%rbp), %rsi
leaq 0x27b41(%rip), %rdi     # 0x305e6
movb $0x0, %al
callq 0x3200
movq -0x1048(%rbp), %rsi
movq -0x1040(%rbp), %rdx
leaq 0x27b64(%rip), %rdi     # 0x30625
movb $0x0, %al
callq 0x3200
movq -0x1030(%rbp), %rsi
leaq 0x27b83(%rip), %rdi     # 0x30659
movb $0x0, %al
callq 0x3200
movl $0x1, %eax
addq $0x10b0, %rsp           # imm = 0x10B0
popq %rbp
retq
nopl (%rax,%rax)