// File: /workspace/github_c_makefile_star2_O3/georgeweigt[P]eigenmath/eigenmath.c
// Function: nfindroot
void
nfindroot(double cr[], double ci[], int n, double *par, double *pai)
{
	int i, j;
	double d;
	double ar, br, dfr, dxr, far, fbr, xr, yr;
	double ai, bi, dfi, dxi, fai, fbi, xi, yi;

	// if const term is zero then root is zero

	// note: use exact zero, not "close to zero"

	// term will be exactly zero from coeffs(), no need for arbitrary cutoff

	if (cr[0] == 0.0 && ci[0] == 0.0) {
		*par = 0.0;
		*pai = 0.0;
		return;
	}

	// secant method

	for (i = 0; i < 100; i++) {

		ar = urandom();
		ai = urandom();

		fata(cr, ci, n, ar, ai, &far, &fai);

		br = ar;
		bi = ai;

		fbr = far;
		fbi = fai;

		ar = urandom();
		ai = urandom();

		for (j = 0; j < 1000; j++) {

			fata(cr, ci, n, ar, ai, &far, &fai);

			if (zabs(far, fai) < EPSILON) {
				*par = ar;
				*pai = ai;
				return;
			}

			if (zabs(far, fai) < zabs(fbr, fbi)) {

				xr = ar;
				xi = ai;

				ar = br;
				ai = bi;

				br = xr;
				bi = xi;

				xr = far;
				xi = fai;

				far = fbr;
				fai = fbi;

				fbr = xr;
				fbi = xi;
			}

			// dx = b - a

			dxr = br - ar;
			dxi = bi - ai;

			// df = fb - fa

			dfr = fbr - far;
			dfi = fbi - fai;

			// y = dx / df

			d = dfr * dfr + dfi * dfi;

			if (d == 0.0)
				break;

			yr = (dxr * dfr + dxi * dfi) / d;
			yi = (dxi * dfr - dxr * dfi) / d;

			// a = b - y * fb

			ar = br - (yr * fbr - yi * fbi);
			ai = bi - (yr * fbi + yi * fbr);
		}
	}

	stopf("nroots: convergence error");
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movl %edx, %r15d
movq %rsi, %r12
movq %rdi, %r13
movsd (%rdi), %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0x1b63f
jp 0x1b63f
movsd (%r12), %xmm0
ucomisd %xmm1, %xmm0
jne 0x1b63f
jp 0x1b63f
xorpd %xmm2, %xmm2
jmp 0x1b967
movq %rcx, -0x40(%rbp)
movq %r8, -0x48(%rbp)
leal -0x1(%r15), %eax
movslq %eax, %rbx
leal -0x2(%r15), %eax
movq %rax, -0x38(%rbp)
movl $0x0, -0x2c(%rbp)
callq 0x7380
movl %eax, %r14d
callq 0x7380
xorps %xmm5, %xmm5
cvtsi2sd %r14d, %xmm5
movsd 0xeae6(%rip), %xmm3     # 0x2a160
divsd %xmm3, %xmm5
movsd 0xeaea(%rip), %xmm1     # 0x2a170
mulsd %xmm1, %xmm5
movsd 0xeae6(%rip), %xmm0     # 0x2a178
addsd %xmm0, %xmm5
xorps %xmm2, %xmm2
cvtsi2sd %eax, %xmm2
divsd %xmm3, %xmm2
mulsd %xmm1, %xmm2
addsd %xmm0, %xmm2
movsd (%r13,%rbx,8), %xmm0
movsd (%r12,%rbx,8), %xmm3
cmpl $0x2, %r15d
jl 0x1b6fc
movq -0x38(%rbp), %r14
movq %r14, %rax
movapd %xmm2, %xmm1
mulsd %xmm3, %xmm1
movapd %xmm5, %xmm4
mulsd %xmm0, %xmm4
subsd %xmm1, %xmm4
mulsd %xmm2, %xmm0
mulsd %xmm5, %xmm3
addsd %xmm0, %xmm3
addsd (%r13,%rax,8), %xmm4
addsd (%r12,%rax,8), %xmm3
movapd %xmm4, %xmm0
addq $-0x1, %rax
jb 0x1b6c3
jmp 0x1b704
movapd %xmm0, %xmm4
movq -0x38(%rbp), %r14
movsd %xmm4, -0x60(%rbp)
movsd %xmm3, -0x68(%rbp)
movsd %xmm2, -0x70(%rbp)
movsd %xmm5, -0x58(%rbp)
callq 0x7380
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
movsd 0xea34(%rip), %xmm0     # 0x2a160
divsd %xmm0, %xmm1
movsd 0xea38(%rip), %xmm0     # 0x2a170
mulsd %xmm0, %xmm1
movsd 0xea34(%rip), %xmm0     # 0x2a178
addsd %xmm0, %xmm1
movsd %xmm1, -0x50(%rbp)
callq 0x7380
movsd -0x50(%rbp), %xmm15
movsd -0x58(%rbp), %xmm14
movsd -0x60(%rbp), %xmm13
movsd -0x68(%rbp), %xmm12
movsd -0x70(%rbp), %xmm11
xorps %xmm2, %xmm2
cvtsi2sd %eax, %xmm2
divsd 0xe9e1(%rip), %xmm2     # 0x2a160
mulsd 0xe9e9(%rip), %xmm2     # 0x2a170
addsd 0xe9e9(%rip), %xmm2     # 0x2a178
movsd (%r13,%rbx,8), %xmm0
movsd (%r12,%rbx,8), %xmm1
xorl %eax, %eax
xorpd %xmm9, %xmm9
movsd 0xe9d4(%rip), %xmm10    # 0x2a180
movapd %xmm0, %xmm3
movapd %xmm1, %xmm4
cmpl $0x2, %r15d
jl 0x1b7fe
movq %r14, %rcx
movapd %xmm1, %xmm4
movapd %xmm0, %xmm5
movapd %xmm2, %xmm6
mulsd %xmm4, %xmm6
movapd %xmm15, %xmm3
mulsd %xmm5, %xmm3
subsd %xmm6, %xmm3
mulsd %xmm2, %xmm5
mulsd %xmm15, %xmm4
addsd %xmm5, %xmm4
addsd (%r13,%rcx,8), %xmm3
addsd (%r12,%rcx,8), %xmm4
movapd %xmm3, %xmm5
addq $-0x1, %rcx
jb 0x1b7c5
movapd %xmm4, %xmm5
mulsd %xmm4, %xmm5
movapd %xmm3, %xmm6
mulsd %xmm3, %xmm6
addsd %xmm5, %xmm6
xorps %xmm5, %xmm5
sqrtsd %xmm6, %xmm5
ucomisd %xmm5, %xmm10
ja 0x1b95a
movapd %xmm12, %xmm6
mulsd %xmm12, %xmm6
movapd %xmm13, %xmm7
mulsd %xmm13, %xmm7
addsd %xmm6, %xmm7
xorps %xmm6, %xmm6
sqrtsd %xmm7, %xmm6
ucomisd %xmm5, %xmm6
jbe 0x1b86e
movapd %xmm13, %xmm7
movapd %xmm12, %xmm8
movapd %xmm3, %xmm13
movapd %xmm11, %xmm5
movapd %xmm2, %xmm11
movapd %xmm15, %xmm3
movapd %xmm4, %xmm12
jmp 0x1b885
movapd %xmm3, %xmm7
movapd %xmm4, %xmm8
movapd %xmm2, %xmm5
movapd %xmm14, %xmm3
movapd %xmm15, %xmm14
movapd %xmm13, %xmm6
subsd %xmm7, %xmm6
movapd %xmm12, %xmm7
subsd %xmm8, %xmm7
movapd %xmm7, %xmm2
mulsd %xmm7, %xmm2
movapd %xmm6, %xmm8
mulsd %xmm6, %xmm8
addsd %xmm2, %xmm8
ucomisd %xmm9, %xmm8
jne 0x1b8bc
jnp 0x1b93d
movapd %xmm11, %xmm2
subsd %xmm5, %xmm2
movapd %xmm3, %xmm4
subsd %xmm14, %xmm4
movapd %xmm2, %xmm5
mulsd %xmm7, %xmm5
mulsd %xmm4, %xmm7
mulsd %xmm6, %xmm4
addsd %xmm5, %xmm4
divsd %xmm8, %xmm4
mulsd %xmm6, %xmm2
subsd %xmm7, %xmm2
divsd %xmm8, %xmm2
movapd %xmm12, %xmm5
mulsd %xmm2, %xmm5
movapd %xmm4, %xmm6
mulsd %xmm13, %xmm6
subsd %xmm5, %xmm6
movapd %xmm3, %xmm15
subsd %xmm6, %xmm15
mulsd %xmm13, %xmm2
mulsd %xmm12, %xmm4
addsd %xmm2, %xmm4
movapd %xmm11, %xmm2
subsd %xmm4, %xmm2
incl %eax
movapd %xmm3, %xmm14
cmpl $0x3e8, %eax            # imm = 0x3E8
jne 0x1b7ac
movl -0x2c(%rbp), %eax
incl %eax
movl %eax, -0x2c(%rbp)
cmpl $0x64, %eax
jne 0x1b65d
leaq 0x12f7d(%rip), %rdi     # 0x2e8d2
callq 0x780b
movapd %xmm15, %xmm1
movq -0x48(%rbp), %r8
movq -0x40(%rbp), %rcx
movsd %xmm1, (%rcx)
movsd %xmm2, (%r8)
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq