// File: /workspace/github_c_makefile_star2_O0/georgeweigt[P]eigenmath/eigenmath.c
// Function: power
void
power(void)
{
	int h, i, n;
	struct atom *p1, *BASE, *EXPO;

	EXPO = pop();
	BASE = pop();

	if (istensor(BASE) && istensor(EXPO)) {
		push_symbol(POWER);
		push(BASE);
		push(EXPO);
		list(3);
		return;
	}

	if (istensor(EXPO)) {
		p1 = copy_tensor(EXPO);
		n = p1->u.tensor->nelem;
		for (i = 0; i < n; i++) {
			push(BASE);
			push(p1->u.tensor->elem[i]);
			power();
			p1->u.tensor->elem[i] = pop();
		}
		push(p1);
		return;
	}

	if (istensor(BASE)) {
		p1 = copy_tensor(BASE);
		n = p1->u.tensor->nelem;
		for (i = 0; i < n; i++) {
			push(p1->u.tensor->elem[i]);
			push(EXPO);
			power();
			p1->u.tensor->elem[i] = pop();
		}
		push(p1);
		return;
	}

	if (BASE == symbol(EXP1) && isdouble(EXPO)) {
		push_double(M_E);
		BASE = pop();
	}

	if (BASE == symbol(PI) && isdouble(EXPO)) {
		push_double(M_PI);
		BASE = pop();
	}

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

	// expr^0

	if (iszero(EXPO)) {
		push_integer(1);
		return;
	}

	// 0^expr

	if (iszero(BASE)) {
		push_symbol(POWER);
		push(BASE);
		push(EXPO);
		list(3);
		return;
	}

	// 1^expr

	if (isplusone(BASE)) {
		push_integer(1);
		return;
	}

	// expr^1

	if (isplusone(EXPO)) {
		push(BASE);
		return;
	}

	// BASE is an integer?

	if (isinteger(BASE)) {
		// raise each factor in BASE to power EXPO
		// EXPO is not numerical, that case was handled by power_numbers() above
		h = tos;
		push(BASE);
		factor_factor();
		n = tos - h;
		for (i = 0; i < n; i++) {
			p1 = stack[h + i];
			if (car(p1) == symbol(POWER)) {
				push_symbol(POWER);
				push(cadr(p1)); // base
				push(caddr(p1)); // expo
				push(EXPO);
				multiply();
				list(3);
			} else {
				push_symbol(POWER);
				push(p1);
				push(EXPO);
				list(3);
			}
			stack[h + i] = pop();
		}
		if (n > 1) {
			sort_factors(n);
			list(n);
			push_symbol(MULTIPLY);
			swap();
			cons(); // prepend MULTIPLY to list
		}
		return;
	}

	// BASE is a numerical fraction?

	if (isfraction(BASE)) {
		// power numerator, power denominator
		// EXPO is not numerical, that case was handled by power_numbers() above
		push(BASE);
		numerator();
		push(EXPO);
		power();
		push(BASE);
		denominator();
		push(EXPO);
		negate();
		power();
		multiply();
		return;
	}

	// BASE = e ?

	if (BASE == symbol(EXP1)) {
		power_natural_number(EXPO);
		return;
	}

	// (a + b) ^ c

	if (car(BASE) == symbol(ADD)) {
		power_sum(BASE, EXPO);
		return;
	}

	// (a b) ^ c  -->  (a ^ c) (b ^ c)

	if (car(BASE) == symbol(MULTIPLY)) {
		h = tos;
		p1 = cdr(BASE);
		while (iscons(p1)) {
			push(car(p1));
			push(EXPO);
			power();
			p1 = cdr(p1);
		}
		multiply_factors(tos - h);
		return;
	}

	// (a ^ b) ^ c  -->  a ^ (b c)

	if (car(BASE) == symbol(POWER)) {
		push(cadr(BASE));
		push(caddr(BASE));
		push(EXPO);
		multiply_expand(); // always expand products of exponents
		power();
		return;
	}

	// none of the above

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

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x100, %rsp            # imm = 0x100
callq 0x91d0
movq %rax, -0x28(%rbp)
callq 0x91d0
movq %rax, -0x20(%rbp)
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x7, %eax
jne 0xb914
movq -0x28(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x7, %eax
jne 0xb914
movq 0x559c2(%rip), %rdi     # 0x612b0
callq 0x9110
movq -0x20(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
movl $0x3, %edi
callq 0x90d0
jmp 0xc271
movq -0x28(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x7, %eax
jne 0xb99e
movq -0x28(%rbp), %rdi
callq 0xa550
movq %rax, -0x18(%rbp)
movq -0x18(%rbp), %rax
movq (%rax), %rax
movl 0x64(%rax), %eax
movl %eax, -0xc(%rbp)
movl $0x0, -0x8(%rbp)
movl -0x8(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0xb990
movq -0x20(%rbp), %rdi
callq 0x9110
movq -0x18(%rbp), %rax
movq (%rax), %rax
movslq -0x8(%rbp), %rcx
movq 0x68(%rax,%rcx,8), %rdi
callq 0x9110
callq 0xb8b0
callq 0x91d0
movq %rax, %rdx
movq -0x18(%rbp), %rax
movq (%rax), %rax
movslq -0x8(%rbp), %rcx
movq %rdx, 0x68(%rax,%rcx,8)
movl -0x8(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x8(%rbp)
jmp 0xb942
movq -0x18(%rbp), %rdi
callq 0x9110
jmp 0xc271
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x7, %eax
jne 0xba28
movq -0x20(%rbp), %rdi
callq 0xa550
movq %rax, -0x18(%rbp)
movq -0x18(%rbp), %rax
movq (%rax), %rax
movl 0x64(%rax), %eax
movl %eax, -0xc(%rbp)
movl $0x0, -0x8(%rbp)
movl -0x8(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0xba1a
movq -0x18(%rbp), %rax
movq (%rax), %rax
movslq -0x8(%rbp), %rcx
movq 0x68(%rax,%rcx,8), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
callq 0xb8b0
callq 0x91d0
movq %rax, %rdx
movq -0x18(%rbp), %rax
movq (%rax), %rax
movslq -0x8(%rbp), %rcx
movq %rdx, 0x68(%rax,%rcx,8)
movl -0x8(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x8(%rbp)
jmp 0xb9cc
movq -0x18(%rbp), %rdi
callq 0x9110
jmp 0xc271
movq -0x20(%rbp), %rax
cmpq 0x55895(%rip), %rax     # 0x612c8
jne 0xba58
movq -0x28(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x5, %eax
jne 0xba58
movsd 0x3964e(%rip), %xmm0    # 0x45098
callq 0xe5e0
callq 0x91d0
movq %rax, -0x20(%rbp)
movq -0x20(%rbp), %rax
cmpq 0x535ed(%rip), %rax     # 0x5f050
jne 0xba88
movq -0x28(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x5, %eax
jne 0xba88
movsd 0x39626(%rip), %xmm0    # 0x450a0
callq 0xe5e0
callq 0x91d0
movq %rax, -0x20(%rbp)
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x4, %eax
je 0xbaa2
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x5, %eax
jne 0xbace
movq -0x28(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x4, %eax
je 0xbabc
movq -0x28(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x5, %eax
jne 0xbace
movq -0x20(%rbp), %rdi
movq -0x28(%rbp), %rsi
callq 0x297c0
jmp 0xc271
movq -0x28(%rbp), %rdi
callq 0xcdc0
cmpl $0x0, %eax
je 0xbaeb
movl $0x1, %edi
callq 0xca50
jmp 0xc271
movq -0x20(%rbp), %rdi
callq 0xcdc0
cmpl $0x0, %eax
je 0xbb26
movq 0x557b0(%rip), %rdi     # 0x612b0
callq 0x9110
movq -0x20(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
movl $0x3, %edi
callq 0x90d0
jmp 0xc271
movq -0x20(%rbp), %rdi
callq 0xd4c0
cmpl $0x0, %eax
je 0xbb43
movl $0x1, %edi
callq 0xca50
jmp 0xc271
movq -0x28(%rbp), %rdi
callq 0xd4c0
cmpl $0x0, %eax
je 0xbb5f
movq -0x20(%rbp), %rdi
callq 0x9110
jmp 0xc271
movq -0x20(%rbp), %rdi
callq 0x9bc0
cmpl $0x0, %eax
je 0xbe2a
movl 0x4c741(%rip), %eax     # 0x582b8
movl %eax, -0x4(%rbp)
movq -0x20(%rbp), %rdi
callq 0x9110
callq 0x247f0
movl 0x4c72a(%rip), %eax     # 0x582b8
subl -0x4(%rbp), %eax
movl %eax, -0xc(%rbp)
movl $0x0, -0x8(%rbp)
movl -0x8(%rbp), %eax
cmpl -0xc(%rbp), %eax
jge 0xbdf9
movl -0x4(%rbp), %eax
addl -0x8(%rbp), %eax
movslq %eax, %rcx
leaq 0x55a09(%rip), %rax     # 0x615c0
movq (%rax,%rcx,8), %rax
movq %rax, -0x18(%rbp)
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbbd9
movq -0x18(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x30(%rbp)
jmp 0xbbe4
movq 0x52e20(%rip), %rax     # 0x5ea00
movq %rax, -0x30(%rbp)
movq -0x30(%rbp), %rax
cmpq 0x556c1(%rip), %rax     # 0x612b0
jne 0xbda7
movq 0x556b4(%rip), %rdi     # 0x612b0
callq 0x9110
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbc1c
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x38(%rbp)
jmp 0xbc27
movq 0x52ddd(%rip), %rax     # 0x5ea00
movq %rax, -0x38(%rbp)
movq -0x38(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbc67
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbc4f
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x40(%rbp)
jmp 0xbc5a
movq 0x52daa(%rip), %rax     # 0x5ea00
movq %rax, -0x40(%rbp)
movq -0x40(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x48(%rbp)
jmp 0xbc72
movq 0x52d92(%rip), %rax     # 0x5ea00
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rdi
callq 0x9110
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbc96
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x50(%rbp)
jmp 0xbca1
movq 0x52d63(%rip), %rax     # 0x5ea00
movq %rax, -0x50(%rbp)
movq -0x50(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbce2
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbcc9
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x58(%rbp)
jmp 0xbcd4
movq 0x52d30(%rip), %rax     # 0x5ea00
movq %rax, -0x58(%rbp)
movq -0x58(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x60(%rbp)
jmp 0xbced
movq 0x52d17(%rip), %rax     # 0x5ea00
movq %rax, -0x60(%rbp)
movq -0x60(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbd79
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbd15
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x68(%rbp)
jmp 0xbd20
movq 0x52ce4(%rip), %rax     # 0x5ea00
movq %rax, -0x68(%rbp)
movq -0x68(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbd61
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbd48
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x70(%rbp)
jmp 0xbd53
movq 0x52cb1(%rip), %rax     # 0x5ea00
movq %rax, -0x70(%rbp)
movq -0x70(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x78(%rbp)
jmp 0xbd6c
movq 0x52c98(%rip), %rax     # 0x5ea00
movq %rax, -0x78(%rbp)
movq -0x78(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x80(%rbp)
jmp 0xbd84
movq 0x52c80(%rip), %rax     # 0x5ea00
movq %rax, -0x80(%rbp)
movq -0x80(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
callq 0xc280
movl $0x3, %edi
callq 0x90d0
jmp 0xbdcf
movq 0x55502(%rip), %rdi     # 0x612b0
callq 0x9110
movq -0x18(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
movl $0x3, %edi
callq 0x90d0
callq 0x91d0
movq %rax, %rdx
movl -0x4(%rbp), %eax
addl -0x8(%rbp), %eax
movslq %eax, %rcx
leaq 0x557d9(%rip), %rax     # 0x615c0
movq %rdx, (%rax,%rcx,8)
movl -0x8(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x8(%rbp)
jmp 0xbb9b
cmpl $0x1, -0xc(%rbp)
jle 0xbe25
movl -0xc(%rbp), %edi
callq 0x26b20
movl -0xc(%rbp), %edi
callq 0x90d0
movq 0x55492(%rip), %rdi     # 0x612a8
callq 0x9110
callq 0xcb10
callq 0x9180
jmp 0xc271
movq -0x20(%rbp), %rdi
callq 0x29d30
cmpl $0x0, %eax
je 0xbe7f
movq -0x20(%rbp), %rdi
callq 0x9110
callq 0x117b0
movq -0x28(%rbp), %rdi
callq 0x9110
callq 0xb8b0
movq -0x20(%rbp), %rdi
callq 0x9110
callq 0x117d0
movq -0x28(%rbp), %rdi
callq 0x9110
callq 0xb4b0
callq 0xb8b0
callq 0xc280
jmp 0xc271
movq -0x20(%rbp), %rax
cmpq 0x5543e(%rip), %rax     # 0x612c8
jne 0xbe9a
movq -0x28(%rbp), %rdi
callq 0x29d80
jmp 0xc271
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbeb7
movq -0x20(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x88(%rbp)
jmp 0xbec5
movq 0x52b42(%rip), %rax     # 0x5ea00
movq %rax, -0x88(%rbp)
movq -0x88(%rbp), %rax
cmpq 0x553cd(%rip), %rax     # 0x612a0
jne 0xbee7
movq -0x20(%rbp), %rdi
movq -0x28(%rbp), %rsi
callq 0x2a5f0
jmp 0xc271
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbf04
movq -0x20(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x90(%rbp)
jmp 0xbf12
movq 0x52af5(%rip), %rax     # 0x5ea00
movq %rax, -0x90(%rbp)
movq -0x90(%rbp), %rax
cmpq 0x55388(%rip), %rax     # 0x612a8
jne 0xc00b
movl 0x4c38c(%rip), %eax     # 0x582b8
movl %eax, -0x4(%rbp)
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbf4d
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0x98(%rbp)
jmp 0xbf5b
movq 0x52aac(%rip), %rax     # 0x5ea00
movq %rax, -0x98(%rbp)
movq -0x98(%rbp), %rax
movq %rax, -0x18(%rbp)
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbff8
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbf94
movq -0x18(%rbp), %rax
movq (%rax), %rax
movq %rax, -0xa0(%rbp)
jmp 0xbfa2
movq 0x52a65(%rip), %rax     # 0x5ea00
movq %rax, -0xa0(%rbp)
movq -0xa0(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
callq 0xb8b0
movq -0x18(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xbfda
movq -0x18(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xa8(%rbp)
jmp 0xbfe8
movq 0x52a1f(%rip), %rax     # 0x5ea00
movq %rax, -0xa8(%rbp)
movq -0xa8(%rbp), %rax
movq %rax, -0x18(%rbp)
jmp 0xbf66
movl 0x4c2ba(%rip), %edi     # 0x582b8
subl -0x4(%rbp), %edi
callq 0xc3d0
jmp 0xc271
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc028
movq -0x20(%rbp), %rax
movq (%rax), %rax
movq %rax, -0xb0(%rbp)
jmp 0xc036
movq 0x529d1(%rip), %rax     # 0x5ea00
movq %rax, -0xb0(%rbp)
movq -0xb0(%rbp), %rax
cmpq 0x5526c(%rip), %rax     # 0x612b0
jne 0xc249
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc068
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xb8(%rbp)
jmp 0xc076
movq 0x52991(%rip), %rax     # 0x5ea00
movq %rax, -0xb8(%rbp)
movq -0xb8(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc0c5
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc0a4
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xc0(%rbp)
jmp 0xc0b2
movq 0x52955(%rip), %rax     # 0x5ea00
movq %rax, -0xc0(%rbp)
movq -0xc0(%rbp), %rax
movq (%rax), %rax
movq %rax, -0xc8(%rbp)
jmp 0xc0d3
movq 0x52934(%rip), %rax     # 0x5ea00
movq %rax, -0xc8(%rbp)
movq -0xc8(%rbp), %rdi
callq 0x9110
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc0fd
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xd0(%rbp)
jmp 0xc10b
movq 0x528fc(%rip), %rax     # 0x5ea00
movq %rax, -0xd0(%rbp)
movq -0xd0(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc15b
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc139
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xd8(%rbp)
jmp 0xc147
movq 0x528c0(%rip), %rax     # 0x5ea00
movq %rax, -0xd8(%rbp)
movq -0xd8(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xe0(%rbp)
jmp 0xc169
movq 0x5289e(%rip), %rax     # 0x5ea00
movq %rax, -0xe0(%rbp)
movq -0xe0(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc21a
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc19b
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xe8(%rbp)
jmp 0xc1a9
movq 0x5285e(%rip), %rax     # 0x5ea00
movq %rax, -0xe8(%rbp)
movq -0xe8(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc1f9
movq -0x20(%rbp), %rax
movzbl 0x10(%rax), %eax
cmpl $0x1, %eax
jne 0xc1d7
movq -0x20(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xf0(%rbp)
jmp 0xc1e5
movq 0x52822(%rip), %rax     # 0x5ea00
movq %rax, -0xf0(%rbp)
movq -0xf0(%rbp), %rax
movq 0x8(%rax), %rax
movq %rax, -0xf8(%rbp)
jmp 0xc207
movq 0x52800(%rip), %rax     # 0x5ea00
movq %rax, -0xf8(%rbp)
movq -0xf8(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x100(%rbp)
jmp 0xc228
movq 0x527df(%rip), %rax     # 0x5ea00
movq %rax, -0x100(%rbp)
movq -0x100(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
callq 0x27fb0
callq 0xb8b0
jmp 0xc271
movq 0x55060(%rip), %rdi     # 0x612b0
callq 0x9110
movq -0x20(%rbp), %rdi
callq 0x9110
movq -0x28(%rbp), %rdi
callq 0x9110
movl $0x3, %edi
callq 0x90d0
addq $0x100, %rsp            # imm = 0x100
popq %rbp
retq
nopw (%rax,%rax)