#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 volTriangle2D(double *, double *, double *);
void Jac2D(double *volu, const double *yc, const double *m, const int n) {
  double yA[2], yB[2], yC[2], yD[2], yM[2], vol;
  int i1, i2, r, id, idn, comp;
  const int m1 = m[0], m2 = m[1];
  const int nn = (m1 + 1) * (m2 + 1);
  const int inc_x = 1;
  const int inc_y = m1 + 1;
  const double *ptr;
#pragma omp parallel for default(shared) private(                              \
    ptr, i1, i2, r, id, idn, comp, vol, yA, yB, yC, yD, yM)
  for (i1 = 0; i1 < m1; i1++) {
    for (i2 = 0; i2 < m2; i2++) {
      comp = 0;
      idn = i1 + i2 * inc_y;
      for (r = 0; r < 2; r++) {
        ptr = yc + idn + comp;
        yA[r] = *ptr;
        ptr = ptr + inc_x;
        yB[r] = *ptr;
        ptr = ptr + inc_y;
        yD[r] = *ptr;
        ptr = ptr - inc_x;
        yC[r] = *ptr;
        yM[r] = 0.25 * (yA[r] + yB[r] + yC[r] + yD[r]);
        comp += nn;
      }
      vol = volTriangle2D(yA, yB, yM);
      vol += volTriangle2D(yB, yD, yM);
      vol += volTriangle2D(yD, yC, yM);
      vol += volTriangle2D(yC, yA, yM);
      id = i1 + i2 * m1;
      volu[id] = vol;
    }
  }
  return;
}


