#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 */
int /*<<< orphan*/ nrerror(char *);
scalar_t__ strchr(char *, char);
int /*<<< orphan*/ strcpy(char *, char *);
double veclen(double *);
void get_hbond_ij(long i, long j, double *HB_UPPER, long **seidx,
                  char **AtomName, char *HB_ATOM, double **xyz, long *nh,
                  char **hb_atom1, char **hb_atom2, double *hb_dist) {
  double dd, dtmp[4];
  long k, m, n, num_hbonds = 0;
  for (m = seidx[i][1]; m <= seidx[i][2]; m++) {
    if (strchr("P", AtomName[m][1]) || strchr("P", AtomName[m][3]))
      continue;
    for (n = seidx[j][1]; n <= seidx[j][2]; n++) {
      if (strchr("P", AtomName[n][1]) || strchr("P", AtomName[n][3]))
        continue;
      for (k = 1; k <= 3; k++) {
        dtmp[k] = xyz[m][k] - xyz[n][k];
      }
      if ((dd = veclen(dtmp)) < HB_UPPER[0]) {
        if (++num_hbonds > 512)
          nrerror("Too many possible H-bonds between two bases");
        strcpy(hb_atom1[num_hbonds], AtomName[m]);
        strcpy(hb_atom2[num_hbonds], AtomName[n]);
        hb_dist[num_hbonds] = dd;
      }
    }
  }
  *nh = num_hbonds;
}


