void reg_wt(double *sum, double *wt, float *data, int i, int j, int nx, int ny,
            double wt0) {
  float t;
  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):
<reg_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   $0x0,-0x1c(%rbp)
jl     114f <reg_wt+0x4f>
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    114f <reg_wt+0x4f>
cmpl   $0x0,-0x20(%rbp)
jl     114f <reg_wt+0x4f>
mov    -0x20(%rbp),%eax
cmp    0x10(%rbp),%eax
jl     1154 <reg_wt+0x54>
jmp    11e7 <reg_wt+0xe7>
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  0xe7e(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     11ac <reg_wt+0xac>
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe6b(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     11ac <reg_wt+0xac>
jmp    11e7 <reg_wt+0xe7>
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):
<reg_wt>:
test   %ecx,%ecx
js     1140 <reg_wt+0x40>
test   %r8d,%r8d
js     1140 <reg_wt+0x40>
cmp    %r9d,%ecx
jge    1140 <reg_wt+0x40>
cmp    0x8(%rsp),%r8d
jge    1140 <reg_wt+0x40>
imul   %r8d,%r9d
add    %ecx,%r9d
movslq %r9d,%rax
movss  (%rdx,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
ucomisd 0xed0(%rip),%xmm1
jb     1141 <reg_wt+0x41>
movsd  0xece(%rip),%xmm2
ucomisd %xmm1,%xmm2
jb     1141 <reg_wt+0x41>
ret
mulsd  %xmm0,%xmm1
addsd  (%rsi),%xmm0
movsd  %xmm0,(%rsi)
addsd  (%rdi),%xmm1
movsd  %xmm1,(%rdi)
ret

LLVM IR (-O0):
define void @reg_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 slt i32 %18, 0
br i1 %19, label %31, label %20
20: ; preds = %8
%21 = load i32, i32* %12
%22 = load i32, i32* %14
%23 = icmp sge i32 %21, %22
br i1 %23, label %31, label %24
24: ; preds = %20
%25 = load i32, i32* %13
%26 = icmp slt i32 %25, 0
br i1 %26, label %31, label %27
27: ; preds = %24
%28 = load i32, i32* %13
%29 = load i32, i32* %15
%30 = icmp sge i32 %28, %29
br i1 %30, label %31, label %32
31: ; preds = %27, %24, %20, %8
br label %63
32: ; preds = %27
%33 = load float*, float** %11
%34 = load i32, i32* %12
%35 = load i32, i32* %13
%36 = load i32, i32* %14
%37 = mul nsw i32 %35, %36
%38 = add nsw i32 %34, %37
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds float, float* %33, i64 %39
%41 = load float, float* %40
store float %41, float* %17
%42 = load float, float* %17
%43 = fpext float %42 to double
%44 = fcmp oge double %43, 0x444B1A2170AD88AA
br i1 %44, label %45, label %50
45: ; preds = %32
%46 = load float, float* %17
%47 = fpext float %46 to double
%48 = fcmp ole double %47, 9.999100e+20
br i1 %48, label %49, label %50
49: ; preds = %45
br label %63
50: ; preds = %45, %32
%51 = load double*, double** %10
%52 = load double, double* %51
%53 = load double, double* %16
%54 = fadd double %52, %53
%55 = load double*, double** %10
store double %54, double* %55
%56 = load double*, double** %9
%57 = load double, double* %56
%58 = load float, float* %17
%59 = fpext float %58 to double
%60 = load double, double* %16
%61 = call double @llvm.fmuladd.f64(double %59, double %60, double %57)
%62 = load double*, double** %9
store double %61, double* %62
br label %63
63: ; preds = %50, %49, %31
ret void
}

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