// File: /workspace/github_c_makefile_star2_O1/icebert[P]pblat/htslib/cram/cram_decode.c
// Function: cram_dependent_data_series
int cram_dependent_data_series(cram_fd *fd,
			       cram_block_compression_hdr *hdr,
			       cram_slice *s) {
    int *block_used;
    int core_used = 0;
    int i;
    static int i_to_id[] = {
	DS_BF, DS_AP, DS_FP, DS_RL, DS_DL, DS_NF, DS_BA, DS_QS,
	DS_FC, DS_FN, DS_BS, DS_IN, DS_RG, DS_MQ, DS_TL, DS_RN,
	DS_NS, DS_NP, DS_TS, DS_MF, DS_CF, DS_RI, DS_RS, DS_PD,
	DS_HC, DS_SC, DS_BB, DS_QQ,
    };
    uint32_t orig_ds;

    /*
     * Set the data_series bit field based on fd->required_fields
     * contents.
     */
    if (fd->required_fields && fd->required_fields != INT_MAX) {
	hdr->data_series = 0;

	if (fd->required_fields & SAM_QNAME)
	    hdr->data_series |= CRAM_RN;

	if (fd->required_fields & SAM_FLAG)
	    hdr->data_series |= CRAM_BF;

	if (fd->required_fields & SAM_RNAME)
	    hdr->data_series |= CRAM_RI | CRAM_BF;

	if (fd->required_fields & SAM_POS)
	    hdr->data_series |= CRAM_AP | CRAM_BF;

	if (fd->required_fields & SAM_MAPQ)
	    hdr->data_series |= CRAM_MQ;

	if (fd->required_fields & SAM_CIGAR)
	    hdr->data_series |= CRAM_CIGAR;

	if (fd->required_fields & SAM_RNEXT)
	    hdr->data_series |= CRAM_CF | CRAM_NF | CRAM_RI | CRAM_NS |CRAM_BF;

	if (fd->required_fields & SAM_PNEXT)
	    hdr->data_series |= CRAM_CF | CRAM_NF | CRAM_AP | CRAM_NP | CRAM_BF;

	if (fd->required_fields & SAM_TLEN)
	    hdr->data_series |= CRAM_CF | CRAM_NF | CRAM_AP | CRAM_TS |
		CRAM_BF | CRAM_MF | CRAM_RI | CRAM_CIGAR;

	if (fd->required_fields & SAM_SEQ)
	    hdr->data_series |= CRAM_SEQ;

	if (!(fd->required_fields & SAM_AUX))
	    // No easy way to get MD/NM without other tags at present
	    fd->decode_md = 0;

	if (fd->required_fields & SAM_QUAL)
	    hdr->data_series |= CRAM_QUAL;

	if (fd->required_fields & SAM_AUX)
	    hdr->data_series |= CRAM_RG | CRAM_TL | CRAM_aux;

	if (fd->required_fields & SAM_RGAUX)
	    hdr->data_series |= CRAM_RG | CRAM_BF;

	// Always uncompress CORE block
	if (cram_uncompress_block(s->block[0]))
	    return -1;
    } else {
	hdr->data_series = CRAM_ALL;

	for (i = 0; i < s->hdr->num_blocks; i++) {
	    if (cram_uncompress_block(s->block[i]))
		return -1;
	}

	return 0;
    }

    block_used = calloc(s->hdr->num_blocks+1, sizeof(int));
    if (!block_used)
	return -1;

    do {
	/*
	 * Also set data_series based on code prerequisites. Eg if we need
	 * CRAM_QS then we also need to know CRAM_RL so we know how long it
	 * is, or if we need FC/FP then we also need FN (number of features).
	 *
	 * It's not reciprocal though. We may be needing to decode FN
	 * but have no need to decode FC, FP and cigar ops.
	 */
	if (hdr->data_series & CRAM_RS)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_PD)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_HC)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_QS)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_IN)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_SC)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_BS)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_DL)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_BA)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_BB)    hdr->data_series |= CRAM_FC|CRAM_FP;
	if (hdr->data_series & CRAM_QQ)    hdr->data_series |= CRAM_FC|CRAM_FP;

	// cram_decode_seq() needs seq[] array
	if (hdr->data_series & (CRAM_SEQ|CRAM_CIGAR)) hdr->data_series |= CRAM_RL;

	if (hdr->data_series & CRAM_FP)    hdr->data_series |= CRAM_FC;
	if (hdr->data_series & CRAM_FC)    hdr->data_series |= CRAM_FN;
	if (hdr->data_series & CRAM_aux)   hdr->data_series |= CRAM_TL;
	if (hdr->data_series & CRAM_MF)    hdr->data_series |= CRAM_CF;
	if (hdr->data_series & CRAM_MQ)    hdr->data_series |= CRAM_BF;
	if (hdr->data_series & CRAM_BS)    hdr->data_series |= CRAM_RI;
	if (hdr->data_series & (CRAM_MF |CRAM_NS |CRAM_NP |CRAM_TS |CRAM_NF))
	    hdr->data_series |= CRAM_CF;
	if (!hdr->read_names_included && hdr->data_series & CRAM_RN)
	    hdr->data_series |= CRAM_CF | CRAM_NF;
	if (hdr->data_series & (CRAM_BA | CRAM_QS | CRAM_BB | CRAM_QQ))
	    hdr->data_series |= CRAM_BF | CRAM_CF | CRAM_RL;

	orig_ds = hdr->data_series;

	// Find which blocks are in use.
	for (i = 0; i < sizeof(i_to_id)/sizeof(*i_to_id); i++) {
	    int bnum1, bnum2, j;
	    cram_codec *c = hdr->codecs[i_to_id[i]];

	    if (!(hdr->data_series & (1<<i)))
		continue;

	    if (!c)
		continue;

	    bnum1 = cram_codec_to_id(c, &bnum2);

	    for (;;) {
		switch (bnum1) {
		case -2:
		    break;

		case -1:
		    core_used = 1;
		    break;

		default:
		    for (j = 0; j < s->hdr->num_blocks; j++) {
			if (s->block[j]->content_type == EXTERNAL &&
			    s->block[j]->content_id == bnum1) {
			    block_used[j] = 1;
			    if (cram_uncompress_block(s->block[j])) {
				free(block_used);
				return -1;
			    }
			}
		    }
		    break;
		}

		if (bnum2 == -2 || bnum1 == bnum2)
		    break;

		bnum1 = bnum2; // 2nd pass
	    }
	}

	// Tags too
	if ((fd->required_fields & SAM_AUX) ||
	    (hdr->data_series & CRAM_aux)) {
	    for (i = 0; i < CRAM_MAP_HASH; i++) {
		int bnum1, bnum2, j;
		cram_map *m = hdr->tag_encoding_map[i];

		while (m) {
		    cram_codec *c = m->codec;
		    if (!c)
			continue;

		    bnum1 = cram_codec_to_id(c, &bnum2);

		    for (;;) {
			switch (bnum1) {
			case -2:
			    break;

			case -1:
			    core_used = 1;
			    break;

			default:
			    for (j = 0; j < s->hdr->num_blocks; j++) {
				if (s->block[j]->content_type == EXTERNAL &&
				    s->block[j]->content_id == bnum1) {
				    block_used[j] = 1;
				    if (cram_uncompress_block(s->block[j])) {
					free(block_used);
					return -1;
				    }
				}
			    }
			    break;
			}

			if (bnum2 == -2 || bnum1 == bnum2)
			    break;

			bnum1 = bnum2; // 2nd pass
		    }

		    m = m->next;
		}
	    }
	}

	// We now know which blocks are in used, so repeat and find
	// which other data series need to be added.
	for (i = 0; i < sizeof(i_to_id)/sizeof(*i_to_id); i++) {
	    int bnum1, bnum2, j;
	    cram_codec *c = hdr->codecs[i_to_id[i]];

	    if (!c)
		continue;

	    bnum1 = cram_codec_to_id(c, &bnum2);

	    for (;;) {
		switch (bnum1) {
		case -2:
		    break;

		case -1:
		    if (core_used) {
			//printf(" + data series %08x:\n", 1<<i);
			hdr->data_series |= 1<<i;
		    }
		    break;

		default:
		    for (j = 0; j < s->hdr->num_blocks; j++) {
			if (s->block[j]->content_type == EXTERNAL &&
			    s->block[j]->content_id == bnum1) {
			    if (block_used[j]) {
				//printf(" + data series %08x:\n", 1<<i);
				hdr->data_series |= 1<<i;
			    }
			}
		    }
		    break;
		}

		if (bnum2 == -2 || bnum1 == bnum2)
		    break;

		bnum1 = bnum2; // 2nd pass
	    }
	}

	// Tags too
	for (i = 0; i < CRAM_MAP_HASH; i++) {
	    int bnum1, bnum2, j;
	    cram_map *m = hdr->tag_encoding_map[i];

	    while (m) {
		cram_codec *c = m->codec;
		if (!c)
		    continue;

		bnum1 = cram_codec_to_id(c, &bnum2);
		
		for (;;) {
		    switch (bnum1) {
		    case -2:
			break;

		    case -1:
			//printf(" + data series %08x:\n", CRAM_aux);
			hdr->data_series |= CRAM_aux;
			break;

		    default:
			for (j = 0; j < s->hdr->num_blocks; j++) {
			    if (s->block[j]->content_type &&
				s->block[j]->content_id == bnum1) {
				if (block_used[j]) {
				    //printf(" + data series %08x:\n",
				    //       CRAM_aux);
				    hdr->data_series |= CRAM_aux;
				}
			    }
			}
			break;
		    }

		    if (bnum2 == -2 || bnum1 == bnum2)
			break;

		    bnum1 = bnum2; // 2nd pass
		}

		m = m->next;
	    }
	}
    } while (orig_ds != hdr->data_series);

    free(block_used);
    return 0;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %rdx, %rbx
