// File: /workspace/github_c_makefile_star2_O2/nh13[P]ksw/src/parasail/src/nw_striped_sse41_128_32.c
// Function: parasail_nw_rowcol_striped_profile_sse41_128_32
parasail_result_t* PNAME(
        const parasail_profile_t * const restrict profile,
        const char * const restrict s2, const int s2Len,
        const int open, const int gap)
{
    /* declare local variables */
    int32_t i = 0;
    int32_t j = 0;
    int32_t k = 0;
    int s1Len = 0;
    int32_t end_query = 0;
    int32_t end_ref = 0;
    const parasail_matrix_t *matrix = NULL;
    int32_t segWidth = 0;
    int32_t segLen = 0;
    int32_t offset = 0;
    int32_t position = 0;
    __m128i* restrict vProfile = NULL;
    __m128i* restrict pvHStore = NULL;
    __m128i* restrict pvHLoad = NULL;
    __m128i* restrict pvE = NULL;
    int32_t* restrict boundary = NULL;
    __m128i vGapO;
    __m128i vGapE;
    int32_t NEG_LIMIT = 0;
    int32_t POS_LIMIT = 0;
    int32_t score = 0;
    __m128i vNegLimit;
    __m128i vPosLimit;
    __m128i vSaturationCheckMin;
    __m128i vSaturationCheckMax;
    parasail_result_t *result = NULL;

    /* validate inputs */
    PARASAIL_CHECK_NULL(profile);
    PARASAIL_CHECK_NULL(profile->profile32.score);
    PARASAIL_CHECK_NULL(profile->matrix);
    PARASAIL_CHECK_GT0(profile->s1Len);
    PARASAIL_CHECK_NULL(s2);
    PARASAIL_CHECK_GT0(s2Len);
    PARASAIL_CHECK_GE0(open);
    PARASAIL_CHECK_GE0(gap);

    /* initialize stack variables */
    i = 0;
    j = 0;
    k = 0;
    s1Len = profile->s1Len;
    end_query = s1Len-1;
    end_ref = s2Len-1;
    matrix = profile->matrix;
    segWidth = 4; /* number of values in vector unit */
    segLen = (s1Len + segWidth - 1) / segWidth;
    offset = (s1Len - 1) % segLen;
    position = (segWidth - 1) - (s1Len - 1) / segLen;
    vProfile = (__m128i*)profile->profile32.score;
    vGapO = _mm_set1_epi32(open);
    vGapE = _mm_set1_epi32(gap);
    NEG_LIMIT = (-open < matrix->min ? INT32_MIN + open : INT32_MIN - matrix->min) + 1;
    POS_LIMIT = INT32_MAX - matrix->max - 1;
    score = NEG_LIMIT;
    vNegLimit = _mm_set1_epi32(NEG_LIMIT);
    vPosLimit = _mm_set1_epi32(POS_LIMIT);
    vSaturationCheckMin = vPosLimit;
    vSaturationCheckMax = vNegLimit;

    /* initialize result */
#ifdef PARASAIL_TABLE
    result = parasail_result_new_table1(segLen*segWidth, s2Len);
#else
#ifdef PARASAIL_ROWCOL
    result = parasail_result_new_rowcol1(segLen*segWidth, s2Len);
#else
    result = parasail_result_new();
#endif
#endif
    if (!result) return NULL;

    /* set known flags */
    result->flag |= PARASAIL_FLAG_NW | PARASAIL_FLAG_STRIPED
        | PARASAIL_FLAG_BITS_32 | PARASAIL_FLAG_LANES_4;
#ifdef PARASAIL_TABLE
    result->flag |= PARASAIL_FLAG_TABLE;
#endif
#ifdef PARASAIL_ROWCOL
    result->flag |= PARASAIL_FLAG_ROWCOL;
#endif

    /* initialize heap variables */
    pvHStore = parasail_memalign___m128i(16, segLen);
    pvHLoad  = parasail_memalign___m128i(16, segLen);
    pvE      = parasail_memalign___m128i(16, segLen);
    boundary = parasail_memalign_int32_t(16, s2Len+1);

    /* validate heap variables */
    if (!pvHStore) return NULL;
    if (!pvHLoad) return NULL;
    if (!pvE) return NULL;
    if (!boundary) return NULL;

    /* initialize H and E */
    {
        int32_t index = 0;
        for (i=0; i<segLen; ++i) {
            int32_t segNum = 0;
            __m128i_32_t h;
            __m128i_32_t e;
            for (segNum=0; segNum<segWidth; ++segNum) {
                int64_t tmp = -open-gap*(segNum*segLen+i);
                h.v[segNum] = tmp < INT32_MIN ? INT32_MIN : tmp;
                tmp = tmp - open;
                e.v[segNum] = tmp < INT32_MIN ? INT32_MIN : tmp;
            }
            _mm_store_si128(&pvHStore[index], h.m);
            _mm_store_si128(&pvE[index], e.m);
            ++index;
        }
    }

    /* initialize uppder boundary */
    {
        boundary[0] = 0;
        for (i=1; i<=s2Len; ++i) {
            int64_t tmp = -open-gap*(i-1);
            boundary[i] = tmp < INT32_MIN ? INT32_MIN : tmp;
        }
    }

    /* outer loop over database sequence */
    for (j=0; j<s2Len; ++j) {
        __m128i vE;
        /* Initialize F value to -inf.  Any errors to vH values will be
         * corrected in the Lazy_F loop.  */
        __m128i vF = vNegLimit;

        /* load final segment of pvHStore and shift left by 2 bytes */
        __m128i vH = _mm_slli_si128(pvHStore[segLen - 1], 4);

        /* Correct part of the vProfile */
        const __m128i* vP = vProfile + matrix->mapper[(unsigned char)s2[j]] * segLen;

        /* Swap the 2 H buffers. */
        __m128i* pv = pvHLoad;
        pvHLoad = pvHStore;
        pvHStore = pv;

        /* insert upper boundary condition */
        vH = _mm_insert_epi32(vH, boundary[j], 0);

        /* inner loop to process the query sequence */
        for (i=0; i<segLen; ++i) {
            vH = _mm_add_epi32(vH, _mm_load_si128(vP + i));
            vE = _mm_load_si128(pvE + i);

            /* Get max from vH, vE and vF. */
            vH = _mm_max_epi32(vH, vE);
            vH = _mm_max_epi32(vH, vF);
            /* Save vH values. */
            _mm_store_si128(pvHStore + i, vH);
            vSaturationCheckMax = _mm_max_epi32(vSaturationCheckMax, vH);
            vSaturationCheckMin = _mm_min_epi32(vSaturationCheckMin, vH);
            vSaturationCheckMin = _mm_min_epi32(vSaturationCheckMin, vE);
            vSaturationCheckMin = _mm_min_epi32(vSaturationCheckMin, vF);
#ifdef PARASAIL_TABLE
            arr_store_si128(result->tables->score_table, vH, i, segLen, j, s2Len);
#endif

            /* Update vE value. */
            vH = _mm_sub_epi32(vH, vGapO);
            vE = _mm_sub_epi32(vE, vGapE);
            vE = _mm_max_epi32(vE, vH);
            _mm_store_si128(pvE + i, vE);

            /* Update vF value. */
            vF = _mm_sub_epi32(vF, vGapE);
            vF = _mm_max_epi32(vF, vH);

            /* Load the next vH. */
            vH = _mm_load_si128(pvHLoad + i);
        }

        /* Lazy_F loop: has been revised to disallow adjecent insertion and
         * then deletion, so don't update E(i, i), learn from SWPS3 */
        for (k=0; k<segWidth; ++k) {
            int64_t tmp = boundary[j+1]-open;
            int32_t tmp2 = tmp < INT32_MIN ? INT32_MIN : tmp;
            vF = _mm_slli_si128(vF, 4);
            vF = _mm_insert_epi32(vF, tmp2, 0);
            for (i=0; i<segLen; ++i) {
                vH = _mm_load_si128(pvHStore + i);
                vH = _mm_max_epi32(vH,vF);
                _mm_store_si128(pvHStore + i, vH);
                vSaturationCheckMin = _mm_min_epi32(vSaturationCheckMin, vH);
                vSaturationCheckMax = _mm_max_epi32(vSaturationCheckMax, vH);
#ifdef PARASAIL_TABLE
                arr_store_si128(result->tables->score_table, vH, i, segLen, j, s2Len);
#endif
                vH = _mm_sub_epi32(vH, vGapO);
                vF = _mm_sub_epi32(vF, vGapE);
                if (! _mm_movemask_epi8(_mm_cmpgt_epi32(vF, vH))) goto end;
                /*vF = _mm_max_epi32(vF, vH);*/
            }
        }
end:
        {
        }

#ifdef PARASAIL_ROWCOL
        /* extract last value from the column */
        {
            vH = _mm_load_si128(pvHStore + offset);
            for (k=0; k<position; ++k) {
                vH = _mm_slli_si128(vH, 4);
            }
            result->rowcols->score_row[j] = (int32_t) _mm_extract_epi32 (vH, 3);
        }
#endif
    }

#ifdef PARASAIL_ROWCOL
    for (i=0; i<segLen; ++i) {
        __m128i vH = _mm_load_si128(pvHStore+i);
        arr_store_col(result->rowcols->score_col, vH, i, segLen);
    }
#endif

    /* extract last value from the last column */
    {
        __m128i vH = _mm_load_si128(pvHStore + offset);
        for (k=0; k<position; ++k) {
            vH = _mm_slli_si128 (vH, 4);
        }
        score = (int32_t) _mm_extract_epi32 (vH, 3);
    }

    if (_mm_movemask_epi8(_mm_or_si128(
            _mm_cmplt_epi32(vSaturationCheckMin, vNegLimit),
            _mm_cmpgt_epi32(vSaturationCheckMax, vPosLimit)))) {
        result->flag |= PARASAIL_FLAG_SATURATED;
        score = 0;
        end_query = 0;
        end_ref = 0;
    }

    result->score = score;
    result->end_query = end_query;
    result->end_ref = end_ref;

    parasail_free(boundary);
    parasail_free(pvE);
    parasail_free(pvHLoad);
    parasail_free(pvHStore);

    return result;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x88, %rsp
movl %ecx, 0xc(%rsp)
movq %rsi, 0x40(%rsp)
testq %rdi, %rdi
je 0x4010ab
movq 0x48(%rdi), %rax
movq %rax, 0x58(%rsp)
testq %rax, %rax
je 0x4010cf
movq 0x10(%rdi), %rax
movq %rax, 0x38(%rsp)
testq %rax, %rax
je 0x4010f3
movl 0x8(%rdi), %edi
testl %edi, %edi
jle 0x401117
cmpq $0x0, 0x40(%rsp)
je 0x40113b
movl %edx, %ebx
testl %edx, %edx
jle 0x40115c
cmpl $0x0, 0xc(%rsp)
js 0x40117d
movl %r8d, %ebp
testl %r8d, %r8d
js 0x40119e
leal -0x1(%rdi), %eax
addq $0x3, %rdi
movq %rdi, %r13
shrq $0x2, %r13
movl %eax, 0x24(%rsp)
xorl %edx, %edx
divl %r13d
movl %eax, %r12d
movl %edx, %r15d
movl 0xc(%rsp), %ecx
movl %ecx, %r14d
negl %r14d
movq 0x38(%rsp), %rdx
movl 0x20(%rdx), %eax
movl $0x80000000, %esi       # imm = 0x80000000
orl %esi, %ecx
subl %eax, %esi
cmpl %r14d, %eax
cmovgl %ecx, %esi
movl %esi, 0x2c(%rsp)
movl $0x7ffffffe, %eax       # imm = 0x7FFFFFFE
subl 0x1c(%rdx), %eax
movl %eax, 0x28(%rsp)
andl $0x7ffffffc, %edi       # imm = 0x7FFFFFFC
movl %ebx, %esi
callq 0x1ab93b
testq %rax, %rax
je 0x400d80
movl %r12d, 0x1c(%rsp)
movl %r15d, 0x20(%rsp)
movl %ebp, 0x10(%rsp)
movq %rax, 0x30(%rsp)
orl $0x4440801, 0xc(%rax)   # imm = 0x4440801
pushq $0x10
popq %rbp
movq %rbp, %rdi
movq %r13, %rsi
callq 0x297cd8
movq %rbx, %r15
movq %rax, %rbx
movq %rbp, %rdi
movq %r13, %rsi
callq 0x297cd8
movq %rax, %r12
movq %rbp, %rdi
movq %r13, %rsi
callq 0x297cd8
movq %rax, %rbp
leal 0x1(%r15), %esi
pushq $0x10
popq %rdi
movq %rsi, 0x50(%rsp)
callq 0x1ab783
testq %rbx, %rbx
sete %cl
testq %r12, %r12
sete %dl
orb %cl, %dl
testq %rbp, %rbp
sete %cl
testq %rax, %rax
sete %sil
orb %cl, %sil
orb %dl, %sil
je 0x400d87
xorl %ebx, %ebx
jmp 0x4011c6
movq %r15, 0x48(%rsp)
leal -0x1(%r15), %ecx
movl %ecx, 0x18(%rsp)
pushq $0x3
popq %rdx
subl 0x1c(%rsp), %edx
movl 0xc(%rsp), %ecx
movd %ecx, %xmm0
pshufd $0x0, %xmm0, %xmm2      # xmm2 = xmm0[0,0,0,0]
movl 0x10(%rsp), %r11d
movd %r11d, %xmm0
pshufd $0x0, %xmm0, %xmm3      # xmm3 = xmm0[0,0,0,0]
movl 0x2c(%rsp), %esi
incl %esi
movd %esi, %xmm0
pshufd $0x0, %xmm0, %xmm0      # xmm0 = xmm0[0,0,0,0]
movd 0x28(%rsp), %xmm1
pshufd $0x0, %xmm1, %xmm1      # xmm1 = xmm1[0,0,0,0]
movl %ecx, %esi
movl %r11d, %edi
imull %r13d, %edi
xorl %r8d, %r8d
movq $-0x80000000, %r9       # imm = 0x80000000
movl %r14d, %r10d
cmpq %r13, %r8
je 0x400e44
movl %r10d, %r15d
xorl %ecx, %ecx
cmpq $0x4, %rcx
je 0x400e1d
movslq %r15d, %r15
movq %r15, %r11
subq %rsi, %r11
cmpq $-0x7fffffff, %r11      # imm = 0x80000001
cmovlq %r9, %r11
movl %r15d, 0x70(%rsp,%rcx,4)
movl %r11d, 0x60(%rsp,%rcx,4)
incq %rcx
subl %edi, %r15d
jmp 0x400df1
movq %r8, %rcx
shlq $0x4, %rcx
movaps 0x70(%rsp), %xmm4
movaps %xmm4, (%rbx,%rcx)
movaps 0x60(%rsp), %xmm4
movaps %xmm4, (%rbp,%rcx)
incq %r8
movl 0x10(%rsp), %r11d
subl %r11d, %r10d
jmp 0x400de7
andl $0x0, (%rax)
pushq $0x1
popq %rcx
movl 0x20(%rsp), %esi
movq 0x50(%rsp), %rdi
cmpq %rcx, %rdi
je 0x400e64
movl %r14d, (%rax,%rcx,4)
incq %rcx
subl %r11d, %r14d
jmp 0x400e53
leal -0x1(%r13), %edi
movl %esi, %ecx
movl 0x48(%rsp), %r8d
movq %r13, %r9
shlq $0x4, %r9
leaq 0x10(%r9), %r10
xorl %r11d, %r11d
shlq $0x4, %rdi
shlq $0x4, %rcx
movq %rcx, 0x10(%rsp)
movdqa %xmm1, %xmm5
movdqa %xmm0, %xmm4
movq %rbx, %r15
cmpq %r8, %r11
je 0x400fe7
movq 0x38(%rsp), %rcx
movq 0x10(%rcx), %rcx
movq 0x40(%rsp), %rsi
movzbl (%rsi,%r11), %ebx
movl (%rcx,%rbx,4), %ecx
imull %r13d, %ecx
movslq %ecx, %rbx
shlq $0x4, %rbx
addq 0x58(%rsp), %rbx
pshufd $0x90, (%r15,%rdi), %xmm7 # xmm7 = mem[0,0,1,2]
pinsrd $0x0, (%rax,%r11,4), %xmm7
xorl %r14d, %r14d
movdqa %xmm0, %xmm6
cmpq %r14, %r9
je 0x400f47
paddd (%rbx,%r14), %xmm7
movdqa (%rbp,%r14), %xmm8
movdqa %xmm8, %xmm9
pmaxsd %xmm6, %xmm9
pmaxsd %xmm7, %xmm9
movdqa %xmm9, (%r12,%r14)
pmaxsd %xmm9, %xmm4
pminsd %xmm8, %xmm5
pminsd %xmm6, %xmm5
pminsd %xmm9, %xmm5
psubd %xmm2, %xmm9
psubd %xmm3, %xmm8
pmaxsd %xmm9, %xmm8
movdqa %xmm8, (%rbp,%r14)
psubd %xmm3, %xmm6
pmaxsd %xmm9, %xmm6
movdqa (%r15,%r14), %xmm7
addq $0x10, %r14
jmp 0x400eda
xorl %ebx, %ebx
cmpl $0x4, %ebx
je 0x400fab
movl 0x4(%rax,%r11,4), %ecx
subl 0xc(%rsp), %ecx
pshufd $0x90, %xmm6, %xmm6     # xmm6 = xmm6[0,0,1,2]
pinsrd $0x0, %ecx, %xmm6
xorl %ecx, %ecx
leaq 0x10(%rcx), %r14
cmpq %r14, %r10
je 0x400fa7
movdqa (%r12,%rcx), %xmm7
pmaxsd %xmm6, %xmm7
movdqa %xmm7, (%r12,%rcx)
pminsd %xmm7, %xmm5
pmaxsd %xmm7, %xmm4
psubd %xmm2, %xmm7
psubd %xmm3, %xmm6
movdqa %xmm6, %xmm8
pcmpgtd %xmm7, %xmm8
movmskps %xmm8, %esi
movq %r14, %rcx
testl %esi, %esi
jne 0x400f64
jmp 0x400fab
incl %ebx
jmp 0x400f49
movq 0x10(%rsp), %rcx
movdqa (%r12,%rcx), %xmm6
xorl %ecx, %ecx
movq 0x30(%rsp), %rsi
cmpl %edx, %ecx
jge 0x400fca
pslldq $0x4, %xmm6             # xmm6 = zero,zero,zero,zero,xmm6[0,1,2,3,4,5,6,7,8,9,10,11]
incl %ecx
jmp 0x400fbd
movq 0x10(%rsi), %rcx
movq (%rcx), %rcx
pextrd $0x3, %xmm6, (%rcx,%r11,4)
incq %r11
movq %r12, %rbx
movq %r15, %r12
jmp 0x400e92
imulq $0xc, %r13, %rdi
leaq (,%r13,8), %r8
shlq $0x2, %r13
xorl %ecx, %ecx
movq 0x30(%rsp), %rbx
cmpq %rcx, %r13
je 0x401037
movdqa (%r15,%rcx,4), %xmm2
movq 0x10(%rbx), %rsi
movq 0x8(%rsi), %rsi
movd %xmm2, (%rsi,%rcx)
addq %rcx, %rsi
pextrd $0x1, %xmm2, (%r13,%rsi)
pextrd $0x2, %xmm2, (%r8,%rsi)
pextrd $0x3, %xmm2, (%rdi,%rsi)
addq $0x4, %rcx
jmp 0x400ffe
movq 0x10(%rsp), %rcx
movdqa (%r15,%rcx), %xmm2
xorl %ecx, %ecx
movl 0x18(%rsp), %esi
cmpl %edx, %ecx
jge 0x401055
pslldq $0x4, %xmm2             # xmm2 = zero,zero,zero,zero,xmm2[0,1,2,3,4,5,6,7,8,9,10,11]
incl %ecx
jmp 0x401048
pcmpgtd %xmm5, %xmm0
pcmpgtd %xmm1, %xmm4
por %xmm0, %xmm4
movmskps %xmm4, %ecx
testl %ecx, %ecx
je 0x401074
orb $0x40, 0xc(%rbx)
xorl %ecx, %ecx
xorl %esi, %esi
xorl %edx, %edx
jmp 0x40107e
pextrd $0x3, %xmm2, %ecx
movl 0x24(%rsp), %edx
movl %ecx, (%rbx)
movl %edx, 0x4(%rbx)
movl %esi, 0x8(%rbx)
movq %rax, %rdi
callq 0x1ab763
movq %rbp, %rdi
callq 0x1ab763
movq %r12, %rdi
callq 0x1ab763
movq %r15, %rdi
callq 0x1ab763
jmp 0x4011c6
movq 0x339f26(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c5f2f(%rip), %rsi    # 0x5c6feb
leaq 0x1b39a9(%rip), %rdx    # 0x5b4a6c
leaq 0x19c413(%rip), %rcx    # 0x59d4dd
jmp 0x4011bd
movq 0x339f02(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c5f0b(%rip), %rsi    # 0x5c6feb
leaq 0x1b3985(%rip), %rdx    # 0x5b4a6c
leaq 0x1c68b8(%rip), %rcx    # 0x5c79a6
jmp 0x4011bd
movq 0x339ede(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c5ee7(%rip), %rsi    # 0x5c6feb
leaq 0x1b3961(%rip), %rdx    # 0x5b4a6c
leaq 0x1c6827(%rip), %rcx    # 0x5c7939
jmp 0x4011bd
movq 0x339eba(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c65d9(%rip), %rsi    # 0x5c7701
leaq 0x1b393d(%rip), %rdx    # 0x5b4a6c
leaq 0x1c6813(%rip), %rcx    # 0x5c7949
jmp 0x4011bd
movq 0x339e96(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c5e9f(%rip), %rsi    # 0x5c6feb
leaq 0x1b3919(%rip), %rdx    # 0x5b4a6c
leaq 0x1c65a4(%rip), %rcx    # 0x5c76fe
jmp 0x4011bd
movq 0x339e75(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c6594(%rip), %rsi    # 0x5c7701
leaq 0x1b38f8(%rip), %rdx    # 0x5b4a6c
leaq 0x1c659a(%rip), %rcx    # 0x5c7715
jmp 0x4011bd
movq 0x339e54(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c658d(%rip), %rsi    # 0x5c771b
leaq 0x1b38d7(%rip), %rdx    # 0x5b4a6c
leaq 0x1c6594(%rip), %rcx    # 0x5c7730
jmp 0x4011bd
movq 0x339e33(%rip), %rax    # 0x73afd8
movq (%rax), %rdi
leaq 0x1c656c(%rip), %rsi    # 0x5c771b
leaq 0x1b38b6(%rip), %rdx    # 0x5b4a6c
leaq 0x1c6578(%rip), %rcx    # 0x5c7735
xorl %ebx, %ebx
xorl %eax, %eax
callq 0x19e180
movq %rbx, %rax
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop