// File: /workspace/github_c_makefile_star2_O3/icebert[P]pblat/htslib/hts.c
// Function: hts_hopen
htsFile *hts_hopen(struct hFILE *hfile, const char *fn, const char *mode)
{
    htsFile *fp = (htsFile*)calloc(1, sizeof(htsFile));
    char simple_mode[101], *cp, *opts;
    simple_mode[100] = '\0';

    if (fp == NULL) goto error;

    fp->fn = strdup(fn);
    fp->is_be = ed_is_big();

    // Split mode into simple_mode,opts strings
    if ((cp = strchr(mode, ','))) {
        strncpy(simple_mode, mode, cp-mode <= 100 ? cp-mode : 100);
        simple_mode[cp-mode] = '\0';
        opts = cp+1;
    } else {
        strncpy(simple_mode, mode, 100);
        opts = NULL;
    }

    if (strchr(simple_mode, 'r')) {
        if (hts_detect_format(hfile, &fp->format) < 0) goto error;
    }
    else if (strchr(simple_mode, 'w') || strchr(simple_mode, 'a')) {
        htsFormat *fmt = &fp->format;
        fp->is_write = 1;

        if (strchr(simple_mode, 'b')) fmt->format = binary_format;
        else if (strchr(simple_mode, 'c')) fmt->format = cram;
        else fmt->format = text_format;

        if (strchr(simple_mode, 'z')) fmt->compression = bgzf;
        else if (strchr(simple_mode, 'g')) fmt->compression = gzip;
        else if (strchr(simple_mode, 'u')) fmt->compression = no_compression;
        else {
            // No compression mode specified, set to the default for the format
            switch (fmt->format) {
            case binary_format: fmt->compression = bgzf; break;
            case cram: fmt->compression = custom; break;
            case text_format: fmt->compression = no_compression; break;
            default: abort();
            }
        }

        // Fill in category (if determinable; e.g. 'b' could be BAM or BCF)
        fmt->category = format_category(fmt->format);

        fmt->version.major = fmt->version.minor = -1;
        fmt->compression_level = -1;
        fmt->specific = NULL;
    }
    else goto error;

    switch (fp->format.format) {
    case binary_format:
    case bam:
    case bcf:
        fp->fp.bgzf = bgzf_hopen(hfile, simple_mode);
        if (fp->fp.bgzf == NULL) goto error;
        fp->is_bin = 1;
        break;

    case cram:
        fp->fp.cram = cram_dopen(hfile, fn, simple_mode);
        if (fp->fp.cram == NULL) goto error;
        if (!fp->is_write)
            cram_set_option(fp->fp.cram, CRAM_OPT_DECODE_MD, 1);
        fp->is_cram = 1;
        break;

    case text_format:
    case sam:
    case vcf:
        if (!fp->is_write) {
        #if KS_BGZF
            BGZF *gzfp = bgzf_hopen(hfile, simple_mode);
        #else
            // TODO Implement gzip hFILE adaptor
            hclose(hfile); // This won't work, especially for stdin
            gzFile gzfp = strcmp(fn, "-")? gzopen(fn, "rb") : gzdopen(fileno(stdin), "rb");
        #endif
            if (gzfp) fp->fp.voidp = ks_init(gzfp);
            else goto error;
        }
        else if (fp->format.compression != no_compression) {
            fp->fp.bgzf = bgzf_hopen(hfile, simple_mode);
            if (fp->fp.bgzf == NULL) goto error;
        }
        else
            fp->fp.hfile = hfile;
        break;

    default:
        goto error;
    }

    if (opts)
        hts_process_opts(fp, opts);

    return fp;

error:
    if (hts_verbose >= 2)
        fprintf(stderr, "[E::%s] fail to open file '%s'\n", __func__, fn);

    if (fp) {
        free(fp->fn);
        free(fp->fn_aux);
        free(fp);
    }
    return NULL;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x88, %rsp
movq %rdx, %r13
movq %rsi, %r15
movq %rdi, %r12
movl $0x1, %edi
movl $0x60, %esi
callq 0x6470
movq %rax, %rbx
movb $0x0, 0x64(%rsp)
testq %rax, %rax
je 0x4df60
movq %r15, %rdi
callq 0x6dd0
movq %rax, 0x28(%rbx)
movq %r13, %rdi
movl $0x2c, %esi
callq 0x6270
testq %rax, %rax
je 0x4dd3b
movq %rax, %r14
movq %rax, %rbp
subq %r13, %rbp
cmpq $0x64, %rbp
movl $0x64, %edx
cmovlq %rbp, %rdx
movq %rsp, %rdi
movq %r13, %rsi
callq 0x6380
movb $0x0, (%rsp,%rbp)
incq %r14
jmp 0x4dd4e
movq %rsp, %rdi
movl $0x64, %edx
movq %r13, %rsi
callq 0x6380
xorl %r14d, %r14d
movq %rsp, %rdi
movl $0x72, %esi
callq 0x6270
testq %rax, %rax
je 0x4dd7d
leaq 0x40(%rbx), %rsi
movq %r12, %rdi
callq 0x4ce35
testl %eax, %eax
js 0x4df60
movl 0x44(%rbx), %r13d
jmp 0x4de72
movq %rsp, %rdi
movl $0x77, %esi
callq 0x6270
testq %rax, %rax
jne 0x4dda5
movq %rsp, %rdi
movl $0x61, %esi
callq 0x6270
testq %rax, %rax
je 0x4df60
movl $0x2, (%rbx)
movq %rsp, %rdi
movl $0x62, %esi
callq 0x6270
testq %rax, %rax
je 0x4ddcc
movl $0x1, 0x44(%rbx)
movl $0x1, %r13d
jmp 0x4ddfa
movq %rsp, %rdi
movl $0x63, %esi
callq 0x6270
testq %rax, %rax
je 0x4dded
movl $0x6, 0x44(%rbx)
movl $0x6, %r13d
jmp 0x4ddfa
movl $0x2, 0x44(%rbx)
movl $0x2, %r13d
movq %rsp, %rdi
movl $0x7a, %esi
callq 0x6270
movl $0x2, %ecx
testq %rax, %rax
jne 0x4de40
movq %rsp, %rdi
movl $0x67, %esi
callq 0x6270
movl $0x1, %ecx
testq %rax, %rax
jne 0x4de40
movq %rsp, %rdi
movl $0x75, %esi
callq 0x6270
xorl %ecx, %ecx
testq %rax, %rax
je 0x4e002
movl %ecx, 0x4c(%rbx)
leal -0x3(%r13), %ecx
xorl %eax, %eax
cmpl $0x4, %ecx
jae 0x4de5a
movl %ecx, %eax
leaq 0x4a3f9(%rip), %rcx     # 0x98250
movl (%rcx,%rax,4), %eax
movl %eax, 0x40(%rbx)
movl $0xffffffff, 0x48(%rbx) # imm = 0xFFFFFFFF
movw $0xffff, 0x50(%rbx)     # imm = 0xFFFF
movq $0x0, 0x58(%rbx)
cmpl $0x9, %r13d
ja 0x4df60
movl $0x10c, %eax            # imm = 0x10C
btl %r13d, %eax
jb 0x4deb2
movl $0x212, %eax            # imm = 0x212
btl %r13d, %eax
jae 0x4df06
movq %rsp, %rsi
movq %r12, %rdi
callq 0x66e61
movq %rax, 0x38(%rbx)
testq %rax, %rax
je 0x4df60
orb $0x1, (%rbx)
jmp 0x4dfc7
testb $0x2, (%rbx)
jne 0x4df46
movq %rsp, %rsi
movq %r12, %rdi
callq 0x66e61
testq %rax, %rax
je 0x4df60
movq %rax, %r12
movl $0x1, %edi
movl $0x28, %esi
callq 0x6470
movq %rax, %r15
movq %r12, 0x18(%rax)
movl $0x40000, 0x8(%rax)     # imm = 0x40000
movl $0x10000, %edi          # imm = 0x10000
callq 0x69d0
movq %rax, 0x20(%r15)
movq %r15, 0x38(%rbx)
jmp 0x4dfc7
cmpl $0x6, %r13d
jne 0x4df60
movq %rsp, %rdx
movq %r12, %rdi
movq %r15, %rsi
callq 0x5c535
movq %rax, 0x38(%rbx)
testq %rax, %rax
je 0x4df60
movl (%rbx), %ebp
testb $0x2, %bpl
jne 0x4df3c
movq %rax, %rdi
xorl %esi, %esi
movl $0x1, %edx
xorl %eax, %eax
callq 0x5cfe5
orl $0x8, %ebp
movl %ebp, (%rbx)
jmp 0x4dfc7
cmpl $0x0, 0x4c(%rbx)
je 0x4dfc3
movq %rsp, %rsi
movq %r12, %rdi
callq 0x66e61
movq %rax, 0x38(%rbx)
testq %rax, %rax
jne 0x4dfc7
cmpl $0x2, 0x72e61(%rip)     # 0xc0dc8
jge 0x4df9f
testq %rbx, %rbx
je 0x4df88
movq 0x28(%rbx), %rdi
callq 0x68a0
movq 0x30(%rbx), %rdi
callq 0x68a0
movq %rbx, %rdi
callq 0x68a0
xorl %ebx, %ebx
movq %rbx, %rax
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq 0x71032(%rip), %rax     # 0xbefd8
movq (%rax), %rdi
leaq 0x4a554(%rip), %rsi     # 0x98504
leaq 0x4a6f5(%rip), %rdx     # 0x986ac
movq %r15, %rcx
xorl %eax, %eax
callq 0x6a90
jmp 0x4df69
movq %r12, 0x38(%rbx)
testq %r14, %r14
je 0x4df8a
leaq 0x68(%rsp), %rdi
movq $0x0, 0x18(%rdi)
movq %r14, %rsi
callq 0x4e5ad
testl %eax, %eax
jne 0x4df8a
movq 0x80(%rsp), %r14
movq %rbx, %rdi
movq %r14, %rsi
callq 0x4e03d
movq %r14, %rdi
callq 0x4e57d
jmp 0x4df8a
cmpl $0x1, %r13d
je 0x4e029
cmpl $0x2, %r13d
je 0x4e025
cmpl $0x6, %r13d
jne 0x4e038
movl $0x3, 0x4c(%rbx)
movl $0x1, %eax
jmp 0x4de5a
xorl %eax, %eax
jmp 0x4e02e
movl $0x2, %eax
movl %eax, 0x4c(%rbx)
xorl %eax, %eax
jmp 0x4de5a
callq 0x63c0