// File: /workspace/github_c_makefile_star2_O0/mohuangrui[P]ArtraCFD/convective_flux.c
// Function: CharacteristicVariable
static void CharacteristicVariable(const int tn, const int s, const int k, const int j,
        const int i, const int sL, const int sR, const int partn[restrict],
        const Node *const node, Real L[restrict][DIMU], Real W[restrict][DIMU])
{
    const int h[DIMS][DIMS] = {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}}; /* direction indicator */
    int idx = 0; /* linear array index math variable */
    const Real *restrict U = NULL;
    for (int n = sL, m = 0; n <= sR; ++n, ++m) {
        idx = IndexNode(k + n * h[s][Z], j + n * h[s][Y], i + n * h[s][X], partn[Y], partn[X]);
        U = node[idx].U[tn];
        for (int r = 0; r < DIMU; ++r) {
            W[m][r] = 0.0;
            for (int c = 0; c < DIMU; ++c) {
                W[m][r] = W[m][r] + L[r][c] * U[c];
            }
        }
    }
    return;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x60, %rsp
movq 0x30(%rbp), %rax
movq 0x28(%rbp), %rax
movq 0x20(%rbp), %rax
movq 0x18(%rbp), %rax
movl 0x10(%rbp), %eax
movl %edi, -0x4(%rbp)
movl %esi, -0x8(%rbp)
movl %edx, -0xc(%rbp)
movl %ecx, -0x10(%rbp)
movl %r8d, -0x14(%rbp)
movl %r9d, -0x18(%rbp)
leaq -0x40(%rbp), %rdi
xorl %esi, %esi
movl $0x24, %edx
callq 0x20d0
movl $0x1, -0x40(%rbp)
movl $0x1, -0x30(%rbp)
movl $0x1, -0x20(%rbp)
movl $0x0, -0x44(%rbp)
movq $0x0, -0x50(%rbp)
movl -0x18(%rbp), %eax
movl %eax, -0x54(%rbp)
movl $0x0, -0x58(%rbp)
movl -0x54(%rbp), %eax
cmpl 0x10(%rbp), %eax
jg 0x10b41
movl -0xc(%rbp), %edi
movl -0x54(%rbp), %eax
movslq -0x8(%rbp), %rdx
leaq -0x40(%rbp), %rcx
imulq $0xc, %rdx, %rdx
addq %rdx, %rcx
imull 0x8(%rcx), %eax
addl %eax, %edi
movl -0x10(%rbp), %esi
movl -0x54(%rbp), %eax
movslq -0x8(%rbp), %rdx
leaq -0x40(%rbp), %rcx
imulq $0xc, %rdx, %rdx
addq %rdx, %rcx
imull 0x4(%rcx), %eax
addl %eax, %esi
movl -0x14(%rbp), %edx
movl -0x54(%rbp), %eax
movslq -0x8(%rbp), %r8
leaq -0x40(%rbp), %rcx
imulq $0xc, %r8, %r8
addq %r8, %rcx
imull (%rcx), %eax
addl %eax, %edx
movq 0x18(%rbp), %rax
movl 0x4(%rax), %ecx
movq 0x18(%rbp), %rax
movl (%rax), %r8d
callq 0x9500
movl %eax, -0x44(%rbp)
movq 0x20(%rbp), %rax
movslq -0x44(%rbp), %rcx
imulq $0x88, %rcx, %rcx
addq %rcx, %rax
addq $0x10, %rax
movslq -0x4(%rbp), %rcx
imulq $0x28, %rcx, %rcx
addq %rcx, %rax
movq %rax, -0x50(%rbp)
movl $0x0, -0x5c(%rbp)
cmpl $0x5, -0x5c(%rbp)
jge 0x10b28
movq 0x30(%rbp), %rax
movslq -0x58(%rbp), %rcx
imulq $0x28, %rcx, %rcx
addq %rcx, %rax
movslq -0x5c(%rbp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movl $0x0, -0x60(%rbp)
cmpl $0x5, -0x60(%rbp)
jge 0x10b18
movq 0x30(%rbp), %rax
movslq -0x58(%rbp), %rcx
imulq $0x28, %rcx, %rcx
addq %rcx, %rax
movslq -0x5c(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x28(%rbp), %rax
movslq -0x5c(%rbp), %rcx
imulq $0x28, %rcx, %rcx
addq %rcx, %rax
movslq -0x60(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq -0x50(%rbp), %rax
movslq -0x60(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x30(%rbp), %rax
movslq -0x58(%rbp), %rcx
imulq $0x28, %rcx, %rcx
addq %rcx, %rax
movslq -0x5c(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl -0x60(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x60(%rbp)
jmp 0x10aaa
jmp 0x10b1a
movl -0x5c(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x5c(%rbp)
jmp 0x10a7e
jmp 0x10b2a
movl -0x54(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x54(%rbp)
movl -0x58(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x58(%rbp)
jmp 0x109e0
addq $0x60, %rsp
popq %rbp
retq
nopw (%rax,%rax)