// File: /workspace/github_c_makefile_star2_O0/BernardTatin[P]s7sh/s7.c
// Function: g_sqrt
static s7_pointer g_sqrt(s7_scheme *sc, s7_pointer args)
    {
#define H_sqrt "(sqrt z) returns the square root of z"
#define Q_sqrt sc->pl_nn

        s7_pointer n;
        s7_double sqx;

        n = car(args);
        switch (type(n))
        {
            case T_INTEGER:
                if (integer(n) >= 0)
                {
                    s7_int ix;
                    sqx = sqrt((s7_double)integer(n));
                    ix = (s7_int)sqx;
                    if ((ix * ix) == integer(n))
                        return(make_integer(sc, ix));
                    return(make_real(sc, sqx));
                    /* Mark Weaver notes that
                     *     (zero? (- (sqrt 9007199136250226) 94906265.0)) -> #t
                     * but (* 94906265 94906265) -> 9007199136250225 -- oops
                     * at least we return a real here, not an incorrect integer and
                     *     (sqrt 9007199136250225) -> 94906265
                     */
                }
                sqx = (s7_double)integer(n); /* we're trying to protect against (sqrt -9223372036854775808) where we can't negate the integer argument */
                return(s7_make_complex(sc, 0.0, sqrt((s7_double)(-sqx))));

            case T_RATIO:
                sqx = (s7_double)fraction(n);
                if (sqx > 0.0) /* else it's complex, so it can't be a ratio */
                {
                    s7_int nm = 0, dn = 1;
                    if (c_rationalize(sqx, 1.0e-16, &nm, &dn)) /* 1e-16 so that (sqrt 1/1099511627776) returns 1/1048576 */
                    {
#if HAVE_OVERFLOW_CHECKS
                        s7_int nm2, dn2;
                        if ((multiply_overflow(nm, nm, &nm2)) ||
                                (multiply_overflow(dn, dn, &dn2)))
                            return(make_real(sc, sqrt(sqx)));
                        if ((nm2 == numerator(n)) &&
                                (dn2 == denominator(n)))
                            return(s7_make_ratio(sc, nm, dn));
#else
                        if ((nm * nm == numerator(n)) &&
                                (dn * dn == denominator(n)))
                            return(s7_make_ratio(sc, nm, dn));
#endif
                    }
                    return(make_real(sc, sqrt(sqx)));
                }
                return(s7_make_complex(sc, 0.0, sqrt(-sqx)));

            case T_REAL:
                if (is_NaN(real(n)))
                    return(real_NaN);
                if (real(n) >= 0.0)
                    return(make_real(sc, sqrt(real(n))));
                return(s7_make_complex(sc, 0.0, sqrt(-real(n))));

            case T_COMPLEX:
                /* (* inf.0 (sqrt -1)) -> -nan+infi, but (sqrt -inf.0) -> 0+infi */
#if HAVE_COMPLEX_NUMBERS
                return(s7_from_c_complex(sc, csqrt(as_c_complex(n))));
#else
                return(out_of_range(sc, sc->sqrt_symbol, small_int(1), n, no_complex_numbers_string));
#endif

            default:
                return(method_or_bust_with_type_one_arg(sc, n, sc->sqrt_symbol, args, a_number_string));
        }
    }

