void PT_div(double *sum, double *p1, double *p2, int n) {
  int i;
  for (i = 0; i < n; i++) {
    p2[i] = 1.0 / p2[i];
    sum[i] = p1[i] * p2[i];
  }
}


x86 assembly (clang, -O0):
<PT_div>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    117d <PT_div+0x7d>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  0xeca(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    111a <PT_div+0x1a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<PT_div>:
push   %rbp
push   %r15
push   %r14
push   %rbx
test   %ecx,%ecx
jle    12b1 <PT_div+0x1b1>
mov    %ecx,%r8d
cmp    $0x1,%ecx
jne    11a2 <PT_div+0xa2>
xor    %ecx,%ecx
mov    %rcx,%rax
not    %rax
test   $0x1,%r8b
je     1148 <PT_div+0x48>
movsd  0xee0(%rip),%xmm0
divsd  (%rdx,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
mulsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
or     $0x1,%rcx
add    %r8,%rax
je     12b1 <PT_div+0x1b1>
movsd  0xeb7(%rip),%xmm0
nopl   0x0(%rax)
movapd %xmm0,%xmm1
divsd  (%rdx,%rcx,8),%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
mulsd  (%rsi,%rcx,8),%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
movapd %xmm0,%xmm1
divsd  0x8(%rdx,%rcx,8),%xmm1
movsd  %xmm1,0x8(%rdx,%rcx,8)
mulsd  0x8(%rsi,%rcx,8),%xmm1
movsd  %xmm1,0x8(%rdi,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r8
jne    1160 <PT_div+0x60>
jmp    12b1 <PT_div+0x1b1>
lea    (%rdi,%r8,8),%rcx
cmp    %rdx,%rcx
seta   %r15b
lea    (%rdx,%r8,8),%rbx
cmp    %rdi,%rbx
seta   %r14b
lea    (%rsi,%r8,8),%rbp
cmp    %rdx,%rbp
seta   %al
cmp    %rsi,%rbx
seta   %r10b
cmp    %rdi,%rbp
seta   %r9b
cmp    %rsi,%rcx
seta   %r11b
xor    %ecx,%ecx
test   %r14b,%r15b
jne    111c <PT_div+0x1c>
and    %r10b,%al
jne    111c <PT_div+0x1c>
and    %r11b,%r9b
jne    111c <PT_div+0x1c>
mov    %r8d,%ecx
and    $0xfffffffe,%ecx
lea    -0x2(%rcx),%rax
mov    %rax,%r9
shr    %r9
add    $0x1,%r9
test   %rax,%rax
je     12b8 <PT_div+0x1b8>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %eax,%eax
movapd 0xdde(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rdx,%rax,8),%xmm1
movapd %xmm0,%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdx,%rax,8)
movupd (%rsi,%rax,8),%xmm1
mulpd  %xmm2,%xmm1
movupd %xmm1,(%rdi,%rax,8)
movupd 0x10(%rdx,%rax,8),%xmm1
movapd %xmm0,%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,0x10(%rdx,%rax,8)
movupd 0x10(%rsi,%rax,8),%xmm1
mulpd  %xmm2,%xmm1
movupd %xmm1,0x10(%rdi,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%r10
jne    1230 <PT_div+0x130>
test   $0x1,%r9b
je     12a8 <PT_div+0x1a8>
movupd (%rdx,%rax,8),%xmm0
movapd 0xd6f(%rip),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%rdx,%rax,8)
movupd (%rsi,%rax,8),%xmm0
mulpd  %xmm1,%xmm0
movupd %xmm0,(%rdi,%rax,8)
cmp    %r8,%rcx
jne    111c <PT_div+0x1c>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %eax,%eax
test   $0x1,%r9b
jne    1284 <PT_div+0x184>
jmp    12a8 <PT_div+0x1a8>

LLVM IR (-O0):
define void @PT_div(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 = %40, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %8
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %43
14: ; preds = %10
%15 = load double*, double** %7
%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 = fdiv double 1.000000e+00, %19
%21 = load double*, double** %7
%22 = load i32, i32* %9
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %21, i64 %23
store double %20, double* %24
%25 = load double*, double** %6
%26 = load i32, i32* %9
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load double*, double** %7
%31 = load i32, i32* %9
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
%35 = fmul double %29, %34
%36 = load double*, double** %5
%37 = load i32, i32* %9
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %36, i64 %38
store double %35, double* %39
br label %40
40: ; preds = %14
%41 = load i32, i32* %9
%42 = add nsw i32 %41, 1
store i32 %42, i32* %9
br label %10
43: ; preds = %10
ret void
}

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