void radialDiffRecurrence(double x, double *y1, double *yd, int n) {
  int k;
  for (k = 0; k < n; k++) {
    yd[k] = y1[k] / x * n - y1[k + 1];
  }
}


x86 assembly (clang, -O0):
<radialDiffRecurrence>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    118a <radialDiffRecurrence+0x8a>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x8(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    111b <radialDiffRecurrence+0x1b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<radialDiffRecurrence>:
test   %edx,%edx
jle    126a <radialDiffRecurrence+0x16a>
cvtsi2sd %edx,%xmm1
mov    %edx,%r9d
cmp    $0x1,%edx
je     1132 <radialDiffRecurrence+0x32>
lea    (%rdi,%r9,8),%rax
add    $0x8,%rax
cmp    %rsi,%rax
jbe    11b3 <radialDiffRecurrence+0xb3>
lea    (%rsi,%r9,8),%rax
cmp    %rdi,%rax
jbe    11b3 <radialDiffRecurrence+0xb3>
xor    %ecx,%ecx
mov    %rcx,%rax
not    %rax
test   $0x1,%r9b
je     1161 <radialDiffRecurrence+0x61>
movsd  (%rdi,%rcx,8),%xmm2
divsd  %xmm0,%xmm2
mov    %rcx,%rdx
or     $0x1,%rdx
mulsd  %xmm1,%xmm2
subsd  (%rdi,%rdx,8),%xmm2
movsd  %xmm2,(%rsi,%rcx,8)
mov    %rdx,%rcx
add    %r9,%rax
je     126a <radialDiffRecurrence+0x16a>
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rcx,8),%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm2
subsd  0x8(%rdi,%rcx,8),%xmm2
movsd  %xmm2,(%rsi,%rcx,8)
movsd  0x8(%rdi,%rcx,8),%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm2
subsd  0x10(%rdi,%rcx,8),%xmm2
movsd  %xmm2,0x8(%rsi,%rcx,8)
lea    0x2(%rcx),%rax
mov    %rax,%rcx
cmp    %rax,%r9
jne    1170 <radialDiffRecurrence+0x70>
jmp    126a <radialDiffRecurrence+0x16a>
mov    %r9d,%ecx
and    $0xfffffffe,%ecx
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
lea    -0x2(%rcx),%rax
mov    %rax,%r8
shr    %r8
add    $0x1,%r8
test   %rax,%rax
je     126b <radialDiffRecurrence+0x16b>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
xor    %edx,%edx
nopl   0x0(%rax)
movupd (%rdi,%rdx,8),%xmm4
movupd 0x8(%rdi,%rdx,8),%xmm5
divpd  %xmm3,%xmm4
mulpd  %xmm2,%xmm4
subpd  %xmm5,%xmm4
movupd %xmm4,(%rsi,%rdx,8)
movupd 0x10(%rdi,%rdx,8),%xmm4
movupd 0x18(%rdi,%rdx,8),%xmm5
divpd  %xmm3,%xmm4
mulpd  %xmm2,%xmm4
subpd  %xmm5,%xmm4
movupd %xmm4,0x10(%rsi,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rax
jne    11f0 <radialDiffRecurrence+0xf0>
test   $0x1,%r8b
je     1261 <radialDiffRecurrence+0x161>
lea    0x0(,%rdx,8),%rax
movupd (%rdi,%rdx,8),%xmm4
divpd  %xmm3,%xmm4
or     $0x8,%rax
movupd (%rdi,%rax,1),%xmm3
mulpd  %xmm2,%xmm4
subpd  %xmm3,%xmm4
movupd %xmm4,(%rsi,%rdx,8)
cmp    %r9,%rcx
jne    1134 <radialDiffRecurrence+0x34>
ret
xor    %edx,%edx
test   $0x1,%r8b
jne    123a <radialDiffRecurrence+0x13a>
jmp    1261 <radialDiffRecurrence+0x161>

LLVM IR (-O0):
define void @radialDiffRecurrence(double %0, double* %1, double* %2, i32 %3) {
%5 = alloca double
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
store double %0, double* %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
store i32 0, i32* %9
br label %10
10: ; preds = %36, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %8
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %39
14: ; preds = %10
%15 = load double*, double** %6
%16 = load i32, i32* %9
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds double, double* %15, i64 %17
%19 = load double, double* %18
%20 = load double, double* %5
%21 = fdiv double %19, %20
%22 = load i32, i32* %8
%23 = sitofp i32 %22 to double
%24 = load double*, double** %6
%25 = load i32, i32* %9
%26 = add nsw i32 %25, 1
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %24, i64 %27
%29 = load double, double* %28
%30 = fneg double %29
%31 = call double @llvm.fmuladd.f64(double %21, double %23, double %30)
%32 = load double*, double** %7
%33 = load i32, i32* %9
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
store double %31, double* %35
br label %36
36: ; preds = %14
%37 = load i32, i32* %9
%38 = add nsw i32 %37, 1
store i32 %38, i32* %9
br label %10
39: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @radialDiffRecurrence(double %0, double* nocapture %1, double* nocapture writeonly %2, i32 %3) {
%5 = icmp sgt i32 %3, 0
br i1 %5, label %6, label %119
6: ; preds = %4
%7 = sitofp i32 %3 to double
%8 = zext i32 %3 to i64
%9 = icmp eq i32 %3, 1
br i1 %9, label %79, label %10
10: ; preds = %6
%11 = getelementptr double, double* %2, i64 %8
%12 = add nuw nsw i64 %8, 1
%13 = getelementptr double, double* %1, i64 %12
%14 = icmp ugt double* %13, %2
%15 = icmp ugt double* %11, %1
%16 = and i1 %14, %15
br i1 %16, label %79, label %17
17: ; preds = %10
%18 = and i64 %8, 4294967294
%19 = insertelement <2 x double> poison, double %0, i64 0
%20 = shufflevector <2 x double> %19, <2 x double> poison, <2 x i32> zeroinitializer
%21 = insertelement <2 x double> poison, double %7, i64 0
%22 = shufflevector <2 x double> %21, <2 x double> poison, <2 x i32> zeroinitializer
%23 = add nsw i64 %18, -2
%24 = lshr exact i64 %23, 1
%25 = add nuw i64 %24, 1
%26 = and i64 %25, 1
%27 = icmp eq i64 %23, 0
br i1 %27, label %61, label %28
28: ; preds = %17
%29 = and i64 %25, -2
br label %30
30: ; preds = %30, %28
%31 = phi i64 [ 0, %28 ], [ %58, %30 ]
%32 = phi i64 [ 0, %28 ], [ %59, %30 ]
%33 = getelementptr inbounds double, double* %1, i64 %31
%34 = bitcast double* %33 to <2 x double>*
%35 = load <2 x double>, <2 x double>* %34
%36 = fdiv <2 x double> %35, %20
%37 = or i64 %31, 1
%38 = getelementptr inbounds double, double* %1, i64 %37
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = fneg <2 x double> %40
%42 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %36, <2 x double> %22, <2 x double> %41)
%43 = getelementptr inbounds double, double* %2, i64 %31
%44 = bitcast double* %43 to <2 x double>*
store <2 x double> %42, <2 x double>* %44
%45 = or i64 %31, 2
%46 = getelementptr inbounds double, double* %1, i64 %45
%47 = bitcast double* %46 to <2 x double>*
%48 = load <2 x double>, <2 x double>* %47
%49 = fdiv <2 x double> %48, %20
%50 = or i64 %31, 3
%51 = getelementptr inbounds double, double* %1, i64 %50
%52 = bitcast double* %51 to <2 x double>*
%53 = load <2 x double>, <2 x double>* %52
%54 = fneg <2 x double> %53
%55 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %49, <2 x double> %22, <2 x double> %54)
%56 = getelementptr inbounds double, double* %2, i64 %45
%57 = bitcast double* %56 to <2 x double>*
store <2 x double> %55, <2 x double>* %57
%58 = add nuw i64 %31, 4
%59 = add nuw i64 %32, 2
%60 = icmp eq i64 %59, %29
br i1 %60, label %61, label %30
61: ; preds = %30, %17
%62 = phi i64 [ 0, %17 ], [ %58, %30 ]
%63 = icmp eq i64 %26, 0
br i1 %63, label %77, label %64
64: ; preds = %61
%65 = getelementptr inbounds double, double* %1, i64 %62
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = fdiv <2 x double> %67, %20
%69 = or i64 %62, 1
%70 = getelementptr inbounds double, double* %1, i64 %69
%71 = bitcast double* %70 to <2 x double>*
%72 = load <2 x double>, <2 x double>* %71
%73 = fneg <2 x double> %72
%74 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %68, <2 x double> %22, <2 x double> %73)
%75 = getelementptr inbounds double, double* %2, i64 %62
%76 = bitcast double* %75 to <2 x double>*
store <2 x double> %74, <2 x double>* %76
br label %77
77: ; preds = %61, %64
%78 = icmp eq i64 %18, %8
br i1 %78, label %119, label %79
79: ; preds = %10, %6, %77
%80 = phi i64 [ 0, %10 ], [ 0, %6 ], [ %18, %77 ]
%81 = xor i64 %80, -1
%82 = and i64 %8, 1
%83 = icmp eq i64 %82, 0
br i1 %83, label %94, label %84
84: ; preds = %79
%85 = getelementptr inbounds double, double* %1, i64 %80
%86 = load double, double* %85
%87 = fdiv double %86, %0
%88 = or i64 %80, 1
%89 = getelementptr inbounds double, double* %1, i64 %88
%90 = load double, double* %89
%91 = fneg double %90
%92 = tail call double @llvm.fmuladd.f64(double %87, double %7, double %91)
%93 = getelementptr inbounds double, double* %2, i64 %80
store double %92, double* %93
br label %94
94: ; preds = %84, %79
%95 = phi i64 [ %80, %79 ], [ %88, %84 ]
%96 = sub nsw i64 0, %8
%97 = icmp eq i64 %81, %96
br i1 %97, label %119, label %98
98: ; preds = %94, %98
%99 = phi i64 [ %112, %98 ], [ %95, %94 ]
%100 = getelementptr inbounds double, double* %1, i64 %99
%101 = load double, double* %100
%102 = fdiv double %101, %0
%103 = add nuw nsw i64 %99, 1
%104 = getelementptr inbounds double, double* %1, i64 %103
%105 = load double, double* %104
%106 = fneg double %105
%107 = tail call double @llvm.fmuladd.f64(double %102, double %7, double %106)
%108 = getelementptr inbounds double, double* %2, i64 %99
store double %107, double* %108
%109 = getelementptr inbounds double, double* %1, i64 %103
%110 = load double, double* %109
%111 = fdiv double %110, %0
%112 = add nuw nsw i64 %99, 2
%113 = getelementptr inbounds double, double* %1, i64 %112
%114 = load double, double* %113
%115 = fneg double %114
%116 = tail call double @llvm.fmuladd.f64(double %111, double %7, double %115)
%117 = getelementptr inbounds double, double* %2, i64 %103
store double %116, double* %117
%118 = icmp eq i64 %112, %8
br i1 %118, label %119, label %98
119: ; preds = %94, %98, %77, %4
ret void
}