// File: /workspace/github_cpp_makefile_star2_O3/dannycalegari[P]circle_pack/./Trigonometry_routine.cc
// Function: ROT(double)
Matrix ROT(double angle){	// rotation in xy-plane about origin
	// returns matrix C=cos(angle), S=sin(angle)
	// (C -S  0)
	// (S  C  0)
	// (0  0  1)
	Matrix M;
	M.e.clear();
	vector<double> row;
	row.clear();
	row.push_back(cos(angle));
	row.push_back(-sin(angle));
	row.push_back(0.0);
	M.e.push_back(row);
	row.clear();
	row.push_back(sin(angle));
	row.push_back(cos(angle));
	row.push_back(0.0);
	M.e.push_back(row);
	row.clear();
	row.push_back(0.0);
	row.push_back(0.0);
	row.push_back(1.0);
	M.e.push_back(row);
	return(M);
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
subq $0x30, %rsp
movsd %xmm0, 0x28(%rsp)
movq %rdi, %rbx
xorl %r15d, %r15d
movq %r15, 0x10(%rdi)
xorps %xmm0, %xmm0
movups %xmm0, (%rdi)
xorl %esi, %esi
callq 0x986a
leaq 0x10(%rsp), %r14
movq %r15, 0x10(%r14)
xorps %xmm0, %xmm0
movaps %xmm0, (%r14)
movsd 0x28(%rsp), %xmm0
callq 0x23d0
leaq 0x8(%rsp), %rdx
movsd %xmm0, (%rdx)
movq %r14, %rdi
xorl %esi, %esi
callq 0x94a6
movsd 0x28(%rsp), %xmm0
callq 0x2060
xorps 0x7761(%rip), %xmm0     # 0xa020
movlps %xmm0, 0x8(%rsp)
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %rax
cmpq %rax, %rsi
je 0x28e2
movsd %xmm0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x28fb
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %rax
movq $0x0, 0x8(%rsp)
cmpq %rax, %rsi
je 0x291b
movq $0x0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x292a
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
leaq 0x10(%rsp), %rsi
movq %rbx, %rdi
callq 0x8dc8
movq 0x10(%rsp), %rax
cmpq %rax, 0x18(%rsp)
je 0x2948
movq %rax, 0x18(%rsp)
movsd 0x28(%rsp), %xmm0
callq 0x2060
movsd %xmm0, 0x8(%rsp)
movq 0x18(%rsp), %rsi
cmpq 0x20(%rsp), %rsi
je 0x2974
movsd %xmm0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x2983
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
movsd 0x28(%rsp), %xmm0
callq 0x23d0
movsd %xmm0, 0x8(%rsp)
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %rax
cmpq %rax, %rsi
je 0x29b2
movsd %xmm0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x29cb
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %rax
movq $0x0, 0x8(%rsp)
cmpq %rax, %rsi
je 0x29eb
movq $0x0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x29fa
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
leaq 0x10(%rsp), %rsi
movq %rbx, %rdi
callq 0x8dc8
movq 0x10(%rsp), %rax
movq 0x18(%rsp), %rsi
cmpq %rax, %rsi
je 0x2a1e
movq %rax, 0x18(%rsp)
movq %rax, %rsi
movq $0x0, 0x8(%rsp)
movq 0x20(%rsp), %rax
cmpq %rax, %rsi
je 0x2a43
movq $0x0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x2a5c
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %rax
movq $0x0, 0x8(%rsp)
cmpq %rax, %rsi
je 0x2a7c
movq $0x0, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x2a95
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %rax
movabsq $0x3ff0000000000000, %rcx # imm = 0x3FF0000000000000
movq %rcx, 0x8(%rsp)
cmpq %rax, %rsi
je 0x2ab7
movq %rcx, (%rsi)
addq $0x8, %rsi
movq %rsi, 0x18(%rsp)
jmp 0x2ac6
leaq 0x10(%rsp), %rdi
leaq 0x8(%rsp), %rdx
callq 0x94a6
leaq 0x10(%rsp), %rsi
movq %rbx, %rdi
callq 0x8dc8
movq 0x10(%rsp), %rdi
testq %rdi, %rdi
je 0x2aea
movq 0x20(%rsp), %rsi
subq %rdi, %rsi
callq 0x21f0
movq %rbx, %rax
addq $0x30, %rsp
popq %rbx
popq %r14
popq %r15
retq
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
jmp 0x2b09
movq %rax, %r14
movq 0x10(%rsp), %rdi
testq %rdi, %rdi
je 0x2b23
movq 0x20(%rsp), %rsi
subq %rdi, %rsi
callq 0x21f0
movq %rbx, %rdi
callq 0x9106
movq %r14, %rdi
callq 0x23b0