// File: /workspace/github_c_makefile_star2_O2/icebert[P]pblat/lib/chain.c
// Function: chainFastSubsetOnT
void chainFastSubsetOnT(struct chain *chain, struct cBlock *firstBlock,
	int subStart, int subEnd, struct chain **retSubChain,  struct chain **retChainToFree)
/* Get subchain as in chainSubsetOnT. Pass in initial block that may
 * be known from some index to speed things up. */
{
struct chain *sub = NULL;
struct cBlock *oldB, *b, *bList = NULL;
int qStart = BIGNUM, qEnd = -BIGNUM;
int tStart = BIGNUM, tEnd = -BIGNUM;

/* Check for easy case. */
if (subStart <= chain->tStart && subEnd >= chain->tEnd)
    {
    *retSubChain = chain;
    *retChainToFree = NULL;
    return;
    }
/* Build new block list and calculate bounds. */
for (oldB = firstBlock; oldB != NULL; oldB = oldB->next)
    {
    if (oldB->tStart >= subEnd)
        break;
    b = CloneVar(oldB);
    if (b->tStart < subStart)
        {
	b->qStart += subStart - b->tStart;
	b->tStart = subStart;
	}
    if (b->tEnd > subEnd)
        {
	b->qEnd -= b->tEnd - subEnd;
	b->tEnd = subEnd;
	}
    slAddHead(&bList, b);
    if (qStart > b->qStart)
        qStart = b->qStart;
    if (qEnd < b->qEnd)
        qEnd = b->qEnd;
    if (tStart > b->tStart)
        tStart = b->tStart;
    if (tEnd < b->tEnd)
        tEnd = b->tEnd;
    }
slReverse(&bList);

/* Make new chain based on old. */
if (bList != NULL)
    {
    double sizeRatio;
    AllocVar(sub);
    sub->blockList = bList;
    sub->qName = cloneString(chain->qName);
    sub->qSize = chain->qSize;
    sub->qStrand = chain->qStrand;
    sub->qStart = qStart;
    sub->qEnd = qEnd;
    sub->tName = cloneString(chain->tName);
    sub->tSize = chain->tSize;
    sub->tStart = tStart;
    sub->tEnd = tEnd;
    sub->id = chain->id;

    /* Fake new score. */
    sizeRatio = (sub->tEnd - sub->tStart);
    sizeRatio /= (chain->tEnd - chain->tStart);
    sub->score = sizeRatio * chain->score;
    }
*retSubChain = *retChainToFree = sub;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movl %ecx, %ebp
andq $0x0, 0x8(%rsp)
movq %rsi, %r13
movl %edx, 0x4(%rsp)
cmpl %edx, 0x24(%rdi)
jl 0x1b72d
cmpl %ebp, 0x28(%rdi)
jle 0x1b886
movq %rdi, 0x10(%rsp)
movq %r8, 0x18(%rsp)
movq %r9, 0x20(%rsp)
movl $0xc0000001, %ebx       # imm = 0xC0000001
movl $0x3fffffff, %r12d      # imm = 0x3FFFFFFF
xorl %r14d, %r14d
movl $0x3fffffff, (%rsp)     # imm = 0x3FFFFFFF
movl $0xc0000001, %r15d      # imm = 0xC0000001
testq %r13, %r13
je 0x1b7cd
cmpl %ebp, 0x8(%r13)
jge 0x1b7cd
movq %r13, %rdi
pushq $0x28
popq %rsi
callq 0x1c27c
movl 0x8(%rax), %ecx
movl 0x4(%rsp), %edx
subl %ecx, %edx
jle 0x1b782
addl %edx, 0x10(%rax)
movl 0x4(%rsp), %ecx
movl %ecx, 0x8(%rax)
movl 0xc(%rax), %edx
movl %ebp, %esi
subl %edx, %esi
jge 0x1b798
addl 0x14(%rax), %esi
movl %esi, 0x14(%rax)
movl %ebp, 0xc(%rax)
movl %ebp, %edx
jmp 0x1b79b
movl 0x14(%rax), %esi
movq %r14, (%rax)
movq %rax, 0x8(%rsp)
movl 0x10(%rax), %edi
cmpl %edi, %r12d
cmovgel %edi, %r12d
cmpl %esi, %ebx
cmovlel %esi, %ebx
movl (%rsp), %esi
cmpl %ecx, %esi
cmovgel %ecx, %esi
movl %esi, (%rsp)
cmpl %edx, %r15d
cmovlel %edx, %r15d
movq (%r13), %r13
movq %rax, %r14
jmp 0x1b757
leaq 0x8(%rsp), %r13
movq %r13, %rdi
callq 0x1c567
cmpq $0x0, (%r13)
je 0x1b871
pushq $0x50
popq %rdi
callq 0x28bb0
movq %rax, %r13
movq 0x8(%rsp), %rax
movq %rax, 0x8(%r13)
movq 0x10(%rsp), %r14
movq 0x30(%r14), %rdi
callq 0x1c306
movq %rax, 0x30(%r13)
movl 0x38(%r14), %eax
movl %eax, 0x38(%r13)
movb 0x3c(%r14), %al
movb %al, 0x3c(%r13)
movl %r12d, 0x40(%r13)
movl %ebx, 0x44(%r13)
movq 0x18(%r14), %rdi
callq 0x1c306
movq %rax, 0x18(%r13)
movl 0x20(%r14), %eax
movl %eax, 0x20(%r13)
movl (%rsp), %ecx
movl %ecx, 0x24(%r13)
movl %r15d, 0x28(%r13)
movl 0x48(%r14), %eax
movl %eax, 0x48(%r13)
subl %ecx, %r15d
cvtsi2sd %r15d, %xmm0
movl 0x28(%r14), %eax
subl 0x24(%r14), %eax
cvtsi2sd %eax, %xmm1
divsd %xmm1, %xmm0
mulsd 0x10(%r14), %xmm0
movsd %xmm0, 0x10(%r13)
jmp 0x1b874
xorl %r13d, %r13d
movq 0x20(%rsp), %rax
movq %r13, (%rax)
movq 0x18(%rsp), %rax
movq %r13, (%rax)
jmp 0x1b88d
movq %rdi, (%r8)
andq $0x0, (%r9)
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq