double PolygonArea_old(int np, double *xp, double *yp) {
  int i, j;
  double area = 0;
  for (i = 0; i < np; i++) {
    j = (i + 1) % np;
    area += xp[i] * yp[j];
    area -= yp[i] * xp[j];
  }
  area /= 2;
  return (area);
}


x86 assembly (clang, -O0):
<PolygonArea_old>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b4 <PolygonArea_old+0xb4>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
cltd
idivl  -0x4(%rbp)
mov    %edx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    111e <PolygonArea_old+0x1e>
movsd  -0x28(%rbp),%xmm0
movsd  0xe3f(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<PolygonArea_old>:
test   %edi,%edi
jle    1115 <PolygonArea_old+0x15>
push   %rbx
mov    %edi,%r8d
cmp    $0x1,%edi
jne    1119 <PolygonArea_old+0x19>
xorpd  %xmm0,%xmm0
xor    %edi,%edi
jmp    1190 <PolygonArea_old+0x90>
xorps  %xmm0,%xmm0
ret
mov    %r8d,%ecx
and    $0xfffffffe,%ecx
lea    -0x1(%r8),%r9
lea    -0x2(%r8),%r10
xorpd  %xmm1,%xmm1
xor    %r11d,%r11d
xor    %eax,%eax
lea    0x1(%rax),%edi
cmp    %rax,%r9
cmove  %r11,%rdi
movsd  (%rsi,%rax,8),%xmm2
movsd  0x8(%rsi,%rax,8),%xmm0
mulsd  (%rdx,%rdi,8),%xmm2
addsd  %xmm1,%xmm2
movsd  (%rdx,%rax,8),%xmm1
movsd  0x8(%rdx,%rax,8),%xmm3
mulsd  (%rsi,%rdi,8),%xmm1
lea    0x2(%rax),%rdi
mov    %edi,%ebx
and    $0xfffffffe,%ebx
cmp    %rax,%r10
cmove  %r11,%rbx
subsd  %xmm1,%xmm2
mulsd  (%rdx,%rbx,8),%xmm0
addsd  %xmm2,%xmm0
mulsd  (%rsi,%rbx,8),%xmm3
subsd  %xmm3,%xmm0
mov    %rdi,%rax
movapd %xmm0,%xmm1
cmp    %rdi,%rcx
jne    1130 <PolygonArea_old+0x30>
test   $0x1,%r8b
je     11c5 <PolygonArea_old+0xc5>
lea    0x1(%rdi),%rax
xor    %ecx,%ecx
cmp    %r8,%rax
mov    %eax,%eax
cmovne %rax,%rcx
movsd  (%rsi,%rdi,8),%xmm1
mulsd  (%rdx,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  (%rdx,%rdi,8),%xmm0
mulsd  (%rsi,%rcx,8),%xmm0
subsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
mulsd  0xe33(%rip),%xmm0
pop    %rbx
ret

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

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