#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*/ crossVectors(float *, float *, float *);
float sqrt(float);
float computeAreaWGradient(const float v0[3], const float v1[3],
                           const float v2[3], float grad[3]) {
  float res[3], area;
  float v1v2[3], v1v0[3], tmp[3];
  float sc;
  v1v2[0] = v2[0] - v1[0];
  v1v2[1] = v2[1] - v1[1];
  v1v2[2] = v2[2] - v1[2];
  v1v0[0] = v0[0] - v1[0];
  v1v0[1] = v0[1] - v1[1];
  v1v0[2] = v0[2] - v1[2];
  crossVectors(v1v2, v1v0, res);
  area =
      (float)(sqrt(res[0] * res[0] + res[1] * res[1] + res[2] * res[2]) / 2.0f);
  crossVectors(v1v2, v1v0, tmp);
  if (area >= 1.0 - 0.9999999)
    sc = (float)(0.25f / area);
  else
    sc = 0.0f;
  grad[0] = (float)sc * (tmp[1] * (v2[2] - v1[2]) + tmp[2] * (v1[1] - v2[1]));
  grad[1] = (float)sc * (tmp[0] * (v1[2] - v2[2]) + tmp[2] * (v2[0] - v1[0]));
  grad[2] = (float)sc * (tmp[0] * (v2[1] - v1[1]) + tmp[1] * (v1[0] - v2[0]));
  return area;
}