x86 assembly (clang, -O0):
<Jac2D>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x18(%rbp),%rax
cvttsd2si (%rax),%eax
mov    %eax,-0x94(%rbp)
mov    -0x18(%rbp),%rax
cvttsd2si 0x8(%rax),%eax
mov    %eax,-0x98(%rbp)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    -0x98(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x9c(%rbp)
movl   $0x1,-0xa0(%rbp)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
movl   $0x0,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x94(%rbp),%eax
jge    13c9 <Jac2D+0x2b9>
movl   $0x0,-0x80(%rbp)
mov    -0x80(%rbp),%eax
cmp    -0x98(%rbp),%eax
jge    13b6 <Jac2D+0x2a6>
movl   $0x0,-0x90(%rbp)
mov    -0x7c(%rbp),%eax
mov    -0x80(%rbp),%ecx
imul   -0xa4(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x8c(%rbp)
movl   $0x0,-0x84(%rbp)
cmpl   $0x2,-0x84(%rbp)
jge    1317 <Jac2D+0x207>
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x90(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0xb0(%rbp)
mov    -0xb0(%rbp),%rax
movsd  (%rax),%xmm0
movslq -0x84(%rbp),%rax
movsd  %xmm0,-0x30(%rbp,%rax,8)
mov    -0xb0(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0xb0(%rbp)
mov    -0xb0(%rbp),%rax
movsd  (%rax),%xmm0
movslq -0x84(%rbp),%rax
movsd  %xmm0,-0x40(%rbp,%rax,8)
mov    -0xb0(%rbp),%rax
movslq -0xa4(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0xb0(%rbp)
mov    -0xb0(%rbp),%rax
movsd  (%rax),%xmm0
movslq -0x84(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0xb0(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0xb0(%rbp)
mov    -0xb0(%rbp),%rax
movsd  (%rax),%xmm0
movslq -0x84(%rbp),%rax
movsd  %xmm0,-0x50(%rbp,%rax,8)
movslq -0x84(%rbp),%rax
movsd  -0x30(%rbp,%rax,8),%xmm1
movslq -0x84(%rbp),%rax
addsd  -0x40(%rbp,%rax,8),%xmm1
movslq -0x84(%rbp),%rax
addsd  -0x50(%rbp,%rax,8),%xmm1
movslq -0x84(%rbp),%rax
addsd  -0x60(%rbp,%rax,8),%xmm1
movsd  0xd20(%rip),%xmm0
mulsd  %xmm1,%xmm0
movslq -0x84(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
mov    -0x9c(%rbp),%eax
add    -0x90(%rbp),%eax
mov    %eax,-0x90(%rbp)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    11d0 <Jac2D+0xc0>
lea    -0x30(%rbp),%rdi
lea    -0x40(%rbp),%rsi
lea    -0x70(%rbp),%rdx
call   1030 <volTriangle2D@plt>
movsd  %xmm0,-0x78(%rbp)
lea    -0x40(%rbp),%rdi
lea    -0x60(%rbp),%rsi
lea    -0x70(%rbp),%rdx
call   1030 <volTriangle2D@plt>
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
lea    -0x60(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x70(%rbp),%rdx
call   1030 <volTriangle2D@plt>
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
lea    -0x50(%rbp),%rdi
lea    -0x30(%rbp),%rsi
lea    -0x70(%rbp),%rdx
call   1030 <volTriangle2D@plt>
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x7c(%rbp),%eax
mov    -0x80(%rbp),%ecx
imul   -0x94(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x88(%rbp)
movsd  -0x78(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x88(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    1198 <Jac2D+0x88>
jmp    13bb <Jac2D+0x2ab>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    1182 <Jac2D+0x72>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Jac2D>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
mov    %rdi,0x68(%rsp)
cvttsd2si (%rdx),%eax
test   %eax,%eax
jle    137e <Jac2D+0x26e>
cvttsd2si 0x8(%rdx),%edi
lea    0x1(%rax),%ecx
lea    0x1(%rdi),%edx
imul   %ecx,%edx
mov    %ecx,%ebp
movslq %edx,%rbx
mov    %eax,%eax
mov    %rdi,0x78(%rsp)
mov    %edi,%edi
mov    %rdi,0x60(%rsp)
mov    %rax,0x70(%rsp)
lea    0x0(,%rax,8),%r8
lea    0x8(%rsi),%rax
mov    %rax,0x58(%rsp)
lea    0x0(,%rbp,8),%r9
lea    (%rsi,%rbp,8),%rax
add    $0x8,%rax
mov    %rax,0x50(%rsp)
add    %rbx,%rbp
lea    (%rsi,%rbp,8),%rax
add    $0x8,%rax
mov    %rax,0x48(%rsp)
lea    (%rsi,%rbx,8),%rax
add    $0x8,%rax
mov    %rax,0x40(%rsp)
xor    %esi,%esi
movapd 0xe5a(%rip),%xmm8
lea    0x10(%rsp),%rdi
lea    0xb0(%rsp),%r12
lea    0xc0(%rsp),%r15
xor    %eax,%eax
mov    %r8,0x90(%rsp)
mov    %r9,0x88(%rsp)
jmp    11eb <Jac2D+0xdb>
nop
mov    0x80(%rsp),%rax
add    $0x1,%rax
add    $0x8,%rsi
cmp    0x70(%rsp),%rax
je     137e <Jac2D+0x26e>
mov    %rax,0x80(%rsp)
cmpl   $0x0,0x78(%rsp)
jle    11d0 <Jac2D+0xc0>
mov    0x60(%rsp),%rax
mov    0x40(%rsp),%r13
mov    0x48(%rsp),%rbp
mov    0x50(%rsp),%rbx
mov    0x58(%rsp),%r14
mov    0x68(%rsp),%rcx
mov    %rsi,0x98(%rsp)
mov    %rcx,0xa0(%rsp)
mov    %rax,0xa8(%rsp)
movsd  -0x8(%r14,%rsi,1),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  (%r14,%rsi,1),%xmm1
movsd  %xmm1,0xb0(%rsp)
movsd  (%rbx,%rsi,1),%xmm2
movsd  %xmm2,0x20(%rsp)
movsd  -0x8(%rbx,%rsi,1),%xmm3
movsd  %xmm3,0x30(%rsp)
movsd  -0x8(%r13,%rsi,1),%xmm4
movsd  %xmm4,0x18(%rsp)
movsd  0x0(%r13,%rsi,1),%xmm5
movsd  %xmm5,0xb8(%rsp)
movsd  0x0(%rbp,%rsi,1),%xmm6
movsd  %xmm6,0x28(%rsp)
movsd  -0x8(%rbp,%rsi,1),%xmm7
movsd  %xmm7,0x38(%rsp)
unpcklpd %xmm4,%xmm0
unpcklpd %xmm5,%xmm1
addpd  %xmm0,%xmm1
unpcklpd %xmm7,%xmm3
addpd  %xmm1,%xmm3
unpcklpd %xmm6,%xmm2
addpd  %xmm3,%xmm2
mulpd  %xmm8,%xmm2
movapd %xmm2,0xc0(%rsp)
mov    %r12,%rsi
mov    %r15,%rdx
call   1030 <volTriangle2D@plt>
movsd  %xmm0,0x8(%rsp)
mov    %r12,%rdi
lea    0x20(%rsp),%rsi
mov    %r15,%rdx
call   1030 <volTriangle2D@plt>
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0x20(%rsp),%rdi
lea    0x30(%rsp),%rsi
mov    %r15,%rdx
call   1030 <volTriangle2D@plt>
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0x30(%rsp),%rdi
lea    0x10(%rsp),%rsi
mov    %r15,%rdx
call   1030 <volTriangle2D@plt>
mov    0xa0(%rsp),%rcx
mov    0xa8(%rsp),%rax
mov    0x98(%rsp),%rsi
lea    0x10(%rsp),%rdi
movapd 0xcbb(%rip),%xmm8
mov    0x88(%rsp),%r9
mov    0x90(%rsp),%r8
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%rcx,%rsi,1)
add    %r8,%rcx
add    %r9,%r14
add    %r9,%rbx
add    %r9,%rbp
add    %r9,%r13
add    $0xffffffffffffffff,%rax
jne    1220 <Jac2D+0x110>
jmp    11d0 <Jac2D+0xc0>
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @Jac2D(double* nocapture writeonly %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = alloca [2 x double]
%6 = alloca [2 x double]
%7 = alloca [2 x double]
%8 = alloca [2 x double]
%9 = alloca <2 x double>
%10 = bitcast [2 x double]* %5 to i8*
%11 = bitcast [2 x double]* %6 to i8*
%12 = bitcast [2 x double]* %7 to i8*
%13 = bitcast [2 x double]* %8 to i8*
%14 = bitcast <2 x double>* %9 to i8*
%15 = load double, double* %2
%16 = fptosi double %15 to i32
%17 = getelementptr inbounds double, double* %2, i64 1
%18 = load double, double* %17
%19 = fptosi double %18 to i32
%20 = add nsw i32 %16, 1
%21 = icmp sgt i32 %16, 0
br i1 %21, label %22, label %89
22: ; preds = %4
%23 = add nsw i32 %19, 1
%24 = mul i32 %23, %20
%25 = icmp sgt i32 %19, 0
%26 = zext i32 %20 to i64
%27 = getelementptr inbounds [2 x double], [2 x double]* %5, i64 0, i64 0
%28 = getelementptr inbounds [2 x double], [2 x double]* %6, i64 0, i64 0
%29 = getelementptr inbounds <2 x double>, <2 x double>* %9, i64 0, i64 0
%30 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
%31 = getelementptr inbounds [2 x double], [2 x double]* %7, i64 0, i64 0
%32 = sext i32 %24 to i64
%33 = zext i32 %16 to i64
%34 = zext i32 %16 to i64
%35 = zext i32 %19 to i64
%36 = getelementptr inbounds [2 x double], [2 x double]* %5, i64 0, i64 1
%37 = getelementptr inbounds [2 x double], [2 x double]* %6, i64 0, i64 1
%38 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 1
%39 = getelementptr inbounds [2 x double], [2 x double]* %7, i64 0, i64 1
br label %40
40: ; preds = %22, %86
%41 = phi i64 [ 0, %22 ], [ %87, %86 ]
br i1 %25, label %42, label %86
42: ; preds = %40, %42
%43 = phi i64 [ %84, %42 ], [ 0, %40 ]
%44 = mul nsw i64 %43, %26
%45 = add nuw nsw i64 %44, %41
%46 = getelementptr inbounds double, double* %1, i64 %45
%47 = load double, double* %46
store double %47, double* %27
%48 = getelementptr inbounds double, double* %46, i64 1
%49 = load double, double* %48
store double %49, double* %28
%50 = getelementptr inbounds double, double* %48, i64 %26
%51 = load double, double* %50
store double %51, double* %30
%52 = getelementptr inbounds double, double* %50, i64 -1
%53 = load double, double* %52
store double %53, double* %31
%54 = getelementptr inbounds double, double* %46, i64 %32
%55 = load double, double* %54
store double %55, double* %36
%56 = getelementptr inbounds double, double* %54, i64 1
%57 = load double, double* %56
store double %57, double* %37
%58 = getelementptr inbounds double, double* %56, i64 %26
%59 = load double, double* %58
store double %59, double* %38
%60 = getelementptr inbounds double, double* %58, i64 -1
%61 = load double, double* %60
store double %61, double* %39
%62 = insertelement <2 x double> poison, double %47, i64 0
%63 = insertelement <2 x double> %62, double %55, i64 1
%64 = insertelement <2 x double> poison, double %49, i64 0
%65 = insertelement <2 x double> %64, double %57, i64 1
%66 = fadd <2 x double> %63, %65
%67 = insertelement <2 x double> poison, double %53, i64 0
%68 = insertelement <2 x double> %67, double %61, i64 1
%69 = fadd <2 x double> %66, %68
%70 = insertelement <2 x double> poison, double %51, i64 0
%71 = insertelement <2 x double> %70, double %59, i64 1
%72 = fadd <2 x double> %71, %69
%73 = fmul <2 x double> %72, <double 2.500000e-01, double 2.500000e-01>
store <2 x double> %73, <2 x double>* %9
%74 = call double @volTriangle2D(double* %27, double* %28, double* %29)
%75 = call double @volTriangle2D(double* %28, double* %30, double* %29)
%76 = fadd double %74, %75
%77 = call double @volTriangle2D(double* %30, double* %31, double* %29)
%78 = fadd double %76, %77
%79 = call double @volTriangle2D(double* %31, double* %27, double* %29)
%80 = fadd double %78, %79
%81 = mul nsw i64 %43, %33
%82 = add nuw nsw i64 %81, %41
%83 = getelementptr inbounds double, double* %0, i64 %82
store double %80, double* %83
%84 = add nuw nsw i64 %43, 1
%85 = icmp eq i64 %84, %35
br i1 %85, label %86, label %42
86: ; preds = %42, %40
%87 = add nuw nsw i64 %41, 1
%88 = icmp eq i64 %87, %34
br i1 %88, label %89, label %40
89: ; preds = %86, %4
ret void
}