#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*/ assert(int);
void l2projectScalar(const int nSimplices, const int nPointsPerSimplex,
                     double *dV, double *r) {
  int eN, k;
  double ravg, vol;
  for (eN = 0; eN < nSimplices; eN++) {
    ravg = 0.0;
    vol = 0.0;
    for (k = 0; k < nPointsPerSimplex; k++) {
      vol += dV[eN * nPointsPerSimplex + k];
      ravg += r[eN * nPointsPerSimplex + k] * dV[eN * nPointsPerSimplex + k];
    }
    assert(vol > 0.0);
    ravg = ravg / vol;
    for (k = 0; k < nPointsPerSimplex; k++) {
      r[eN * nPointsPerSimplex + k] = ravg;
    }
  }
}


x86 assembly (clang, -O0):
<l2projectScalar>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    123f <l2projectScalar+0x12f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11c8 <l2projectScalar+0xb8>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1150 <l2projectScalar+0x40>
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
movsd  -0x28(%rbp),%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    122c <l2projectScalar+0x11c>
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11f7 <l2projectScalar+0xe7>
jmp    1231 <l2projectScalar+0x121>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112d <l2projectScalar+0x1d>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<l2projectScalar>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %rcx,0x10(%rsp)
mov    %rdx,0x8(%rsp)
test   %edi,%edi
jle    1394 <l2projectScalar+0x284>
movslq %esi,%rax
mov    %edi,%ecx
mov    %rcx,0x58(%rsp)
mov    %eax,%r8d
mov    %r8d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0x38(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x30(%rsp)
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %r8d,%edx
and    $0xfffffffe,%edx
mov    %rdx,0x28(%rsp)
mov    %ecx,%ebx
and    $0x3,%ebx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x18(%rsp)
mov    0x10(%rsp),%r14
lea    0x8(%r14),%r12
mov    %rax,0x40(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,0x50(%rsp)
mov    0x8(%rsp),%rax
lea    0x8(%rax),%r13
lea    0x70(%r14),%rbp
lea    0x10(%r14),%rdi
mov    %rbx,0x20(%rsp)
shl    $0x5,%rbx
xor    %r15d,%r15d
mov    %r8,0x48(%rsp)
mov    %esi,0x4(%rsp)
jmp    11e3 <l2projectScalar+0xd3>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r15
mov    0x50(%rsp),%rax
add    %rax,%r12
add    %rax,%r13
add    %rax,%rbp
add    %rax,%rdi
add    %rax,%r14
cmp    0x58(%rsp),%r15
je     1394 <l2projectScalar+0x284>
mov    %rdi,0x68(%rsp)
xorpd  %xmm2,%xmm2
xorpd  %xmm3,%xmm3
test   %esi,%esi
jle    1286 <l2projectScalar+0x176>
xorpd  %xmm3,%xmm3
xor    %eax,%eax
cmp    $0x1,%esi
jne    1210 <l2projectScalar+0x100>
xorpd  %xmm2,%xmm2
jmp    1258 <l2projectScalar+0x148>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
mov    0x28(%rsp),%rcx
nopl   0x0(%rax)
movsd  -0x8(%r13,%rax,8),%xmm1
addsd  %xmm1,%xmm3
mulsd  -0x8(%r12,%rax,8),%xmm1
movsd  0x0(%r13,%rax,8),%xmm2
addsd  %xmm2,%xmm3
mulsd  (%r12,%rax,8),%xmm2
addsd  %xmm0,%xmm1
addsd  %xmm1,%xmm2
add    $0x2,%rax
movapd %xmm2,%xmm0
cmp    %rax,%rcx
jne    1220 <l2projectScalar+0x110>
test   $0x1,%r8b
je     1286 <l2projectScalar+0x176>
mov    %r15,%rcx
imul   0x40(%rsp),%rcx
add    %rcx,%rax
mov    0x8(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
addsd  %xmm0,%xmm3
mov    0x10(%rsp),%rcx
mulsd  (%rcx,%rax,8),%xmm0
addsd  %xmm0,%xmm2
movsd  %xmm3,0x60(%rsp)
movapd %xmm2,0x70(%rsp)
xor    %edi,%edi
ucomisd 0xd64(%rip),%xmm3
seta   %dil
call   1030 <assert@plt>
movapd 0x70(%rsp),%xmm1
mov    0x48(%rsp),%r8
mov    0x4(%rsp),%esi
test   %esi,%esi
mov    0x68(%rsp),%rdi
jle    11c0 <l2projectScalar+0xb0>
divsd  0x60(%rsp),%xmm1
cmp    $0x4,%esi
jae    12e0 <l2projectScalar+0x1d0>
xor    %eax,%eax
jmp    1380 <l2projectScalar+0x270>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
cmpq   $0xc,0x30(%rsp)
jae    12f4 <l2projectScalar+0x1e4>
xor    %eax,%eax
jmp    133a <l2projectScalar+0x22a>
mov    0x18(%rsp),%rcx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x70(%rbp,%rax,8)
movupd %xmm0,-0x60(%rbp,%rax,8)
movupd %xmm0,-0x50(%rbp,%rax,8)
movupd %xmm0,-0x40(%rbp,%rax,8)
movupd %xmm0,-0x30(%rbp,%rax,8)
movupd %xmm0,-0x20(%rbp,%rax,8)
movupd %xmm0,-0x10(%rbp,%rax,8)
movupd %xmm0,0x0(%rbp,%rax,8)
add    $0x10,%rax
add    $0xfffffffffffffffc,%rcx
jne    1300 <l2projectScalar+0x1f0>
cmpq   $0x0,0x20(%rsp)
je     1364 <l2projectScalar+0x254>
lea    (%rdi,%rax,8),%rax
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%rax,%rcx,1)
movupd %xmm0,(%rax,%rcx,1)
add    $0x20,%rcx
cmp    %rcx,%rbx
jne    1350 <l2projectScalar+0x240>
mov    0x38(%rsp),%rcx
mov    %rcx,%rax
cmp    %r8,%rcx
je     11c0 <l2projectScalar+0xb0>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm1,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rax,%r8
jne    1380 <l2projectScalar+0x270>
jmp    11c0 <l2projectScalar+0xb0>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @l2projectScalar(i32 %0, i32 %1, double* %2, double* %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
store i32 %0, i32* %5
store i32 %1, i32* %6
store double* %2, double** %7
store double* %3, double** %8
store i32 0, i32* %9
br label %13
13: ; preds = %83, %4
%14 = load i32, i32* %9
%15 = load i32, i32* %5
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %86
17: ; preds = %13
store double 0.000000e+00, double* %11
store double 0.000000e+00, double* %12
store i32 0, i32* %10
br label %18
18: ; preds = %54, %17
%19 = load i32, i32* %10
%20 = load i32, i32* %6
%21 = icmp slt i32 %19, %20
br i1 %21, label %22, label %57
22: ; preds = %18
%23 = load double*, double** %7
%24 = load i32, i32* %9
%25 = load i32, i32* %6
%26 = mul nsw i32 %24, %25
%27 = load i32, i32* %10
%28 = add nsw i32 %26, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %23, i64 %29
%31 = load double, double* %30
%32 = load double, double* %12
%33 = fadd double %32, %31
store double %33, double* %12
%34 = load double*, double** %8
%35 = load i32, i32* %9
%36 = load i32, i32* %6
%37 = mul nsw i32 %35, %36
%38 = load i32, i32* %10
%39 = add nsw i32 %37, %38
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %34, i64 %40
%42 = load double, double* %41
%43 = load double*, double** %7
%44 = load i32, i32* %9
%45 = load i32, i32* %6
%46 = mul nsw i32 %44, %45
%47 = load i32, i32* %10
%48 = add nsw i32 %46, %47
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %43, i64 %49
%51 = load double, double* %50
%52 = load double, double* %11
%53 = call double @llvm.fmuladd.f64(double %42, double %51, double %52)
store double %53, double* %11
br label %54
54: ; preds = %22
%55 = load i32, i32* %10
%56 = add nsw i32 %55, 1
store i32 %56, i32* %10
br label %18
57: ; preds = %18
%58 = load double, double* %12
%59 = fcmp ogt double %58, 0.000000e+00
%60 = zext i1 %59 to i32
%61 = call i32 @assert(i32 %60)
%62 = load double, double* %11
%63 = load double, double* %12
%64 = fdiv double %62, %63
store double %64, double* %11
store i32 0, i32* %10
br label %65
65: ; preds = %79, %57
%66 = load i32, i32* %10
%67 = load i32, i32* %6
%68 = icmp slt i32 %66, %67
br i1 %68, label %69, label %82
69: ; preds = %65
%70 = load double, double* %11
%71 = load double*, double** %8
%72 = load i32, i32* %9
%73 = load i32, i32* %6
%74 = mul nsw i32 %72, %73
%75 = load i32, i32* %10
%76 = add nsw i32 %74, %75
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %71, i64 %77
store double %70, double* %78
br label %79
79: ; preds = %69
%80 = load i32, i32* %10
%81 = add nsw i32 %80, 1
store i32 %81, i32* %10
br label %65
82: ; preds = %65
br label %83
83: ; preds = %82
%84 = load i32, i32* %9
%85 = add nsw i32 %84, 1
store i32 %85, i32* %9
br label %13
86: ; preds = %13
ret void
}

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