#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 */
float dotVectors(float const *, float const *);
void computeChainRuleFromProjOntoTri(const float *v0, const float *point,
                                     const float *norm, float *vec_x,
                                     float *vec_y, float *vec_z) {
  float C1, C2;
  C1 = dotVectors(v0, norm) / dotVectors(point, norm);
  C2 = C1 / dotVectors(point, norm);
  vec_x[0] = C1 - C2 * point[0] * norm[0];
  vec_x[1] = 0.0 - C2 * point[1] * norm[0];
  vec_x[2] = 0.0 - C2 * point[2] * norm[0];
  vec_y[0] = 0.0 - C2 * point[0] * norm[1];
  vec_y[1] = C1 - C2 * point[1] * norm[1];
  vec_y[2] = 0.0 - C2 * point[2] * norm[1];
  vec_z[0] = 0.0 - C2 * point[0] * norm[2];
  vec_z[1] = 0.0 - C2 * point[1] * norm[2];
  vec_z[2] = C1 - C2 * point[2] * norm[2];
}


x86 assembly (clang, -O0):
<computeChainRuleFromProjOntoTri>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   1030 <dotVectors@plt>
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   1030 <dotVectors@plt>
movaps %xmm0,%xmm1
movss  -0x40(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
call   1030 <dotVectors@plt>
movaps %xmm0,%xmm1
movss  -0x3c(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x34(%rbp),%xmm1
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  0x4(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x8(%rax)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x34(%rbp),%xmm1
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  0x4(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,0x8(%rax)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  0x4(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x34(%rbp),%xmm1
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mulss  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,0x8(%rax)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<computeChainRuleFromProjOntoTri>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %r9,%r14
mov    %r8,%r12
mov    %rcx,%r13
mov    %rdx,%rbx
mov    %rsi,%r15
mov    %rdx,%rsi
call   1030 <dotVectors@plt>
movss  %xmm0,0xc(%rsp)
mov    %r15,%rdi
mov    %rbx,%rsi
call   1030 <dotVectors@plt>
movss  0xc(%rsp),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,0xc(%rsp)
mov    %r15,%rdi
mov    %rbx,%rsi
call   1030 <dotVectors@plt>
movss  0xc(%rsp),%xmm4
movaps %xmm4,%xmm1
divss  %xmm0,%xmm1
movss  (%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  (%rbx),%xmm0
movaps %xmm4,%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x0(%r13)
movss  0x4(%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  (%rbx),%xmm0
xorps  %xmm2,%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x4(%r13)
movss  0x8(%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  (%rbx),%xmm0
xorps  %xmm2,%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x8(%r13)
movss  (%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  0x4(%rbx),%xmm0
xorps  %xmm2,%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,(%r12)
movss  0x4(%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  0x4(%rbx),%xmm0
movaps %xmm4,%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x4(%r12)
movss  0x8(%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  0x4(%rbx),%xmm0
xorps  %xmm2,%xmm2
subss  %xmm0,%xmm2
movss  %xmm2,0x8(%r12)
movss  (%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  0x8(%rbx),%xmm0
xorps  %xmm2,%xmm2
xorps  %xmm3,%xmm3
subss  %xmm0,%xmm3
movss  %xmm3,(%r14)
movss  0x4(%r15),%xmm0
mulss  %xmm1,%xmm0
mulss  0x8(%rbx),%xmm0
subss  %xmm0,%xmm2
movss  %xmm2,0x4(%r14)
mulss  0x8(%r15),%xmm1
mulss  0x8(%rbx),%xmm1
subss  %xmm1,%xmm4
movss  %xmm4,0x8(%r14)
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @computeChainRuleFromProjOntoTri(float* %0, float* %1, float* %2, float* nocapture writeonly %3, float* nocapture writeonly %4, float* nocapture writeonly %5) {
%7 = tail call float @dotVectors(float* %0, float* %2)
%8 = tail call float @dotVectors(float* %1, float* %2)
%9 = fdiv float %7, %8
%10 = tail call float @dotVectors(float* %1, float* %2)
%11 = fdiv float %9, %10
%12 = load float, float* %1
%13 = load float, float* %2
%14 = fneg float %11
%15 = fmul float %12, %14
%16 = tail call float @llvm.fmuladd.f32(float %15, float %13, float %9)
store float %16, float* %3
%17 = getelementptr inbounds float, float* %1, i64 1
%18 = load float, float* %17
%19 = fmul float %11, %18
%20 = load float, float* %2
%21 = fmul float %19, %20
%22 = fsub float 0.000000e+00, %21
%23 = getelementptr inbounds float, float* %3, i64 1
store float %22, float* %23
%24 = getelementptr inbounds float, float* %1, i64 2
%25 = load float, float* %24
%26 = fmul float %11, %25
%27 = load float, float* %2
%28 = fmul float %26, %27
%29 = fsub float 0.000000e+00, %28
%30 = getelementptr inbounds float, float* %3, i64 2
store float %29, float* %30
%31 = load float, float* %1
%32 = fmul float %11, %31
%33 = getelementptr inbounds float, float* %2, i64 1
%34 = load float, float* %33
%35 = fmul float %32, %34
%36 = fsub float 0.000000e+00, %35
store float %36, float* %4
%37 = load float, float* %17
%38 = load float, float* %33
%39 = fmul float %37, %14
%40 = tail call float @llvm.fmuladd.f32(float %39, float %38, float %9)
%41 = getelementptr inbounds float, float* %4, i64 1
store float %40, float* %41
%42 = load float, float* %24
%43 = fmul float %11, %42
%44 = load float, float* %33
%45 = fmul float %43, %44
%46 = fsub float 0.000000e+00, %45
%47 = getelementptr inbounds float, float* %4, i64 2
store float %46, float* %47
%48 = load float, float* %1
%49 = fmul float %11, %48
%50 = getelementptr inbounds float, float* %2, i64 2
%51 = load float, float* %50
%52 = fmul float %49, %51
%53 = fsub float 0.000000e+00, %52
store float %53, float* %5
%54 = load float, float* %17
%55 = fmul float %11, %54
%56 = load float, float* %50
%57 = fmul float %55, %56
%58 = fsub float 0.000000e+00, %57
%59 = getelementptr inbounds float, float* %5, i64 1
store float %58, float* %59
%60 = load float, float* %24
%61 = load float, float* %50
%62 = fmul float %60, %14
%63 = tail call float @llvm.fmuladd.f32(float %62, float %61, float %9)
%64 = getelementptr inbounds float, float* %5, i64 2
store float %63, float* %64
ret void
}