double WattersonsTheta(int nsam, int s) {
  int i;
  double a1, tw;
  if (s == 0) {
    tw = 0.0;
  } else {
    a1 = 0.0;
    i = nsam - 1;
    while (i >= 1) {
      a1 += 1.0 / i;
      i--;
    }
    tw = s / a1;
  }
  return tw;
}


x86 assembly (clang, -O0):
<WattersonsTheta>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
cmpl   $0x0,-0x8(%rbp)
jne    1121 <WattersonsTheta+0x21>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1174 <WattersonsTheta+0x74>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0xc(%rbp)
cmpl   $0x1,-0xc(%rbp)
jl     1165 <WattersonsTheta+0x65>
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  0xeb7(%rip),%xmm0
divsd  %xmm1,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0xc(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0xc(%rbp)
jmp    1132 <WattersonsTheta+0x32>
cvtsi2sdl -0x8(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<WattersonsTheta>:
test   %esi,%esi
je     114a <WattersonsTheta+0x4a>
xorpd  %xmm1,%xmm1
cmp    $0x2,%edi
jl     113e <WattersonsTheta+0x3e>
add    $0x1,%edi
xorpd  %xmm1,%xmm1
movsd  0xee4(%rip),%xmm0
nopl   0x0(%rax)
lea    -0x2(%rdi),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movapd %xmm0,%xmm3
divsd  %xmm2,%xmm3
addsd  %xmm3,%xmm1
add    $0xffffffff,%edi
cmp    $0x2,%edi
ja     1120 <WattersonsTheta+0x20>
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
divsd  %xmm1,%xmm0
ret
xorps  %xmm0,%xmm0
ret

LLVM IR (-O0):
define double @WattersonsTheta(i32 %0, i32 %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
%7 = alloca double
store i32 %0, i32* %3
store i32 %1, i32* %4
%8 = load i32, i32* %4
%9 = icmp eq i32 %8, 0
br i1 %9, label %10, label %11
10: ; preds = %2
store double 0.000000e+00, double* %7
br label %30
11: ; preds = %2
store double 0.000000e+00, double* %6
%12 = load i32, i32* %3
%13 = sub nsw i32 %12, 1
store i32 %13, i32* %5
br label %14
14: ; preds = %17, %11
%15 = load i32, i32* %5
%16 = icmp sge i32 %15, 1
br i1 %16, label %17, label %25
17: ; preds = %14
%18 = load i32, i32* %5
%19 = sitofp i32 %18 to double
%20 = fdiv double 1.000000e+00, %19
%21 = load double, double* %6
%22 = fadd double %21, %20
store double %22, double* %6
%23 = load i32, i32* %5
%24 = add nsw i32 %23, -1
store i32 %24, i32* %5
br label %14
25: ; preds = %14
%26 = load i32, i32* %4
%27 = sitofp i32 %26 to double
%28 = load double, double* %6
%29 = fdiv double %27, %28
store double %29, double* %7
br label %30
30: ; preds = %25, %10
%31 = load double, double* %7
ret double %31
}

LLVM IR (-O2):
define double @WattersonsTheta(i32 %0, i32 %1) {
%3 = icmp eq i32 %1, 0
br i1 %3, label %18, label %4
4: ; preds = %2
%5 = icmp sgt i32 %0, 1
br i1 %5, label %6, label %14
6: ; preds = %4, %6
%7 = phi i32 [ %9, %6 ], [ %0, %4 ]
%8 = phi double [ %12, %6 ], [ 0.000000e+00, %4 ]
%9 = add nsw i32 %7, -1
%10 = sitofp i32 %9 to double
%11 = fdiv double 1.000000e+00, %10
%12 = fadd double %8, %11
%13 = icmp ugt i32 %7, 2
br i1 %13, label %6, label %14
14: ; preds = %6, %4
%15 = phi double [ 0.000000e+00, %4 ], [ %12, %6 ]
%16 = sitofp i32 %1 to double
%17 = fdiv double %16, %15
br label %18
18: ; preds = %2, %14
%19 = phi double [ %17, %14 ], [ 0.000000e+00, %2 ]
ret double %19
}