x86 assembly (clang, -O0):
<get_hbond_ij>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movq   $0x0,-0x80(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x8(%rax),%rax
mov    %rax,-0x70(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x8(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
cmp    0x10(%rcx),%rax
jg     1394 <get_hbond_ij+0x254>
mov    -0x28(%rbp),%rax
mov    -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
lea    0xe40(%rip),%rdi
movsbl 0x1(%rax),%esi
call   1050 <strchr@plt>
cmp    $0x0,%rax
jne    11f9 <get_hbond_ij+0xb9>
mov    -0x28(%rbp),%rax
mov    -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
lea    0xe1a(%rip),%rdi
movsbl 0x3(%rax),%esi
call   1050 <strchr@plt>
cmp    $0x0,%rax
je     11fe <get_hbond_ij+0xbe>
jmp    1383 <get_hbond_ij+0x243>
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x8(%rax),%rax
mov    %rax,-0x78(%rbp)
mov    -0x78(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x10(%rbp),%rdx
mov    (%rcx,%rdx,8),%rcx
cmp    0x10(%rcx),%rax
jg     137e <get_hbond_ij+0x23e>
mov    -0x28(%rbp),%rax
mov    -0x78(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
lea    0xdc1(%rip),%rdi
movsbl 0x1(%rax),%esi
call   1050 <strchr@plt>
cmp    $0x0,%rax
jne    1278 <get_hbond_ij+0x138>
mov    -0x28(%rbp),%rax
mov    -0x78(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
lea    0xd9b(%rip),%rdi
movsbl 0x3(%rax),%esi
call   1050 <strchr@plt>
cmp    $0x0,%rax
je     127d <get_hbond_ij+0x13d>
jmp    136d <get_hbond_ij+0x22d>
movq   $0x1,-0x68(%rbp)
cmpq   $0x3,-0x68(%rbp)
jg     12d5 <get_hbond_ij+0x195>
mov    0x10(%rbp),%rax
mov    -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
mov    -0x78(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x68(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x68(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    1285 <get_hbond_ij+0x145>
lea    -0x60(%rbp),%rdi
call   1060 <veclen@plt>
movaps %xmm0,%xmm1
movsd  %xmm1,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    1368 <get_hbond_ij+0x228>
mov    -0x80(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x80(%rbp)
cmp    $0x200,%rax
jle    131c <get_hbond_ij+0x1dc>
lea    0xceb(%rip),%rdi
call   1040 <nrerror@plt>
mov    0x20(%rbp),%rax
mov    -0x80(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x28(%rbp),%rax
mov    -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
call   1030 <strcpy@plt>
mov    0x28(%rbp),%rax
mov    -0x80(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x28(%rbp),%rax
mov    -0x78(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
call   1030 <strcpy@plt>
movsd  -0x38(%rbp),%xmm0
mov    0x30(%rbp),%rax
mov    -0x80(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    136d <get_hbond_ij+0x22d>
mov    -0x78(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x78(%rbp)
jmp    1212 <get_hbond_ij+0xd2>
jmp    1383 <get_hbond_ij+0x243>
mov    -0x70(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
jmp    1193 <get_hbond_ij+0x53>
mov    -0x80(%rbp),%rcx
mov    0x18(%rbp),%rax
mov    %rcx,(%rax)
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_hbond_ij>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rdx,0x20(%rsp)
mov    %rdi,0x28(%rsp)
mov    (%rcx,%rdi,8),%rax
mov    0x8(%rax),%r15
cmp    0x10(%rax),%r15
jle    118c <get_hbond_ij+0x4c>
xor    %eax,%eax
mov    %rax,0x8(%rsp)
mov    0x98(%rsp),%rax
mov    0x8(%rsp),%rcx
mov    %rcx,(%rax)
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r8,%rbx
mov    %rcx,%r14
mov    %rsi,%r12
xor    %eax,%eax
mov    %rax,0x8(%rsp)
lea    0xe5d(%rip),%rbp
jmp    11cd <get_hbond_ij+0x8d>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0x1(%r15),%rax
mov    0x28(%rsp),%rcx
mov    (%r14,%rcx,8),%rcx
cmp    0x10(%rcx),%r15
mov    %rax,%r15
lea    0xe35(%rip),%rbp
jge    116d <get_hbond_ij+0x2d>
mov    (%rbx,%r15,8),%rax
movsbl 0x1(%rax),%esi
mov    %rbp,%rdi
call   1050 <strchr@plt>
test   %rax,%rax
jne    11b0 <get_hbond_ij+0x70>
mov    (%rbx,%r15,8),%rax
movsbl 0x3(%rax),%esi
mov    %rbp,%rdi
call   1050 <strchr@plt>
test   %rax,%rax
jne    11b0 <get_hbond_ij+0x70>
mov    (%r14,%r12,8),%rax
mov    0x8(%rax),%r13
cmp    0x10(%rax),%r13
jg     11b0 <get_hbond_ij+0x70>
add    $0xffffffffffffffff,%r13
jmp    1282 <get_hbond_ij+0x142>
mov    0x8(%rsp),%rbp
lea    0x1(%rbp),%rax
mov    %rax,0x10(%rsp)
mov    0xa0(%rsp),%rax
mov    0x8(%rax,%rbp,8),%rdi
mov    (%rbx,%r15,8),%rsi
call   1030 <strcpy@plt>
mov    0xa8(%rsp),%rax
mov    0x8(%rax,%rbp,8),%rdi
mov    0x8(%rbx,%r13,8),%rsi
call   1030 <strcpy@plt>
mov    0xb0(%rsp),%rax
movsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x8(%rax,%rbp,8)
mov    0x10(%rsp),%rax
mov    %rax,0x8(%rsp)
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    (%r14,%r12,8),%rax
add    $0x1,%r13
cmp    0x10(%rax),%r13
jge    11b0 <get_hbond_ij+0x70>
mov    0x8(%rbx,%r13,8),%rax
movsbl 0x1(%rax),%esi
lea    0xd6e(%rip),%rbp
mov    %rbp,%rdi
call   1050 <strchr@plt>
test   %rax,%rax
jne    1270 <get_hbond_ij+0x130>
mov    0x8(%rbx,%r13,8),%rax
movsbl 0x3(%rax),%esi
mov    %rbp,%rdi
call   1050 <strchr@plt>
test   %rax,%rax
jne    1270 <get_hbond_ij+0x130>
mov    0x90(%rsp),%rcx
mov    (%rcx,%r15,8),%rax
mov    0x8(%rcx,%r13,8),%rcx
movupd 0x8(%rax),%xmm0
movupd 0x8(%rcx),%xmm1
subpd  %xmm1,%xmm0
movupd %xmm0,0x38(%rsp)
movsd  0x18(%rax),%xmm0
subsd  0x18(%rcx),%xmm0
movsd  %xmm0,0x48(%rsp)
lea    0x30(%rsp),%rdi
call   1060 <veclen@plt>
mov    0x20(%rsp),%rax
movsd  (%rax),%xmm1
ucomisd %xmm0,%xmm1
jbe    1270 <get_hbond_ij+0x130>
movsd  %xmm0,0x18(%rsp)
cmpq   $0x200,0x8(%rsp)
jl     120b <get_hbond_ij+0xcb>
lea    0xcdf(%rip),%rdi
call   1040 <nrerror@plt>
jmp    120b <get_hbond_ij+0xcb>

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

LLVM IR (-O2):
@.str = [2 x i8] c"P\00"
@.str.1 = [44 x i8] c"Too many possible H-bonds between two bases\00"
define void @get_hbond_ij(i64 %0, i64 %1, double* nocapture %2, i64** nocapture %3, i8** nocapture %4, i8* nocapture readnone %5, double** nocapture %6, i64* nocapture writeonly %7, i8** nocapture %8, i8** nocapture %9, double* nocapture writeonly %10) {
%12 = alloca [4 x double]
%13 = bitcast [4 x double]* %12 to i8*
%14 = getelementptr inbounds i64*, i64** %3, i64 %0
%15 = load i64*, i64** %14
%16 = getelementptr inbounds i64, i64* %15, i64 1
%17 = load i64, i64* %16
%18 = getelementptr inbounds i64, i64* %15, i64 2
%19 = load i64, i64* %18
%20 = icmp sgt i64 %17, %19
br i1 %20, label %114, label %21
21: ; preds = %11
%22 = getelementptr inbounds i64*, i64** %3, i64 %1
%23 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 0
%24 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 1
%25 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 0, i64 3
%26 = bitcast double* %24 to <2 x double>*
br label %27
27: ; preds = %21, %107
%28 = phi i64 [ 0, %21 ], [ %108, %107 ]
%29 = phi i64 [ %17, %21 ], [ %109, %107 ]
%30 = getelementptr inbounds i8*, i8** %4, i64 %29
%31 = load i8*, i8** %30
%32 = getelementptr inbounds i8, i8* %31, i64 1
%33 = load i8, i8* %32
%34 = call i64 @strchr(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str, i64 0, i64 0), i8 signext %33)
%35 = icmp eq i64 %34, 0
br i1 %35, label %36, label %107
36: ; preds = %27
%37 = load i8*, i8** %30
%38 = getelementptr inbounds i8, i8* %37, i64 3
%39 = load i8, i8* %38
%40 = call i64 @strchr(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str, i64 0, i64 0), i8 signext %39)
%41 = icmp eq i64 %40, 0
br i1 %41, label %42, label %107
42: ; preds = %36
%43 = load i64*, i64** %22
%44 = getelementptr inbounds i64, i64* %43, i64 1
%45 = load i64, i64* %44
%46 = getelementptr inbounds i64, i64* %43, i64 2
%47 = load i64, i64* %46
%48 = icmp sgt i64 %45, %47
br i1 %48, label %107, label %49
49: ; preds = %42
%50 = getelementptr inbounds double*, double** %6, i64 %29
br label %51
51: ; preds = %49, %100
%52 = phi i64 [ %28, %49 ], [ %101, %100 ]
%53 = phi i64 [ %45, %49 ], [ %102, %100 ]
%54 = getelementptr inbounds i8*, i8** %4, i64 %53
%55 = load i8*, i8** %54
%56 = getelementptr inbounds i8, i8* %55, i64 1
%57 = load i8, i8* %56
%58 = call i64 @strchr(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str, i64 0, i64 0), i8 signext %57)
%59 = icmp eq i64 %58, 0
br i1 %59, label %60, label %100
60: ; preds = %51
%61 = load i8*, i8** %54
%62 = getelementptr inbounds i8, i8* %61, i64 3
%63 = load i8, i8* %62
%64 = call i64 @strchr(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str, i64 0, i64 0), i8 signext %63)
%65 = icmp eq i64 %64, 0
br i1 %65, label %66, label %100
66: ; preds = %60
%67 = load double*, double** %50
%68 = getelementptr inbounds double*, double** %6, i64 %53
%69 = load double*, double** %68
%70 = getelementptr inbounds double, double* %67, i64 1
%71 = getelementptr inbounds double, double* %69, i64 1
%72 = bitcast double* %70 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = bitcast double* %71 to <2 x double>*
%75 = load <2 x double>, <2 x double>* %74
%76 = fsub <2 x double> %73, %75
store <2 x double> %76, <2 x double>* %26
%77 = getelementptr inbounds double, double* %67, i64 3
%78 = load double, double* %77
%79 = getelementptr inbounds double, double* %69, i64 3
%80 = load double, double* %79
%81 = fsub double %78, %80
store double %81, double* %25
%82 = call double @veclen(double* %23)
%83 = load double, double* %2
%84 = fcmp olt double %82, %83
br i1 %84, label %85, label %100
85: ; preds = %66
%86 = add nsw i64 %52, 1
%87 = icmp sgt i64 %52, 511
br i1 %87, label %88, label %90
88: ; preds = %85
%89 = call i32 @nrerror(i8* getelementptr inbounds ([44 x i8], [44 x i8]* @.str.1, i64 0, i64 0))
br label %90
90: ; preds = %88, %85
%91 = getelementptr inbounds i8*, i8** %8, i64 %86
%92 = load i8*, i8** %91
%93 = load i8*, i8** %30
%94 = call i32 @strcpy(i8* %92, i8* %93)
%95 = getelementptr inbounds i8*, i8** %9, i64 %86
%96 = load i8*, i8** %95
%97 = load i8*, i8** %54
%98 = call i32 @strcpy(i8* %96, i8* %97)
%99 = getelementptr inbounds double, double* %10, i64 %86
store double %82, double* %99
br label %100
100: ; preds = %66, %90, %51, %60
%101 = phi i64 [ %52, %51 ], [ %52, %60 ], [ %86, %90 ], [ %52, %66 ]
%102 = add nsw i64 %53, 1
%103 = load i64*, i64** %22
%104 = getelementptr inbounds i64, i64* %103, i64 2
%105 = load i64, i64* %104
%106 = icmp slt i64 %53, %105
br i1 %106, label %51, label %107
107: ; preds = %100, %42, %27, %36
%108 = phi i64 [ %28, %27 ], [ %28, %36 ], [ %28, %42 ], [ %101, %100 ]
%109 = add nsw i64 %29, 1
%110 = load i64*, i64** %14
%111 = getelementptr inbounds i64, i64* %110, i64 2
%112 = load i64, i64* %111
%113 = icmp slt i64 %29, %112
br i1 %113, label %27, label %114
114: ; preds = %107, %11
%115 = phi i64 [ 0, %11 ], [ %108, %107 ]
store i64 %115, i64* %7
ret void
}