LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: normalize_power_factors
void
normalize_power_factors(int h)
{
	int i, k;
	struct atom *p1;
	k = tos;
	for (i = h; i < k; i++) {
		p1 = stack[i];
		if (car(p1) == symbol(POWER)) {
			push(cadr(p1));
			push(caddr(p1));
			power();
			p1 = pop();
			if (car(p1) == symbol(MULTIPLY)) {
				stack[i] = cadr(p1);
				p1 = cddr(p1);
				while (iscons(p1)) {
					push(car(p1));
					p1 = cdr(p1);
				}
			} else
				stack[i] = p1;
		}
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movl 0x21fab(%rip), %ebx     # 0x3a2b8
cmpl %edi, %ebx
jle 0x1843c
movslq %edi, %r14
leaq 0x2b2a1(%rip), %r15     # 0x435c0
leaq 0x286da(%rip), %r12     # 0x40a00
movl %ebx, -0x2c(%rbp)
movq (%r15,%r14,8), %r13
movb 0x10(%r13), %al
cmpb $0x1, %al
movq %r12, %rcx
cmoveq %r13, %rcx
movq (%rcx), %rcx
cmpq 0x2af6c(%rip), %rcx     # 0x432b0
jne 0x18430
movq %r15, %rbx
leaq 0x8(%r13), %r15
cmpb $0x1, %al
movq %r12, %rax
cmoveq %r15, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmovneq %r12, %rax
movq (%rax), %rdi
callq 0x8b01
cmpb $0x1, 0x10(%r13)
cmovneq %r12, %r15
movq (%r15), %rax
movb 0x10(%rax), %cl
addq $0x8, %rax
cmpb $0x1, %cl
movq %r12, %rdx
cmoveq %rax, %rdx
movq (%rdx), %rdx
cmpb $0x1, 0x10(%rdx)
movq %r12, %rdx
jne 0x183a0
cmpb $0x1, %cl
cmovneq %r12, %rax
movq (%rax), %rdx
movq (%rdx), %rdi
callq 0x8b01
callq 0xa234
callq 0x8b78
movb 0x10(%rax), %cl
cmpb $0x1, %cl
movq %r12, %rdx
cmoveq %rax, %rdx
movq (%rdx), %rdx
cmpq 0x2aedf(%rip), %rdx     # 0x432a8
movq %rbx, %r15
je 0x183d7
movq %rax, (%r15,%r14,8)
movl -0x2c(%rbp), %ebx
jmp 0x18430
addq $0x8, %rax
cmpb $0x1, %cl
cmovneq %r12, %rax
movq (%rax), %rcx
cmpb $0x1, 0x10(%rcx)
cmovneq %r12, %rcx
movq (%rcx), %rcx
movq %rcx, (%r15,%r14,8)
movq (%rax), %rax
leaq 0x8(%rax), %rcx
cmpb $0x1, 0x10(%rax)
cmovneq %r12, %rcx
movq (%rcx), %r13
cmpb $0x1, 0x10(%r13)
movl -0x2c(%rbp), %ebx
jne 0x18430
movq (%r13), %rdi
callq 0x8b01
leaq 0x8(%r13), %rax
cmpb $0x1, 0x10(%r13)
cmovneq %r12, %rax
movq (%rax), %r13
cmpb $0x1, 0x10(%r13)
je 0x18410
incq %r14
cmpl %r14d, %ebx
jne 0x18329
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: reduce_radical_factors
struct atom *
reduce_radical_factors(int h, struct atom *COEF)
{
	if (!any_radical_factors(h))
		return COEF;

	if (isrational(COEF))
		return reduce_radical_rational(h, COEF);
	else
		return reduce_radical_double(h, COEF);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r14
pushq %rbx
movq %rsi, %rbx
movl %edi, %r14d
callq 0x18b11
testl %eax, %eax
je 0x18476
cmpb $0x4, 0x10(%rbx)
jne 0x1847e
movl %r14d, %edi
movq %rbx, %rsi
popq %rbx
popq %r14
popq %rbp
jmp 0x18b62
movq %rbx, %rax
popq %rbx
popq %r14
popq %rbp
retq
movl %r14d, %edi
movq %rbx, %rsi
popq %rbx
popq %r14
popq %rbp
jmp 0x18e7e
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: expand_sum_factors
void
expand_sum_factors(int h)
{
	int i, n;
	struct atom *p1, *p2;

	if (tos - h < 2)
		return;

	// search for a sum factor

	for (i = h; i < tos; i++) {
		p2 = stack[i];
		if (car(p2) == symbol(ADD))
			break;
	}

	if (i == tos)
		return; // no sum factors

	// remove the sum factor

	slice(i, 1);

	n = tos - h;

	if (n > 1) {
		sort_factors(n);
		list(n);
		push_symbol(MULTIPLY);
		swap();
		cons(); // prepend MULTIPLY to list
	}

	p1 = pop(); // p1 is the multiplier

	p2 = cdr(p2); // p2 is the sum

	while (iscons(p2)) {
		push(p1);
		push(car(p2));
		multiply();
		p2 = cdr(p2);
	}

	add_terms(tos - h);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movl 0x21e1a(%rip), %eax     # 0x3a2b8
movl %eax, %ecx
subl %edi, %ecx
cmpl $0x2, %ecx
jl 0x184fc
movl %edi, %ebx
cmpl %edi, %eax
jle 0x184f8
movslq %eax, %rcx
movq 0x2ade9(%rip), %rdx     # 0x432a0
movslq %ebx, %rsi
subq %rsi, %rcx
shlq $0x3, %rsi
xorl %edi, %edi
leaq 0x2b0f6(%rip), %r8      # 0x435c0
leaq 0x2852f(%rip), %r9      # 0x40a00
leaq (%rsi,%rdi,8), %r10
movq (%r8,%r10), %r15
cmpb $0x1, 0x10(%r15)
movq %r9, %r10
cmoveq %r15, %r10
cmpq %rdx, (%r10)
je 0x184f6
incq %rdi
cmpq %rdi, %rcx
jne 0x184d1
movl %eax, %edi
jmp 0x184f8
addl %ebx, %edi
cmpl %eax, %edi
jne 0x18505
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
movl $0x1, %esi
callq 0xadfc
movslq 0x21da2(%rip), %rax     # 0x3a2b8
movl %eax, %r14d
subl %ebx, %r14d
cmpl $0x2, %r14d
jl 0x1856a
leaq 0x2b097(%rip), %rcx     # 0x435c0
leaq (%rcx,%rax,8), %rdi
movl %r14d, %esi
leaq (,%rsi,8), %rax
subq %rax, %rdi
leaq 0x368(%rip), %rcx       # 0x188aa
movl $0x8, %edx
callq 0x70a0
movl %r14d, %edi
callq 0x8ad9
movq 0x2ad4d(%rip), %rdi     # 0x432a8
callq 0x8b01
callq 0xad38
callq 0x8b47
callq 0x8b78
cmpb $0x1, 0x10(%r15)
leaq 0x8(%r15), %rcx
leaq 0x28481(%rip), %r15     # 0x40a00
cmovneq %r15, %rcx
movq (%rcx), %r12
cmpb $0x1, 0x10(%r12)
jne 0x185d2
movq %rax, %r14
movq %r14, %rdi
callq 0x8b01
cmpb $0x1, 0x10(%r12)
movq %r15, %rax
cmoveq %r12, %rax
movq (%rax), %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
leaq 0x8(%r12), %rax
cmpb $0x1, 0x10(%r12)
cmovneq %r15, %rax
movq (%rax), %r12
cmpb $0x1, 0x10(%r12)
je 0x18591
movl 0x21ce0(%rip), %edi     # 0x3a2b8
subl %ebx, %edi
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
jmp 0xa975
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: multiply_numbers
void
multiply_numbers(struct atom *p1, struct atom *p2)
{
	double d1, d2;

	if (isrational(p1) && isrational(p2)) {
		multiply_rationals(p1, p2);
		return;
	}

	push(p1);
	d1 = pop_double();

	push(p2);
	d2 = pop_double();

	push_double(d1 * d2);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
pushq %rax
movq %rsi, %rbx
cmpb $0x4, 0x10(%rdi)
jne 0x18643
cmpb $0x4, 0x10(%rbx)
jne 0x18643
movq %rbx, %rsi
addq $0x8, %rsp
popq %rbx
popq %rbp
jmp 0x189bc
callq 0x8b01
callq 0x9086
movsd %xmm0, -0x10(%rbp)
movq %rbx, %rdi
callq 0x8b01
callq 0x9086
mulsd -0x10(%rbp), %xmm0
movsd %xmm0, -0x10(%rbp)
callq 0x74ec
movb $0x5, 0x10(%rax)
movsd -0x10(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
addq $0x8, %rsp
popq %rbx
popq %rbp
jmp 0x8b01
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: combine_factors_nib
int
combine_factors_nib(int i, int j)
{
	struct atom *p1, *p2, *BASE1, *EXPO1, *BASE2, *EXPO2;

	p1 = stack[i];
	p2 = stack[j];

	if (car(p1) == symbol(POWER)) {
		BASE1 = cadr(p1);
		EXPO1 = caddr(p1);
	} else {
		BASE1 = p1;
		EXPO1 = one;
	}

	if (car(p2) == symbol(POWER)) {
		BASE2 = cadr(p2);
		EXPO2 = caddr(p2);
	} else {
		BASE2 = p2;
		EXPO2 = one;
	}

	if (!equal(BASE1, BASE2))
		return 0;

	if (isdouble(BASE2))
		BASE1 = BASE2; // if mixed rational and double, use double

	push_symbol(POWER);
	push(BASE1);
	push(EXPO1);
	push(EXPO2);
	add();
	list(3);

	stack[i] = pop();

	return 1;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movslq %edi, %r13
leaq 0x2aef1(%rip), %rcx     # 0x435c0
movq (%rcx,%r13,8), %rbx
movslq %esi, %rax
movq (%rcx,%rax,8), %r14
movb 0x10(%rbx), %cl
cmpb $0x1, %cl
leaq 0x28319(%rip), %rax     # 0x40a00
movq %rax, %rdi
cmoveq %rbx, %rdi
movq 0x2abbb(%rip), %rsi     # 0x432b0
leaq 0xee3dc(%rip), %rdx     # 0x106ad8
cmpq %rsi, (%rdi)
je 0x18706
movq %rdx, %rcx
jmp 0x18751
addq $0x8, %rbx
cmpb $0x1, %cl
leaq 0x282ec(%rip), %rcx     # 0x40a00
cmovneq %rcx, %rbx
movq (%rbx), %r9
movb 0x10(%r9), %r8b
leaq 0x8(%r9), %rdi
cmpb $0x1, %r8b
cmovneq %rcx, %r9
movq %rcx, %r10
cmoveq %rdi, %r10
movq (%r9), %rbx
movq (%r10), %r9
cmpb $0x1, 0x10(%r9)
jne 0x18751
cmpb $0x1, %r8b
leaq 0x282b6(%rip), %rcx     # 0x40a00
cmoveq %rdi, %rcx
movq (%rcx), %rcx
movb 0x10(%r14), %dil
cmpb $0x1, %dil
movq %rax, %r8
cmoveq %r14, %r8
cmpq %rsi, (%r8)
jne 0x187ae
addq $0x8, %r14
cmpb $0x1, %dil
cmovneq %rax, %r14
movq (%r14), %rdi
movb 0x10(%rdi), %sil
leaq 0x8(%rdi), %rdx
cmpb $0x1, %sil
cmovneq %rax, %rdi
movq %rax, %r8
cmoveq %rdx, %r8
movq (%rdi), %r14
movq (%r8), %rdi
cmpb $0x1, 0x10(%rdi)
jne 0x187ab
cmpb $0x1, %sil
leaq 0x2825e(%rip), %rax     # 0x40a00
cmoveq %rdx, %rax
movq (%rax), %rdx
jmp 0x187ae
movq %rax, %rdx
movq (%rcx), %r12
movq (%rdx), %r15
movq %rbx, %rdi
movq %r14, %rsi
callq 0x8c79
movl %eax, %ecx
xorl %eax, %eax
testl %ecx, %ecx
jne 0x1881d
cmpb $0x5, 0x10(%r14)
cmoveq %r14, %rbx
movq 0x2aad9(%rip), %rdi     # 0x432b0
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movq %r12, %rdi
callq 0x8b01
movq %r15, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa975
movl $0x3, %edi
callq 0x8ad9
callq 0x8b78
leaq 0x2adac(%rip), %rcx     # 0x435c0
movq %rax, (%rcx,%r13,8)
movl $0x1, %eax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: cmp_factors_provisional
int
cmp_factors_provisional(struct atom *p1, struct atom *p2)
{
	if (car(p1) == symbol(POWER))
		p1 = cadr(p1); // p1 = base

	if (car(p2) == symbol(POWER))
		p2 = cadr(p2); // p2 = base

	return cmp(p1, p2);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
movb 0x10(%rdi), %dl
cmpb $0x1, %dl
leaq 0x281b3(%rip), %rax     # 0x40a00
movq %rax, %r8
cmoveq %rdi, %r8
movq 0x2aa55(%rip), %rcx     # 0x432b0
cmpq %rcx, (%r8)
jne 0x18879
addq $0x8, %rdi
cmpb $0x1, %dl
cmovneq %rax, %rdi
movq (%rdi), %rdx
cmpb $0x1, 0x10(%rdx)
cmovneq %rax, %rdx
movq (%rdx), %rdi
movb 0x10(%rsi), %dl
cmpb $0x1, %dl
movq %rax, %r8
cmoveq %rsi, %r8
cmpq %rcx, (%r8)
jne 0x188a4
addq $0x8, %rsi
cmpb $0x1, %dl
cmovneq %rax, %rsi
movq (%rsi), %rcx
cmpb $0x1, 0x10(%rcx)
cmoveq %rcx, %rax
movq (%rax), %rsi
popq %rbp
jmp 0x8c79
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: order_factor
int
order_factor(struct atom *p)
{
	if (isnum(p))
		return 1;

	if (p == symbol(EXP1))
		return 5;

	if (car(p) == symbol(DERIVATIVE) || car(p) == symbol(D_LOWER))
		return 6;

	if (car(p) == symbol(POWER)) {

		p = cadr(p); // p = base

		if (isminusone(p))
			return 3;

		if (isnum(p))
			return 2;

		if (p == symbol(EXP1))
			return 5;

		if (car(p) == symbol(DERIVATIVE) || car(p) == symbol(D_LOWER))
			return 6;
	}

	return 4;
}

// Assembly:
movb 0x10(%rdi), %cl
movl %ecx, %edx
andb $-0x2, %dl
movl $0x1, %eax
cmpb $0x4, %dl
je 0x189bb
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq 0x2a9e3(%rip), %r15     # 0x432c8
movl $0x5, %eax
cmpq %rdi, %r15
je 0x189ad
cmpb $0x1, %cl
leaq 0x28103(%rip), %r14     # 0x40a00
movq %r14, %rax
cmoveq %rdi, %rax
movq (%rax), %rdx
movq 0x261d2(%rip), %r12     # 0x3eae0
movl $0x6, %eax
cmpq %r12, %rdx
je 0x189ad
movq 0x261a5(%rip), %r13     # 0x3eac8
cmpq %r13, %rdx
je 0x189ad
cmpq 0x2a97d(%rip), %rdx     # 0x432b0
je 0x1893c
movl $0x4, %eax
jmp 0x189ad
addq $0x8, %rdi
cmpb $0x1, %cl
cmovneq %r14, %rdi
movq (%rdi), %rax
cmpb $0x1, 0x10(%rax)
cmovneq %r14, %rax
movq (%rax), %rbx
movq %rbx, %rdi
movl $0xffffffff, %esi       # imm = 0xFFFFFFFF
movl $0x1, %edx
callq 0xcb82
movl %eax, %ecx
movl $0x3, %eax
testl %ecx, %ecx
jne 0x189ad
movb 0x10(%rbx), %cl
movl %ecx, %edx
andb $-0x2, %dl
movl $0x2, %eax
cmpb $0x4, %dl
je 0x189ad
movl $0x5, %eax
cmpq %r15, %rbx
je 0x189ad
cmpb $0x1, %cl
cmoveq %rbx, %r14
movq (%r14), %rax
cmpq %r12, %rax
sete %cl
cmpq %r13, %rax
sete %al
orb %cl, %al
movl $0x6, %eax
je 0x18935
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: multiply_rationals
void
multiply_rationals(struct atom *p1, struct atom *p2)
{
	int sign;
	uint32_t *a, *b, *c;

	if (iszero(p1) || iszero(p2)) {
		push_integer(0);
		return;
	}

	if (p1->sign == p2->sign)
		sign = MPLUS;
	else
		sign = MMINUS;

	if (isinteger(p1) && isinteger(p2)) {
		push_bignum(sign, mmul(p1->u.q.a, p2->u.q.a), mint(1));
		return;
	}

	a = mmul(p1->u.q.a, p2->u.q.a);
	b = mmul(p1->u.q.b, p2->u.q.b);
	c = mgcd(a, b);
	push_bignum(sign, mdiv(a, c), mdiv(b, c));

	mfree(a);
	mfree(b);
	mfree(c);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rsi, %r14
movq %rdi, %r15
callq 0xaece
testl %eax, %eax
jne 0x189e5
movq %r14, %rdi
callq 0xaece
testl %eax, %eax
je 0x189ff
movq 0x25754(%rip), %rdi     # 0x3e140
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x8b01
movb 0x12(%r15), %al
xorl %ebx, %ebx
cmpb 0x12(%r14), %al
setne %bl
cmpb $0x4, 0x10(%r15)
jne 0x18a6b
movq 0x8(%r15), %rax
cmpl $0x1, -0x4(%rax)
jne 0x18a6b
cmpl $0x1, (%rax)
jne 0x18a6b
cmpb $0x4, 0x10(%r14)
jne 0x18a6b
movq 0x8(%r14), %rax
cmpl $0x1, -0x4(%rax)
jne 0x18a6b
cmpl $0x1, (%rax)
jne 0x18a6b
movq (%r15), %rdi
movq (%r14), %rsi
callq 0x7b16
movq %rax, %r14
movl $0x1, %edi
callq 0x7ad5
movl %ebx, %edi
movq %r14, %rsi
movq %rax, %rdx
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xb9f7
movq (%r15), %rdi
movq (%r14), %rsi
callq 0x7b16
movq %rax, %r12
movq 0x8(%r15), %rdi
movq 0x8(%r14), %rsi
callq 0x7b16
movq %rax, %r15
movq %r12, %rdi
movq %rax, %rsi
callq 0x8713
movq %rax, %r14
movq %r12, %rdi
movq %rax, %rsi
callq 0x8049
movq %rax, %r13
movq %r15, %rdi
movq %r14, %rsi
callq 0x8049
movl %ebx, %edi
movq %r13, %rsi
movq %rax, %rdx
callq 0xb9f7
addq $-0x4, %r12
movq %r12, %rdi
callq 0x7040
movl 0x2568d(%rip), %ebx     # 0x3e15c
leal -0x1(%rbx), %eax
movl %eax, 0x25684(%rip)     # 0x3e15c
addq $-0x4, %r15
movq %r15, %rdi
callq 0x7040
leal -0x2(%rbx), %eax
movl %eax, 0x2566f(%rip)     # 0x3e15c
addq $-0x4, %r14
movq %r14, %rdi
callq 0x7040
addl $-0x3, %ebx
movl %ebx, 0x2565a(%rip)     # 0x3e15c
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: any_radical_factors
int
any_radical_factors(int h)
{
	int i;
	for (i = h; i < tos; i++)
		if (isradical(stack[i]))
			return 1;
	return 0;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movl 0x21796(%rip), %ebx     # 0x3a2b8
subl %edi, %ebx
jle 0x18b50
movslq %edi, %r14
shlq $0x3, %r14
xorl %r15d, %r15d
leaq 0x2aa89(%rip), %r12     # 0x435c0
leaq (%r14,%r15,8), %rax
movq (%r12,%rax), %rdi
callq 0x18f9e
testl %eax, %eax
jne 0x18b54
incq %r15
cmpl %r15d, %ebx
jne 0x18b37
xorl %eax, %eax
jmp 0x18b59
movl $0x1, %eax
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: reduce_radical_double
struct atom *
reduce_radical_double(int h, struct atom *COEF)
{
	int i;
	double a, b, c;
	struct atom *p1;

	c = COEF->u.d;

	for (i = h; i < tos; i++) {

		p1 = stack[i];

		if (isradical(p1)) {

			push(cadr(p1)); // base
			a = pop_double();

			push(caddr(p1)); // exponent
			b = pop_double();

			c = c * pow(a, b); // a > 0 by isradical above

			slice(i, 1); // remove factor

			i--; // use same index again
		}
	}

	push_double(c);
	COEF = pop();

	return COEF;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movsd (%rsi), %xmm0
movsd %xmm0, -0x30(%rbp)
cmpl %edi, 0x2141a(%rip)     # 0x3a2b8
jle 0x18f71
movl %edi, %ebx
leaq 0x2a713(%rip), %r15     # 0x435c0
leaq 0x27b4c(%rip), %r12     # 0x40a00
movslq %ebx, %rax
movq (%r15,%rax,8), %r14
movq %r14, %rdi
callq 0x18f9e
testl %eax, %eax
je 0x18f63
cmpb $0x1, 0x10(%r14)
leaq 0x8(%r14), %r13
movq %r12, %rax
cmoveq %r13, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmovneq %r12, %rax
movq (%rax), %rdi
callq 0x8b01
callq 0x9086
cmpb $0x1, 0x10(%r14)
cmovneq %r12, %r13
movsd %xmm0, -0x38(%rbp)
movq (%r13), %rax
movb 0x10(%rax), %cl
addq $0x8, %rax
cmpb $0x1, %cl
movq %r12, %rdx
cmoveq %rax, %rdx
movq (%rdx), %rdx
cmpb $0x1, 0x10(%rdx)
movq %r12, %rdx
jne 0x18f2c
cmpb $0x1, %cl
cmovneq %r12, %rax
movq (%rax), %rdx
movq (%rdx), %rdi
callq 0x8b01
callq 0x9086
movapd %xmm0, %xmm1
movsd -0x38(%rbp), %xmm0
callq 0x70d0
movsd -0x30(%rbp), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, -0x30(%rbp)
movl %ebx, %edi
movl $0x1, %esi
callq 0xadfc
decl %ebx
incl %ebx
cmpl 0x2134d(%rip), %ebx     # 0x3a2b8
jl 0x18eb4
callq 0x74ec
movb $0x5, 0x10(%rax)
movsd -0x30(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x8b78
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: multiply_noexpand
void
multiply_noexpand(void)
{
	int t;
	t = expanding;
	expanding = 0;
	multiply();
	expanding = t;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
pushq %rax
movl 0xeda5f(%rip), %ebx     # 0x106ad0
movl $0x0, 0xeda55(%rip)     # 0x106ad0
movl $0x2, %edi
callq 0xa8bb
movl %ebx, 0xeda45(%rip)     # 0x106ad0
addq $0x8, %rsp
popq %rbx
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: eval_noexpand
void
eval_noexpand(struct atom *p1)
{
	int t;

	t = expanding;
	expanding = 0;

	push(cadr(p1));
	evalf();

	expanding = t;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
pushq %rax
movl 0xeda21(%rip), %ebx     # 0x106ad0
movl $0x0, 0xeda17(%rip)     # 0x106ad0
cmpb $0x1, 0x10(%rdi)
leaq 0x8(%rdi), %rax
leaq 0x27938(%rip), %rcx     # 0x40a00
cmovneq %rcx, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmoveq %rax, %rcx
movq (%rcx), %rdi
callq 0x8b01
callq 0x95e2
movl %ebx, 0xed9e6(%rip)     # 0x106ad0
addq $0x8, %rsp
popq %rbx
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: eval_not
void
eval_not(struct atom *p1)
{
	push(cadr(p1));
	evalp();
	p1 = pop();
	if (iszero(p1))
		push_integer(1);
	else
		push_integer(0);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
cmpb $0x1, 0x10(%rdi)
leaq 0x8(%rdi), %rax
leaq 0x278fc(%rip), %rcx     # 0x40a00
cmovneq %rcx, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmoveq %rax, %rcx
movq (%rcx), %rdi
callq 0x8b01
callq 0x9ae6
callq 0x8b78
movq %rax, %rdi
callq 0xaece
testl %eax, %eax
leaq 0x2500a(%rip), %rax     # 0x3e140
leaq 0xed99b(%rip), %rcx     # 0x106ad8
cmoveq %rax, %rcx
movq (%rcx), %rdi
popq %rbp
jmp 0x8b01
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: eval_nroots
void
eval_nroots(struct atom *p1)
{
	push(cadr(p1));
	evalf();

	p1 = cddr(p1);

	if (iscons(p1)) {
		push(car(p1));
		evalf();
	} else
		push_symbol(X_LOWER);

	nroots();
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
leaq 0x8(%rdi), %r14
cmpb $0x1, 0x10(%rdi)
leaq 0x2789a(%rip), %r15     # 0x40a00
movq %r15, %rax
cmoveq %r14, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmovneq %r15, %rax
movq (%rax), %rdi
callq 0x8b01
callq 0x95e2
cmpb $0x1, 0x10(%rbx)
cmovneq %r15, %r14
movq (%r14), %rax
leaq 0x8(%rax), %rcx
cmpb $0x1, 0x10(%rax)
cmovneq %r15, %rcx
movq (%rcx), %rax
cmpb $0x1, 0x10(%rax)
jne 0x191b4
movq (%rax), %rdi
callq 0x8b01
callq 0x95e2
jmp 0x191c0
movq 0x2978d(%rip), %rdi     # 0x42948
callq 0x8b01
addq $0x8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
jmp 0x191cf
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: coeffs
void
coeffs(struct atom *P, struct atom *X)
{
	struct atom *C;

	for (;;) {

		push(P);
		push(X);
		push_integer(0);
		subst();
		evalf();
		C = pop();

		push(C);

		push(P);
		push(C);
		subtract();
		P = pop();

		if (iszero(P))
			break;

		push(P);
		push(X);
		divide();
		P = pop();
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r14
movq %r14, %rdi
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movq 0x24b9a(%rip), %rdi     # 0x3e140
callq 0x8b01
callq 0x9402
callq 0x95e2
callq 0x8b78
movq %rax, %r15
movq %rax, %rdi
callq 0x8b01
movq %r14, %rdi
callq 0x8b01
movq %r15, %rdi
callq 0x8b01
callq 0xba6d
callq 0x8b78
movq %rax, %r14
movq %rax, %rdi
callq 0xaece
testl %eax, %eax
jne 0x1960d
movq %r14, %rdi
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
callq 0xd9a4
callq 0x8b78
movq %rax, %r14
jmp 0x1958f
addq $0x8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: nreduce
void
nreduce(double cr[], double ci[], int n, double ar, double ai)
{
	int k;

	// divide

	for (k = n - 1; k > 0; k--) {
		cr[k - 1] += cr[k] * ar - ci[k] * ai;
		ci[k - 1] += ci[k] * ar + cr[k] * ai;
	}

	if (zabs(cr[0], ci[0]) > DELTA)
		stopf("nroots: residual error"); // not a root

	// shift

	for (k = 0; k < n - 1; k++) {
		cr[k] = cr[k + 1];
		ci[k] = ci[k + 1];
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movl %edx, %r15d
movq %rsi, %rbx
movq %rdi, %r14
leal -0x1(%r15), %r12d
cmpl $0x2, %edx
jl 0x19ab9
movl %r12d, %eax
incq %rax
movsd -0x8(%r14,%rax,8), %xmm2
movsd -0x8(%rbx,%rax,8), %xmm3
mulsd %xmm1, %xmm3
movapd %xmm2, %xmm4
mulsd %xmm0, %xmm4
subsd %xmm3, %xmm4
addsd -0x10(%r14,%rax,8), %xmm4
movsd %xmm4, -0x10(%r14,%rax,8)
mulsd %xmm1, %xmm2
movsd -0x8(%rbx,%rax,8), %xmm3
mulsd %xmm0, %xmm3
addsd %xmm2, %xmm3
addsd -0x10(%rbx,%rax,8), %xmm3
movsd %xmm3, -0x10(%rbx,%rax,8)
decq %rax
cmpq $0x1, %rax
jg 0x19a67
movsd (%r14), %xmm0
movsd (%rbx), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
addsd %xmm1, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x19ade
sqrtsd %xmm0, %xmm0
jmp 0x19ae3
callq 0x7350
ucomisd 0xc65d(%rip), %xmm0     # 0x26148
ja 0x19b25
cmpl $0x2, %r15d
jl 0x19b1c
movl %r12d, %eax
xorl %ecx, %ecx
movsd 0x8(%r14,%rcx,8), %xmm0
movsd %xmm0, (%r14,%rcx,8)
movsd 0x8(%rbx,%rcx,8), %xmm0
movsd %xmm0, (%rbx,%rcx,8)
leaq 0x1(%rcx), %rdx
movq %rdx, %rcx
cmpq %rdx, %rax
jne 0x19af8
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
leaq 0x10da0(%rip), %rdi     # 0x2a8cc
callq 0x77c8
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: fata
void
fata(double cr[], double ci[], int n, double ar, double ai, double *far, double *fai)
{
	int k;
	double xr, xi, yr, yi;

	yr = cr[n - 1];
	yi = ci[n - 1];

	for (k = n - 2; k >= 0; k--) {

		// x = a * y

		xr = ar * yr - ai * yi;
		xi = ar * yi + ai * yr;

		// y = x + c

		yr = xr + cr[k];
		yi = xi + ci[k];
	}

	*far = yr;
	*fai = yi;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
movslq %edx, %rax
movsd -0x8(%rdi,%rax,8), %xmm3
movsd -0x8(%rsi,%rax,8), %xmm2
cmpl $0x2, %eax
jl 0x19bb1
addl $-0x2, %edx
incq %rdx
movapd %xmm0, %xmm5
mulsd %xmm2, %xmm5
mulsd %xmm1, %xmm2
movapd %xmm0, %xmm4
mulsd %xmm3, %xmm4
subsd %xmm2, %xmm4
mulsd %xmm1, %xmm3
addsd %xmm5, %xmm3
addsd -0x8(%rdi,%rdx,8), %xmm4
movapd %xmm3, %xmm2
addsd -0x8(%rsi,%rdx,8), %xmm2
movapd %xmm4, %xmm3
decq %rdx
jg 0x19b76
jmp 0x19bb5
movapd %xmm3, %xmm4
movsd %xmm4, (%rcx)
movsd %xmm2, (%r8)
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: eval_number
void
eval_number(struct atom *p1)
{
	push(cadr(p1));
	evalf();

	p1 = pop();

	if (isnum(p1))
		push_integer(1);
	else
		push_integer(0);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
cmpb $0x1, 0x10(%rdi)
leaq 0x8(%rdi), %rax
leaq 0x26e07(%rip), %rcx     # 0x40a00
cmovneq %rcx, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmoveq %rax, %rcx
movq (%rcx), %rdi
callq 0x8b01
callq 0x95e2
callq 0x8b78
movb 0x10(%rax), %al
andb $-0x2, %al
cmpb $0x4, %al
leaq 0xeceb0(%rip), %rax     # 0x106ad8
leaq 0x24511(%rip), %rcx     # 0x3e140
cmoveq %rax, %rcx
movq (%rcx), %rdi
popq %rbp
jmp 0x8b01
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: eval_outer
void
eval_outer(struct atom *p1)
{
	push(cadr(p1));
	evalf();
	p1 = cddr(p1);
	while (iscons(p1)) {
		push(car(p1));
		evalf();
		outer();
		p1 = cdr(p1);
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
leaq 0x8(%rdi), %r15
cmpb $0x1, 0x10(%rdi)
leaq 0x26cd8(%rip), %r14     # 0x40a00
movq %r14, %rax
cmoveq %r15, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmovneq %r14, %rax
movq (%rax), %rdi
callq 0x8b01
callq 0x95e2
cmpb $0x1, 0x10(%rbx)
cmovneq %r14, %r15
movq (%r15), %rax
leaq 0x8(%rax), %rcx
cmpb $0x1, 0x10(%rax)
cmovneq %r14, %rcx
movq (%rcx), %rbx
cmpb $0x1, 0x10(%rbx)
jne 0x19d8a
movq (%rbx), %rdi
callq 0x8b01
callq 0x95e2
callq 0x19d95
leaq 0x8(%rbx), %rax
cmpb $0x1, 0x10(%rbx)
cmovneq %r14, %rax
movq (%rax), %rbx
jmp 0x19d61
addq $0x8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: outer
void
outer(void)
{
	int i, j, k, n, ncol, nrow;
	struct atom *p1, *p2, *p3;

	p2 = pop();
	p1 = pop();

	if (!istensor(p1) || !istensor(p2)) {
		push(p1);
		push(p2);
		multiply();
		return;
	}

	if (p1->u.tensor->ndim + p2->u.tensor->ndim > MAXDIM)
		stopf("outer: rank exceeds max");

	// sync diffs

	nrow = p1->u.tensor->nelem;
	ncol = p2->u.tensor->nelem;

	p3 = alloc_tensor(nrow * ncol);

	for (i = 0; i < nrow; i++)
		for (j = 0; j < ncol; j++) {
			push(p1->u.tensor->elem[i]);
			push(p2->u.tensor->elem[j]);
			multiply();
			p3->u.tensor->elem[i * ncol + j] = pop();
		}

	// dim info

	p3->u.tensor->ndim = p1->u.tensor->ndim + p2->u.tensor->ndim;

	k = 0;

	n = p1->u.tensor->ndim;

	for (i = 0; i < n; i++)
		p3->u.tensor->dim[k++] = p1->u.tensor->dim[i];

	n = p2->u.tensor->ndim;

	for (i = 0; i < n; i++)
		p3->u.tensor->dim[k++] = p2->u.tensor->dim[i];

	push(p3);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
callq 0x8b78
movq %rax, %rbx
callq 0x8b78
movq %rax, %rdi
cmpb $0x7, 0x10(%rax)
jne 0x19f01
cmpb $0x7, 0x10(%rbx)
jne 0x19f01
movq (%rdi), %rax
movq (%rbx), %rcx
movl (%rcx), %ecx
addl (%rax), %ecx
cmpl $0x19, %ecx
jge 0x19f26
movq %rdi, -0x30(%rbp)
movq (%rdi), %rax
movslq 0x64(%rax), %r14
movq %rbx, -0x38(%rbp)
movq (%rbx), %rax
movl 0x64(%rax), %ebx
movl %ebx, %edi
imull %r14d, %edi
callq 0x7850
movq %rax, %r9
movq %r14, -0x48(%rbp)
testq %r14, %r14
jle 0x19e84
leaq (,%rbx,8), %rax
movq %rax, -0x40(%rbp)
movl $0x68, %r15d
xorl %ecx, %ecx
movq -0x30(%rbp), %r13
testl %ebx, %ebx
jle 0x19e77
xorl %r14d, %r14d
movq %rcx, -0x50(%rbp)
movq (%r13), %rax
movq 0x68(%rax,%rcx,8), %rdi
movq %r9, %r12
callq 0x8b01
movq -0x38(%rbp), %rax
movq (%rax), %rax
movq 0x68(%rax,%r14,8), %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
callq 0x8b78
movq %r12, %r9
movq (%r12), %rcx
addq %r15, %rcx
movq %rax, (%rcx,%r14,8)
movq -0x50(%rbp), %rcx
incq %r14
cmpq %r14, %rbx
jne 0x19e2c
incq %rcx
addq -0x40(%rbp), %r15
cmpq -0x48(%rbp), %rcx
jne 0x19e21
movq -0x30(%rbp), %rdi
movq (%rdi), %rax
movq -0x38(%rbp), %r8
movq (%r8), %rcx
movl (%rcx), %ecx
addl (%rax), %ecx
movq (%r9), %rax
movl %ecx, (%rax)
movq (%rdi), %rax
movslq (%rax), %rcx
xorl %eax, %eax
testq %rcx, %rcx
jle 0x19ebe
movq (%rdi), %rdx
movl 0x4(%rdx,%rax,4), %edx
movq (%r9), %rsi
movl %edx, 0x4(%rsi,%rax,4)
incq %rax
cmpq %rax, %rcx
jne 0x19ea8
movq (%r8), %rcx
movslq (%rcx), %rcx
testq %rcx, %rcx
jle 0x19eeb
leaq 0x4(,%rax,4), %rax
xorl %edx, %edx
movq (%r8), %rsi
movl 0x4(%rsi,%rdx,4), %esi
movq (%r9), %rdi
addq %rax, %rdi
movl %esi, (%rdi,%rdx,4)
incq %rdx
cmpq %rdx, %rcx
jne 0x19ed3
movq %r9, %rdi
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x8b01
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movl $0x2, %edi
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xa8bb
leaq 0x109b6(%rip), %rdi     # 0x2a8e3
callq 0x77c8
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: eval_power
void
eval_power(struct atom *p1)
{
	int t;
	struct atom *p2;

	expanding--; // undo expanding++ in evalf

	// base

	push(cadr(p1));

	// exponent

	push(caddr(p1));
	evalg();
	dupl();
	p2 = pop();

	// if exponent is negative then evaluate base without expanding

	swap();
	if (isnegativenumber(p2)) {
		t = expanding;
		expanding = 0;
		evalg();
		expanding = t;
	} else
		evalg();
	swap();

	power();

	expanding++;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %rbx
pushq %rax
decl 0xecb59(%rip)           # 0x106ad0
movq %rdi, %rbx
leaq 0x8(%rdi), %r15
cmpb $0x1, 0x10(%rdi)
leaq 0x26a77(%rip), %r14     # 0x40a00
movq %r14, %rax
cmoveq %r15, %rax
movq (%rax), %rax
cmpb $0x1, 0x10(%rax)
cmovneq %r14, %rax
movq (%rax), %rdi
callq 0x8b01
cmpb $0x1, 0x10(%rbx)
cmovneq %r14, %r15
movq (%r15), %rax
movb 0x10(%rax), %cl
addq $0x8, %rax
cmpb $0x1, %cl
movq %r14, %rdx
cmoveq %rax, %rdx
movq (%rdx), %rdx
cmpb $0x1, 0x10(%rdx)
jne 0x19fd9
cmpb $0x1, %cl
leaq 0x26a2e(%rip), %rcx     # 0x40a00
cmoveq %rax, %rcx
movq (%rcx), %r14
movq (%r14), %rdi
callq 0x8b01
callq 0x95ac
callq 0x8b78
movq %rax, %rbx
movq %rax, %rdi
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
callq 0x8b78
movq %rax, %rbx
callq 0xad38
movzbl 0x10(%rbx), %eax
cmpl $0x5, %eax
je 0x1a022
cmpl $0x4, %eax
jne 0x1a04e
cmpb $0x1, 0x12(%rbx)
sete %al
jmp 0x1a02d
xorpd %xmm0, %xmm0
ucomisd (%rbx), %xmm0
seta %al
testb %al, %al
je 0x1a04e
movl 0xeca99(%rip), %ebx     # 0x106ad0
movl $0x0, 0xeca8f(%rip)     # 0x106ad0
callq 0x95ac
movl %ebx, 0xeca84(%rip)     # 0x106ad0
jmp 0x1a053
callq 0x95ac
callq 0xad38
callq 0xa234
incl 0xeca6d(%rip)           # 0x106ad0
addq $0x8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_natural_number
void
power_natural_number(struct atom *EXPO)
{
	double x, y;

	// exp(x + i y) = exp(x) (cos(y) + i sin(y))

	if (isdoublez(EXPO)) {
		if (car(EXPO) == symbol(ADD)) {
			x = cadr(EXPO)->u.d;
			y = cadaddr(EXPO)->u.d;
		} else {
			x = 0.0;
			y = cadr(EXPO)->u.d;
		}
		push_double(exp(x));
		push_double(y);
		cosfunc();
		push(imaginaryunit);
		push_double(y);
		sinfunc();
		multiply();
		add();
		multiply();
		return;
	}

	// e^log(expr) = expr

	if (car(EXPO) == symbol(LOG)) {
		push(cadr(EXPO));
		return;
	}

	if (isdenormalpolar(EXPO)) {
		normalize_polar(EXPO);
		return;
	}

	push_symbol(POWER);
	push_symbol(EXP1);
	push(EXPO);
	list(3);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
subq $0x18, %rsp
movq %rdi, %rbx
callq 0xc3a9
movb 0x10(%rbx), %dl
cmpb $0x1, %dl
leaq 0x2652f(%rip), %rcx     # 0x40a00
movq %rcx, %rsi
cmoveq %rbx, %rsi
movq (%rsi), %rsi
testl %eax, %eax
je 0x1a510
addq $0x8, %rbx
cmpb $0x1, %dl
cmovneq %rcx, %rbx
movq (%rbx), %rax
movb 0x10(%rax), %dil
cmpb $0x1, %dil
movq %rcx, %rdx
cmoveq %rax, %rdx
xorps %xmm0, %xmm0
cmpq 0x28d9a(%rip), %rsi     # 0x432a0
je 0x1a53b
movq %rdx, %rcx
jmp 0x1a5ba
cmpq 0x25eb9(%rip), %rsi     # 0x403d0
je 0x1a657
movq %rbx, %rdi
callq 0x1b13d
testl %eax, %eax
je 0x1a682
movq %rbx, %rdi
addq $0x18, %rsp
popq %rbx
popq %rbp
jmp 0x1b1c7
addq $0x8, %rax
cmpb $0x1, %dil
movq (%rdx), %rsi
cmovneq %rcx, %rax
movq (%rax), %rax
movb 0x10(%rax), %dl
cmpb $0x1, %dl
movq %rcx, %rdi
cmoveq %rax, %rdi
movq (%rdi), %rdi
cmpb $0x1, 0x10(%rdi)
movq %rcx, %rdi
jne 0x1a57b
cmpb $0x1, %dl
leaq 0x26490(%rip), %rdi     # 0x40a00
cmoveq %rax, %rdi
movq (%rdi), %rdi
addq $0x8, %rdi
movsd (%rsi), %xmm0
movq (%rdi), %rsi
cmpb $0x1, 0x10(%rsi)
jne 0x1a5ba
cmpb $0x1, %dl
leaq 0x2646e(%rip), %rcx     # 0x40a00
movq %rcx, %rsi
cmoveq %rax, %rsi
movq (%rsi), %rsi
cmpb $0x1, 0x10(%rsi)
jne 0x1a5b7
cmpb $0x1, %dl
leaq 0x26454(%rip), %rcx     # 0x40a00
cmoveq %rax, %rcx
movq (%rcx), %rcx
addq $0x8, %rcx
movq (%rcx), %rcx
movq (%rcx), %rax
movsd (%rax), %xmm1
movsd %xmm1, -0x10(%rbp)
callq 0x7380
movsd %xmm0, -0x18(%rbp)
callq 0x74ec
movb $0x5, %bl
movb %bl, 0x10(%rax)
movsd -0x18(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
callq 0x74ec
movb %bl, 0x10(%rax)
movsd -0x10(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
callq 0xf2ab
movq 0xec4d0(%rip), %rdi     # 0x106ae0
callq 0x8b01
callq 0x74ec
movb %bl, 0x10(%rax)
movsd -0x10(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
callq 0xfade
movl $0x2, %edi
callq 0xa8bb
movl $0x2, %edi
callq 0xa975
movl $0x2, %edi
addq $0x18, %rsp
popq %rbx
popq %rbp
jmp 0xa8bb
addq $0x8, %rbx
cmpb $0x1, %dl
leaq 0x2639b(%rip), %rax     # 0x40a00
cmovneq %rax, %rbx
movq (%rbx), %rcx
cmpb $0x1, 0x10(%rcx)
cmoveq %rcx, %rax
movq (%rax), %rdi
addq $0x18, %rsp
popq %rbx
popq %rbp
jmp 0x8b01
movq 0x28c27(%rip), %rdi     # 0x432b0
callq 0x8b01
movq 0x28c33(%rip), %rdi     # 0x432c8
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movl $0x3, %edi
addq $0x18, %rsp
popq %rbx
popq %rbp
jmp 0x8ad9
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_sum
void
power_sum(struct atom *BASE, struct atom *EXPO)
{
	int h, i, n;
	struct atom *p1, *p2;

	if (iscomplexnumber(BASE) && isnum(EXPO)) {
		power_complex_number(BASE, EXPO);
		return;
	}

	if (expanding == 0 || !issmallinteger(EXPO) || isnegativenumber(EXPO)) {
		push_symbol(POWER);
		push(BASE);
		push(EXPO);
		list(3);
		return;
	}

	push(EXPO);
	n = pop_integer();

	// square the sum first (prevents infinite loop through multiply)

	h = tos;

	p1 = cdr(BASE);

	while (iscons(p1)) {
		p2 = cdr(BASE);
		while (iscons(p2)) {
			push(car(p1));
			push(car(p2));
			multiply();
			p2 = cdr(p2);
		}
		p1 = cdr(p1);
	}

	add_terms(tos - h);

	// continue up to power n

	for (i = 2; i < n; i++) {
		push(BASE);
		multiply();
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x18, %rsp
movq %rsi, %r14
movq %rdi, %rbx
callq 0x1bb01
testl %eax, %eax
je 0x1a6f5
movb 0x10(%r14), %al
andb $-0x2, %al
cmpb $0x4, %al
jne 0x1a6f5
movq %rbx, %rdi
movq %r14, %rsi
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x1bbd2
cmpl $0x0, 0xec3d4(%rip)     # 0x106ad0
je 0x1a733
movq %r14, %rdi
callq 0x1433e
testl %eax, %eax
je 0x1a733
movzbl 0x10(%r14), %eax
cmpl $0x5, %eax
je 0x1a723
cmpl $0x4, %eax
jne 0x1a767
cmpb $0x1, 0x12(%r14)
sete %al
jmp 0x1a72f
xorpd %xmm0, %xmm0
ucomisd (%r14), %xmm0
seta %al
testb %al, %al
je 0x1a767
movq 0x28b76(%rip), %rdi     # 0x432b0
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movq %r14, %rdi
callq 0x8b01
movl $0x3, %edi
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x8ad9
movq %r14, %rdi
callq 0x8b01
callq 0xeb5b
movl %eax, %r14d
leaq 0x8(%rbx), %rcx
cmpb $0x1, 0x10(%rbx)
leaq 0x2627a(%rip), %r13     # 0x40a00
movq %r13, %rax
movq %rcx, -0x40(%rbp)
cmoveq %rcx, %rax
movl 0x1fb21(%rip), %ecx     # 0x3a2b8
movl %ecx, -0x2c(%rbp)
movq (%rax), %r12
cmpb $0x1, 0x10(%r12)
jne 0x1a822
leaq 0x10(%r12), %r15
movq %rbx, -0x38(%rbp)
cmpb $0x1, 0x10(%rbx)
movq %r13, %rax
cmoveq -0x40(%rbp), %rax
jmp 0x1a7f8
cmpb $0x1, (%r15)
movq %r13, %rax
cmoveq %r12, %rax
movq (%rax), %rdi
callq 0x8b01
cmpb $0x1, 0x10(%rbx)
movq %r13, %rax
cmoveq %rbx, %rax
movq (%rax), %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
leaq 0x8(%rbx), %rax
cmpb $0x1, 0x10(%rbx)
cmovneq %r13, %rax
movq (%rax), %rbx
cmpb $0x1, 0x10(%rbx)
je 0x1a7bc
addq $0x8, %r12
cmpb $0x1, (%r15)
cmovneq %r13, %r12
movq (%r12), %r12
leaq 0x10(%r12), %r15
cmpb $0x1, 0x10(%r12)
movq -0x38(%rbp), %rbx
je 0x1a7ae
movl 0x1fa90(%rip), %edi     # 0x3a2b8
subl -0x2c(%rbp), %edi
callq 0xa975
cmpl $0x3, %r14d
jl 0x1a851
addl $-0x2, %r14d
movq %rbx, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
decl %r14d
jne 0x1a83a
addq $0x18, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_double
void
power_double(struct atom *BASE, struct atom *EXPO)
{
	double base, d, expo;

	push(BASE);
	base = pop_double();

	push(EXPO);
	expo = pop_double();

	if (base > 0.0 || expo == floor(expo)) {
		d = pow(base, expo);
		push_double(d);
		return;
	}

	// BASE is negative and EXPO is fractional

	power_minusone(EXPO);

	if (base == -1.0)
		return;

	d = pow(-base, expo);
	push_double(d);

	multiply();
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
subq $0x28, %rsp
movq %rsi, %rbx
callq 0x8b01
callq 0x9086
movapd %xmm0, -0x20(%rbp)
movq %rbx, %rdi
callq 0x8b01
callq 0x9086
movsd %xmm0, -0x28(%rbp)
callq 0x7300
movsd -0x28(%rbp), %xmm1
xorpd %xmm2, %xmm2
movapd -0x20(%rbp), %xmm3
ucomisd %xmm2, %xmm3
movapd %xmm3, %xmm2
ja 0x1a8d2
ucomisd %xmm0, %xmm1
jne 0x1a8b2
jnp 0x1a8d2
movq %rbx, %rdi
callq 0x1abf4
movapd -0x20(%rbp), %xmm1
ucomisd 0xb8a1(%rip), %xmm1     # 0x26168
jne 0x1a900
jp 0x1a900
addq $0x28, %rsp
popq %rbx
popq %rbp
retq
movapd %xmm2, %xmm0
callq 0x70d0
movsd %xmm0, -0x20(%rbp)
callq 0x74ec
movb $0x5, 0x10(%rax)
movsd -0x20(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
addq $0x28, %rsp
popq %rbx
popq %rbp
jmp 0x8b01
movapd 0xb7e8(%rip), %xmm0     # 0x260f0
xorpd %xmm1, %xmm0
movsd -0x28(%rbp), %xmm1
callq 0x70d0
movsd %xmm0, -0x20(%rbp)
callq 0x74ec
movb $0x5, 0x10(%rax)
movsd -0x20(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
movl $0x2, %edi
addq $0x28, %rsp
popq %rbx
popq %rbp
jmp 0xa8bb
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_minusone
void
power_minusone(struct atom *EXPO)
{
	// optimization for i

	if (isequalq(EXPO, 1, 2)) {
		push(imaginaryunit);
		return;
	}

	// root is an odd number?

	if (isrational(EXPO) && EXPO->u.q.b[0] & 1) {
		if (EXPO->u.q.a[0] & 1)
			push_integer(-1);
		else
			push_integer(1);
		return;
	}

	if (isrational(EXPO)) {
		normalize_clock_rational(EXPO);
		return;
	}

	if (isdouble(EXPO)) {
		normalize_clock_double(EXPO);
		rect();
		return;
	}

	push_symbol(POWER);
	push_integer(-1);
	push(EXPO);
	list(3);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
pushq %rax
movq %rdi, %rbx
movl $0x1, %esi
movl $0x2, %edx
callq 0xcb82
testl %eax, %eax
je 0x1ac22
movq 0xebec9(%rip), %rdi     # 0x106ae0
addq $0x8, %rsp
popq %rbx
popq %rbp
jmp 0x8b01
movzbl 0x10(%rbx), %eax
cmpl $0x5, %eax
je 0x1ac51
cmpl $0x4, %eax
jne 0x1ac64
movq 0x8(%rbx), %rcx
testb $0x1, (%rcx)
jne 0x1ac94
cmpl $0x5, %eax
je 0x1ac51
cmpl $0x4, %eax
jne 0x1ac64
movq %rbx, %rdi
addq $0x8, %rsp
popq %rbx
popq %rbp
jmp 0x1acb4
movq %rbx, %rdi
callq 0x1af0c
addq $0x8, %rsp
popq %rbx
popq %rbp
jmp 0x145aa
movq 0x28645(%rip), %rdi     # 0x432b0
callq 0x8b01
movq 0xf6779(%rip), %rdi     # 0x1113f0
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movl $0x3, %edi
addq $0x8, %rsp
popq %rbx
popq %rbp
jmp 0x8ad9
movq (%rbx), %rax
testb $0x1, (%rax)
jne 0x1aca8
movq 0xebe35(%rip), %rdi     # 0x106ad8
jmp 0x1ac17
movq 0xf6741(%rip), %rdi     # 0x1113f0
jmp 0x1ac17
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: isdenormalpolar
int
isdenormalpolar(struct atom *p)
{
	if (car(p) == symbol(ADD)) {
		p = cdr(p);
		while (iscons(p)) {
			if (isdenormalpolarterm(car(p)))
				return 1;
			p = cdr(p);
		}
		return 0;
	}

	return isdenormalpolarterm(p);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movb 0x10(%rdi), %al
cmpb $0x1, %al
leaq 0x258ac(%rip), %r14     # 0x40a00
movq %r14, %rcx
cmoveq %rdi, %rcx
movq (%rcx), %rcx
cmpq 0x2813b(%rip), %rcx     # 0x432a0
je 0x1b174
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
jmp 0x1b2a6
addq $0x8, %rdi
cmpb $0x1, %al
cmovneq %r14, %rdi
movq (%rdi), %r15
xorl %ebx, %ebx
cmpb $0x1, 0x10(%r15)
jne 0x1b1bc
leaq 0x10(%r15), %r12
movq (%r15), %rdi
callq 0x1b2a6
testl %eax, %eax
jne 0x1b1b7
addq $0x8, %r15
cmpb $0x1, (%r12)
cmovneq %r14, %r15
movq (%r15), %r15
leaq 0x10(%r15), %r12
cmpb $0x1, 0x10(%r15)
je 0x1b18e
jmp 0x1b1bc
movl $0x1, %ebx
movl %ebx, %eax
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: normalize_polar
void
normalize_polar(struct atom *EXPO)
{
	int h;
	struct atom *p1;
	if (car(EXPO) == symbol(ADD)) {
		h = tos;
		p1 = cdr(EXPO);
		while (iscons(p1)) {
			EXPO = car(p1);
			if (isdenormalpolarterm(EXPO))
				normalize_polar_term(EXPO);
			else {
				push_symbol(POWER);
				push_symbol(EXP1);
				push(EXPO);
				list(3);
			}
			p1 = cdr(p1);
		}
		multiply_factors(tos - h);
	} else
		normalize_polar_term(EXPO);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movb 0x10(%rdi), %al
cmpb $0x1, %al
leaq 0x2581f(%rip), %r14     # 0x40a00
movq %r14, %rcx
cmoveq %rdi, %rcx
movq (%rcx), %rcx
cmpq 0x280ae(%rip), %rcx     # 0x432a0
je 0x1b207
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x1b52e
addq $0x8, %rdi
cmpb $0x1, %al
cmovneq %r14, %rdi
movl 0x1f0a0(%rip), %r15d    # 0x3a2b8
movq (%rdi), %r12
cmpb $0x1, 0x10(%r12)
jne 0x1b28a
leaq 0x10(%r12), %r13
movq (%r12), %rbx
movq %rbx, %rdi
callq 0x1b2a6
testl %eax, %eax
je 0x1b242
movq %rbx, %rdi
callq 0x1b52e
jmp 0x1b26c
movq 0x28067(%rip), %rdi     # 0x432b0
callq 0x8b01
movq 0x28073(%rip), %rdi     # 0x432c8
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movl $0x3, %edi
callq 0x8ad9
addq $0x8, %r12
cmpb $0x1, (%r13)
cmovneq %r14, %r12
movq (%r12), %r12
leaq 0x10(%r12), %r13
cmpb $0x1, 0x10(%r12)
je 0x1b228
movl 0x1f028(%rip), %edi     # 0x3a2b8
subl %r15d, %edi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xa8bb
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: normalize_polar_term
void
normalize_polar_term(struct atom *EXPO)
{
	struct atom *R;

	// exp(i pi) = -1

	if (lengthf(EXPO) == 3) {
		push_integer(-1);
		return;
	}

	R = cadr(EXPO); // R = coeff of term

	if (isrational(R))
		normalize_polar_term_rational(R);
	else
		normalize_polar_term_double(R);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
movb 0x10(%rdi), %al
cmpb $0x1, %al
jne 0x1b551
movl $0x3, %ecx
movq %rdi, %rdx
movq 0x8(%rdx), %rdx
decl %ecx
cmpb $0x1, 0x10(%rdx)
je 0x1b541
testl %ecx, %ecx
je 0x1b582
addq $0x8, %rdi
cmpb $0x1, %al
leaq 0x254a2(%rip), %rax     # 0x40a00
cmovneq %rax, %rdi
movq (%rdi), %rcx
cmpb $0x1, 0x10(%rcx)
cmoveq %rcx, %rax
movq (%rax), %rdi
cmpb $0x4, 0x10(%rdi)
jne 0x1b57c
popq %rbp
jmp 0x1b58f
popq %rbp
jmp 0x1b83c
movq 0xf5e67(%rip), %rdi     # 0x1113f0
popq %rbp
jmp 0x8b01
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_complex_double
void
power_complex_double(struct atom *X, struct atom *Y, struct atom *EXPO)
{
	double expo, r, theta, x, y;

	push(EXPO);
	expo = pop_double();

	push(X);
	x = pop_double();

	push(Y);
	y = pop_double();

	r = hypot(x, y);
	theta = atan2(y, x);

	r = pow(r, expo);
	theta = expo * theta;

	x = r * cos(theta);
	y = r * sin(theta);

	push_double(x);
	push_double(y);
	push(imaginaryunit);
	multiply();
	add();
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r14
pushq %rbx
subq $0x20, %rsp
movq %rsi, %rbx
movq %rdi, %r14
movq %rdx, %rdi
callq 0x8b01
callq 0x9086
movsd %xmm0, -0x28(%rbp)
movq %r14, %rdi
callq 0x8b01
callq 0x9086
movsd %xmm0, -0x18(%rbp)
movq %rbx, %rdi
callq 0x8b01
callq 0x9086
movaps %xmm0, %xmm1
movsd %xmm0, -0x30(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x72b0
movsd %xmm0, -0x20(%rbp)
movsd -0x30(%rbp), %xmm0
movsd -0x18(%rbp), %xmm1
callq 0x72a0
movsd %xmm0, -0x18(%rbp)
movsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
callq 0x70d0
movsd %xmm0, -0x20(%rbp)
movsd -0x18(%rbp), %xmm0
mulsd -0x28(%rbp), %xmm0
movsd %xmm0, -0x18(%rbp)
callq 0x7170
mulsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x28(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x72f0
mulsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x18(%rbp)
callq 0x74ec
movb $0x5, %bl
movb %bl, 0x10(%rax)
movsd -0x28(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
callq 0x74ec
movb %bl, 0x10(%rax)
movsd -0x18(%rbp), %xmm0
movsd %xmm0, (%rax)
movq %rax, %rdi
callq 0x8b01
movq 0xeab5c(%rip), %rdi     # 0x106ae0
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movl $0x2, %edi
addq $0x20, %rsp
popq %rbx
popq %r14
popq %rbp
jmp 0xa975
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_complex_rational
void
power_complex_rational(struct atom *X, struct atom *Y, struct atom *EXPO)
{
	// calculate sqrt(X^2 + Y^2) ^ (1/2 * EXPO)

	push(X);
	push(X);
	multiply();
	push(Y);
	push(Y);
	multiply();
	add();
	push_rational(1, 2);
	push(EXPO);
	multiply();
	power();

	// calculate (-1) ^ (EXPO * arctan(Y, X) / pi)

	push(Y);
	push(X);
	arctan();
	push_symbol(PI);
	divide();
	push(EXPO);
	multiply();
	EXPO = pop();
	power_minusone(EXPO);

	// result = sqrt(X^2 + Y^2) ^ (1/2 * EXPO) * (-1) ^ (EXPO * arctan(Y, X) / pi)

	multiply();
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %rbx
pushq %rax
movq %rdx, %rbx
movq %rsi, %r15
movq %rdi, %r14
callq 0x8b01
movq %r14, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movq %r15, %rdi
callq 0x8b01
movq %r15, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movl $0x2, %edi
callq 0xa975
movl $0x1, %edi
movl $0x2, %esi
callq 0xa1f4
movq %rbx, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
callq 0xa234
movq %r15, %rdi
callq 0x8b01
movq %r14, %rdi
callq 0x8b01
callq 0xd34a
movq 0x2501b(%rip), %rdi     # 0x41050
callq 0x8b01
movq 0xf53af(%rip), %rdi     # 0x1113f0
callq 0x8b01
callq 0xa234
movl $0x2, %edi
callq 0xa8bb
movq %rbx, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
callq 0x8b78
movq %rax, %rdi
callq 0x1abf4
movl $0x2, %edi
addq $0x8, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
jmp 0xa8bb
// File: /workspace/github_c_makefile_star2_O1/georgeweigt[P]eigenmath/eigenmath.c
// Function: power_complex_plus
void
power_complex_plus(struct atom *X, struct atom *Y, int n)
{
	int i;
	struct atom *PX, *PY;

	PX = X;
	PY = Y;

	for (i = 1; i < n; i++) {

		push(PX);
		push(X);
		multiply();
		push(PY);
		push(Y);
		multiply();
		subtract();

		push(PX);
		push(Y);
		multiply();
		push(PY);
		push(X);
		multiply();
		add();

		PY = pop();
		PX = pop();
	}

	// X + i*Y

	push(PX);
	push(imaginaryunit);
	push(PY);
	multiply();
	add();
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rsi, %rbx
movq %rdi, %r15
cmpl $0x2, %edx
jl 0x1c15e
movl %edx, %r14d
decl %r14d
movq %rbx, %r12
movq %r15, %r13
movq %r13, %rdi
callq 0x8b01
movq %r15, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movq %r12, %rdi
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movq 0xf5304(%rip), %rdi     # 0x1113f0
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movl $0x2, %edi
callq 0xa975
movq %r13, %rdi
callq 0x8b01
movq %rbx, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movq %r12, %rdi
callq 0x8b01
movq %r15, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movl $0x2, %edi
callq 0xa975
callq 0x8b78
movq %rax, %r12
callq 0x8b78
movq %rax, %r13
decl %r14d
jne 0x1c0b1
jmp 0x1c164
movq %r15, %r13
movq %rbx, %r12
movq %r13, %rdi
callq 0x8b01
movq 0xea96d(%rip), %rdi     # 0x106ae0
callq 0x8b01
movq %r12, %rdi
callq 0x8b01
movl $0x2, %edi
callq 0xa8bb
movl $0x2, %edi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0xa975