// File: /workspace/github_cpp_makefile_star2_O3/karel-brinda[P]ococo/ext/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 0x4360
movq %rax, %rbx
movb $0x0, 0x64(%rsp)
testq %rax, %rax
je 0xf4bd
movq %r15, %rdi
callq 0x49c0
movq %rax, 0x28(%rbx)
movq %r13, %rdi
movl $0x2c, %esi
callq 0x4200
testq %rax, %rax
je 0xf298
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 0x42b0
movb $0x0, (%rsp,%rbp)
incq %r14
jmp 0xf2ab
movq %rsp, %rdi
movl $0x64, %edx
movq %r13, %rsi
callq 0x42b0
xorl %r14d, %r14d
movq %rsp, %rdi
movl $0x72, %esi
callq 0x4200
testq %rax, %rax
je 0xf2da
leaq 0x40(%rbx), %rsi
movq %r12, %rdi
callq 0xe2bd
testl %eax, %eax
js 0xf4bd
movl 0x44(%rbx), %r13d
jmp 0xf3cf
movq %rsp, %rdi
movl $0x77, %esi
callq 0x4200
testq %rax, %rax
jne 0xf302
movq %rsp, %rdi
movl $0x61, %esi
callq 0x4200
testq %rax, %rax
je 0xf4bd
movl $0x2, (%rbx)
movq %rsp, %rdi
movl $0x62, %esi
callq 0x4200
testq %rax, %rax
je 0xf329
movl $0x1, 0x44(%rbx)
movl $0x1, %r13d
jmp 0xf357
movq %rsp, %rdi
movl $0x63, %esi
callq 0x4200
testq %rax, %rax
je 0xf34a
movl $0x6, 0x44(%rbx)
movl $0x6, %r13d
jmp 0xf357
movl $0x2, 0x44(%rbx)
movl $0x2, %r13d
movq %rsp, %rdi
movl $0x7a, %esi
callq 0x4200
movl $0x2, %ecx
testq %rax, %rax
jne 0xf39d
movq %rsp, %rdi
movl $0x67, %esi
callq 0x4200
movl $0x1, %ecx
testq %rax, %rax
jne 0xf39d
movq %rsp, %rdi
movl $0x75, %esi
callq 0x4200
xorl %ecx, %ecx
testq %rax, %rax
je 0xf55f
movl %ecx, 0x4c(%rbx)
leal -0x3(%r13), %ecx
xorl %eax, %eax
cmpl $0x4, %ecx
jae 0xf3b7
movl %ecx, %eax
leaq 0x415ec(%rip), %rcx     # 0x509a0
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 0xf4bd
movl $0x10c, %eax            # imm = 0x10C
btl %r13d, %eax
jb 0xf40f
movl $0x212, %eax            # imm = 0x212
btl %r13d, %eax
jae 0xf463
movq %rsp, %rsi
movq %r12, %rdi
callq 0x42b9b
movq %rax, 0x38(%rbx)
testq %rax, %rax
je 0xf4bd
orb $0x1, (%rbx)
jmp 0xf524
testb $0x2, (%rbx)
jne 0xf4a3
movq %rsp, %rsi
movq %r12, %rdi
callq 0x42b9b
testq %rax, %rax
je 0xf4bd
movq %rax, %r12
movl $0x1, %edi
movl $0x28, %esi
callq 0x4360
movq %rax, %r15
movq %r12, 0x18(%rax)
movl $0x40000, 0x8(%rax)     # imm = 0x40000
movl $0x10000, %edi          # imm = 0x10000
callq 0x4700
movq %rax, 0x20(%r15)
movq %r15, 0x38(%rbx)
jmp 0xf524
cmpl $0x6, %r13d
jne 0xf4bd
movq %rsp, %rdx
movq %r12, %rdi
movq %r15, %rsi
callq 0x37e32
movq %rax, 0x38(%rbx)
testq %rax, %rax
je 0xf4bd
movl (%rbx), %ebp
testb $0x2, %bpl
jne 0xf499
movq %rax, %rdi
xorl %esi, %esi
movl $0x1, %edx
xorl %eax, %eax
callq 0x388fc
orl $0x8, %ebp
movl %ebp, (%rbx)
jmp 0xf524
cmpl $0x0, 0x4c(%rbx)
je 0xf520
movq %rsp, %rsi
movq %r12, %rdi
callq 0x42b9b
movq %rax, 0x38(%rbx)
testq %rax, %rax
jne 0xf524
cmpl $0x2, 0x5308c(%rip)     # 0x62550
jge 0xf4fc
testq %rbx, %rbx
je 0xf4e5
movq 0x28(%rbx), %rdi
callq 0x4610
movq 0x30(%rbx), %rdi
callq 0x4610
movq %rbx, %rdi
callq 0x4610
xorl %ebx, %ebx
movq %rbx, %rax
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq 0x52ac5(%rip), %rax     # 0x61fc8
movq (%rax), %rdi
leaq 0x41761(%rip), %rsi     # 0x50c6e
leaq 0x41932(%rip), %rdx     # 0x50e46
movq %r15, %rcx
xorl %eax, %eax
callq 0x4770
jmp 0xf4c6
movq %r12, 0x38(%rbx)
testq %r14, %r14
je 0xf4e7
leaq 0x68(%rsp), %rdi
movq $0x0, 0x18(%rdi)
movq %r14, %rsi
callq 0xfb7b
testl %eax, %eax
jne 0xf4e7
movq 0x80(%rsp), %r14
movq %rbx, %rdi
movq %r14, %rsi
callq 0xf59a
movq %r14, %rdi
callq 0xfb4b
jmp 0xf4e7
cmpl $0x1, %r13d
je 0xf586
cmpl $0x2, %r13d
je 0xf582
cmpl $0x6, %r13d
jne 0xf595
movl $0x3, 0x4c(%rbx)
movl $0x1, %eax
jmp 0xf3b7
xorl %eax, %eax
jmp 0xf58b
movl $0x2, %eax
movl %eax, 0x4c(%rbx)
xorl %eax, %eax
jmp 0xf3b7
callq 0x42e0