void get_face_2d_upmulti(int iX, int iY, int nX, int nY, int X_nbr, int Y_nbr,
                         int nrm, int tid, double delta_t, double u_face,
                         double u_face_t, double *sc_face, double *sc_old) {
  int offset = 1;
  int iwind_t = 1;
  int nbrid;
  double wind_t = 1.;
  double abs_crnt_t;
  double sc_face_tmp;
  if (u_face * nrm > 0.)
    offset = 0;
  if (u_face_t < 0.) {
    wind_t = -1.;
    iwind_t = -1;
  }
  abs_crnt_t = wind_t * u_face_t * delta_t * 1.;
  nbrid = (iY + offset * Y_nbr) * nX + iX + offset * X_nbr;
  sc_face_tmp = *(sc_old + nbrid) * (1. - 0.5 * abs_crnt_t);
  nbrid = (iY + offset * Y_nbr - iwind_t * (1 - Y_nbr * nrm)) * nX + iX +
          offset * X_nbr - iwind_t * (1 - X_nbr * nrm);
  sc_face_tmp += *(sc_old + nbrid) * (0.5 * abs_crnt_t);
  *sc_face = sc_face_tmp;
}


x86 assembly (clang, -O0):
<get_face_2d_upmulti>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9d,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  %xmm2,-0x30(%rbp)
movl   $0x1,-0x34(%rbp)
movl   $0x1,-0x38(%rbp)
movsd  0xeb5(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x28(%rbp),%xmm0
cvtsi2sdl 0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1172 <get_face_2d_upmulti+0x72>
movl   $0x0,-0x34(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    1194 <get_face_2d_upmulti+0x94>
movsd  0xe80(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0xffffffff,-0x38(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xe55(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x34(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x3c(%rbp)
mov    0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm1
movsd  0xe2b(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xe07(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x34(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
mov    -0x38(%rbp),%ecx
mov    -0x18(%rbp),%esi
imul   0x10(%rbp),%esi
mov    $0x1,%edx
sub    %esi,%edx
imul   %edx,%ecx
sub    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    -0x38(%rbp),%ecx
mov    -0x14(%rbp),%esi
imul   0x10(%rbp),%esi
mov    $0x1,%edx
sub    %esi,%edx
imul   %edx,%ecx
sub    %ecx,%eax
mov    %eax,-0x3c(%rbp)
mov    0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xdaa(%rip),%xmm2
mulsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_face_2d_upmulti>:
push   %rbp
push   %rbx
mov    0x28(%rsp),%r10
mov    0x18(%rsp),%ecx
cvtsi2sd %ecx,%xmm3
xorpd  %xmm4,%xmm4
xor    %ebx,%ebx
ucomisd %xmm2,%xmm4
setbe  %bl
mulsd  %xmm1,%xmm3
xor    %eax,%eax
ucomisd %xmm4,%xmm3
mov    0x30(%rsp),%r11
lea    -0x1(%rbx,%rbx,1),%ebx
movapd 0xec9(%rip),%xmm3
xorpd  %xmm2,%xmm3
movapd %xmm2,%xmm1
cmpltsd %xmm4,%xmm1
andpd  %xmm1,%xmm3
andnpd %xmm2,%xmm1
orpd   %xmm3,%xmm1
mov    $0x0,%ebp
cmovbe %r9d,%ebp
mulsd  %xmm0,%xmm1
cmovbe %r8d,%eax
add    %esi,%ebp
imul   %ecx,%r9d
add    $0xffffffff,%r9d
imul   %ebx,%r9d
add    %ebp,%r9d
imul   %edx,%ebp
add    %edi,%eax
add    %eax,%ebp
movsd  0xe8f(%rip),%xmm0
mulsd  %xmm1,%xmm0
addsd  0xe8b(%rip),%xmm0
movslq %ebp,%rsi
mulsd  (%r11,%rsi,8),%xmm0
imul   %edx,%r9d
imul   %r8d,%ecx
add    $0xffffffff,%ecx
imul   %ebx,%ecx
add    %eax,%ecx
add    %r9d,%ecx
movslq %ecx,%rax
mulsd  0xe6c(%rip),%xmm1
mulsd  (%r11,%rax,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r10)
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define void @get_face_2d_upmulti(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, double %8, double %9, double %10, double* %11, double* %12) {
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double*
%26 = alloca double*
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca double
%31 = alloca double
%32 = alloca double
store i32 %0, i32* %14
store i32 %1, i32* %15
store i32 %2, i32* %16
store i32 %3, i32* %17
store i32 %4, i32* %18
store i32 %5, i32* %19
store i32 %6, i32* %20
store i32 %7, i32* %21
store double %8, double* %22
store double %9, double* %23
store double %10, double* %24
store double* %11, double** %25
store double* %12, double** %26
store i32 1, i32* %27
store i32 1, i32* %28
store double 1.000000e+00, double* %30
%33 = load double, double* %23
%34 = load i32, i32* %20
%35 = sitofp i32 %34 to double
%36 = fmul double %33, %35
%37 = fcmp ogt double %36, 0.000000e+00
br i1 %37, label %38, label %39
38: ; preds = %13
store i32 0, i32* %27
br label %39
39: ; preds = %38, %13
%40 = load double, double* %24
%41 = fcmp olt double %40, 0.000000e+00
br i1 %41, label %42, label %43
42: ; preds = %39
store double -1.000000e+00, double* %30
store i32 -1, i32* %28
br label %43
43: ; preds = %42, %39
%44 = load double, double* %30
%45 = load double, double* %24
%46 = fmul double %44, %45
%47 = load double, double* %22
%48 = fmul double %46, %47
%49 = fmul double %48, 1.000000e+00
store double %49, double* %31
%50 = load i32, i32* %15
%51 = load i32, i32* %27
%52 = load i32, i32* %19
%53 = mul nsw i32 %51, %52
%54 = add nsw i32 %50, %53
%55 = load i32, i32* %16
%56 = mul nsw i32 %54, %55
%57 = load i32, i32* %14
%58 = add nsw i32 %56, %57
%59 = load i32, i32* %27
%60 = load i32, i32* %18
%61 = mul nsw i32 %59, %60
%62 = add nsw i32 %58, %61
store i32 %62, i32* %29
%63 = load double*, double** %26
%64 = load i32, i32* %29
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = load double, double* %31
%69 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %68, double 1.000000e+00)
%70 = fmul double %67, %69
store double %70, double* %32
%71 = load i32, i32* %15
%72 = load i32, i32* %27
%73 = load i32, i32* %19
%74 = mul nsw i32 %72, %73
%75 = add nsw i32 %71, %74
%76 = load i32, i32* %28
%77 = load i32, i32* %19
%78 = load i32, i32* %20
%79 = mul nsw i32 %77, %78
%80 = sub nsw i32 1, %79
%81 = mul nsw i32 %76, %80
%82 = sub nsw i32 %75, %81
%83 = load i32, i32* %16
%84 = mul nsw i32 %82, %83
%85 = load i32, i32* %14
%86 = add nsw i32 %84, %85
%87 = load i32, i32* %27
%88 = load i32, i32* %18
%89 = mul nsw i32 %87, %88
%90 = add nsw i32 %86, %89
%91 = load i32, i32* %28
%92 = load i32, i32* %18
%93 = load i32, i32* %20
%94 = mul nsw i32 %92, %93
%95 = sub nsw i32 1, %94
%96 = mul nsw i32 %91, %95
%97 = sub nsw i32 %90, %96
store i32 %97, i32* %29
%98 = load double*, double** %26
%99 = load i32, i32* %29
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %98, i64 %100
%102 = load double, double* %101
%103 = load double, double* %31
%104 = fmul double 5.000000e-01, %103
%105 = load double, double* %32
%106 = call double @llvm.fmuladd.f64(double %102, double %104, double %105)
store double %106, double* %32
%107 = load double, double* %32
%108 = load double*, double** %25
store double %107, double* %108
ret void
}

LLVM IR (-O2):
define void @get_face_2d_upmulti(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, double %8, double %9, double %10, double* nocapture writeonly %11, double* nocapture %12) {
%14 = sitofp i32 %6 to double
%15 = fmul double %14, %9
%16 = fcmp ule double %15, 0.000000e+00
%17 = fcmp olt double %10, 0.000000e+00
%18 = select i1 %17, i32 -1, i32 1
%19 = fneg double %10
%20 = select i1 %17, double %19, double %10
%21 = fmul double %20, %8
%22 = select i1 %16, i32 %5, i32 0
%23 = add nsw i32 %22, %1
%24 = mul nsw i32 %23, %2
%25 = select i1 %16, i32 %4, i32 0
%26 = add i32 %25, %0
%27 = add i32 %26, %24
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %12, i64 %28
%30 = load double, double* %29
%31 = tail call double @llvm.fmuladd.f64(double %21, double -5.000000e-01, double 1.000000e+00)
%32 = fmul double %31, %30
%33 = mul nsw i32 %6, %5
%34 = add i32 %33, -1
%35 = mul i32 %18, %34
%36 = add i32 %35, %23
%37 = mul nsw i32 %36, %2
%38 = mul nsw i32 %6, %4
%39 = add i32 %38, -1
%40 = mul i32 %18, %39
%41 = add i32 %26, %40
%42 = add i32 %41, %37
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %12, i64 %43
%45 = load double, double* %44
%46 = fmul double %21, 5.000000e-01
%47 = tail call double @llvm.fmuladd.f64(double %45, double %46, double %32)
store double %47, double* %11
ret void
}