feature. See also
. The project being documented here (as the example) is the Zig library itself.
hexagon.__hexagon_sqrtdf2
fn __hexagon_sqrtdf2() align(32) callconv(.naked) noreturn
File
Code
fn __hexagon_sqrtdf2() align(32) callconv(.naked) noreturn {
asm volatile (
\\ {
\\ r15:14 = extractu(r1:0,#23 +1,#52 -23)
\\ r28 = extractu(r1,#11,#52 -32)
\\ r5:4 = combine(##0x3f000004,#1)
\\ }
\\ {
\\ p2 = dfclass(r1:0,#0x02)
\\ p2 = cmp.gt(r1,#-1)
\\ if (!p2.new) jump:nt .Lsqrt_abnormal
\\ r9 = or(r5,r14)
\\ }
\\
\\ .Ldenormal_restart:
\\ {
\\ r11:10 = r1:0
\\ r7,p0 = sfinvsqrta(r9)
\\ r5 = and(r5,#-16)
\\ r3:2 = #0
\\ }
\\ {
\\ r3 += sfmpy(r7,r9):lib
\\ r2 += sfmpy(r7,r5):lib
\\ r6 = r5
\\
\\
\\ r9 = and(r28,#1)
\\ }
\\ {
\\ r6 -= sfmpy(r3,r2):lib
\\ r11 = insert(r4,#11 +1,#52 -32)
\\ p1 = cmp.gtu(r9,#0)
\\ }
\\ {
\\ r3 += sfmpy(r3,r6):lib
\\ r2 += sfmpy(r2,r6):lib
\\ r6 = r5
\\ r9 = mux(p1,#8,#9)
\\ }
\\ {
\\ r6 -= sfmpy(r3,r2):lib
\\ r11:10 = asl(r11:10,r9)
\\ r9 = mux(p1,#3,#2)
\\ }
\\ {
\\ r2 += sfmpy(r2,r6):lib
\\
\\ r15:14 = asl(r11:10,r9)
\\ }
\\ {
\\ r2 = and(r2,##0x007fffff)
\\ }
\\ {
\\ r2 = add(r2,##0x00800000 - 3)
\\ r9 = mux(p1,#7,#8)
\\ }
\\ {
\\ r8 = asl(r2,r9)
\\ r9 = mux(p1,#15-(1+1),#15-(1+0))
\\ }
\\ {
\\ r13:12 = mpyu(r8,r15)
\\ }
\\ {
\\ r1:0 = asl(r11:10,#15)
\\ r15:14 = mpyu(r13,r13)
\\ p1 = cmp.eq(r0,r0)
\\ }
\\ {
\\ r1:0 -= asl(r15:14,#15)
\\ r15:14 = mpyu(r13,r12)
\\ p2 = cmp.eq(r0,r0)
\\ }
\\ {
\\ r1:0 -= lsr(r15:14,#16)
\\ p3 = cmp.eq(r0,r0)
\\ }
\\ {
\\ r1:0 = mpyu(r1,r8)
\\ }
\\ {
\\ r13:12 += lsr(r1:0,r9)
\\ r9 = add(r9,#16)
\\ r1:0 = asl(r11:10,#31)
\\ }
\\
\\ {
\\ r15:14 = mpyu(r13,r13)
\\ r1:0 -= mpyu(r13,r12)
\\ }
\\ {
\\ r1:0 -= asl(r15:14,#31)
\\ r15:14 = mpyu(r12,r12)
\\ }
\\ {
\\ r1:0 -= lsr(r15:14,#33)
\\ }
\\ {
\\ r1:0 = mpyu(r1,r8)
\\ }
\\ {
\\ r13:12 += lsr(r1:0,r9)
\\ r9 = add(r9,#16)
\\ r1:0 = asl(r11:10,#47)
\\ }
\\
\\ {
\\ r15:14 = mpyu(r13,r13)
\\ }
\\ {
\\ r1:0 -= asl(r15:14,#47)
\\ r15:14 = mpyu(r13,r12)
\\ }
\\ {
\\ r1:0 -= asl(r15:14,#16)
\\ r15:14 = mpyu(r12,r12)
\\ }
\\ {
\\ r1:0 -= lsr(r15:14,#17)
\\ }
\\ {
\\ r1:0 = mpyu(r1,r8)
\\ }
\\ {
\\ r13:12 += lsr(r1:0,r9)
\\ }
\\ {
\\ r3:2 = mpyu(r13,r12)
\\ r5:4 = mpyu(r12,r12)
\\ r15:14 = #0
\\ r1:0 = #0
\\ }
\\ {
\\ r3:2 += lsr(r5:4,#33)
\\ r5:4 += asl(r3:2,#33)
\\ p1 = cmp.eq(r0,r0)
\\ }
\\ {
\\ r7:6 = mpyu(r13,r13)
\\ r1:0 = sub(r1:0,r5:4,p1):carry
\\ r9:8 = #1
\\ }
\\ {
\\ r7:6 += lsr(r3:2,#31)
\\ r9:8 += asl(r13:12,#1)
\\ }
\\
\\
\\
\\
\\
\\ {
\\ r15:14 = sub(r11:10,r7:6,p1):carry
\\ r5:4 = sub(r1:0,r9:8,p2):carry
\\
\\
\\
\\
\\ r7:6 = #1
\\ r11:10 = #0
\\ }
\\ {
\\ r3:2 = sub(r15:14,r11:10,p2):carry
\\ r7:6 = add(r13:12,r7:6)
\\ r28 = add(r28,#-0x3ff)
\\ }
\\ {
\\
\\ if (p2) r13:12 = r7:6
\\ if (p2) r1:0 = r5:4
\\ if (p2) r15:14 = r3:2
\\ }
\\ {
\\ r5:4 = sub(r1:0,r9:8,p3):carry
\\ r7:6 = #1
\\ r28 = asr(r28,#1)
\\ }
\\ {
\\ r3:2 = sub(r15:14,r11:10,p3):carry
\\ r7:6 = add(r13:12,r7:6)
\\ }
\\ {
\\ if (p3) r13:12 = r7:6
\\ if (p3) r1:0 = r5:4
\\
\\
\\
\\
\\
\\ r2 = #1
\\ }
\\ {
\\ p0 = cmp.eq(r1:0,r11:10)
\\ if (!p0.new) r12 = or(r12,r2)
\\ r3 = cl0(r13:12)
\\ r28 = add(r28,#-63)
\\ }
\\
\\
\\
\\ {
\\ r1:0 = convert_ud2df(r13:12)
\\ r28 = add(r28,r3)
\\ }
\\ {
\\ r1 += asl(r28,#52 -32)
\\ jumpr r31
\\ }
\\ .Lsqrt_abnormal:
\\ {
\\ p0 = dfclass(r1:0,#0x01)
\\ if (p0.new) jumpr:t r31
\\ }
\\ {
\\ p0 = dfclass(r1:0,#0x10)
\\ if (p0.new) jump:nt .Lsqrt_nan
\\ }
\\ {
\\ p0 = cmp.gt(r1,#-1)
\\ if (!p0.new) jump:nt .Lsqrt_invalid_neg
\\ if (!p0.new) r28 = ##0x7F800001
\\ }
\\ {
\\ p0 = dfclass(r1:0,#0x08)
\\ if (p0.new) jumpr:nt r31
\\ }
\\
\\
\\ {
\\ r1:0 = extractu(r1:0,#52,#0)
\\ }
\\ {
\\ r28 = add(clb(r1:0),#-11)
\\ }
\\ {
\\ r1:0 = asl(r1:0,r28)
\\ r28 = sub(#1,r28)
\\ }
\\ {
\\ r1 = insert(r28,#1,#52 -32)
\\ }
\\ {
\\ r3:2 = extractu(r1:0,#23 +1,#52 -23)
\\ r5 = ##0x3f000004
\\ }
\\ {
\\ r9 = or(r5,r2)
\\ r5 = and(r5,#-16)
\\ jump .Ldenormal_restart
\\ }
\\ .Lsqrt_nan:
\\ {
\\ r28 = convert_df2sf(r1:0)
\\ r1:0 = #-1
\\ jumpr r31
\\ }
\\ .Lsqrt_invalid_neg:
\\ {
\\ r1:0 = convert_sf2df(r28)
\\ jumpr r31
\\ }
);
}
