int ti_ad(int size, double const *const *inputs, double const *options,
          double *const *outputs) {
  const double *high = inputs[0];
  const double *low = inputs[1];
  const double *close = inputs[2];
  const double *volume = inputs[3];
  (void)options;
  double *output = outputs[0];
  double sum = 0;
  int i;
  for (i = 0; i < size; ++i) {
    const double hl = (high[i] - low[i]);
    if (hl != 0.0) {
      sum += (close[i] - low[i] - high[i] + close[i]) / hl * volume[i];
    }
    output[i] = sum;
  }
  return 0;
}


x86 assembly (clang, -O0):
<ti_ad>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    0x10(%rax),%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    0x18(%rax),%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    121c <ti_ad+0x11c>
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11a4 <ti_ad+0xa4>
jp     11a4 <ti_ad+0xa4>
jmp    11fc <ti_ad+0xfc>
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
divsd  -0x60(%rbp),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    115c <ti_ad+0x5c>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<ti_ad>:
test   %edi,%edi
jle    117e <ti_ad+0x7e>
mov    (%rsi),%rax
mov    0x8(%rsi),%rdx
mov    0x10(%rsi),%r8
mov    0x18(%rsi),%r9
mov    (%rcx),%rcx
mov    %edi,%edi
xorpd  %xmm0,%xmm0
xor    %esi,%esi
xorpd  %xmm1,%xmm1
jmp    113e <ti_ad+0x3e>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm1,(%rcx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rdi
je     117e <ti_ad+0x7e>
movsd  (%rax,%rsi,8),%xmm3
movsd  (%rdx,%rsi,8),%xmm4
movapd %xmm3,%xmm2
subsd  %xmm4,%xmm2
ucomisd %xmm0,%xmm2
jne    1158 <ti_ad+0x58>
jnp    1130 <ti_ad+0x30>
movsd  (%r8,%rsi,8),%xmm5
movapd %xmm5,%xmm6
subsd  %xmm4,%xmm6
subsd  %xmm3,%xmm6
addsd  %xmm5,%xmm6
divsd  %xmm2,%xmm6
mulsd  (%r9,%rsi,8),%xmm6
addsd  %xmm6,%xmm1
jmp    1130 <ti_ad+0x30>
xor    %eax,%eax
ret

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

LLVM IR (-O2):
define i32 @ti_ad(i32 %0, double** nocapture %1, double* nocapture readnone %2, double** nocapture %3) {
%5 = load double*, double** %1
%6 = getelementptr inbounds double*, double** %1, i64 1
%7 = load double*, double** %6
%8 = getelementptr inbounds double*, double** %1, i64 2
%9 = load double*, double** %8
%10 = getelementptr inbounds double*, double** %1, i64 3
%11 = load double*, double** %10
%12 = load double*, double** %3
%13 = icmp sgt i32 %0, 0
br i1 %13, label %14, label %40
14: ; preds = %4
%15 = zext i32 %0 to i64
br label %16
16: ; preds = %14, %35
%17 = phi i64 [ 0, %14 ], [ %38, %35 ]
%18 = phi double [ 0.000000e+00, %14 ], [ %36, %35 ]
%19 = getelementptr inbounds double, double* %5, i64 %17
%20 = load double, double* %19
%21 = getelementptr inbounds double, double* %7, i64 %17
%22 = load double, double* %21
%23 = fsub double %20, %22
%24 = fcmp une double %23, 0.000000e+00
br i1 %24, label %25, label %35
25: ; preds = %16
%26 = getelementptr inbounds double, double* %9, i64 %17
%27 = load double, double* %26
%28 = fsub double %27, %22
%29 = fsub double %28, %20
%30 = fadd double %27, %29
%31 = fdiv double %30, %23
%32 = getelementptr inbounds double, double* %11, i64 %17
%33 = load double, double* %32
%34 = tail call double @llvm.fmuladd.f64(double %31, double %33, double %18)
br label %35
35: ; preds = %25, %16
%36 = phi double [ %34, %25 ], [ %18, %16 ]
%37 = getelementptr inbounds double, double* %12, i64 %17
store double %36, double* %37
%38 = add nuw nsw i64 %17, 1
%39 = icmp eq i64 %38, %15
br i1 %39, label %40, label %16
40: ; preds = %35, %4
ret i32 0
}