x86 assembly (clang, -O0):
<computeAreaWGradient>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  0x4(%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  0x8(%rax),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  0x4(%rax),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
subss  0x8(%rax),%xmm0
movss  %xmm0,-0x40(%rbp)
lea    -0x3c(%rbp),%rdi
lea    -0x48(%rbp),%rsi
lea    -0x2c(%rbp),%rdx
call   1030 <crossVectors@plt>
movss  -0x2c(%rbp),%xmm1
movss  -0x2c(%rbp),%xmm2
movss  -0x28(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x24(%rbp),%xmm0
movss  -0x24(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
call   1040 <sqrt@plt>
movss  0xdfe(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
lea    -0x3c(%rbp),%rdi
lea    -0x48(%rbp),%rsi
lea    -0x54(%rbp),%rdx
call   1030 <crossVectors@plt>
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xdcb(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1256 <computeAreaWGradient+0x136>
movss  0xdc5(%rip),%xmm0
divss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
jmp    125e <computeAreaWGradient+0x13e>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x58(%rbp),%xmm0
movss  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
mov    -0x10(%rbp),%rax
subss  0x8(%rax),%xmm3
movss  -0x4c(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm4
mov    -0x18(%rbp),%rax
subss  0x4(%rax),%xmm4
mulss  %xmm4,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x58(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm3
mov    -0x18(%rbp),%rax
subss  0x8(%rax),%xmm3
movss  -0x4c(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm4
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm4
mulss  %xmm4,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x58(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x10(%rbp),%rax
subss  0x4(%rax),%xmm3
movss  -0x50(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm4
mov    -0x18(%rbp),%rax
subss  (%rax),%xmm4
mulss  %xmm4,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x8(%rax)
movss  -0x30(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<computeAreaWGradient>:
push   %r15
push   %r14
push   %rbx
sub    $0x40,%rsp
mov    %rcx,%r14
mov    %rdx,%r15
mov    %rsi,%rbx
movss  (%rdx),%xmm0
movss  (%rsi),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,0x1c(%rsp)
movss  0x4(%rdx),%xmm0
movss  0x4(%rsi),%xmm2
subss  %xmm2,%xmm0
movss  %xmm0,0x20(%rsp)
movss  0x8(%rdx),%xmm0
movss  0x8(%rsi),%xmm3
subss  %xmm3,%xmm0
movss  %xmm0,0x24(%rsp)
movss  (%rdi),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,0x10(%rsp)
movss  0x4(%rdi),%xmm0
subss  %xmm2,%xmm0
movss  %xmm0,0x14(%rsp)
movss  0x8(%rdi),%xmm0
subss  %xmm3,%xmm0
movss  %xmm0,0x18(%rsp)
lea    0x1c(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x34(%rsp),%rdx
call   1030 <crossVectors@plt>
movss  0x34(%rsp),%xmm1
movss  0x38(%rsp),%xmm0
mulss  %xmm0,%xmm0
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  0x3c(%rsp),%xmm0
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     11e0 <computeAreaWGradient+0xc0>
sqrtss %xmm0,%xmm0
jmp    11e5 <computeAreaWGradient+0xc5>
call   1040 <sqrt@plt>
mulss  0xe13(%rip),%xmm0
movss  %xmm0,0xc(%rsp)
lea    0x1c(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x28(%rsp),%rdx
call   1030 <crossVectors@plt>
movss  0xc(%rsp),%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
ucomisd 0xdec(%rip),%xmm1
xorps  %xmm5,%xmm5
jb     122d <computeAreaWGradient+0x10d>
movss  0xddb(%rip),%xmm5
divss  %xmm0,%xmm5
movss  0x8(%r15),%xmm1
subss  0x8(%rbx),%xmm1
movss  0x2c(%rsp),%xmm6
movss  0x30(%rsp),%xmm2
movss  0x4(%rbx),%xmm3
subss  0x4(%r15),%xmm3
mulss  %xmm2,%xmm3
mulss  %xmm6,%xmm1
addss  %xmm3,%xmm1
mulss  %xmm5,%xmm1
movss  %xmm1,(%r14)
movss  0x8(%rbx),%xmm1
subss  0x8(%r15),%xmm1
movss  0x28(%rsp),%xmm3
movss  (%r15),%xmm4
subss  (%rbx),%xmm4
mulss  %xmm2,%xmm4
mulss  %xmm3,%xmm1
addss  %xmm4,%xmm1
mulss  %xmm5,%xmm1
movss  %xmm1,0x4(%r14)
movss  0x4(%r15),%xmm1
subss  0x4(%rbx),%xmm1
movss  (%rbx),%xmm2
subss  (%r15),%xmm2
mulss  %xmm6,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm5,%xmm1
movss  %xmm1,0x8(%r14)
add    $0x40,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define float @computeAreaWGradient(float* nocapture %0, float* nocapture %1, float* nocapture %2, float* nocapture writeonly %3) {
%5 = alloca [3 x float]
%6 = alloca [3 x float]
%7 = alloca [3 x float]
%8 = alloca [3 x float]
%9 = bitcast [3 x float]* %5 to i8*
%10 = bitcast [3 x float]* %6 to i8*
%11 = bitcast [3 x float]* %7 to i8*
%12 = bitcast [3 x float]* %8 to i8*
%13 = load float, float* %2
%14 = load float, float* %1
%15 = fsub float %13, %14
%16 = getelementptr inbounds [3 x float], [3 x float]* %6, i64 0, i64 0
store float %15, float* %16
%17 = getelementptr inbounds float, float* %2, i64 1
%18 = load float, float* %17
%19 = getelementptr inbounds float, float* %1, i64 1
%20 = load float, float* %19
%21 = fsub float %18, %20
%22 = getelementptr inbounds [3 x float], [3 x float]* %6, i64 0, i64 1
store float %21, float* %22
%23 = getelementptr inbounds float, float* %2, i64 2
%24 = load float, float* %23
%25 = getelementptr inbounds float, float* %1, i64 2
%26 = load float, float* %25
%27 = fsub float %24, %26
%28 = getelementptr inbounds [3 x float], [3 x float]* %6, i64 0, i64 2
store float %27, float* %28
%29 = load float, float* %0
%30 = fsub float %29, %14
%31 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 0, i64 0
store float %30, float* %31
%32 = getelementptr inbounds float, float* %0, i64 1
%33 = load float, float* %32
%34 = fsub float %33, %20
%35 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 0, i64 1
store float %34, float* %35
%36 = getelementptr inbounds float, float* %0, i64 2
%37 = load float, float* %36
%38 = fsub float %37, %26
%39 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 0, i64 2
store float %38, float* %39
%40 = getelementptr inbounds [3 x float], [3 x float]* %5, i64 0, i64 0
%41 = call i32 @crossVectors(float* %16, float* %31, float* %40)
%42 = load float, float* %40
%43 = getelementptr inbounds [3 x float], [3 x float]* %5, i64 0, i64 1
%44 = load float, float* %43
%45 = fmul float %44, %44
%46 = call float @llvm.fmuladd.f32(float %42, float %42, float %45)
%47 = getelementptr inbounds [3 x float], [3 x float]* %5, i64 0, i64 2
%48 = load float, float* %47
%49 = call float @llvm.fmuladd.f32(float %48, float %48, float %46)
%50 = call float @sqrt(float %49)
%51 = fmul float %50, 5.000000e-01
%52 = getelementptr inbounds [3 x float], [3 x float]* %8, i64 0, i64 0
%53 = call i32 @crossVectors(float* %16, float* %31, float* %52)
%54 = fpext float %51 to double
%55 = fcmp ult double %54, 0x3E7AD7F29A800000
%56 = fdiv float 2.500000e-01, %51
%57 = select i1 %55, float 0.000000e+00, float %56
%58 = getelementptr inbounds [3 x float], [3 x float]* %8, i64 0, i64 1
%59 = load float, float* %58
%60 = load float, float* %23
%61 = load float, float* %25
%62 = fsub float %60, %61
%63 = getelementptr inbounds [3 x float], [3 x float]* %8, i64 0, i64 2
%64 = load float, float* %63
%65 = load float, float* %19
%66 = load float, float* %17
%67 = fsub float %65, %66
%68 = fmul float %64, %67
%69 = call float @llvm.fmuladd.f32(float %59, float %62, float %68)
%70 = fmul float %57, %69
store float %70, float* %3
%71 = load float, float* %52
%72 = load float, float* %25
%73 = load float, float* %23
%74 = fsub float %72, %73
%75 = load float, float* %2
%76 = load float, float* %1
%77 = fsub float %75, %76
%78 = fmul float %64, %77
%79 = call float @llvm.fmuladd.f32(float %71, float %74, float %78)
%80 = fmul float %57, %79
%81 = getelementptr inbounds float, float* %3, i64 1
store float %80, float* %81
%82 = load float, float* %17
%83 = load float, float* %19
%84 = fsub float %82, %83
%85 = load float, float* %1
%86 = load float, float* %2
%87 = fsub float %85, %86
%88 = fmul float %59, %87
%89 = call float @llvm.fmuladd.f32(float %71, float %84, float %88)
%90 = fmul float %57, %89
%91 = getelementptr inbounds float, float* %3, i64 2
store float %90, float* %91
ret float %51
}