void gbl_wt(double *sum, double *wt, float *data, int i, int j, int nx, int ny,
            double wt0) {
  float t;
  i = (i == -1) ? nx - 1 : i;
  i = (i == nx) ? 0 : i;
  if (i < 0 || i >= nx || j < 0 || j >= ny)
    return;
  t = data[i + j * nx];
  if (((t) >= 9.9989e20 && (t) <= 9.9991e20))
    return;
  *wt = *wt + wt0;
  *sum = *sum + t * wt0;
  return;
}


x86 assembly (clang, -O0):
<gbl_wt>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movsd  %xmm0,-0x30(%rbp)
cmpl   $0xffffffff,-0x1c(%rbp)
jne    113b <gbl_wt+0x3b>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1141 <gbl_wt+0x41>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jne    115d <gbl_wt+0x5d>
xor    %eax,%eax
mov    %eax,-0x3c(%rbp)
jmp    1163 <gbl_wt+0x63>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jl     1195 <gbl_wt+0x95>
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1195 <gbl_wt+0x95>
cmpl   $0x0,-0x20(%rbp)
jl     1195 <gbl_wt+0x95>
mov    -0x20(%rbp),%eax
cmp    0x10(%rbp),%eax
jl     119a <gbl_wt+0x9a>
jmp    122d <gbl_wt+0x12d>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe38(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     11f2 <gbl_wt+0xf2>
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe25(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     11f2 <gbl_wt+0xf2>
jmp    122d <gbl_wt+0x12d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
addsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm1
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<gbl_wt>:
lea    -0x1(%r9),%eax
cmp    $0xffffffff,%ecx
cmovne %ecx,%eax
xor    %ecx,%ecx
cmp    %r9d,%eax
cmovne %eax,%ecx
test   %ecx,%ecx
js     1152 <gbl_wt+0x52>
cmp    0x8(%rsp),%r8d
jge    1152 <gbl_wt+0x52>
test   %r8d,%r8d
js     1152 <gbl_wt+0x52>
cmp    %r9d,%ecx
jge    1152 <gbl_wt+0x52>
imul   %r8d,%r9d
add    %r9d,%ecx
movslq %ecx,%rax
movss  (%rdx,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
ucomisd 0xebe(%rip),%xmm1
jb     1153 <gbl_wt+0x53>
movsd  0xebc(%rip),%xmm2
ucomisd %xmm1,%xmm2
jb     1153 <gbl_wt+0x53>
ret
mulsd  %xmm0,%xmm1
addsd  (%rsi),%xmm0
movsd  %xmm0,(%rsi)
addsd  (%rdi),%xmm1
movsd  %xmm1,(%rdi)
ret

LLVM IR (-O0):
define void @gbl_wt(double* %0, double* %1, float* %2, i32 %3, i32 %4, i32 %5, i32 %6, double %7) {
%9 = alloca double*
%10 = alloca double*
%11 = alloca float*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca float
store double* %0, double** %9
store double* %1, double** %10
store float* %2, float** %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store double %7, double* %16
%18 = load i32, i32* %12
%19 = icmp eq i32 %18, -1
br i1 %19, label %20, label %23
20: ; preds = %8
%21 = load i32, i32* %14
%22 = sub nsw i32 %21, 1
br label %25
23: ; preds = %8
%24 = load i32, i32* %12
br label %25
25: ; preds = %23, %20
%26 = phi i32 [ %22, %20 ], [ %24, %23 ]
store i32 %26, i32* %12
%27 = load i32, i32* %12
%28 = load i32, i32* %14
%29 = icmp eq i32 %27, %28
br i1 %29, label %30, label %31
30: ; preds = %25
br label %33
31: ; preds = %25
%32 = load i32, i32* %12
br label %33
33: ; preds = %31, %30
%34 = phi i32 [ 0, %30 ], [ %32, %31 ]
store i32 %34, i32* %12
%35 = load i32, i32* %12
%36 = icmp slt i32 %35, 0
br i1 %36, label %48, label %37
37: ; preds = %33
%38 = load i32, i32* %12
%39 = load i32, i32* %14
%40 = icmp sge i32 %38, %39
br i1 %40, label %48, label %41
41: ; preds = %37
%42 = load i32, i32* %13
%43 = icmp slt i32 %42, 0
br i1 %43, label %48, label %44
44: ; preds = %41
%45 = load i32, i32* %13
%46 = load i32, i32* %15
%47 = icmp sge i32 %45, %46
br i1 %47, label %48, label %49
48: ; preds = %44, %41, %37, %33
br label %80
49: ; preds = %44
%50 = load float*, float** %11
%51 = load i32, i32* %12
%52 = load i32, i32* %13
%53 = load i32, i32* %14
%54 = mul nsw i32 %52, %53
%55 = add nsw i32 %51, %54
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds float, float* %50, i64 %56
%58 = load float, float* %57
store float %58, float* %17
%59 = load float, float* %17
%60 = fpext float %59 to double
%61 = fcmp oge double %60, 0x444B1A2170AD88AA
br i1 %61, label %62, label %67
62: ; preds = %49
%63 = load float, float* %17
%64 = fpext float %63 to double
%65 = fcmp ole double %64, 9.999100e+20
br i1 %65, label %66, label %67
66: ; preds = %62
br label %80
67: ; preds = %62, %49
%68 = load double*, double** %10
%69 = load double, double* %68
%70 = load double, double* %16
%71 = fadd double %69, %70
%72 = load double*, double** %10
store double %71, double* %72
%73 = load double*, double** %9
%74 = load double, double* %73
%75 = load float, float* %17
%76 = fpext float %75 to double
%77 = load double, double* %16
%78 = call double @llvm.fmuladd.f64(double %76, double %77, double %74)
%79 = load double*, double** %9
store double %78, double* %79
br label %80
80: ; preds = %67, %66, %48
ret void
}

LLVM IR (-O2):
define void @gbl_wt(double* nocapture %0, double* nocapture %1, float* nocapture %2, i32 %3, i32 %4, i32 %5, i32 %6, double %7) {
%9 = icmp eq i32 %3, -1
%10 = add nsw i32 %5, -1
%11 = select i1 %9, i32 %10, i32 %3
%12 = icmp eq i32 %11, %5
%13 = select i1 %12, i32 0, i32 %11
%14 = icmp slt i32 %13, 0
br i1 %14, label %36, label %15
15: ; preds = %8
%16 = icmp slt i32 %13, %5
%17 = icmp sgt i32 %4, -1
%18 = and i1 %17, %16
%19 = icmp slt i32 %4, %6
%20 = and i1 %19, %18
br i1 %20, label %21, label %36
21: ; preds = %15
%22 = mul nsw i32 %5, %4
%23 = add nsw i32 %13, %22
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds float, float* %2, i64 %24
%26 = load float, float* %25
%27 = fpext float %26 to double
%28 = fcmp ult double %27, 0x444B1A2170AD88AA
%29 = fcmp ugt double %27, 9.999100e+20
%30 = or i1 %28, %29
br i1 %30, label %31, label %36
31: ; preds = %21
%32 = load double, double* %1
%33 = fadd double %32, %7
store double %33, double* %1
%34 = load double, double* %0
%35 = tail call double @llvm.fmuladd.f64(double %27, double %7, double %34)
store double %35, double* %0
br label %36
36: ; preds = %21, %8, %15, %31
ret void
}