void corr_hum(long gridsize, double *q, double q_min) {
  long i;
  for (i = 0; i < gridsize; ++i) {
    if (q[i] < q_min)
      q[i] = q_min;
  }
}


x86 assembly (clang, -O0):
<corr_hum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    116b <corr_hum+0x6b>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1155 <corr_hum+0x55>
movsd  -0x18(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    115a <corr_hum+0x5a>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1119 <corr_hum+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<corr_hum>:
test   %rdi,%rdi
jle    111f <corr_hum+0x1f>
cmp    $0x1,%rdi
jne    1120 <corr_hum+0x20>
xor    %eax,%eax
test   $0x1,%dil
je     111f <corr_hum+0x1f>
ucomisd (%rsi,%rax,8),%xmm0
jbe    111f <corr_hum+0x1f>
movsd  %xmm0,(%rsi,%rax,8)
ret
mov    %rdi,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %eax,%eax
jmp    1139 <corr_hum+0x39>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rax
cmp    %rax,%rcx
je     110d <corr_hum+0xd>
ucomisd (%rsi,%rax,8),%xmm0
jbe    1145 <corr_hum+0x45>
movsd  %xmm0,(%rsi,%rax,8)
ucomisd 0x8(%rsi,%rax,8),%xmm0
jbe    1130 <corr_hum+0x30>
movsd  %xmm0,0x8(%rsi,%rax,8)
jmp    1130 <corr_hum+0x30>

LLVM IR (-O0):
define void @corr_hum(i64 %0, double* %1, double %2) {
%4 = alloca i64
%5 = alloca double*
%6 = alloca double
%7 = alloca i64
store i64 %0, i64* %4
store double* %1, double** %5
store double %2, double* %6
store i64 0, i64* %7
br label %8
8: ; preds = %25, %3
%9 = load i64, i64* %7
%10 = load i64, i64* %4
%11 = icmp slt i64 %9, %10
br i1 %11, label %12, label %28
12: ; preds = %8
%13 = load double*, double** %5
%14 = load i64, i64* %7
%15 = getelementptr inbounds double, double* %13, i64 %14
%16 = load double, double* %15
%17 = load double, double* %6
%18 = fcmp olt double %16, %17
br i1 %18, label %19, label %24
19: ; preds = %12
%20 = load double, double* %6
%21 = load double*, double** %5
%22 = load i64, i64* %7
%23 = getelementptr inbounds double, double* %21, i64 %22
store double %20, double* %23
br label %24
24: ; preds = %19, %12
br label %25
25: ; preds = %24
%26 = load i64, i64* %7
%27 = add nsw i64 %26, 1
store i64 %27, i64* %7
br label %8
28: ; preds = %8
ret void
}

LLVM IR (-O2):
define void @corr_hum(i64 %0, double* nocapture %1, double %2) {
%4 = icmp sgt i64 %0, 0
br i1 %4, label %5, label %35
5: ; preds = %3
%6 = and i64 %0, 1
%7 = icmp eq i64 %0, 1
br i1 %7, label %27, label %8
8: ; preds = %5
%9 = and i64 %0, -2
br label %10
10: ; preds = %23, %8
%11 = phi i64 [ 0, %8 ], [ %24, %23 ]
%12 = phi i64 [ 0, %8 ], [ %25, %23 ]
%13 = getelementptr inbounds double, double* %1, i64 %11
%14 = load double, double* %13
%15 = fcmp olt double %14, %2
br i1 %15, label %16, label %17
16: ; preds = %10
store double %2, double* %13
br label %17
17: ; preds = %10, %16
%18 = or i64 %11, 1
%19 = getelementptr inbounds double, double* %1, i64 %18
%20 = load double, double* %19
%21 = fcmp olt double %20, %2
br i1 %21, label %22, label %23
22: ; preds = %17
store double %2, double* %19
br label %23
23: ; preds = %22, %17
%24 = add nuw nsw i64 %11, 2
%25 = add i64 %12, 2
%26 = icmp eq i64 %25, %9
br i1 %26, label %27, label %10
27: ; preds = %23, %5
%28 = phi i64 [ 0, %5 ], [ %24, %23 ]
%29 = icmp eq i64 %6, 0
br i1 %29, label %35, label %30
30: ; preds = %27
%31 = getelementptr inbounds double, double* %1, i64 %28
%32 = load double, double* %31
%33 = fcmp olt double %32, %2
br i1 %33, label %34, label %35
34: ; preds = %30
store double %2, double* %31
br label %35
35: ; preds = %27, %34, %30, %3
ret void
}