void compute_dfactor(double *dx, int N, double *dfactor) {
  int ii;
  for (ii = 1; ii < N - 1; ii++)
    dfactor[ii] = 2. / (dx[ii] + dx[ii - 1]);
  dfactor[0] = 2. / dx[0];
  dfactor[N - 1] = 2. / dx[N - 2];
}


x86 assembly (clang, -O0):
<compute_dfactor>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    116d <compute_dfactor+0x6d>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
movsd  0xeb2(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1116 <compute_dfactor+0x16>
mov    -0x8(%rbp),%rax
movsd  0xe87(%rip),%xmm0
divsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  0xe66(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_dfactor>:
push   %rbx
lea    -0x1(%rsi),%r9d
cmp    $0x3,%esi
jl     126a <compute_dfactor+0x16a>
mov    %r9d,%ecx
lea    -0x1(%rcx),%r8
mov    $0x1,%eax
cmp    $0x2,%r8
jb     11f2 <compute_dfactor+0xf2>
lea    0x8(%rdx),%r10
lea    (%rdi,%rcx,8),%rbx
cmp    %rbx,%r10
jae    113e <compute_dfactor+0x3e>
lea    (%rdx,%rcx,8),%rbx
cmp    %rdi,%rbx
ja     11f2 <compute_dfactor+0xf2>
mov    %r8,%r10
and    $0xfffffffffffffffe,%r10
lea    -0x2(%r10),%rax
mov    %rax,%r11
shr    %r11
add    $0x1,%r11
test   %rax,%rax
je     1291 <compute_dfactor+0x191>
mov    %r11,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %eax,%eax
movapd 0xe93(%rip),%xmm0
nopl   (%rax)
movupd (%rdi,%rax,8),%xmm1
movupd 0x8(%rdi,%rax,8),%xmm2
addpd  %xmm1,%xmm2
movapd %xmm0,%xmm1
divpd  %xmm2,%xmm1
movupd %xmm1,0x8(%rdx,%rax,8)
movupd 0x10(%rdi,%rax,8),%xmm1
movupd 0x18(%rdi,%rax,8),%xmm2
addpd  %xmm1,%xmm2
movapd %xmm0,%xmm1
divpd  %xmm2,%xmm1
movupd %xmm1,0x18(%rdx,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rbx
jne    1170 <compute_dfactor+0x70>
test   $0x1,%r11b
je     11e6 <compute_dfactor+0xe6>
lea    0x0(,%rax,8),%rbx
or     $0x8,%rbx
movupd (%rdi,%rbx,1),%xmm0
movupd (%rdi,%rax,8),%xmm1
addpd  %xmm0,%xmm1
movapd 0xe23(%rip),%xmm0
divpd  %xmm1,%xmm0
movupd %xmm0,(%rdx,%rbx,1)
cmp    %r10,%r8
je     126a <compute_dfactor+0x16a>
or     $0x1,%r8
mov    %r8,%rax
mov    %rax,%rbx
not    %rbx
test   $0x1,%cl
jne    121d <compute_dfactor+0x11d>
movsd  (%rdi,%rax,8),%xmm0
addsd  -0x8(%rdi,%rax,8),%xmm0
movsd  0xe00(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rax,8)
add    $0x1,%rax
add    %rcx,%rbx
je     126a <compute_dfactor+0x16a>
movsd  0xde6(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm1
addsd  -0x8(%rdi,%rax,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rdx,%rax,8)
movsd  0x8(%rdi,%rax,8),%xmm1
addsd  (%rdi,%rax,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%rdx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rcx
jne    1230 <compute_dfactor+0x130>
movsd  0xd9e(%rip),%xmm0
movapd %xmm0,%xmm1
divsd  (%rdi),%xmm1
movsd  %xmm1,(%rdx)
movslq %esi,%rax
divsd  -0x10(%rdi,%rax,8),%xmm0
movslq %r9d,%rax
movsd  %xmm0,(%rdx,%rax,8)
pop    %rbx
ret
xor    %eax,%eax
test   $0x1,%r11b
jne    11bb <compute_dfactor+0xbb>
jmp    11e6 <compute_dfactor+0xe6>

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

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