movq %rsi, %r14
movl 0x224(%rdi), %eax
testl %eax, %eax
je 0x76c9f
cmpl $0x7fffffff, %eax       # imm = 0x7FFFFFFF
jne 0x76ce1
movl $0x7fffffff, 0x410(%r14) # imm = 0x7FFFFFFF
movq (%rbx), %rax
cmpl $0x0, 0x20(%rax)
jle 0x76cda
xorl %r14d, %r14d
movq 0x10(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x567e1
testl %eax, %eax
jne 0x77351
incq %r14
movq (%rbx), %rax
movslq 0x20(%rax), %rax
cmpq %rax, %r14
jl 0x76cb6
xorl %ecx, %ecx
jmp 0x77356
movq %rdi, %r15
shll $0xf, %eax
movzwl %ax, %ecx
movl %ecx, 0x410(%r14)
movl 0x224(%rdi), %eax
testb $0x2, %al
je 0x76d05
orl $0x1, %ecx
movl %ecx, 0x410(%r14)
testb $0x4, %al
je 0x76d14
orl $0x200001, 0x410(%r14) # imm = 0x200001
testb $0x8, %al
je 0x76d20
orb $0x3, 0x410(%r14)
testb $0x10, %al
je 0x76d2c
orb $0x20, 0x411(%r14)
testb $0x20, %al
je 0x76d3b
orl $0x3c00b1d, 0x410(%r14) # imm = 0x3C00B1D
testb $0x40, %al
je 0x76d4a
orl $0x310021, 0x410(%r14) # imm = 0x310021
testb %al, %al
jns 0x76d59
orl $0x120023, 0x410(%r14) # imm = 0x120023
btl $0x8, %eax
jae 0x76d6a
orl $0x3fc0b3f, 0x410(%r14) # imm = 0x3FC0B3F
btl $0x9, %eax
jae 0x76d7b
orl $0x7c00f5f, 0x410(%r14) # imm = 0x7C00F5F
btl $0xb, %eax
jb 0x76d8c
movl $0x0, 0x1f8(%r15)
btl $0xa, %eax
jae 0x76d9d
orl $0x13c00b9f, 0x410(%r14) # imm = 0x13C00B9F
btl $0xb, %eax
jae 0x76dae
orl $0x40005000, 0x410(%r14) # imm = 0x40005000
btl $0xc, %eax
jae 0x76dbf
orl $0x1001, 0x410(%r14) # imm = 0x1001
movq 0x10(%rbx), %rax
movq (%rax), %rdi
callq 0x567e1
movl $0xffffffff, %ecx       # imm = 0xFFFFFFFF
testl %eax, %eax
jne 0x77356
movq %r15, 0x18(%rsp)
movq (%rbx), %rax
movslq 0x20(%rax), %rdi
incq %rdi
movl $0x4, %esi
callq 0x6470
testq %rax, %rax
je 0x77351
movq %rax, %r12
xorl %ecx, %ecx
movl 0x410(%r14), %eax
btl $0x16, %eax
jae 0x76e18
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x17, %eax
jae 0x76e31
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x18, %eax
jae 0x76e4a
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
testb %al, %al
jns 0x76e61
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0xb, %eax
jae 0x76e7a
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x19, %eax
jae 0x76e93
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0xa, %eax
jae 0x76eac
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
testb $0x10, %al
je 0x76ec3
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
testb $0x40, %al
je 0x76eda
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x1a, %eax
jae 0x76ef3
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x1c, %eax
jae 0x76f0c
orl $0x104, %eax            # imm = 0x104
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
testl $0x7c00f5f, %eax        # imm = 0x7C00F5F
je 0x76f24
orl $0x8, %eax
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
testb $0x4, %al
je 0x76f3b
orl $0x100, %eax            # imm = 0x100
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x8, %eax
jae 0x76f54
orl $0x200, %eax            # imm = 0x200
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x1e, %eax
jae 0x76f6d
orl $0x4000, %eax           # imm = 0x4000
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0x13, %eax
jae 0x76f86
orl $0x100000, %eax         # imm = 0x100000
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0xd, %eax
jae 0x76f9d
orl $0x1, %eax
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
btl $0xa, %eax
jae 0x76fb6
orl $0x200000, %eax         # imm = 0x200000
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
testl $0xf0020, %eax          # imm = 0xF0020
je 0x76fd0
orl $0x100000, %eax         # imm = 0x100000
movl %eax, 0x410(%r14)
cmpl $0x0, 0x30(%r14)
jne 0x76fef
movl 0x410(%r14), %eax
testw %ax, %ax
jns 0x76fef
orl $0x100020, %eax         # imm = 0x100020
movl %eax, 0x410(%r14)
movq %rcx, 0x8(%rsp)
movl 0x410(%r14), %eax
testl $0x140000c0, %eax       # imm = 0x140000C0
je 0x7700e
orl $0x100009, %eax         # imm = 0x100009
movl %eax, 0x410(%r14)
movl 0x410(%r14), %eax
movl %eax, 0x14(%rsp)
xorl %r13d, %r13d
movl 0x410(%r14), %eax
btl %r13d, %eax
movl $0x9, %eax
jae 0x770cf
leaq 0x23ac7(%rip), %rcx     # 0x9ab00
movslq (%rcx,%r13,4), %rcx
movq 0x280(%r14,%rcx,8), %rdi
testq %rdi, %rdi
je 0x770cf
leaq 0x4(%rsp), %rsi
callq 0x8bc3d
movl %eax, %ebp
cmpl $-0x2, %ebp
je 0x770af
cmpl $-0x1, %ebp
jne 0x77070
movl $0x1, %eax
movq %rax, 0x8(%rsp)
jmp 0x770af
movq (%rbx), %rax
cmpl $0x0, 0x20(%rax)
jle 0x770af
xorl %r15d, %r15d
movq 0x10(%rbx), %rax
movq (%rax,%r15,8), %rdi
cmpl $0x4, 0x8(%rdi)
jne 0x770a0
cmpl %ebp, 0xc(%rdi)
jne 0x770a0
movl $0x1, (%r12,%r15,4)
callq 0x567e1
testl %eax, %eax
jne 0x770c2
incq %r15
movq (%rbx), %rax
movslq 0x20(%rax), %rax
cmpq %rax, %r15
jl 0x7707c
movl 0x4(%rsp), %ecx
xorl %eax, %eax
cmpl $-0x2, %ecx
je 0x770cf
cmpl %ecx, %ebp
movl %ecx, %ebp
jne 0x7705a
jmp 0x770cf
movq %r12, %rdi
callq 0x68a0
movl $0x1, %eax
cmpl $0x9, %eax
je 0x770dc
testl %eax, %eax
jne 0x77351
incq %r13
cmpq $0x1c, %r13
jne 0x7701c
movq 0x18(%rsp), %rax
testb $0x8, 0x225(%rax)
jne 0x7710b
testb $0x40, 0x413(%r14)
jne 0x7710b
movq 0x8(%rsp), %r13
jmp 0x771ff
xorl %eax, %eax
movq 0x8(%rsp), %r13
movq %rax, 0x20(%rsp)
movq 0x180(%r14,%rax,8), %rcx
testq %rcx, %rcx
je 0x771e3
movq 0x10(%rcx), %rdi
testq %rdi, %rdi
je 0x771af
movq %rcx, 0x8(%rsp)
leaq 0x4(%rsp), %rsi
callq 0x8bc3d
movl %eax, %ebp
cmpl $-0x2, %ebp
je 0x77193
cmpl $-0x1, %ebp
jne 0x77154
movl $0x1, %r13d
jmp 0x77193
movq (%rbx), %rax
cmpl $0x0, 0x20(%rax)
jle 0x77193
xorl %r15d, %r15d
movq 0x10(%rbx), %rax
movq (%rax,%r15,8), %rdi
cmpl $0x4, 0x8(%rdi)
jne 0x77184
cmpl %ebp, 0xc(%rdi)
jne 0x77184
movl $0x1, (%r12,%r15,4)
callq 0x567e1
testl %eax, %eax
jne 0x771b6
incq %r15
movq (%rbx), %rax
movslq 0x20(%rax), %rax
cmpq %rax, %r15
jl 0x77160
movl 0x4(%rsp), %eax
cmpl $-0x2, %eax
je 0x771a2
cmpl %eax, %ebp
movl %eax, %ebp
jne 0x77142
movq 0x8(%rsp), %rcx
movq 0x18(%rcx), %rcx
xorl %eax, %eax
jmp 0x771c8
movl $0x13, %eax
jmp 0x771c8
movq %r12, %rdi
callq 0x68a0
movl $0x1, %eax
movq 0x8(%rsp), %rcx
cmpl $0x13, %eax
je 0x771d1
testl %eax, %eax
jne 0x771dc
testq %rcx, %rcx
jne 0x77128
xorl %eax, %eax
testl %eax, %eax
sete %al
jmp 0x771e5
movb $0x1, %al
testb %al, %al
je 0x77351
movq 0x20(%rsp), %rax
incq %rax
cmpq $0x20, %rax
jne 0x77112
xorl %ebp, %ebp
leaq 0x238f8(%rip), %rax     # 0x9ab00
movslq (%rax,%rbp,4), %rax
movq 0x280(%r14,%rax,8), %rdi
testq %rdi, %rdi
je 0x7728d
leaq 0x4(%rsp), %rsi
callq 0x8bc3d
movl 0x4(%rsp), %edx
movl $0x1, %esi
movl %ebp, %ecx
shll %cl, %esi
cmpl $-0x2, %eax
je 0x77282
cmpl $-0x1, %eax
jne 0x77248
testl %r13d, %r13d
je 0x77282
orl %esi, 0x410(%r14)
jmp 0x77282
movq (%rbx), %rcx
movslq 0x20(%rcx), %rcx
testq %rcx, %rcx
jle 0x77282
movq 0x10(%rbx), %rdi
xorl %r8d, %r8d
movq (%rdi,%r8,8), %r9
cmpl $0x4, 0x8(%r9)
jne 0x7727a
cmpl %eax, 0xc(%r9)
jne 0x7727a
cmpl $0x0, (%r12,%r8,4)
je 0x7727a
orl %esi, 0x410(%r14)
incq %r8
cmpq %r8, %rcx
jne 0x7725b
cmpl $-0x2, %edx
je 0x7728d
cmpl %edx, %eax
movl %edx, %eax
jne 0x77230
incq %rbp
cmpq $0x1c, %rbp
jne 0x77201
movq %r13, %r15
xorl %r13d, %r13d
movq 0x180(%r14,%r13,8), %rbp
testq %rbp, %rbp
je 0x77323
movq 0x10(%rbp), %rdi
testq %rdi, %rdi
je 0x772a8
leaq 0x4(%rsp), %rsi
callq 0x8bc3d
movl 0x4(%rsp), %ecx
cmpl $-0x2, %eax
je 0x77312
cmpl $-0x1, %eax
jne 0x772d8
orb $0x40, 0x413(%r14)
jmp 0x77312
movq (%rbx), %rdx
movslq 0x20(%rdx), %rdx
testq %rdx, %rdx
jle 0x77312
movq 0x10(%rbx), %rsi
xorl %edi, %edi
movq (%rsi,%rdi,8), %r8
cmpl $0x0, 0x8(%r8)
je 0x7730a
cmpl %eax, 0xc(%r8)
jne 0x7730a
cmpl $0x0, (%r12,%rdi,4)
je 0x7730a
orb $0x40, 0x413(%r14)
incq %rdi
cmpq %rdi, %rdx
jne 0x772ea
cmpl $-0x2, %ecx
je 0x7731d
cmpl %ecx, %eax
movl %ecx, %eax
jne 0x772c4
movq 0x18(%rbp), %rbp
jmp 0x772a8
incq %r13
cmpq $0x20, %r13
jne 0x772a0
movl 0x14(%rsp), %eax
cmpl 0x410(%r14), %eax
movq %r15, %rcx
jne 0x76dff
movq %r12, %rdi
callq 0x68a0
jmp 0x76cda
movl $0xffffffff, %ecx       # imm = 0xFFFFFFFF
movl %ecx, %eax
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq