void integrate(s, nt, dt, iv) float *s, *dt, *iv;
int nt;
{
  float dt3, dt8, s0, s1, s2, s3, s4;
  float s_start, s_end;
  int it;
  float four = 4.0;
  float six = 6.0;
  dt3 = (*dt) / 3.0;
  dt8 = (*dt) / 8.0;
  s0 = *iv;
  for (it = 0; it < nt; it++) {
    s0 = s[it] * (*dt) + s0;
    s[it] = s0;
  }
}


x86 assembly (clang, -O0):
<integrate>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  0xef9(%rip),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  0xee8(%rip),%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xec7(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe9e(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d6 <integrate+0xd6>
mov    -0x8(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1183 <integrate+0x83>
pop    %rbp
ret

x86 assembly (clang, -O2):
<integrate>:
test   %esi,%esi
jle    1168 <integrate+0x68>
movss  (%rcx),%xmm0
mov    %esi,%eax
cmp    $0x1,%esi
jne    1113 <integrate+0x13>
xor    %ecx,%ecx
jmp    1152 <integrate+0x52>
mov    %eax,%esi
and    $0xfffffffe,%esi
xor    %ecx,%ecx
movaps %xmm0,%xmm1
nopl   (%rax)
movss  (%rdi,%rcx,4),%xmm2
movss  0x4(%rdi,%rcx,4),%xmm0
mulss  (%rdx),%xmm2
addss  %xmm1,%xmm2
movss  %xmm2,(%rdi,%rcx,4)
mulss  (%rdx),%xmm0
addss  %xmm2,%xmm0
movss  %xmm0,0x4(%rdi,%rcx,4)
add    $0x2,%rcx
movaps %xmm0,%xmm1
cmp    %rcx,%rsi
jne    1120 <integrate+0x20>
test   $0x1,%al
je     1168 <integrate+0x68>
movss  (%rdi,%rcx,4),%xmm1
mulss  (%rdx),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rdi,%rcx,4)
ret

LLVM IR (-O0):
define void @integrate(float* %0, i32 %1, float* %2, float* %3) {
%5 = alloca float*
%6 = alloca i32
%7 = alloca float*
%8 = alloca float*
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca i32
%19 = alloca float
%20 = alloca float
store float* %0, float** %5
store i32 %1, i32* %6
store float* %2, float** %7
store float* %3, float** %8
store float 4.000000e+00, float* %19
store float 6.000000e+00, float* %20
%21 = load float*, float** %7
%22 = load float, float* %21
%23 = fpext float %22 to double
%24 = fdiv double %23, 3.000000e+00
%25 = fptrunc double %24 to float
store float %25, float* %9
%26 = load float*, float** %7
%27 = load float, float* %26
%28 = fpext float %27 to double
%29 = fdiv double %28, 8.000000e+00
%30 = fptrunc double %29 to float
store float %30, float* %10
%31 = load float*, float** %8
%32 = load float, float* %31
store float %32, float* %11
store i32 0, i32* %18
br label %33
33: ; preds = %52, %4
%34 = load i32, i32* %18
%35 = load i32, i32* %6
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %55
37: ; preds = %33
%38 = load float*, float** %5
%39 = load i32, i32* %18
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
%42 = load float, float* %41
%43 = load float*, float** %7
%44 = load float, float* %43
%45 = load float, float* %11
%46 = call float @llvm.fmuladd.f32(float %42, float %44, float %45)
store float %46, float* %11
%47 = load float, float* %11
%48 = load float*, float** %5
%49 = load i32, i32* %18
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %48, i64 %50
store float %47, float* %51
br label %52
52: ; preds = %37
%53 = load i32, i32* %18
%54 = add nsw i32 %53, 1
store i32 %54, i32* %18
br label %33
55: ; preds = %33
ret void
}

LLVM IR (-O2):
define void @integrate(float* nocapture %0, i32 %1, float* nocapture %2, float* nocapture %3) {
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %38
6: ; preds = %4
%7 = load float, float* %3
%8 = zext i32 %1 to i64
%9 = and i64 %8, 1
%10 = icmp eq i32 %1, 1
br i1 %10, label %29, label %11
11: ; preds = %6
%12 = and i64 %8, 4294967294
br label %13
13: ; preds = %13, %11
%14 = phi i64 [ 0, %11 ], [ %26, %13 ]
%15 = phi float [ %7, %11 ], [ %25, %13 ]
%16 = phi i64 [ 0, %11 ], [ %27, %13 ]
%17 = getelementptr inbounds float, float* %0, i64 %14
%18 = load float, float* %17
%19 = load float, float* %2
%20 = tail call float @llvm.fmuladd.f32(float %18, float %19, float %15)
store float %20, float* %17
%21 = or i64 %14, 1
%22 = getelementptr inbounds float, float* %0, i64 %21
%23 = load float, float* %22
%24 = load float, float* %2
%25 = tail call float @llvm.fmuladd.f32(float %23, float %24, float %20)
store float %25, float* %22
%26 = add nuw nsw i64 %14, 2
%27 = add i64 %16, 2
%28 = icmp eq i64 %27, %12
br i1 %28, label %29, label %13
29: ; preds = %13, %6
%30 = phi i64 [ 0, %6 ], [ %26, %13 ]
%31 = phi float [ %7, %6 ], [ %25, %13 ]
%32 = icmp eq i64 %9, 0
br i1 %32, label %38, label %33
33: ; preds = %29
%34 = getelementptr inbounds float, float* %0, i64 %30
%35 = load float, float* %34
%36 = load float, float* %2
%37 = tail call float @llvm.fmuladd.f32(float %35, float %36, float %31)
store float %37, float* %34
br label %38
38: ; preds = %33, %29, %4
ret void
}