void update(int n, float a, int str, float *x) {
  for (int i = 2; i < n - 1; i += 2)
    x[i * str] += a * (x[(i - 1) * str] + x[(i + 1) * str]);
  if (0 == n % 2)
    x[0] += a * (x[(n - 1) * str] + x[1 * str]);
  else {
    x[0] += 2. * a * x[1 * str];
    x[(n - 1) * str] += 2. * a * x[(n - 2) * str];
  }
}


x86 assembly (clang, -O0):
<update>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movss  %xmm0,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x2,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    118a <update+0x8a>
movss  -0x8(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    111a <update+0x1a>
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
xor    %eax,%eax
cmp    %edx,%eax
jne    11e5 <update+0xe5>
movss  -0x8(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
shl    $0x0,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
jmp    1286 <update+0x186>
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe0a(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
shl    $0x0,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax)
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdc3(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
pop    %rbp
ret

x86 assembly (clang, -O2):
<update>:
lea    -0x1(%rdi),%r8d
movslq %esi,%rcx
cmp    $0x4,%edi
jl     1155 <update+0x55>
mov    %r8d,%r9d
lea    0x0(,%rcx,4),%rax
lea    (%rax,%rax,2),%r10
lea    0x0(,%rcx,8),%r11
mov    $0x2,%eax
mov    %rdx,%rsi
nopl   0x0(%rax,%rax,1)
movss  (%rsi,%rcx,4),%xmm1
addss  (%rsi,%r10,1),%xmm1
mulss  %xmm0,%xmm1
addss  (%rsi,%rcx,8),%xmm1
movss  %xmm1,(%rsi,%rcx,8)
add    $0x2,%rax
add    %r11,%rsi
cmp    %r9,%rax
jb     1130 <update+0x30>
test   $0x1,%dil
jne    1179 <update+0x79>
imul   %ecx,%r8d
movslq %r8d,%rax
movss  (%rdx,%rax,4),%xmm1
addss  (%rdx,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
addss  (%rdx),%xmm0
movss  %xmm0,(%rdx)
ret
cvtss2sd %xmm0,%xmm0
addsd  %xmm0,%xmm0
movss  (%rdx,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  (%rdx),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rdx)
add    $0xfffffffe,%edi
imul   %ecx,%edi
movslq %edi,%rax
movss  (%rdx,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
imul   %ecx,%r8d
movslq %r8d,%rax
movss  (%rdx,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,(%rdx,%rax,4)
ret

LLVM IR (-O0):
define void @update(i32 %0, float %1, i32 %2, float* %3) {
%5 = alloca i32
%6 = alloca float
%7 = alloca i32
%8 = alloca float*
%9 = alloca i32
store i32 %0, i32* %5
store float %1, float* %6
store i32 %2, i32* %7
store float* %3, float** %8
store i32 2, i32* %9
br label %10
10: ; preds = %42, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %5
%13 = sub nsw i32 %12, 1
%14 = icmp slt i32 %11, %13
br i1 %14, label %15, label %45
15: ; preds = %10
%16 = load float, float* %6
%17 = load float*, float** %8
%18 = load i32, i32* %9
%19 = sub nsw i32 %18, 1
%20 = load i32, i32* %7
%21 = mul nsw i32 %19, %20
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds float, float* %17, i64 %22
%24 = load float, float* %23
%25 = load float*, float** %8
%26 = load i32, i32* %9
%27 = add nsw i32 %26, 1
%28 = load i32, i32* %7
%29 = mul nsw i32 %27, %28
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %25, i64 %30
%32 = load float, float* %31
%33 = fadd float %24, %32
%34 = load float*, float** %8
%35 = load i32, i32* %9
%36 = load i32, i32* %7
%37 = mul nsw i32 %35, %36
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds float, float* %34, i64 %38
%40 = load float, float* %39
%41 = call float @llvm.fmuladd.f32(float %16, float %33, float %40)
store float %41, float* %39
br label %42
42: ; preds = %15
%43 = load i32, i32* %9
%44 = add nsw i32 %43, 2
store i32 %44, i32* %9
br label %10
45: ; preds = %10
%46 = load i32, i32* %5
%47 = srem i32 %46, 2
%48 = icmp eq i32 0, %47
br i1 %48, label %49, label %70
49: ; preds = %45
%50 = load float, float* %6
%51 = load float*, float** %8
%52 = load i32, i32* %5
%53 = sub nsw i32 %52, 1
%54 = load i32, i32* %7
%55 = mul nsw i32 %53, %54
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds float, float* %51, i64 %56
%58 = load float, float* %57
%59 = load float*, float** %8
%60 = load i32, i32* %7
%61 = mul nsw i32 1, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %59, i64 %62
%64 = load float, float* %63
%65 = fadd float %58, %64
%66 = load float*, float** %8
%67 = getelementptr inbounds float, float* %66, i64 0
%68 = load float, float* %67
%69 = call float @llvm.fmuladd.f32(float %50, float %65, float %68)
store float %69, float* %67
br label %110
70: ; preds = %45
%71 = load float, float* %6
%72 = fpext float %71 to double
%73 = fmul double 2.000000e+00, %72
%74 = load float*, float** %8
%75 = load i32, i32* %7
%76 = mul nsw i32 1, %75
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %74, i64 %77
%79 = load float, float* %78
%80 = fpext float %79 to double
%81 = load float*, float** %8
%82 = getelementptr inbounds float, float* %81, i64 0
%83 = load float, float* %82
%84 = fpext float %83 to double
%85 = call double @llvm.fmuladd.f64(double %73, double %80, double %84)
%86 = fptrunc double %85 to float
store float %86, float* %82
%87 = load float, float* %6
%88 = fpext float %87 to double
%89 = fmul double 2.000000e+00, %88
%90 = load float*, float** %8
%91 = load i32, i32* %5
%92 = sub nsw i32 %91, 2
%93 = load i32, i32* %7
%94 = mul nsw i32 %92, %93
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds float, float* %90, i64 %95
%97 = load float, float* %96
%98 = fpext float %97 to double
%99 = load float*, float** %8
%100 = load i32, i32* %5
%101 = sub nsw i32 %100, 1
%102 = load i32, i32* %7
%103 = mul nsw i32 %101, %102
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds float, float* %99, i64 %104
%106 = load float, float* %105
%107 = fpext float %106 to double
%108 = call double @llvm.fmuladd.f64(double %89, double %98, double %107)
%109 = fptrunc double %108 to float
store float %109, float* %105
br label %110
110: ; preds = %70, %49
ret void
}

LLVM IR (-O2):
define void @update(i32 %0, float %1, i32 %2, float* nocapture %3) {
%5 = add nsw i32 %0, -1
%6 = icmp sgt i32 %0, 3
br i1 %6, label %7, label %10
7: ; preds = %4
%8 = sext i32 %2 to i64
%9 = zext i32 %5 to i64
br label %13
10: ; preds = %13, %4
%11 = and i32 %0, 1
%12 = icmp eq i32 %11, 0
br i1 %12, label %30, label %41
13: ; preds = %7, %13
%14 = phi i64 [ 2, %7 ], [ %28, %13 ]
%15 = add nsw i64 %14, -1
%16 = mul nsw i64 %15, %8
%17 = getelementptr inbounds float, float* %3, i64 %16
%18 = load float, float* %17
%19 = or i64 %14, 1
%20 = mul nsw i64 %19, %8
%21 = getelementptr inbounds float, float* %3, i64 %20
%22 = load float, float* %21
%23 = fadd float %18, %22
%24 = mul nsw i64 %14, %8
%25 = getelementptr inbounds float, float* %3, i64 %24
%26 = load float, float* %25
%27 = tail call float @llvm.fmuladd.f32(float %1, float %23, float %26)
store float %27, float* %25
%28 = add nuw nsw i64 %14, 2
%29 = icmp ult i64 %28, %9
br i1 %29, label %13, label %10
30: ; preds = %10
%31 = mul nsw i32 %5, %2
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %3, i64 %32
%34 = load float, float* %33
%35 = sext i32 %2 to i64
%36 = getelementptr inbounds float, float* %3, i64 %35
%37 = load float, float* %36
%38 = fadd float %34, %37
%39 = load float, float* %3
%40 = tail call float @llvm.fmuladd.f32(float %1, float %38, float %39)
store float %40, float* %3
br label %65
41: ; preds = %10
%42 = fpext float %1 to double
%43 = fmul double %42, 2.000000e+00
%44 = sext i32 %2 to i64
%45 = getelementptr inbounds float, float* %3, i64 %44
%46 = load float, float* %45
%47 = fpext float %46 to double
%48 = load float, float* %3
%49 = fpext float %48 to double
%50 = tail call double @llvm.fmuladd.f64(double %43, double %47, double %49)
%51 = fptrunc double %50 to float
store float %51, float* %3
%52 = add nsw i32 %0, -2
%53 = mul nsw i32 %52, %2
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %3, i64 %54
%56 = load float, float* %55
%57 = fpext float %56 to double
%58 = mul nsw i32 %5, %2
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %3, i64 %59
%61 = load float, float* %60
%62 = fpext float %61 to double
%63 = tail call double @llvm.fmuladd.f64(double %43, double %57, double %62)
%64 = fptrunc double %63 to float
store float %64, float* %60
br label %65
65: ; preds = %41, %30
ret void
}