double lnLTaylor(int nb, double z[], double g[], double H[]) {
  int i, j;
  double lnL = 0;
  for (i = 0; i < nb; i++)
    lnL += z[i] * g[i];
  for (i = 0; i < nb; i++) {
    lnL += z[i] * z[i] * H[i * nb + i] / 2;
    for (j = 0; j < i; j++)
      lnL += z[i] * z[j] * H[i * nb + j];
  }
  return (lnL);
}


x86 assembly (clang, -O0):
<lnLTaylor>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1168 <lnLTaylor+0x68>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1122 <lnLTaylor+0x22>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1237 <lnLTaylor+0x137>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xe4d(%rip),%xmm1
divsd  %xmm1,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1224 <lnLTaylor+0x124>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11c8 <lnLTaylor+0xc8>
jmp    1229 <lnLTaylor+0x129>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    116f <lnLTaylor+0x6f>
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<lnLTaylor>:
test   %edi,%edi
jle    1120 <lnLTaylor+0x20>
mov    %edi,%r8d
lea    -0x1(%r8),%rax
mov    %r8d,%r9d
and    $0x3,%r9d
cmp    $0x3,%rax
jae    1124 <lnLTaylor+0x24>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    1187 <lnLTaylor+0x87>
xorps  %xmm0,%xmm0
ret
mov    %r8d,%r10d
and    $0xfffffffc,%r10d
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rax,8),%xmm1
movsd  0x8(%rsi,%rax,8),%xmm2
mulsd  (%rdx,%rax,8),%xmm1
mulsd  0x8(%rdx,%rax,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x10(%rsi,%rax,8),%xmm3
mulsd  0x10(%rdx,%rax,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x18(%rsi,%rax,8),%xmm0
mulsd  0x18(%rdx,%rax,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rax
cmp    %rax,%r10
jne    1140 <lnLTaylor+0x40>
test   %r9,%r9
je     11b8 <lnLTaylor+0xb8>
lea    (%rdx,%rax,8),%r10
lea    (%rsi,%rax,8),%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm1
mulsd  (%r10,%rdx,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%r9
jne    11a0 <lnLTaylor+0xa0>
test   %edi,%edi
jle    1297 <lnLTaylor+0x197>
lea    0x8(%rcx),%rdx
lea    0x0(,%r8,8),%r9
xor    %r10d,%r10d
movsd  0xe29(%rip),%xmm1
jmp    11f0 <lnLTaylor+0xf0>
nopl   0x0(%rax)
add    $0x1,%r10
add    %r9,%rdx
cmp    %r8,%r10
je     1297 <lnLTaylor+0x197>
movapd %xmm0,%xmm3
movsd  (%rsi,%r10,8),%xmm2
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mov    %r10d,%eax
imul   %edi,%eax
cltq
add    %r10,%rax
mulsd  (%rcx,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
test   %r10,%r10
je     11e0 <lnLTaylor+0xe0>
cmp    $0x1,%r10
jne    1230 <lnLTaylor+0x130>
xor    %eax,%eax
jmp    126f <lnLTaylor+0x16f>
nopl   0x0(%rax)
mov    %r10,%r11
and    $0xfffffffffffffffe,%r11
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%rsi,%rax,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  -0x8(%rdx,%rax,8),%xmm3
addsd  %xmm0,%xmm3
movsd  0x8(%rsi,%rax,8),%xmm0
mulsd  %xmm2,%xmm0
mulsd  (%rdx,%rax,8),%xmm0
addsd  %xmm3,%xmm0
add    $0x2,%rax
cmp    %rax,%r11
jne    1240 <lnLTaylor+0x140>
test   $0x1,%r10b
je     11e0 <lnLTaylor+0xe0>
mov    %r10,%r11
imul   %r8,%r11
mulsd  (%rsi,%rax,8),%xmm2
add    %rax,%r11
mulsd  (%rcx,%r11,8),%xmm2
addsd  %xmm2,%xmm0
jmp    11e0 <lnLTaylor+0xe0>
ret

LLVM IR (-O0):
define double @lnLTaylor(i32 %0, double* %1, double* %2, double* %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
store i32 %0, i32* %5
store double* %1, double** %6
store double* %2, double** %7
store double* %3, double** %8
store double 0.000000e+00, double* %11
store i32 0, i32* %9
br label %12
12: ; preds = %29, %4
%13 = load i32, i32* %9
%14 = load i32, i32* %5
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %32
16: ; preds = %12
%17 = load double*, double** %6
%18 = load i32, i32* %9
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double, double* %17, i64 %19
%21 = load double, double* %20
%22 = load double*, double** %7
%23 = load i32, i32* %9
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
%26 = load double, double* %25
%27 = load double, double* %11
%28 = call double @llvm.fmuladd.f64(double %21, double %26, double %27)
store double %28, double* %11
br label %29
29: ; preds = %16
%30 = load i32, i32* %9
%31 = add nsw i32 %30, 1
store i32 %31, i32* %9
br label %12
32: ; preds = %12
store i32 0, i32* %9
br label %33
33: ; preds = %93, %32
%34 = load i32, i32* %9
%35 = load i32, i32* %5
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %96
37: ; preds = %33
%38 = load double*, double** %6
%39 = load i32, i32* %9
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = load double*, double** %6
%44 = load i32, i32* %9
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %43, i64 %45
%47 = load double, double* %46
%48 = fmul double %42, %47
%49 = load double*, double** %8
%50 = load i32, i32* %9
%51 = load i32, i32* %5
%52 = mul nsw i32 %50, %51
%53 = load i32, i32* %9
%54 = add nsw i32 %52, %53
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %49, i64 %55
%57 = load double, double* %56
%58 = fmul double %48, %57
%59 = fdiv double %58, 2.000000e+00
%60 = load double, double* %11
%61 = fadd double %60, %59
store double %61, double* %11
store i32 0, i32* %10
br label %62
62: ; preds = %89, %37
%63 = load i32, i32* %10
%64 = load i32, i32* %9
%65 = icmp slt i32 %63, %64
br i1 %65, label %66, label %92
66: ; preds = %62
%67 = load double*, double** %6
%68 = load i32, i32* %9
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
%71 = load double, double* %70
%72 = load double*, double** %6
%73 = load i32, i32* %10
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = fmul double %71, %76
%78 = load double*, double** %8
%79 = load i32, i32* %9
%80 = load i32, i32* %5
%81 = mul nsw i32 %79, %80
%82 = load i32, i32* %10
%83 = add nsw i32 %81, %82
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %78, i64 %84
%86 = load double, double* %85
%87 = load double, double* %11
%88 = call double @llvm.fmuladd.f64(double %77, double %86, double %87)
store double %88, double* %11
br label %89
89: ; preds = %66
%90 = load i32, i32* %10
%91 = add nsw i32 %90, 1
store i32 %91, i32* %10
br label %62
92: ; preds = %62
br label %93
93: ; preds = %92
%94 = load i32, i32* %9
%95 = add nsw i32 %94, 1
store i32 %95, i32* %9
br label %33
96: ; preds = %33
%97 = load double, double* %11
ret double %97
}

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