// Assembly:
subq $0x138, %rsp            # imm = 0x138
movq %rdi, 0x128(%rsp)
movq %rsi, 0x120(%rsp)
movq 0x120(%rsp), %rax
movq 0x8(%rax), %rax
movq %rax, 0x118(%rsp)
movq 0x118(%rsp), %rax
movzbl (%rax), %eax
addl $-0xb, %eax
movl %eax, %ecx
movq %rcx, 0x28(%rsp)
subl $0x3, %eax
ja 0x769d7
movq 0x28(%rsp), %rax
leaq 0x1687bb(%rip), %rcx    # 0x1dea60
movslq (%rcx,%rax,4), %rax
addq %rcx, %rax
jmpq *%rax
movq 0x118(%rsp), %rax
cmpq $0x0, 0x8(%rax)
jl 0x764bf
movq 0x118(%rsp), %rax
cvtsi2sdq 0x8(%rax), %xmm0
callq 0x9690
movsd %xmm0, 0x110(%rsp)
cvttsd2si 0x110(%rsp), %rax
movq %rax, 0x108(%rsp)
movq 0x108(%rsp), %rax
imulq 0x108(%rsp), %rax
movq 0x118(%rsp), %rcx
cmpq 0x8(%rcx), %rax
jne 0x7640d
movq 0x108(%rsp), %rax
movq %rax, 0x100(%rsp)
movq 0x100(%rsp), %rax
andq $-0x800, %rax           # imm = 0xF800
cmpq $0x0, %rax
jne 0x76353
movq 0x100(%rsp), %rcx
leaq 0x1a013b(%rip), %rax    # 0x216480
movq (%rax,%rcx,8), %rax
movq %rax, 0x20(%rsp)
jmp 0x763eb
jmp 0x76355
movq 0x128(%rsp), %rax
movq 0x90(%rax), %rax
movq 0x128(%rsp), %rcx
cmpq 0x98(%rcx), %rax
ja 0x76382
movq 0x128(%rsp), %rdi
callq 0xc670
movq 0x128(%rsp), %rcx
movq 0x90(%rcx), %rax
movq %rax, %rdx
addq $-0x8, %rdx
movq %rdx, 0x90(%rcx)
movq -0x8(%rax), %rax
movq %rax, 0xf0(%rsp)
movq 0xf0(%rsp), %rax
movq $0xb, (%rax)
movq 0x100(%rsp), %rcx
movq 0xf0(%rsp), %rax
movq %rcx, 0x8(%rax)
movq 0xf0(%rsp), %rax
movq %rax, 0xe8(%rsp)
movq 0xe8(%rsp), %rax
movq %rax, 0x20(%rsp)
movq 0x20(%rsp), %rax
movq %rax, 0xf8(%rsp)
movq 0xf8(%rsp), %rax
movq %rax, 0x130(%rsp)
jmp 0x76a12
movsd 0x110(%rsp), %xmm0
movsd %xmm0, 0xd8(%rsp)
movq 0x128(%rsp), %rax
movq 0x90(%rax), %rax
movq 0x128(%rsp), %rcx
cmpq 0x98(%rcx), %rax
ja 0x7644c
movq 0x128(%rsp), %rdi
callq 0xc670
movq 0x128(%rsp), %rcx
movq 0x90(%rcx), %rax
movq %rax, %rdx
addq $-0x8, %rdx
movq %rdx, 0x90(%rcx)
movq -0x8(%rax), %rax
movq %rax, 0xe0(%rsp)
movq 0xe0(%rsp), %rax
movq $0xd, (%rax)
movsd 0xd8(%rsp), %xmm0
movq 0xe0(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movq 0xe0(%rsp), %rax
movq %rax, 0xd0(%rsp)
movq 0xd0(%rsp), %rax
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x118(%rsp), %rax
cvtsi2sdq 0x8(%rax), %xmm0
movsd %xmm0, 0x110(%rsp)
movq 0x128(%rsp), %rax
movq %rax, 0x18(%rsp)
movsd 0x110(%rsp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
callq 0x9690
movq 0x18(%rsp), %rdi
movaps %xmm0, %xmm1
xorps %xmm0, %xmm0
callq 0x122a0
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x118(%rsp), %rax
movq 0x8(%rax), %rcx
movq 0x10(%rax), %rax
movq %rcx, 0x40(%rsp)
fildll 0x40(%rsp)
movq %rax, 0x38(%rsp)
fildll 0x38(%rsp)
fdivrp %st, %st(1)
fstpl 0x30(%rsp)
movsd 0x30(%rsp), %xmm0
movsd %xmm0, 0x110(%rsp)
movsd 0x110(%rsp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jbe 0x767c4
movq $0x0, 0xc8(%rsp)
movq $0x1, 0xc0(%rsp)
movsd 0x110(%rsp), %xmm0
movsd 0x16b47d(%rip), %xmm1   # 0x1e1a18
leaq 0xc8(%rsp), %rdi
leaq 0xc0(%rsp), %rsi
callq 0x114e0
testb $0x1, %al
jne 0x765b9
jmp 0x7670d
movq 0xc8(%rsp), %rcx
imulq 0xc8(%rsp), %rcx
seto %al
movq %rcx, 0xb8(%rsp)
testb $0x1, %al
jne 0x765fe
movq 0xc0(%rsp), %rcx
imulq 0xc0(%rsp), %rcx
seto %al
movq %rcx, 0xb0(%rsp)
testb $0x1, %al
jne 0x765fe
jmp 0x766b5
movsd 0x110(%rsp), %xmm0
callq 0x9690
movsd %xmm0, 0xa0(%rsp)
movq 0x128(%rsp), %rax
movq 0x90(%rax), %rax
movq 0x128(%rsp), %rcx
cmpq 0x98(%rcx), %rax
ja 0x76642
movq 0x128(%rsp), %rdi
callq 0xc670
movq 0x128(%rsp), %rcx
movq 0x90(%rcx), %rax
movq %rax, %rdx
addq $-0x8, %rdx
movq %rdx, 0x90(%rcx)
movq -0x8(%rax), %rax
movq %rax, 0xa8(%rsp)
movq 0xa8(%rsp), %rax
movq $0xd, (%rax)
movsd 0xa0(%rsp), %xmm0
movq 0xa8(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movq 0xa8(%rsp), %rax
movq %rax, 0x98(%rsp)
movq 0x98(%rsp), %rax
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0xb8(%rsp), %rax
movq 0x118(%rsp), %rcx
cmpq 0x8(%rcx), %rax
jne 0x7670b
movq 0xb0(%rsp), %rax
movq 0x118(%rsp), %rcx
cmpq 0x10(%rcx), %rax
jne 0x7670b
movq 0x128(%rsp), %rdi
movq 0xc8(%rsp), %rsi
movq 0xc0(%rsp), %rdx
callq 0x11a80
movq %rax, 0x130(%rsp)
jmp 0x76a12
jmp 0x7670d
movsd 0x110(%rsp), %xmm0
callq 0x9690
movsd %xmm0, 0x88(%rsp)
movq 0x128(%rsp), %rax
movq 0x90(%rax), %rax
movq 0x128(%rsp), %rcx
cmpq 0x98(%rcx), %rax
ja 0x76751
movq 0x128(%rsp), %rdi
callq 0xc670
movq 0x128(%rsp), %rcx
movq 0x90(%rcx), %rax
movq %rax, %rdx
addq $-0x8, %rdx
movq %rdx, 0x90(%rcx)
movq -0x8(%rax), %rax
movq %rax, 0x90(%rsp)
movq 0x90(%rsp), %rax
movq $0xd, (%rax)
movsd 0x88(%rsp), %xmm0
movq 0x90(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movq 0x90(%rsp), %rax
movq %rax, 0x80(%rsp)
movq 0x80(%rsp), %rax
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x128(%rsp), %rax
movq %rax, 0x10(%rsp)
movsd 0x110(%rsp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
callq 0x9690
movq 0x10(%rsp), %rdi
movaps %xmm0, %xmm1
xorps %xmm0, %xmm0
callq 0x122a0
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x118(%rsp), %rax
movsd 0x8(%rax), %xmm0
callq 0x13860
testb $0x1, %al
jne 0x7682b
jmp 0x7683f
movq 0x1a719e(%rip), %rax    # 0x21d9d0
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x118(%rsp), %rax
movsd 0x8(%rax), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x768fc
movq 0x118(%rsp), %rax
movsd 0x8(%rax), %xmm0
callq 0x9690
movsd %xmm0, 0x70(%rsp)
movq 0x128(%rsp), %rax
movq 0x90(%rax), %rax
movq 0x128(%rsp), %rcx
cmpq 0x98(%rcx), %rax
ja 0x7689e
movq 0x128(%rsp), %rdi
callq 0xc670
movq 0x128(%rsp), %rcx
movq 0x90(%rcx), %rax
movq %rax, %rdx
addq $-0x8, %rdx
movq %rdx, 0x90(%rcx)
movq -0x8(%rax), %rax
movq %rax, 0x78(%rsp)
movq 0x78(%rsp), %rax
movq $0xd, (%rax)
movsd 0x70(%rsp), %xmm0
movq 0x78(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movq 0x78(%rsp), %rax
movq %rax, 0x68(%rsp)
movq 0x68(%rsp), %rax
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x128(%rsp), %rax
movq %rax, 0x8(%rsp)
movq 0x118(%rsp), %rax
movsd 0x8(%rax), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
callq 0x9690
movq 0x8(%rsp), %rdi
movaps %xmm0, %xmm1
xorps %xmm0, %xmm0
callq 0x122a0
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x128(%rsp), %rax
movq %rax, (%rsp)
movq 0x118(%rsp), %rax
movsd 0x8(%rax), %xmm1
movq 0x118(%rsp), %rax
movsd 0x10(%rax), %xmm0
xorps %xmm3, %xmm3
movaps %xmm0, %xmm2
mulsd %xmm3, %xmm2
movsd 0x16afe1(%rip), %xmm3   # 0x1e1968
mulsd %xmm3, %xmm0
addsd %xmm2, %xmm1
movsd %xmm1, 0x58(%rsp)
movsd %xmm0, 0x60(%rsp)
movsd 0x58(%rsp), %xmm0
movsd 0x60(%rsp), %xmm1
callq 0x9560
movq (%rsp), %rdi
movsd %xmm0, 0x48(%rsp)
movsd %xmm1, 0x50(%rsp)
movsd 0x48(%rsp), %xmm0
movsd 0x50(%rsp), %xmm1
callq 0x1bb560
movq %rax, 0x130(%rsp)
jmp 0x76a12
movq 0x128(%rsp), %rdi
movq 0x118(%rsp), %rsi
movq 0x128(%rsp), %rax
movq 0x14f0(%rax), %rdx
movq 0x120(%rsp), %rcx
movq 0x1ae4eb(%rip), %r8     # 0x224ef0
callq 0x1ba7e0
movq %rax, 0x130(%rsp)
movq 0x130(%rsp), %rax
addq $0x138, %rsp            # imm = 0x138
retq
nopw %cs:(%rax,%rax)