#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double const fabs(double);
double cityblock(int n, double **data1, double **data2, int **mask1,
                 int **mask2, const double weight[], int index1, int index2,
                 int transpose) {
  double result = 0.;
  double tweight = 0;
  int i;
  if (transpose == 0) {
    for (i = 0; i < n; i++) {
      if (mask1[index1][i] && mask2[index2][i]) {
        double term = data1[index1][i] - data2[index2][i];
        result = result + weight[i] * fabs(term);
        tweight += weight[i];
      }
    }
  } else {
    for (i = 0; i < n; i++) {
      if (mask1[i][index1] && mask2[i][index2]) {
        double term = data1[i][index1] - data2[i][index2];
        result = result + weight[i] * fabs(term);
        tweight += weight[i];
      }
    }
  }
  if (!tweight)
    return 0;
  result /= tweight;
  return result;
}


x86 assembly (clang, -O0):
<cityblock>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
cmpl   $0x0,0x20(%rbp)
jne    123a <cityblock+0x12a>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1235 <cityblock+0x125>
mov    -0x28(%rbp),%rax
movslq 0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1222 <cityblock+0x112>
mov    -0x30(%rbp),%rax
movslq 0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1222 <cityblock+0x112>
mov    -0x18(%rbp),%rax
movslq 0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq 0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
call   1030 <fabs@plt>
movsd  -0x70(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x68(%rbp),%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    1227 <cityblock+0x117>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    115c <cityblock+0x4c>
jmp    131f <cityblock+0x20f>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    131a <cityblock+0x20a>
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq 0x10(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1307 <cityblock+0x1f7>
mov    -0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq 0x18(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1307 <cityblock+0x1f7>
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq 0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq 0x18(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x60(%rbp),%xmm0
call   1030 <fabs@plt>
movsd  -0x80(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x78(%rbp),%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    130c <cityblock+0x1fc>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1241 <cityblock+0x131>
jmp    131f <cityblock+0x20f>
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1344 <cityblock+0x234>
jp     1344 <cityblock+0x234>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1361 <cityblock+0x251>
movsd  -0x48(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cityblock>:
push   %rbx
mov    0x18(%rsp),%ebx
mov    0x10(%rsp),%eax
cmpl   $0x0,0x20(%rsp)
je     118d <cityblock+0x8d>
test   %edi,%edi
jle    1204 <cityblock+0x104>
cltq
movslq %ebx,%r10
mov    %edi,%r11d
xorpd  %xmm1,%xmm1
xor    %edi,%edi
movapd 0xed2(%rip),%xmm0
xorpd  %xmm2,%xmm2
jmp    114d <cityblock+0x4d>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdi
cmp    %rdi,%r11
je     120c <cityblock+0x10c>
mov    (%rcx,%rdi,8),%rbx
cmpl   $0x0,(%rbx,%rax,4)
je     1140 <cityblock+0x40>
mov    (%r8,%rdi,8),%rbx
cmpl   $0x0,(%rbx,%r10,4)
je     1140 <cityblock+0x40>
mov    (%rsi,%rdi,8),%rbx
movsd  (%rbx,%rax,8),%xmm3
mov    (%rdx,%rdi,8),%rbx
subsd  (%rbx,%r10,8),%xmm3
andpd  %xmm0,%xmm3
movsd  (%r9,%rdi,8),%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm3,%xmm1
addsd  %xmm4,%xmm2
jmp    1140 <cityblock+0x40>
test   %edi,%edi
jle    1204 <cityblock+0x104>
movslq %eax,%r10
mov    (%rcx,%r10,8),%r11
movslq %ebx,%rax
mov    %edi,%edi
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
movapd 0xe55(%rip),%xmm0
xorpd  %xmm2,%xmm2
jmp    11c9 <cityblock+0xc9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%rcx
cmp    %rcx,%rdi
je     120c <cityblock+0x10c>
cmpl   $0x0,(%r11,%rcx,4)
je     11c0 <cityblock+0xc0>
mov    (%r8,%rax,8),%rbx
cmpl   $0x0,(%rbx,%rcx,4)
je     11c0 <cityblock+0xc0>
mov    (%rsi,%r10,8),%rbx
movsd  (%rbx,%rcx,8),%xmm3
mov    (%rdx,%rax,8),%rbx
subsd  (%rbx,%rcx,8),%xmm3
andpd  %xmm0,%xmm3
movsd  (%r9,%rcx,8),%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm3,%xmm1
addsd  %xmm4,%xmm2
jmp    11c0 <cityblock+0xc0>
xorpd  %xmm2,%xmm2
xorpd  %xmm1,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm2
jne    1218 <cityblock+0x118>
jnp    1220 <cityblock+0x120>
divsd  %xmm2,%xmm1
movapd %xmm1,%xmm0
pop    %rbx
ret

LLVM IR (-O0):
define double @cityblock(i32 %0, double** %1, double** %2, i32** %3, i32** %4, double* %5, i32 %6, i32 %7, i32 %8) {
%10 = alloca double
%11 = alloca i32
%12 = alloca double**
%13 = alloca double**
%14 = alloca i32**
%15 = alloca i32**
%16 = alloca double*
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca double
%21 = alloca double
%22 = alloca i32
%23 = alloca double
%24 = alloca double
store i32 %0, i32* %11
store double** %1, double*** %12
store double** %2, double*** %13
store i32** %3, i32*** %14
store i32** %4, i32*** %15
store double* %5, double** %16
store i32 %6, i32* %17
store i32 %7, i32* %18
store i32 %8, i32* %19
store double 0.000000e+00, double* %20
store double 0.000000e+00, double* %21
%25 = load i32, i32* %19
%26 = icmp eq i32 %25, 0
br i1 %26, label %27, label %95
27: ; preds = %9
store i32 0, i32* %22
br label %28
28: ; preds = %91, %27
%29 = load i32, i32* %22
%30 = load i32, i32* %11
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %94
32: ; preds = %28
%33 = load i32**, i32*** %14
%34 = load i32, i32* %17
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds i32*, i32** %33, i64 %35
%37 = load i32*, i32** %36
%38 = load i32, i32* %22
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds i32, i32* %37, i64 %39
%41 = load i32, i32* %40
%42 = icmp ne i32 %41, 0
br i1 %42, label %43, label %90
43: ; preds = %32
%44 = load i32**, i32*** %15
%45 = load i32, i32* %18
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds i32*, i32** %44, i64 %46
%48 = load i32*, i32** %47
%49 = load i32, i32* %22
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds i32, i32* %48, i64 %50
%52 = load i32, i32* %51
%53 = icmp ne i32 %52, 0
br i1 %53, label %54, label %90
54: ; preds = %43
%55 = load double**, double*** %12
%56 = load i32, i32* %17
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double*, double** %55, i64 %57
%59 = load double*, double** %58
%60 = load i32, i32* %22
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %59, i64 %61
%63 = load double, double* %62
%64 = load double**, double*** %13
%65 = load i32, i32* %18
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double*, double** %64, i64 %66
%68 = load double*, double** %67
%69 = load i32, i32* %22
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %68, i64 %70
%72 = load double, double* %71
%73 = fsub double %63, %72
store double %73, double* %23
%74 = load double, double* %20
%75 = load double*, double** %16
%76 = load i32, i32* %22
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = load double, double* %23
%81 = call double @fabs(double %80)
%82 = call double @llvm.fmuladd.f64(double %79, double %81, double %74)
store double %82, double* %20
%83 = load double*, double** %16
%84 = load i32, i32* %22
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = load double, double* %21
%89 = fadd double %88, %87
store double %89, double* %21
br label %90
90: ; preds = %54, %43, %32
br label %91
91: ; preds = %90
%92 = load i32, i32* %22
%93 = add nsw i32 %92, 1
store i32 %93, i32* %22
br label %28
94: ; preds = %28
br label %163
95: ; preds = %9
store i32 0, i32* %22
br label %96
96: ; preds = %159, %95
%97 = load i32, i32* %22
%98 = load i32, i32* %11
%99 = icmp slt i32 %97, %98
br i1 %99, label %100, label %162
100: ; preds = %96
%101 = load i32**, i32*** %14
%102 = load i32, i32* %22
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds i32*, i32** %101, i64 %103
%105 = load i32*, i32** %104
%106 = load i32, i32* %17
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds i32, i32* %105, i64 %107
%109 = load i32, i32* %108
%110 = icmp ne i32 %109, 0
br i1 %110, label %111, label %158
111: ; preds = %100
%112 = load i32**, i32*** %15
%113 = load i32, i32* %22
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds i32*, i32** %112, i64 %114
%116 = load i32*, i32** %115
%117 = load i32, i32* %18
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds i32, i32* %116, i64 %118
%120 = load i32, i32* %119
%121 = icmp ne i32 %120, 0
br i1 %121, label %122, label %158
122: ; preds = %111
%123 = load double**, double*** %12
%124 = load i32, i32* %22
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double*, double** %123, i64 %125
%127 = load double*, double** %126
%128 = load i32, i32* %17
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %127, i64 %129
%131 = load double, double* %130
%132 = load double**, double*** %13
%133 = load i32, i32* %22
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double*, double** %132, i64 %134
%136 = load double*, double** %135
%137 = load i32, i32* %18
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds double, double* %136, i64 %138
%140 = load double, double* %139
%141 = fsub double %131, %140
store double %141, double* %24
%142 = load double, double* %20
%143 = load double*, double** %16
%144 = load i32, i32* %22
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds double, double* %143, i64 %145
%147 = load double, double* %146
%148 = load double, double* %24
%149 = call double @fabs(double %148)
%150 = call double @llvm.fmuladd.f64(double %147, double %149, double %142)
store double %150, double* %20
%151 = load double*, double** %16
%152 = load i32, i32* %22
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds double, double* %151, i64 %153
%155 = load double, double* %154
%156 = load double, double* %21
%157 = fadd double %156, %155
store double %157, double* %21
br label %158
158: ; preds = %122, %111, %100
br label %159
159: ; preds = %158
%160 = load i32, i32* %22
%161 = add nsw i32 %160, 1
store i32 %161, i32* %22
br label %96
162: ; preds = %96
br label %163
163: ; preds = %162, %94
%164 = load double, double* %21
%165 = fcmp une double %164, 0.000000e+00
br i1 %165, label %167, label %166
166: ; preds = %163
store double 0.000000e+00, double* %10
br label %172
167: ; preds = %163
%168 = load double, double* %21
%169 = load double, double* %20
%170 = fdiv double %169, %168
store double %170, double* %20
%171 = load double, double* %20
store double %171, double* %10
br label %172
172: ; preds = %167, %166
%173 = load double, double* %10
ret double %173
}

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