#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 LinearAD_DiracIC(int *, double *, int, double, double *, double *);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
int LinearAD_DiracIC_du(int *iwork, double *rwork, int nPoints, double T,
                        double *x, double *f) {
  int i, j, iret;
  double b[3] = {rwork[0], rwork[1], rwork[2]};
  double a = rwork[4], tStart = rwork[5];
  double x0[3] = {rwork[7], rwork[8], rwork[9]};
  double y[3], t = T + tStart;
  double *u = (double *)malloc(nPoints * sizeof(double));
  iret = LinearAD_DiracIC(iwork, rwork, nPoints, T, x, u);
  for (i = 0; i < nPoints; i++) {
    for (j = 0; j < 3; j++) {
      y[j] = x[i * 3 + j] - x0[j] - b[j] * t;
      f[i * 3 + j] = u[i] * 2.0 * y[j] / (4.0 * a * t);
    }
  }
  free(u);
  return iret;
}


x86 assembly (clang, -O0):
<LinearAD_DiracIC_du>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x38(%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x48(%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x20(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movslq -0x14(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0xc0(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%edx
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rcx
mov    -0xc0(%rbp),%r8
call   1050 <LinearAD_DiracIC@plt>
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1307 <LinearAD_DiracIC_du+0x1d7>
movl   $0x0,-0x38(%rbp)
cmpl   $0x3,-0x38(%rbp)
jge    12f4 <LinearAD_DiracIC_du+0x1c4>
mov    -0x28(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movslq -0x38(%rbp),%rax
subsd  -0x90(%rbp,%rax,8),%xmm1
ff ff
movslq -0x38(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm0
movsd  -0xb8(%rbp),%xmm2
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x38(%rbp),%rax
movsd  %xmm0,-0xb0(%rbp,%rax,8)
ff ff
mov    -0xc0(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xd60(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movslq -0x38(%rbp),%rax
mulsd  -0xb0(%rbp,%rax,8),%xmm0
ff ff
movsd  0xd3e(%rip),%xmm1
mulsd  -0x68(%rbp),%xmm1
mulsd  -0xb8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
imul   $0x3,-0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    122d <LinearAD_DiracIC_du+0xfd>
jmp    12f9 <LinearAD_DiracIC_du+0x1c9>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    121a <LinearAD_DiracIC_du+0xea>
mov    -0xc0(%rbp),%rdi
call   1030 <free@plt>
mov    -0x3c(%rbp),%eax
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<LinearAD_DiracIC_du>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %r8,%r14
mov    %rcx,%rbx
movsd  %xmm0,0x8(%rsp)
mov    %edx,%r13d
mov    %rsi,%rbp
mov    %rdi,%r12
movsd  (%rsi),%xmm1
movaps %xmm1,0x80(%rsp)
movsd  0x8(%rsi),%xmm1
movaps %xmm1,0x70(%rsp)
movsd  0x10(%rsi),%xmm1
movaps %xmm1,0x60(%rsp)
movsd  0x20(%rsi),%xmm1
movaps %xmm1,0x50(%rsp)
movsd  0x38(%rsi),%xmm1
movaps %xmm1,0x40(%rsp)
movsd  0x40(%rsi),%xmm1
movaps %xmm1,0x30(%rsp)
movsd  0x48(%rsi),%xmm1
movaps %xmm1,0x20(%rsp)
movsd  0x28(%rsi),%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
lea    0x0(,%r13,8),%edi
call   1040 <malloc@plt>
mov    %rax,%r15
mov    %r12,%rdi
mov    %rbp,%rsi
mov    %r13d,%edx
movsd  0x8(%rsp),%xmm0
mov    %rbx,%rcx
mov    %rax,%r8
call   1050 <LinearAD_DiracIC@plt>
movapd 0x10(%rsp),%xmm1
movapd 0x20(%rsp),%xmm11
movapd 0x30(%rsp),%xmm10
movapd 0x40(%rsp),%xmm15
movapd 0x50(%rsp),%xmm2
movapd 0x60(%rsp),%xmm6
movapd 0x70(%rsp),%xmm14
movapd 0x80(%rsp),%xmm5
mov    %eax,%r12d
test   %r13d,%r13d
jle    1448 <LinearAD_DiracIC_du+0x318>
mulsd  0xdde(%rip),%xmm2
mulsd  %xmm1,%xmm2
mov    %r13d,%eax
movapd 0xddf(%rip),%xmm0
xorpd  %xmm0,%xmm5
xorpd  %xmm0,%xmm14
xorpd  %xmm0,%xmm6
cmp    $0x1,%r13d
jne    124b <LinearAD_DiracIC_du+0x11b>
xor    %ecx,%ecx
jmp    13a5 <LinearAD_DiracIC_du+0x275>
lea    0x0(,%rax,8),%rcx
lea    (%rcx,%rcx,2),%rcx
lea    (%r14,%rcx,1),%rsi
add    %rbx,%rcx
cmp    %r14,%rcx
seta   %dil
cmp    %rbx,%rsi
seta   %r8b
lea    (%r15,%rax,8),%rcx
cmp    %r14,%rcx
seta   %dl
cmp    %rsi,%r15
setb   %sil
xor    %ecx,%ecx
test   %r8b,%dil
jne    13a5 <LinearAD_DiracIC_du+0x275>
and    %sil,%dl
jne    13a5 <LinearAD_DiracIC_du+0x275>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movapd %xmm15,%xmm8
unpcklpd %xmm15,%xmm8
movapd %xmm5,%xmm9
unpcklpd %xmm5,%xmm9
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
movapd %xmm2,%xmm4
unpcklpd %xmm2,%xmm4
unpcklpd %xmm10,%xmm10
movapd %xmm14,%xmm12
unpcklpd %xmm14,%xmm12
unpcklpd %xmm11,%xmm11
movapd %xmm6,%xmm13
unpcklpd %xmm6,%xmm13
xor    %edx,%edx
mulpd  %xmm0,%xmm9
mulpd  %xmm0,%xmm12
mulpd  %xmm0,%xmm13
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movupd (%rbx,%rdx,1),%xmm7
movupd 0x10(%rbx,%rdx,1),%xmm1
movupd 0x20(%rbx,%rdx,1),%xmm3
movapd %xmm1,%xmm0
movsd  %xmm7,%xmm0
subpd  %xmm8,%xmm0
addpd  %xmm9,%xmm0
shufpd $0x1,%xmm3,%xmm7
subpd  %xmm10,%xmm7
addpd  %xmm12,%xmm7
movsd  %xmm1,%xmm3
subpd  %xmm11,%xmm3
addpd  %xmm13,%xmm3
movupd (%r15,%rsi,8),%xmm1
addpd  %xmm1,%xmm1
mulpd  %xmm1,%xmm0
divpd  %xmm4,%xmm0
mulpd  %xmm1,%xmm7
divpd  %xmm4,%xmm7
mulpd  %xmm3,%xmm1
divpd  %xmm4,%xmm1
movapd %xmm0,%xmm3
unpcklpd %xmm7,%xmm3
unpckhpd %xmm1,%xmm7
shufpd $0x2,%xmm0,%xmm1
movupd %xmm1,0x10(%r14,%rdx,1)
movupd %xmm7,0x20(%r14,%rdx,1)
movupd %xmm3,(%r14,%rdx,1)
add    $0x2,%rsi
add    $0x30,%rdx
cmp    %rsi,%rcx
jne    12f0 <LinearAD_DiracIC_du+0x1c0>
cmp    %rax,%rcx
movapd 0x30(%rsp),%xmm10
movapd 0x20(%rsp),%xmm11
movapd 0x10(%rsp),%xmm1
je     1448 <LinearAD_DiracIC_du+0x318>
lea    (%rcx,%rcx,2),%rdx
lea    0x10(,%rdx,8),%rdx
mulsd  %xmm1,%xmm5
mulsd  %xmm1,%xmm14
mulsd  %xmm1,%xmm6
xchg   %ax,%ax
movsd  -0x10(%rbx,%rdx,1),%xmm0
subsd  %xmm15,%xmm0
addsd  %xmm5,%xmm0
movsd  (%r15,%rcx,8),%xmm1
addsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,-0x10(%r14,%rdx,1)
movsd  -0x8(%rbx,%rdx,1),%xmm0
subsd  %xmm10,%xmm0
addsd  %xmm14,%xmm0
movsd  (%r15,%rcx,8),%xmm1
addsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,-0x8(%r14,%rdx,1)
movsd  (%rbx,%rdx,1),%xmm0
subsd  %xmm11,%xmm0
addsd  %xmm6,%xmm0
movsd  (%r15,%rcx,8),%xmm1
addsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,(%r14,%rdx,1)
add    $0x1,%rcx
add    $0x18,%rdx
cmp    %rcx,%rax
jne    13c0 <LinearAD_DiracIC_du+0x290>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r12d,%eax
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @LinearAD_DiracIC_du(i32* %0, double* %1, i32 %2, double %3, double* %4, double* nocapture writeonly %5) {
%7 = load double, double* %1
%8 = getelementptr inbounds double, double* %1, i64 1
%9 = load double, double* %8
%10 = getelementptr inbounds double, double* %1, i64 2
%11 = load double, double* %10
%12 = getelementptr inbounds double, double* %1, i64 4
%13 = load double, double* %12
%14 = getelementptr inbounds double, double* %1, i64 5
%15 = load double, double* %14
%16 = getelementptr inbounds double, double* %1, i64 7
%17 = load double, double* %16
%18 = getelementptr inbounds double, double* %1, i64 8
%19 = load double, double* %18
%20 = getelementptr inbounds double, double* %1, i64 9
%21 = load double, double* %20
%22 = fadd double %15, %3
%23 = shl i32 %2, 3
%24 = tail call i64 @malloc(i32 %23)
%25 = inttoptr i64 %24 to double*
%26 = inttoptr i64 %24 to i8*
%27 = tail call i32 @LinearAD_DiracIC(i32* %0, double* %1, i32 %2, double %3, double* %4, double* %25)
%28 = icmp sgt i32 %2, 0
br i1 %28, label %29, label %145
29: ; preds = %6
%30 = fmul double %13, 4.000000e+00
%31 = fmul double %30, %22
%32 = zext i32 %2 to i64
%33 = fneg double %7
%34 = fneg double %9
%35 = fneg double %11
%36 = icmp eq i32 %2, 1
br i1 %36, label %108, label %37
37: ; preds = %29
%38 = mul nuw nsw i64 %32, 3
%39 = getelementptr double, double* %5, i64 %38
%40 = bitcast double* %39 to i8*
%41 = getelementptr double, double* %4, i64 %38
%42 = getelementptr double, double* %25, i64 %32
%43 = icmp ugt double* %41, %5
%44 = icmp ugt double* %39, %4
%45 = and i1 %43, %44
%46 = icmp ugt double* %42, %5
%47 = icmp ult i8* %26, %40
%48 = and i1 %46, %47
%49 = or i1 %45, %48
br i1 %49, label %108, label %50
50: ; preds = %37
%51 = and i64 %32, 4294967294
%52 = insertelement <2 x double> poison, double %17, i64 0
%53 = shufflevector <2 x double> %52, <2 x double> poison, <2 x i32> zeroinitializer
%54 = insertelement <2 x double> poison, double %33, i64 0
%55 = shufflevector <2 x double> %54, <2 x double> poison, <2 x i32> zeroinitializer
%56 = insertelement <2 x double> poison, double %22, i64 0
%57 = shufflevector <2 x double> %56, <2 x double> poison, <2 x i32> zeroinitializer
%58 = insertelement <2 x double> poison, double %31, i64 0
%59 = shufflevector <2 x double> %58, <2 x double> poison, <2 x i32> zeroinitializer
%60 = insertelement <2 x double> poison, double %19, i64 0
%61 = shufflevector <2 x double> %60, <2 x double> poison, <2 x i32> zeroinitializer
%62 = insertelement <2 x double> poison, double %34, i64 0
%63 = shufflevector <2 x double> %62, <2 x double> poison, <2 x i32> zeroinitializer
%64 = insertelement <2 x double> poison, double %21, i64 0
%65 = shufflevector <2 x double> %64, <2 x double> poison, <2 x i32> zeroinitializer
%66 = insertelement <2 x double> poison, double %35, i64 0
%67 = shufflevector <2 x double> %66, <2 x double> poison, <2 x i32> zeroinitializer
br label %68
68: ; preds = %68, %50
%69 = phi i64 [ 0, %50 ], [ %104, %68 ]
%70 = mul nuw nsw i64 %69, 3
%71 = getelementptr inbounds double, double* %25, i64 %69
%72 = getelementptr inbounds double, double* %4, i64 %70
%73 = bitcast double* %72 to <6 x double>*
%74 = load <6 x double>, <6 x double>* %73
%75 = shufflevector <6 x double> %74, <6 x double> poison, <2 x i32> <i32 0, i32 3>
%76 = shufflevector <6 x double> %74, <6 x double> poison, <2 x i32> <i32 1, i32 4>
%77 = shufflevector <6 x double> %74, <6 x double> poison, <2 x i32> <i32 2, i32 5>
%78 = fsub <2 x double> %75, %53
%79 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %55, <2 x double> %57, <2 x double> %78)
%80 = bitcast double* %71 to <2 x double>*
%81 = load <2 x double>, <2 x double>* %80
%82 = fmul <2 x double> %81, <double 2.000000e+00, double 2.000000e+00>
%83 = fmul <2 x double> %79, %82
%84 = fdiv <2 x double> %83, %59
%85 = fsub <2 x double> %76, %61
%86 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %63, <2 x double> %57, <2 x double> %85)
%87 = bitcast double* %71 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = fmul <2 x double> %88, <double 2.000000e+00, double 2.000000e+00>
%90 = fmul <2 x double> %86, %89
%91 = fdiv <2 x double> %90, %59
%92 = fsub <2 x double> %77, %65
%93 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %67, <2 x double> %57, <2 x double> %92)
%94 = bitcast double* %71 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = fmul <2 x double> %95, <double 2.000000e+00, double 2.000000e+00>
%97 = fmul <2 x double> %93, %96
%98 = fdiv <2 x double> %97, %59
%99 = getelementptr inbounds double, double* %5, i64 %70
%100 = bitcast double* %99 to <6 x double>*
%101 = shufflevector <2 x double> %84, <2 x double> %91, <4 x i32> <i32 0, i32 1, i32 2, i32 3>
%102 = shufflevector <2 x double> %98, <2 x double> poison, <4 x i32> <i32 0, i32 1, i32 undef, i32 undef>
%103 = shufflevector <4 x double> %101, <4 x double> %102, <6 x i32> <i32 0, i32 2, i32 4, i32 1, i32 3, i32 5>
store <6 x double> %103, <6 x double>* %100
%104 = add nuw i64 %69, 2
%105 = icmp eq i64 %104, %51
br i1 %105, label %106, label %68
106: ; preds = %68
%107 = icmp eq i64 %51, %32
br i1 %107, label %145, label %108
108: ; preds = %37, %29, %106
%109 = phi i64 [ 0, %37 ], [ 0, %29 ], [ %51, %106 ]
br label %110
110: ; preds = %108, %110
%111 = phi i64 [ %143, %110 ], [ %109, %108 ]
%112 = mul nuw nsw i64 %111, 3
%113 = getelementptr inbounds double, double* %25, i64 %111
%114 = getelementptr inbounds double, double* %4, i64 %112
%115 = load double, double* %114
%116 = fsub double %115, %17
%117 = tail call double @llvm.fmuladd.f64(double %33, double %22, double %116)
%118 = load double, double* %113
%119 = fmul double %118, 2.000000e+00
%120 = fmul double %117, %119
%121 = fdiv double %120, %31
%122 = getelementptr inbounds double, double* %5, i64 %112
store double %121, double* %122
%123 = add nuw nsw i64 %112, 1
%124 = getelementptr inbounds double, double* %4, i64 %123
%125 = load double, double* %124
%126 = fsub double %125, %19
%127 = tail call double @llvm.fmuladd.f64(double %34, double %22, double %126)
%128 = load double, double* %113
%129 = fmul double %128, 2.000000e+00
%130 = fmul double %127, %129
%131 = fdiv double %130, %31
%132 = getelementptr inbounds double, double* %5, i64 %123
store double %131, double* %132
%133 = add nuw nsw i64 %112, 2
%134 = getelementptr inbounds double, double* %4, i64 %133
%135 = load double, double* %134
%136 = fsub double %135, %21
%137 = tail call double @llvm.fmuladd.f64(double %35, double %22, double %136)
%138 = load double, double* %113
%139 = fmul double %138, 2.000000e+00
%140 = fmul double %137, %139
%141 = fdiv double %140, %31
%142 = getelementptr inbounds double, double* %5, i64 %133
store double %141, double* %142
%143 = add nuw nsw i64 %111, 1
%144 = icmp eq i64 %143, %32
br i1 %144, label %145, label %110
145: ; preds = %110, %106, %6
%146 = tail call i32 @free(double* %25)
ret i32 %27
}