#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 */
typedef int GLdouble;

/* Variables and functions */
int sqrt(int);
void matrixPerspective(int srcdim, int dstdim, GLdouble *campos,
                       GLdouble *vector) {
  int i;
  GLdouble a1 = 1, a2 = 1, b = 1;
  for (i = srcdim - 1; i >= dstdim; i--) {
    a1 *= sqrt(campos[i] * campos[i] - 1);
    a2 *= campos[i];
    b *= campos[i] - vector[i];
  }
  a1 /= b;
  a2 /= b;
  for (i = 0; i < dstdim; i++)
    vector[i] *= a1;
  vector[i] = a2;
}


x86 assembly (clang, -O0):
<matrixPerspective>:
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   $0x1,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jl     11b3 <matrixPerspective+0xa3>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,4),%edi
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   (%rax,%rcx,4),%edi
sub    $0x1,%edi
call   1030 <sqrt@plt>
imul   -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
imul   -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x18(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
sub    (%rcx,%rdx,4),%eax
imul   -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    1144 <matrixPerspective+0x34>
mov    -0x28(%rbp),%ecx
mov    -0x20(%rbp),%eax
cltd
idiv   %ecx
mov    %eax,-0x20(%rbp)
mov    -0x28(%rbp),%ecx
mov    -0x24(%rbp),%eax
cltd
idiv   %ecx
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11fe <matrixPerspective+0xee>
mov    -0x20(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   (%rax,%rcx,4),%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11d2 <matrixPerspective+0xc2>
mov    -0x24(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrixPerspective>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rcx,%rbp
mov    %esi,%r12d
mov    $0x1,%ebx
mov    $0x1,%r13d
mov    $0x1,%r14d
cmp    %esi,%edi
jle    1265 <matrixPerspective+0x155>
movslq %edi,%rsi
mov    %r12d,0x18(%rsp)
movslq %r12d,%rax
mov    %rsi,%rcx
sub    %rax,%rcx
not    %rax
add    %rsi,%rax
mov    %rcx,0x20(%rsp)
mov    %rbp,0x10(%rsp)
jne    1176 <matrixPerspective+0x66>
mov    $0x1,%r14d
mov    $0x1,%r15d
mov    $0x1,%ebx
jmp    1220 <matrixPerspective+0x110>
mov    %rcx,%rax
and    $0xfffffffffffffffe,%rax
mov    $0x1,%r14d
mov    %rsi,%r12
mov    $0x1,%r15d
mov    $0x1,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r14d,0x1c(%rsp)
mov    %rax,0x8(%rsp)
mov    -0x4(%rdx,%r12,4),%edi
imul   %edi,%edi
add    $0xffffffff,%edi
mov    %rdx,%r14
call   1030 <sqrt@plt>
mov    %eax,%ebp
imul   %ebx,%ebp
mov    -0x8(%r14,%r12,4),%edi
mov    -0x4(%r14,%r12,4),%r13d
imul   %r13d,%r15d
mov    0x10(%rsp),%rax
sub    -0x4(%rax,%r12,4),%r13d
imul   0x1c(%rsp),%r13d
lea    -0x2(%r12),%rbx
imul   %edi,%edi
add    $0xffffffff,%edi
call   1030 <sqrt@plt>
mov    %rbx,%rsi
mov    %r14,%rdx
mov    %eax,%ebx
mov    0x8(%rsp),%rax
imul   %ebp,%ebx
mov    -0x8(%r14,%r12,4),%r14d
imul   %r14d,%r15d
mov    0x10(%rsp),%rdi
sub    -0x8(%rdi,%r12,4),%r14d
imul   %r13d,%r14d
mov    %rsi,%r12
add    $0xfffffffffffffffe,%rax
jne    11a0 <matrixPerspective+0x90>
testb  $0x1,0x20(%rsp)
mov    0x10(%rsp),%rbp
mov    %r15d,%r13d
mov    0x18(%rsp),%r12d
je     1265 <matrixPerspective+0x155>
mov    -0x4(%rdx,%rsi,4),%edi
imul   %edi,%edi
add    $0xffffffff,%edi
mov    %rdx,0x8(%rsp)
mov    %rsi,%r15
call   1030 <sqrt@plt>
imul   %eax,%ebx
mov    0x8(%rsp),%rax
mov    -0x4(%rax,%r15,4),%eax
imul   %eax,%r13d
sub    -0x4(%rbp,%r15,4),%eax
imul   %eax,%r14d
mov    %ebx,%eax
cltd
idiv   %r14d
test   %r12d,%r12d
jle    1280 <matrixPerspective+0x170>
mov    %r12d,%ebx
cmp    $0x8,%r12d
jae    1287 <matrixPerspective+0x177>
xor    %edx,%edx
jmp    13d0 <matrixPerspective+0x2c0>
xor    %ebx,%ebx
jmp    13e4 <matrixPerspective+0x2d4>
mov    %ebx,%edx
and    $0xfffffff8,%edx
movd   %eax,%xmm0
pshufd $0x0,%xmm0,%xmm0
lea    -0x8(%rdx),%rsi
mov    %rsi,%r8
shr    $0x3,%r8
add    $0x1,%r8
pshufd $0xf5,%xmm0,%xmm1
test   %rsi,%rsi
je     13fe <matrixPerspective+0x2ee>
mov    %r8,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movdqu 0x0(%rbp,%rsi,4),%xmm3
movdqu 0x10(%rbp,%rsi,4),%xmm4
movdqu 0x20(%rbp,%rsi,4),%xmm5
movdqu 0x30(%rbp,%rsi,4),%xmm2
pshufd $0xf5,%xmm3,%xmm6
pmuludq %xmm0,%xmm3
pshufd $0xe8,%xmm3,%xmm3
pmuludq %xmm1,%xmm6
pshufd $0xe8,%xmm6,%xmm6
punpckldq %xmm6,%xmm3
pshufd $0xf5,%xmm4,%xmm6
pmuludq %xmm0,%xmm4
pshufd $0xe8,%xmm4,%xmm4
pmuludq %xmm1,%xmm6
pshufd $0xe8,%xmm6,%xmm6
punpckldq %xmm6,%xmm4
movdqu %xmm3,0x0(%rbp,%rsi,4)
movdqu %xmm4,0x10(%rbp,%rsi,4)
pshufd $0xf5,%xmm5,%xmm3
pmuludq %xmm0,%xmm5
pshufd $0xe8,%xmm5,%xmm4
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm3
punpckldq %xmm3,%xmm4
pshufd $0xf5,%xmm2,%xmm3
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm3
punpckldq %xmm3,%xmm2
movdqu %xmm4,0x20(%rbp,%rsi,4)
movdqu %xmm2,0x30(%rbp,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rdi
jne    12c0 <matrixPerspective+0x1b0>
test   $0x1,%r8b
je     13be <matrixPerspective+0x2ae>
movdqu 0x0(%rbp,%rsi,4),%xmm2
movdqu 0x10(%rbp,%rsi,4),%xmm3
pshufd $0xf5,%xmm2,%xmm4
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pmuludq %xmm1,%xmm4
pshufd $0xe8,%xmm4,%xmm4
punpckldq %xmm4,%xmm2
pmuludq %xmm3,%xmm0
pshufd $0xe8,%xmm0,%xmm0
pshufd $0xf5,%xmm3,%xmm3
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm1
punpckldq %xmm1,%xmm0
movdqu %xmm2,0x0(%rbp,%rsi,4)
movdqu %xmm0,0x10(%rbp,%rsi,4)
cmp    %rbx,%rdx
je     13e4 <matrixPerspective+0x2d4>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x0(%rbp,%rdx,4),%esi
imul   %eax,%esi
mov    %esi,0x0(%rbp,%rdx,4)
add    $0x1,%rdx
cmp    %rdx,%rbx
jne    13d0 <matrixPerspective+0x2c0>
mov    %r13d,%eax
cltd
idiv   %r14d
mov    %eax,0x0(%rbp,%rbx,4)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %esi,%esi
test   $0x1,%r8b
jne    1370 <matrixPerspective+0x260>
jmp    13be <matrixPerspective+0x2ae>

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

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