double MAD(double *a, int n, int center, double *b, double *tmp) {
  double median_abs(double *, int, double *);
  double median(double *, int, double *);
  int i;
  double med, q;
  if (center > 0) {
    med = median(a, n, b);
  } else {
    med = 0.0;
  };
  for (i = 0; i < n; i++)
    b[i] = a[i] - med;
  q = median_abs(b, n, tmp) * 1.4826;
  return (q);
}


x86 assembly (clang, -O0):
<MAD>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
cmpl   $0x0,-0x10(%rbp)
jle    115e <MAD+0x3e>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x18(%rbp),%rdx
call   1030 <median@plt>
movsd  %xmm0,-0x30(%rbp)
jmp    1166 <MAD+0x46>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11a6 <MAD+0x86>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    116d <MAD+0x4d>
mov    -0x18(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x20(%rbp),%rdx
call   1040 <median_abs@plt>
movsd  0xe42(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<MAD>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %r8,%r14
mov    %rcx,%rbx
mov    %esi,%r15d
mov    %rdi,%rbp
test   %edx,%edx
jle    114f <MAD+0x2f>
mov    %rbp,%rdi
mov    %r15d,%esi
mov    %rbx,%rdx
call   1030 <median@plt>
test   %r15d,%r15d
jg     115c <MAD+0x3c>
jmp    12c7 <MAD+0x1a7>
xorpd  %xmm0,%xmm0
test   %r15d,%r15d
jle    12c7 <MAD+0x1a7>
mov    %r15d,%eax
cmp    $0x4,%r15d
jb     1186 <MAD+0x66>
lea    0x0(,%rax,8),%rcx
add    %rbp,%rcx
cmp    %rbx,%rcx
jbe    121d <MAD+0xfd>
lea    (%rbx,%rax,8),%rcx
cmp    %rbp,%rcx
jbe    121d <MAD+0xfd>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rsi
and    $0x3,%rsi
je     11b9 <MAD+0x99>
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rcx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    11a0 <MAD+0x80>
cmp    $0x3,%rdx
jb     12c7 <MAD+0x1a7>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%rbp,%rcx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
movsd  0x8(%rbp,%rcx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rbx,%rcx,8)
movsd  0x10(%rbp,%rcx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rbx,%rcx,8)
movsd  0x18(%rbp,%rcx,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rbx,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    11d0 <MAD+0xb0>
jmp    12c7 <MAD+0x1a7>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rcx),%rdx
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rdx,%rdx
je     12e8 <MAD+0x1c8>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movupd 0x0(%rbp,%rdx,8),%xmm2
movupd 0x10(%rbp,%rdx,8),%xmm3
subpd  %xmm1,%xmm2
subpd  %xmm1,%xmm3
movupd %xmm2,(%rbx,%rdx,8)
movupd %xmm3,0x10(%rbx,%rdx,8)
movupd 0x20(%rbp,%rdx,8),%xmm2
movupd 0x30(%rbp,%rdx,8),%xmm3
subpd  %xmm1,%xmm2
subpd  %xmm1,%xmm3
movupd %xmm2,0x20(%rbx,%rdx,8)
movupd %xmm3,0x30(%rbx,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1250 <MAD+0x130>
test   $0x1,%sil
je     12be <MAD+0x19e>
movupd 0x0(%rbp,%rdx,8),%xmm2
movupd 0x10(%rbp,%rdx,8),%xmm3
subpd  %xmm1,%xmm2
subpd  %xmm1,%xmm3
movupd %xmm2,(%rbx,%rdx,8)
movupd %xmm3,0x10(%rbx,%rdx,8)
cmp    %rax,%rcx
jne    1188 <MAD+0x68>
mov    %rbx,%rdi
mov    %r15d,%esi
mov    %r14,%rdx
call   1040 <median_abs@plt>
mulsd  0xd23(%rip),%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%sil
jne    129f <MAD+0x17f>
jmp    12be <MAD+0x19e>

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

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