double Gcmp_d(const void *i, const void *j, void *arg) {
  const int *Gi = (int *)i;
  const int *Gj = (int *)j;
  const double *w = (const double *)arg;
  int dv[3] = {Gi[0] * Gi[0] - Gj[0] * Gj[0], Gi[0] * Gi[1] - Gj[0] * Gj[1],
               Gi[1] * Gi[1] - Gj[1] * Gj[1]};
  return dv[0] * w[0] + dv[1] * w[1] + dv[2] * w[2];
}


x86 assembly (clang, -O0):
<Gcmp_d>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    -0x20(%rbp),%rcx
imul   (%rcx),%eax
mov    -0x28(%rbp),%rcx
mov    (%rcx),%ecx
mov    -0x28(%rbp),%rdx
imul   (%rdx),%ecx
sub    %ecx,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    -0x20(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x28(%rbp),%rcx
mov    (%rcx),%ecx
mov    -0x28(%rbp),%rdx
imul   0x4(%rdx),%ecx
sub    %ecx,%eax
mov    %eax,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    0x4(%rax),%eax
mov    -0x20(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x28(%rbp),%rcx
mov    0x4(%rcx),%ecx
mov    -0x28(%rbp),%rdx
imul   0x4(%rdx),%ecx
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
cvtsi2sdl -0x3c(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
cvtsi2sdl -0x38(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
cvtsi2sdl -0x34(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<Gcmp_d>:
mov    (%rdi),%eax
mov    0x4(%rdi),%r8d
mov    %eax,%r9d
imul   %eax,%r9d
mov    (%rsi),%edi
mov    0x4(%rsi),%esi
mov    %edi,%ecx
imul   %edi,%ecx
sub    %ecx,%r9d
imul   %r8d,%eax
imul   %esi,%edi
sub    %edi,%eax
imul   %r8d,%r8d
imul   %esi,%esi
sub    %esi,%r8d
cvtsi2sd %r9d,%xmm1
cvtsi2sd %eax,%xmm2
mulsd  0x8(%rdx),%xmm2
mulsd  (%rdx),%xmm1
cvtsi2sd %r8d,%xmm0
addsd  %xmm2,%xmm1
mulsd  0x10(%rdx),%xmm0
addsd  %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @Gcmp_d(i8* %0, i8* %1, i8* %2) {
%4 = alloca i8*
%5 = alloca i8*
%6 = alloca i8*
%7 = alloca i32*
%8 = alloca i32*
%9 = alloca double*
%10 = alloca [3 x i32]
store i8* %0, i8** %4
store i8* %1, i8** %5
store i8* %2, i8** %6
%11 = load i8*, i8** %4
%12 = bitcast i8* %11 to i32*
store i32* %12, i32** %7
%13 = load i8*, i8** %5
%14 = bitcast i8* %13 to i32*
store i32* %14, i32** %8
%15 = load i8*, i8** %6
%16 = bitcast i8* %15 to double*
store double* %16, double** %9
%17 = getelementptr inbounds [3 x i32], [3 x i32]* %10, i64 0, i64 0
%18 = load i32*, i32** %7
%19 = getelementptr inbounds i32, i32* %18, i64 0
%20 = load i32, i32* %19
%21 = load i32*, i32** %7
%22 = getelementptr inbounds i32, i32* %21, i64 0
%23 = load i32, i32* %22
%24 = mul nsw i32 %20, %23
%25 = load i32*, i32** %8
%26 = getelementptr inbounds i32, i32* %25, i64 0
%27 = load i32, i32* %26
%28 = load i32*, i32** %8
%29 = getelementptr inbounds i32, i32* %28, i64 0
%30 = load i32, i32* %29
%31 = mul nsw i32 %27, %30
%32 = sub nsw i32 %24, %31
store i32 %32, i32* %17
%33 = getelementptr inbounds i32, i32* %17, i64 1
%34 = load i32*, i32** %7
%35 = getelementptr inbounds i32, i32* %34, i64 0
%36 = load i32, i32* %35
%37 = load i32*, i32** %7
%38 = getelementptr inbounds i32, i32* %37, i64 1
%39 = load i32, i32* %38
%40 = mul nsw i32 %36, %39
%41 = load i32*, i32** %8
%42 = getelementptr inbounds i32, i32* %41, i64 0
%43 = load i32, i32* %42
%44 = load i32*, i32** %8
%45 = getelementptr inbounds i32, i32* %44, i64 1
%46 = load i32, i32* %45
%47 = mul nsw i32 %43, %46
%48 = sub nsw i32 %40, %47
store i32 %48, i32* %33
%49 = getelementptr inbounds i32, i32* %33, i64 1
%50 = load i32*, i32** %7
%51 = getelementptr inbounds i32, i32* %50, i64 1
%52 = load i32, i32* %51
%53 = load i32*, i32** %7
%54 = getelementptr inbounds i32, i32* %53, i64 1
%55 = load i32, i32* %54
%56 = mul nsw i32 %52, %55
%57 = load i32*, i32** %8
%58 = getelementptr inbounds i32, i32* %57, i64 1
%59 = load i32, i32* %58
%60 = load i32*, i32** %8
%61 = getelementptr inbounds i32, i32* %60, i64 1
%62 = load i32, i32* %61
%63 = mul nsw i32 %59, %62
%64 = sub nsw i32 %56, %63
store i32 %64, i32* %49
%65 = getelementptr inbounds [3 x i32], [3 x i32]* %10, i64 0, i64 0
%66 = load i32, i32* %65
%67 = sitofp i32 %66 to double
%68 = load double*, double** %9
%69 = getelementptr inbounds double, double* %68, i64 0
%70 = load double, double* %69
%71 = getelementptr inbounds [3 x i32], [3 x i32]* %10, i64 0, i64 1
%72 = load i32, i32* %71
%73 = sitofp i32 %72 to double
%74 = load double*, double** %9
%75 = getelementptr inbounds double, double* %74, i64 1
%76 = load double, double* %75
%77 = fmul double %73, %76
%78 = call double @llvm.fmuladd.f64(double %67, double %70, double %77)
%79 = getelementptr inbounds [3 x i32], [3 x i32]* %10, i64 0, i64 2
%80 = load i32, i32* %79
%81 = sitofp i32 %80 to double
%82 = load double*, double** %9
%83 = getelementptr inbounds double, double* %82, i64 2
%84 = load double, double* %83
%85 = call double @llvm.fmuladd.f64(double %81, double %84, double %78)
ret double %85
}

LLVM IR (-O2):
define double @Gcmp_d(i8* nocapture %0, i8* nocapture %1, i8* nocapture %2) {
%4 = bitcast i8* %0 to i32*
%5 = bitcast i8* %1 to i32*
%6 = bitcast i8* %2 to double*
%7 = load i32, i32* %4
%8 = mul nsw i32 %7, %7
%9 = load i32, i32* %5
%10 = mul nsw i32 %9, %9
%11 = sub nsw i32 %8, %10
%12 = getelementptr inbounds i8, i8* %0, i64 4
%13 = bitcast i8* %12 to i32*
%14 = load i32, i32* %13
%15 = mul nsw i32 %14, %7
%16 = getelementptr inbounds i8, i8* %1, i64 4
%17 = bitcast i8* %16 to i32*
%18 = load i32, i32* %17
%19 = mul nsw i32 %18, %9
%20 = sub nsw i32 %15, %19
%21 = mul nsw i32 %14, %14
%22 = mul nsw i32 %18, %18
%23 = sub nsw i32 %21, %22
%24 = sitofp i32 %11 to double
%25 = load double, double* %6
%26 = sitofp i32 %20 to double
%27 = getelementptr inbounds i8, i8* %2, i64 8
%28 = bitcast i8* %27 to double*
%29 = load double, double* %28
%30 = fmul double %29, %26
%31 = tail call double @llvm.fmuladd.f64(double %24, double %25, double %30)
%32 = sitofp i32 %23 to double
%33 = getelementptr inbounds i8, i8* %2, i64 16
%34 = bitcast i8* %33 to double*
%35 = load double, double* %34
%36 = tail call double @llvm.fmuladd.f64(double %32, double %35, double %31)
ret double %36
}