.file "__vrhypot.S"

#include "libm.h"

RO_DATA

.align 64

.CONST_TBL:

.word 0x7fe00000, 0x7fdfc07f, 0x7fdf81f8, 0x7fdf4465,

.word 0x7fdf07c1, 0x7fdecc07, 0x7fde9131, 0x7fde573a,

.word 0x7fde1e1e, 0x7fdde5d6, 0x7fddae60, 0x7fdd77b6,

.word 0x7fdd41d4, 0x7fdd0cb5, 0x7fdcd856, 0x7fdca4b3,

.word 0x7fdc71c7, 0x7fdc3f8f, 0x7fdc0e07, 0x7fdbdd2b,

.word 0x7fdbacf9, 0x7fdb7d6c, 0x7fdb4e81, 0x7fdb2036,

.word 0x7fdaf286, 0x7fdac570, 0x7fda98ef, 0x7fda6d01,

.word 0x7fda41a4, 0x7fda16d3, 0x7fd9ec8e, 0x7fd9c2d1,

.word 0x7fd99999, 0x7fd970e4, 0x7fd948b0, 0x7fd920fb,

.word 0x7fd8f9c1, 0x7fd8d301, 0x7fd8acb9, 0x7fd886e5,

.word 0x7fd86186, 0x7fd83c97, 0x7fd81818, 0x7fd7f405,

.word 0x7fd7d05f, 0x7fd7ad22, 0x7fd78a4c, 0x7fd767dc,

.word 0x7fd745d1, 0x7fd72428, 0x7fd702e0, 0x7fd6e1f7,

.word 0x7fd6c16c, 0x7fd6a13c, 0x7fd68168, 0x7fd661ec,

.word 0x7fd642c8, 0x7fd623fa, 0x7fd60581, 0x7fd5e75b,

.word 0x7fd5c988, 0x7fd5ac05, 0x7fd58ed2, 0x7fd571ed,

.word 0x7fd55555, 0x7fd53909, 0x7fd51d07, 0x7fd50150,

.word 0x7fd4e5e0, 0x7fd4cab8, 0x7fd4afd6, 0x7fd49539,

.word 0x7fd47ae1, 0x7fd460cb, 0x7fd446f8, 0x7fd42d66,

.word 0x7fd41414, 0x7fd3fb01, 0x7fd3e22c, 0x7fd3c995,

.word 0x7fd3b13b, 0x7fd3991c, 0x7fd38138, 0x7fd3698d,

.word 0x7fd3521c, 0x7fd33ae4, 0x7fd323e3, 0x7fd30d19,

.word 0x7fd2f684, 0x7fd2e025, 0x7fd2c9fb, 0x7fd2b404,

.word 0x7fd29e41, 0x7fd288b0, 0x7fd27350, 0x7fd25e22,

.word 0x7fd24924, 0x7fd23456, 0x7fd21fb7, 0x7fd20b47,

.word 0x7fd1f704, 0x7fd1e2ef, 0x7fd1cf06, 0x7fd1bb4a,

.word 0x7fd1a7b9, 0x7fd19453, 0x7fd18118, 0x7fd16e06,

.word 0x7fd15b1e, 0x7fd1485f, 0x7fd135c8, 0x7fd12358,

.word 0x7fd11111, 0x7fd0fef0, 0x7fd0ecf5, 0x7fd0db20,

.word 0x7fd0c971, 0x7fd0b7e6, 0x7fd0a681, 0x7fd0953f,

.word 0x7fd08421, 0x7fd07326, 0x7fd0624d, 0x7fd05197,

.word 0x7fd04104, 0x7fd03091, 0x7fd02040, 0x7fd01010,

.word 0x42300000, 0 ! D2ON36 = 2**36

.word 0xffffff00, 0 ! DA0

.word 0xfff00000, 0 ! DA1

.word 0x3ff00000, 0 ! DONE = 1.0

.word 0x40000000, 0 ! DTWO = 2.0

.word 0x7fd00000, 0 ! D2ON1022

.word 0x3cb00000, 0 ! D2ONM52

.word 0x43200000, 0 ! D2ON51

.word 0x0007ffff, 0xffffffff ! 0x0007ffffffffffff

#define stridex %l2

#define stridey %l3

#define stridez %l5

#define TBL_SHIFT 512

#define TBL %l1

#define counter %l4

#define _0x7ff00000 %l0

#define _0x00100000 %o5

#define _0x7fffffff %l6

#define D2ON36 %f4

#define DTWO %f6

#define DONE %f8

#define DA0 %f58

#define DA1 %f56

#define dtmp0 STACK_BIAS-0x80

#define dtmp1 STACK_BIAS-0x78

#define dtmp2 STACK_BIAS-0x70

#define dtmp3 STACK_BIAS-0x68

#define dtmp4 STACK_BIAS-0x60

#define dtmp5 STACK_BIAS-0x58

#define dtmp6 STACK_BIAS-0x50

#define dtmp7 STACK_BIAS-0x48

#define dtmp8 STACK_BIAS-0x40

#define dtmp9 STACK_BIAS-0x38

#define dtmp10 STACK_BIAS-0x30

#define dtmp11 STACK_BIAS-0x28

#define dtmp12 STACK_BIAS-0x20

#define dtmp13 STACK_BIAS-0x18

#define dtmp14 STACK_BIAS-0x10

#define dtmp15 STACK_BIAS-0x08

#define ftmp0 STACK_BIAS-0x100

#define tmp_px STACK_BIAS-0x98

#define tmp_py STACK_BIAS-0x90

#define tmp_counter STACK_BIAS-0x88

! sizeof temp storage - must be a multiple of 16 for V9

#define tmps 0x100

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! !!!!! algorithm !!!!!

! hx0 = *(int*)px;

! hy0 = *(int*)py;

!

! ((float*)&x0)[0] = ((float*)px)[0];

! ((float*)&x0)[1] = ((float*)px)[1];

! ((float*)&y0)[0] = ((float*)py)[0];

! ((float*)&y0)[1] = ((float*)py)[1];

!

! hx0 &= 0x7fffffff;

! hy0 &= 0x7fffffff;

!

! diff0 = hy0 - hx0;

! j0 = diff0 >> 31;

! j0 &= diff0;

! j0 = hy0 - j0;

! j0 &= 0x7ff00000;

!

! j0 = 0x7ff00000 - j0;

! ll = (long long)j0 << 32;

! *(long long*)&scl0 = ll;

!

! if ( hx0 >= 0x7ff00000 || hy0 >= 0x7ff00000 )

! {

! lx = ((int*)px)[1];

! ly = ((int*)py)[1];

!

! if ( hx0 == 0x7ff00000 && lx == 0 ) res0 = 0.0;

! else if ( hy0 == 0x7ff00000 && ly == 0 ) res0 = 0.0;

! else res0 = fabs(x0) * fabs(y0);

!

! ((float*)pz)[0] = ((float*)&res0)[0];

! ((float*)pz)[1] = ((float*)&res0)[1];

!

! px += stridex;

! py += stridey;

! pz += stridez;

! continue;

! }

! if ( hx0 < 0x00100000 && hy0 < 0x00100000 )

! {

! lx = ((int*)px)[1];

! ly = ((int*)py)[1];

! ii = hx0 | hy0;

! ii |= lx;

! ii |= ly;

! if ( ii == 0 )

! {

! res0 = 1.0 / 0.0;

! ((float*)pz)[0] = ((float*)&res0)[0];

! ((float*)pz)[1] = ((float*)&res0)[1];

!

! px += stridex;

! py += stridey;

! pz += stridez;

! continue;

! }

! x0 = fabs(x0);

! y0 = fabs(y0);

! if ( hx0 < 0x00080000 )

! {

! x0 = *(long long*)&x0;

! }

! else

! {

! ((long long*)&dtmp0)[0] = 0x0007ffffffffffffULL;

! x0 = vis_fand(x0, dtmp0);

! x0 = *(long long*)&x0;

! x0 += D2ON51;

! }

! x0 *= D2ONM52;

! if ( hy0 < 0x00080000 )

! {

! y0 = *(long long*)&y0;

! }

! else

! {

! ((long long*)&dtmp0)[0] = 0x0007ffffffffffffULL;

! y0 = vis_fand(y0, dtmp0);

! y0 = *(long long*)&y0;

! y0 += D2ON51;

! }

! y0 *= D2ONM52;

! *(long long*)&scl0 = 0x7fd0000000000000ULL;

! }

! else

! {

! x0 *= scl0;

! y0 *= scl0;

! }

!

! x_hi0 = x0 + D2ON36;

! y_hi0 = y0 + D2ON36;

! x_hi0 -= D2ON36;

! y_hi0 -= D2ON36;

! x_lo0 = x0 - x_hi0;

! y_lo0 = y0 - y_hi0;

! res0_hi = x_hi0 * x_hi0;

! dtmp0 = y_hi0 * y_hi0;

! res0_hi += dtmp0;

! res0_lo = x0 + x_hi0;

! res0_lo *= x_lo0;

! dtmp1 = y0 + y_hi0;

! dtmp1 *= y_lo0;

! res0_lo += dtmp1;

!

! dres = res0_hi + res0_lo;

! dexp0 = vis_fand(dres,DA1);

! iarr = ((int*)&dres)[0];

!

! iarr >>= 11;

! iarr &= 0x1fc;

! dtmp0 = ((double*)((char*)dll1 + iarr))[0];

! dd = vis_fpsub32(dtmp0, dexp0);

!

! dtmp0 = dd * dres;

! dtmp0 = DTWO - dtmp0;

! dd *= dtmp0;

! dtmp1 = dd * dres;

! dtmp1 = DTWO - dtmp1;

! dd *= dtmp1;

! dtmp2 = dd * dres;

! dtmp2 = DTWO - dtmp2;

! dres = dd * dtmp2;

!

! res0 = vis_fand(dres,DA0);

!

! dtmp0 = res0_hi * res0;

! dtmp0 = DONE - dtmp0;

! dtmp1 = res0_lo * res0;

! dtmp0 -= dtmp1;

! dtmp0 *= dres;

! res0 += dtmp0;

!

! res0 = sqrt ( res0 );

!

! res0 = scl0 * res0;

!

! ((float*)pz)[0] = ((float*)&res0)[0];

! ((float*)pz)[1] = ((float*)&res0)[1];

!

! px += stridex;

! py += stridey;

! pz += stridez;

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

ENTRY(__vrhypot)

save %sp,-SA(MINFRAME)-tmps,%sp

PIC_SETUP(l7)

PIC_SET(l7,.CONST_TBL,l1)

wr %g0,0x82,%asi

#ifdef __sparcv9

ldx [%fp+STACK_BIAS+176],stridez

#else

ld [%fp+STACK_BIAS+92],stridez

#endif

sll %i2,3,stridex

sethi %hi(0x7ff00000),_0x7ff00000

st %i0,[%fp+tmp_counter]

sll %i4,3,stridey

sethi %hi(0x00100000),_0x00100000

stx %i1,[%fp+tmp_px]

sll stridez,3,stridez

sethi %hi(0x7ffffc00),_0x7fffffff

stx %i3,[%fp+tmp_py]

ldd [TBL+TBL_SHIFT],D2ON36

add _0x7fffffff,1023,_0x7fffffff

ldd [TBL+TBL_SHIFT+8],DA0

ldd [TBL+TBL_SHIFT+16],DA1

ldd [TBL+TBL_SHIFT+24],DONE

ldd [TBL+TBL_SHIFT+32],DTWO

.begin:

ld [%fp+tmp_counter],counter

ldx [%fp+tmp_px],%i4

ldx [%fp+tmp_py],%i3

st %g0,[%fp+tmp_counter]

.begin1:

cmp counter,0

ble,pn %icc,.exit

lda [%i4]0x82,%o1 ! (7_0) hx0 = *(int*)px;

add %i4,stridex,%i1

lda [%i3]0x82,%o4 ! (7_0) hy0 = *(int*)py;

add %i3,stridey,%i0 ! py += stridey

and %o1,_0x7fffffff,%o7 ! (7_0) hx0 &= 0x7fffffff;

cmp %o7,_0x7ff00000 ! (7_0) hx0 ? 0x7ff00000

bge,pn %icc,.spec0 ! (7_0) if ( hx0 >= 0x7ff00000 )

and %o4,_0x7fffffff,%l7 ! (7_0) hy0 &= 0x7fffffff;

cmp %l7,_0x7ff00000 ! (7_0) hy0 ? 0x7ff00000

bge,pn %icc,.spec0 ! (7_0) if ( hy0 >= 0x7ff00000 )

sub %l7,%o7,%o1 ! (7_0) diff0 = hy0 - hx0;

sra %o1,31,%o3 ! (7_0) j0 = diff0 >> 31;

cmp %o7,_0x00100000 ! (7_0) hx0 ? 0x00100000

bl,pn %icc,.spec1 ! (7_0) if ( hx0 < 0x00100000 )

and %o1,%o3,%o1 ! (7_0) j0 &= diff0;

.cont_spec0:

sub %l7,%o1,%o4 ! (7_0) j0 = hy0 - j0;

and %o4,%l0,%o4 ! (7_0) j0 &= 0x7ff00000;

sub %l0,%o4,%g1 ! (7_0) j0 = 0x7ff00000 - j0;

sllx %g1,32,%g1 ! (7_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp15] ! (7_0) *(long long*)&scl0 = ll;

stx %g1,[%fp+dtmp0] ! (7_1) *(long long*)&scl0 = ll;

.cont_spec1:

lda [%i1]0x82,%o1 ! (0_0) hx0 = *(int*)px;

mov %i1,%i2

lda [%i0]0x82,%o4 ! (0_0) hy0 = *(int*)py;

and %o1,_0x7fffffff,%o7 ! (0_0) hx0 &= 0x7fffffff;

mov %i0,%o0

cmp %o7,_0x7ff00000 ! (0_0) hx0 ? 0x7ff00000

bge,pn %icc,.update0 ! (0_0) if ( hx0 >= 0x7ff00000 )

and %o4,_0x7fffffff,%l7 ! (0_0) hy0 &= 0x7fffffff;

cmp %l7,_0x7ff00000 ! (0_0) hy0 ? 0x7ff00000

sub %l7,%o7,%o1 ! (0_0) diff0 = hy0 - hx0;

bge,pn %icc,.update0 ! (0_0) if ( hy0 >= 0x7ff00000 )

sra %o1,31,%o3 ! (0_0) j0 = diff0 >> 31;

cmp %o7,_0x00100000 ! (0_0) hx0 ? 0x00100000

and %o1,%o3,%o1 ! (0_0) j0 &= diff0;

bl,pn %icc,.update1 ! (0_0) if ( hx0 < 0x00100000 )

sub %l7,%o1,%o4 ! (0_0) j0 = hy0 - j0;

.cont0:

and %o4,%l0,%o4 ! (0_0) j0 &= 0x7ff00000;

sub %l0,%o4,%o4 ! (0_0) j0 = 0x7ff00000 - j0;

.cont1:

sllx %o4,32,%o4 ! (0_0) ll = (long long)j0 << 32;

stx %o4,[%fp+dtmp1] ! (0_0) *(long long*)&scl0 = ll;

ldd [%fp+dtmp15],%f62 ! (7_1) *(long long*)&scl0 = ll;

lda [%i4]%asi,%f10 ! (7_1) ((float*)&x0)[0] = ((float*)px)[0];

lda [%i4+4]%asi,%f11 ! (7_1) ((float*)&x0)[1] = ((float*)px)[1];

lda [%i3]%asi,%f12 ! (7_1) ((float*)&y0)[0] = ((float*)py)[0];

add %i1,stridex,%i4 ! px += stridex

lda [%i3+4]%asi,%f13 ! (7_1) ((float*)&y0)[1] = ((float*)py)[1];

fmuld %f10,%f62,%f10 ! (7_1) x0 *= scl0;

add %i4,stridex,%i1 ! px += stridex

fmuld %f12,%f62,%f60 ! (7_1) y0 *= scl0;

lda [%i4]0x82,%o1 ! (1_0) hx0 = *(int*)px;

add %i0,stridey,%i3 ! py += stridey

faddd %f10,D2ON36,%f46 ! (7_1) x_hi0 = x0 + D2ON36;

lda [%i3]0x82,%g1 ! (1_0) hy0 = *(int*)py;

add %i3,stridey,%i0 ! py += stridey

faddd %f60,D2ON36,%f50 ! (7_1) y_hi0 = y0 + D2ON36;

and %o1,_0x7fffffff,%o7 ! (1_0) hx0 &= 0x7fffffff;

cmp %o7,_0x7ff00000 ! (1_0) hx0 ? 0x7ff00000

stx %o4,[%fp+dtmp2] ! (0_0) *(long long*)&scl0 = ll;

and %g1,_0x7fffffff,%l7 ! (1_0) hy0 &= 0x7fffffff;

bge,pn %icc,.update2 ! (1_0) if ( hx0 >= 0x7ff00000 )

fsubd %f46,D2ON36,%f20 ! (7_1) x_hi0 -= D2ON36;

cmp %l7,_0x7ff00000 ! (1_0) hy0 ? 0x7ff00000

sub %l7,%o7,%o1 ! (1_0) diff0 = hy0 - hx0;

bge,pn %icc,.update3 ! (1_0) if ( hy0 >= 0x7ff00000 )

fsubd %f50,D2ON36,%f54 ! (7_1) y_hi0 -= D2ON36;

sra %o1,31,%o3 ! (1_0) j0 = diff0 >> 31;

and %o1,%o3,%o1 ! (1_0) j0 &= diff0;

fmuld %f20,%f20,%f2 ! (7_1) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (1_0) j0 = hy0 - j0;

cmp %o7,_0x00100000 ! (1_0) hx0 ? 0x00100000

fsubd %f10,%f20,%f0 ! (7_1) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (7_1) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (1_0) j0 &= 0x7ff00000;

bl,pn %icc,.update4 ! (1_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (7_1) res0_lo = x0 + x_hi0;

sub %l0,%o4,%o4 ! (1_0) j0 = 0x7ff00000 - j0;

.cont4:

sllx %o4,32,%o4 ! (1_0) ll = (long long)j0 << 32;

stx %o4,[%fp+dtmp3] ! (1_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (7_1) dtmp1 = y0 + y_hi0;

fsubd %f60,%f54,%f12 ! (7_1) y_lo0 = y0 - y_hi0;

fmuld %f62,%f0,%f0 ! (7_1) res0_lo *= x_lo0;

ldd [%fp+dtmp1],%f62 ! (0_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f44 ! (7_1) res0_hi += dtmp0;

lda [%i2]%asi,%f10 ! (0_0) ((float*)&x0)[0] = ((float*)px)[0];

lda [%i2+4]%asi,%f11 ! (0_0) ((float*)&x0)[1] = ((float*)px)[1];

fmuld %f50,%f12,%f26 ! (7_1) dtmp1 *= y_lo0;

lda [%o0]%asi,%f12 ! (0_0) ((float*)&y0)[0] = ((float*)py)[0];

lda [%o0+4]%asi,%f13 ! (0_0) ((float*)&y0)[1] = ((float*)py)[1];

fmuld %f10,%f62,%f10 ! (0_0) x0 *= scl0;

fmuld %f12,%f62,%f60 ! (0_0) y0 *= scl0;

faddd %f0,%f26,%f38 ! (7_1) res0_lo += dtmp1;

lda [%i1]0x82,%o1 ! (2_0) hx0 = *(int*)px;

mov %i1,%i2

faddd %f10,D2ON36,%f46 ! (0_0) x_hi0 = x0 + D2ON36;

lda [%i0]0x82,%g1 ! (2_0) hy0 = *(int*)py;

mov %i0,%o0

faddd %f60,D2ON36,%f12 ! (0_0) y_hi0 = y0 + D2ON36;

faddd %f44,%f38,%f14 ! (7_1) dres = res0_hi + res0_lo;

and %o1,_0x7fffffff,%o7 ! (2_0) hx0 &= 0x7fffffff;

cmp %o7,_0x7ff00000 ! (2_0) hx0 ? 0x7ff00000

bge,pn %icc,.update5 ! (2_0) if ( hx0 >= 0x7ff00000 )

stx %o4,[%fp+dtmp4] ! (1_0) *(long long*)&scl0 = ll;

and %g1,_0x7fffffff,%l7 ! (2_0) hx0 &= 0x7fffffff;

st %f14,[%fp+ftmp0] ! (7_1) iarr = ((int*)&dres)[0];

fsubd %f46,D2ON36,%f20 ! (0_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (2_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (2_0) hy0 ? 0x7ff00000

bge,pn %icc,.update6 ! (2_0) if ( hy0 >= 0x7ff00000 )

fsubd %f12,D2ON36,%f54 ! (0_0) y_hi0 -= D2ON36;

sra %o1,31,%o3 ! (2_0) j0 = diff0 >> 31;

and %o1,%o3,%o1 ! (2_0) j0 &= diff0;

fmuld %f20,%f20,%f2 ! (0_0) res0_hi = x_hi0 * x_hi0;

cmp %o7,_0x00100000 ! (2_0) hx0 ? 0x00100000

sub %l7,%o1,%o4 ! (2_0) j0 = hy0 - j0;

fsubd %f10,%f20,%f0 ! (0_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (0_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (2_0) j0 &= 0x7ff00000;

bl,pn %icc,.update7 ! (2_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (0_0) res0_lo = x0 + x_hi0;

.cont7:

sub %l0,%o4,%g1 ! (2_0) j0 = 0x7ff00000 - j0;

sllx %g1,32,%g1 ! (2_0) ll = (long long)j0 << 32;

.cont8:

stx %g1,[%fp+dtmp5] ! (2_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (0_0) dtmp1 = y0 + y_hi0;

fsubd %f60,%f54,%f12 ! (0_0) y_lo0 = y0 - y_hi0;

fmuld %f62,%f0,%f0 ! (0_0) res0_lo *= x_lo0;

ldd [%fp+dtmp3],%f62 ! (1_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f32 ! (0_0) res0_hi += dtmp0;

lda [%i4]%asi,%f10 ! (1_0) ((float*)&x0)[0] = ((float*)px)[0];

lda [%i4+4]%asi,%f11 ! (1_0) ((float*)&x0)[1] = ((float*)px)[1];

fmuld %f50,%f12,%f28 ! (0_0) dtmp1 *= y_lo0;

lda [%i3]%asi,%f12 ! (1_0) ((float*)&y0)[0] = ((float*)py)[0];

add %i1,stridex,%i4 ! px += stridex

lda [%i3+4]%asi,%f13 ! (1_0) ((float*)&y0)[1] = ((float*)py)[1];

ld [%fp+ftmp0],%o2 ! (7_1) iarr = ((int*)&dres)[0];

add %i4,stridex,%i1 ! px += stridex

fand %f14,DA1,%f2 ! (7_1) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (1_0) x0 *= scl0;

fmuld %f12,%f62,%f60 ! (1_0) y0 *= scl0;

sra %o2,11,%i3 ! (7_1) iarr >>= 11;

faddd %f0,%f28,%f36 ! (0_0) res0_lo += dtmp1;

and %i3,0x1fc,%i3 ! (7_1) iarr &= 0x1fc;

add %i3,TBL,%o4 ! (7_1) (char*)dll1 + iarr

lda [%i4]0x82,%o1 ! (3_0) hx0 = *(int*)px;

add %i0,stridey,%i3 ! py += stridey

ld [%o4],%f26 ! (7_1) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f46 ! (1_0) x_hi0 = x0 + D2ON36;

lda [%i3]0x82,%o4 ! (3_0) hy0 = *(int*)py;

add %i3,stridey,%i0 ! py += stridey

faddd %f60,D2ON36,%f12 ! (1_0) y_hi0 = y0 + D2ON36;

faddd %f32,%f36,%f22 ! (0_0) dres = res0_hi + res0_lo;

and %o1,_0x7fffffff,%o7 ! (3_0) hx0 &= 0x7fffffff;

cmp %o7,_0x7ff00000 ! (3_0) hx0 ? 0x7ff00000

stx %g1,[%fp+dtmp6] ! (2_0) *(long long*)&scl0 = ll;

bge,pn %icc,.update9 ! (3_0) if ( hx0 >= 0x7ff00000 )

fpsub32 %f26,%f2,%f26 ! (7_1) dd = vis_fpsub32(dtmp0, dexp0);

and %o4,_0x7fffffff,%l7 ! (3_0) hy0 &= 0x7fffffff;

st %f22,[%fp+ftmp0] ! (0_0) iarr = ((int*)&dres)[0];

fsubd %f46,D2ON36,%f20 ! (1_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (3_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (3_0) hy0 ? 0x7ff00000

bge,pn %icc,.update10 ! (3_0) if ( hy0 >= 0x7ff00000 )

fsubd %f12,D2ON36,%f54 ! (1_0) y_hi0 -= D2ON36;

fmuld %f26,%f14,%f50 ! (7_1) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (3_0) j0 = diff0 >> 31;

and %o1,%o3,%o1 ! (3_0) j0 &= diff0;

fmuld %f20,%f20,%f2 ! (1_0) res0_hi = x_hi0 * x_hi0;

cmp %o7,_0x00100000 ! (3_0) hx0 ? 0x00100000

sub %l7,%o1,%o4 ! (3_0) j0 = hy0 - j0;

fsubd %f10,%f20,%f0 ! (1_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (1_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (3_0) j0 &= 0x7ff00000;

bl,pn %icc,.update11 ! (3_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (1_0) res0_lo = x0 + x_hi0;

.cont11:

sub %l0,%o4,%g1 ! (3_0) j0 = 0x7ff00000 - j0;

fsubd DTWO,%f50,%f20 ! (7_1) dtmp0 = DTWO - dtmp0;

.cont12:

sllx %g1,32,%g1 ! (3_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp7] ! (3_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (1_0) dtmp1 = y0 + y_hi0;

fsubd %f60,%f54,%f12 ! (1_0) y_lo0 = y0 - y_hi0

fmuld %f62,%f0,%f0 ! (1_0) res0_lo *= x_lo0;

ldd [%fp+dtmp5],%f62 ! (2_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f42 ! (1_0) res0_hi += dtmp0;

lda [%i2]%asi,%f10 ! (2_0) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f26,%f20,%f54 ! (7_1) dd *= dtmp0;

lda [%i2+4]%asi,%f11 ! (2_0) ((float*)&x0)[1] = ((float*)px)[1];

fmuld %f50,%f12,%f26 ! (1_0) dtmp1 *= y_lo0;

lda [%o0]%asi,%f12 ! (2_0) ((float*)&y0)[0] = ((float*)py)[0];

lda [%o0+4]%asi,%f13 ! (2_0) ((float*)&y0)[1] = ((float*)py)[1];

fmuld %f54,%f14,%f50 ! (7_1) dtmp1 = dd * dres;

ld [%fp+ftmp0],%o2 ! (0_0) iarr = ((int*)&dres)[0];

fand %f22,DA1,%f2 ! (0_0) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (2_0) x0 *= scl0;

fmuld %f12,%f62,%f60 ! (2_0) y0 *= scl0;

sra %o2,11,%o4 ! (0_0) iarr >>= 11;

faddd %f0,%f26,%f34 ! (1_0) res0_lo += dtmp1;

and %o4,0x1fc,%o4 ! (0_0) iarr &= 0x1fc;

add %o4,TBL,%o4 ! (0_0) (char*)dll1 + iarr

mov %i1,%i2

lda [%i1]0x82,%o1 ! (4_0) hx0 = *(int*)px;

fsubd DTWO,%f50,%f20 ! (7_1) dtmp1 = DTWO - dtmp1;

ld [%o4],%f28 ! (0_0) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f46 ! (2_0) x_hi0 = x0 + D2ON36;

lda [%i0]0x82,%o4 ! (4_0) hy0 = *(int*)py;

mov %i0,%o0

faddd %f60,D2ON36,%f50 ! (2_0) y_hi0 = y0 + D2ON36;

and %o1,_0x7fffffff,%o7 ! (4_0) hx0 &= 0x7fffffff;

faddd %f42,%f34,%f18 ! (1_0) dres = res0_hi + res0_lo;

fmuld %f54,%f20,%f16 ! (7_1) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (4_0) hx0 ? 0x7ff00000

stx %g1,[%fp+dtmp8] ! (3_0) *(long long*)&scl0 = ll;

fpsub32 %f28,%f2,%f28 ! (0_0) dd = vis_fpsub32(dtmp0, dexp0);

and %o4,_0x7fffffff,%l7 ! (4_0) hy0 &= 0x7fffffff;

bge,pn %icc,.update13 ! (4_0) if ( hx0 >= 0x7ff00000 )

st %f18,[%fp+ftmp0] ! (1_0) iarr = ((int*)&dres)[0];

fsubd %f46,D2ON36,%f20 ! (2_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (4_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (4_0) hy0 ? 0x7ff00000

bge,pn %icc,.update14 ! (4_0) if ( hy0 >= 0x7ff00000 )

fsubd %f50,D2ON36,%f54 ! (2_0) y_hi0 -= D2ON36;

fmuld %f28,%f22,%f50 ! (0_0) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (4_0) j0 = diff0 >> 31;

and %o1,%o3,%o1 ! (4_0) j0 &= diff0;

fmuld %f20,%f20,%f2 ! (2_0) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (4_0) j0 = hy0 - j0;

cmp %o7,_0x00100000 ! (4_0) hx0 ? 0x00100000

fsubd %f10,%f20,%f0 ! (2_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (2_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (4_0) j0 &= 0x7ff00000;

bl,pn %icc,.update15 ! (4_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (2_0) res0_lo = x0 + x_hi0;

.cont15:

sub %l0,%o4,%g1 ! (4_0) j0 = 0x7ff00000 - j0;

fsubd DTWO,%f50,%f20 ! (0_0) dtmp0 = DTWO - dtmp0;

.cont16:

fmuld %f16,%f14,%f14 ! (7_1) dtmp2 = dd * dres;

sllx %g1,32,%g1 ! (4_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp9] ! (4_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (2_0) dtmp1 = y0 + y_hi0;

fsubd %f60,%f54,%f12 ! (2_0) y_lo0 = y0 - y_hi0;

fmuld %f62,%f0,%f0 ! (2_0) res0_lo *= x_lo0;

ldd [%fp+dtmp7],%f62 ! (3_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f30 ! (2_0) res0_hi += dtmp0;

lda [%i4]%asi,%f10 ! (3_0) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f28,%f20,%f54 ! (0_0) dd *= dtmp0;

lda [%i4+4]%asi,%f11 ! (3_0) ((float*)&x0)[1] = ((float*)px)[1];

fmuld %f50,%f12,%f28 ! (2_0) dtmp1 *= y_lo0;

lda [%i3]%asi,%f12 ! (3_0) ((float*)&y0)[0] = ((float*)py)[0];

fsubd DTWO,%f14,%f20 ! (7_1) dtmp2 = DTWO - dtmp2;

lda [%i3+4]%asi,%f13 ! (3_0) ((float*)&y0)[1] = ((float*)py)[1];

add %i1,stridex,%i4 ! px += stridex

fmuld %f54,%f22,%f50 ! (0_0) dtmp1 = dd * dres;

ld [%fp+ftmp0],%o2 ! (1_0) iarr = ((int*)&dres)[0];

add %i4,stridex,%i1 ! px += stridex

fand %f18,DA1,%f2 ! (1_0) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (3_0) x0 *= scl0;

fmuld %f12,%f62,%f60 ! (3_0) y0 *= scl0;

sra %o2,11,%i3 ! (1_0) iarr >>= 11;

faddd %f0,%f28,%f40 ! (2_0) res0_lo += dtmp1;

and %i3,0x1fc,%i3 ! (1_0) iarr &= 0x1fc;

fmuld %f16,%f20,%f28 ! (7_1) dres = dd * dtmp2;

add %i3,TBL,%o4 ! (1_0) (char*)dll1 + iarr

lda [%i4]0x82,%o1 ! (5_0) hx0 = *(int*)px;

fsubd DTWO,%f50,%f20 ! (0_0) dtmp1 = DTWO - dtmp1;

add %i0,stridey,%i3 ! py += stridey

ld [%o4],%f26 ! (1_0) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f46 ! (3_0) x_hi0 = x0 + D2ON36;

lda [%i3]0x82,%o4 ! (5_0) hy0 = *(int*)py;

add %i3,stridey,%i0 ! py += stridey

faddd %f60,D2ON36,%f50 ! (3_0) y_hi0 = y0 + D2ON36;

and %o1,_0x7fffffff,%o7 ! (5_0) hx0 &= 0x7fffffff;

faddd %f30,%f40,%f14 ! (2_0) dres = res0_hi + res0_lo;

fmuld %f54,%f20,%f24 ! (0_0) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (5_0) hx0 ? 0x7ff00000

stx %g1,[%fp+dtmp10] ! (4_0) *(long long*)&scl0 = ll;

fpsub32 %f26,%f2,%f26 ! (1_0) dd = vis_fpsub32(dtmp0, dexp0);

and %o4,_0x7fffffff,%l7 ! (5_0) hy0 &= 0x7fffffff;

st %f14,[%fp+ftmp0] ! (2_0) iarr = ((int*)&dres)[0];

bge,pn %icc,.update17 ! (5_0) if ( hx0 >= 0x7ff00000 )

fsubd %f46,D2ON36,%f20 ! (3_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (5_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (5_0) hy0 ? 0x7ff00000

bge,pn %icc,.update18 ! (5_0) if ( hy0 >= 0x7ff00000 )

fsubd %f50,D2ON36,%f54 ! (3_0) y_hi0 -= D2ON36;

fmuld %f26,%f18,%f50 ! (1_0) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (5_0) j0 = diff0 >> 31;

and %o1,%o3,%o1 ! (5_0) j0 &= diff0;

fand %f28,DA0,%f48 ! (7_1) res0 = vis_fand(dres,DA0);

fmuld %f20,%f20,%f2 ! (3_0) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (5_0) j0 = hy0 - j0;

cmp %o7,_0x00100000 ! (5_0) hx0 ? 0x00100000

fsubd %f10,%f20,%f0 ! (3_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (3_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (5_0) j0 &= 0x7ff00000;

bl,pn %icc,.update19 ! (5_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (3_0) res0_lo = x0 + x_hi0;

.cont19a:

fmuld %f44,%f48,%f10 ! (7_1) dtmp0 = res0_hi * res0;

sub %l0,%o4,%g1 ! (5_0) j0 = 0x7ff00000 - j0;

fsubd DTWO,%f50,%f20 ! (1_0) dtmp0 = DTWO - dtmp0;

.cont19b:

fmuld %f24,%f22,%f22 ! (0_0) dtmp2 = dd * dres;

sllx %g1,32,%g1 ! (5_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp11] ! (5_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (3_0) dtmp1 = y0 + y_hi0;

fmuld %f38,%f48,%f38 ! (7_1) dtmp1 = res0_lo * res0;

fsubd %f60,%f54,%f12 ! (3_0) y_lo0 = y0 - y_hi0;

.cont20:

fmuld %f62,%f0,%f0 ! (3_0) res0_lo *= x_lo0;

ldd [%fp+dtmp9],%f62 ! (4_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f44 ! (3_0) res0_hi += dtmp0;

fsubd DONE,%f10,%f60 ! (7_1) dtmp0 = DONE - dtmp0;

lda [%i2]%asi,%f10 ! (4_0) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f26,%f20,%f54 ! (1_0) dd *= dtmp0;

lda [%i2+4]%asi,%f11 ! (4_0) ((float*)&x0)[1] = ((float*)px)[1];

fmuld %f50,%f12,%f26 ! (3_0) dtmp1 *= y_lo0;

lda [%o0]%asi,%f12 ! (4_0) ((float*)&y0)[0] = ((float*)py)[0];

fsubd DTWO,%f22,%f20 ! (0_0) dtmp2 = DTWO - dtmp2;

lda [%o0+4]%asi,%f13 ! (4_0) ((float*)&y0)[1] = ((float*)py)[1];

fmuld %f54,%f18,%f50 ! (1_0) dtmp1 = dd * dres;

ld [%fp+ftmp0],%o2 ! (2_0) iarr = ((int*)&dres)[0];

fand %f14,DA1,%f2 ! (2_0) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (4_0) x0 *= scl0;

fsubd %f60,%f38,%f46 ! (7_1) dtmp0 -= dtmp1;

fmuld %f12,%f62,%f60 ! (4_0) y0 *= scl0;

sra %o2,11,%o4 ! (2_0) iarr >>= 11;

faddd %f0,%f26,%f38 ! (3_0) res0_lo += dtmp1;

and %o4,0x1fc,%o4 ! (2_0) iarr &= 0x1fc;

fmuld %f24,%f20,%f26 ! (0_0) dres = dd * dtmp2;

add %o4,TBL,%o4 ! (2_0) (char*)dll1 + iarr

mov %i1,%i2

lda [%i1]0x82,%o1 ! (6_0) hx0 = *(int*)px;

fsubd DTWO,%f50,%f52 ! (1_0) dtmp1 = DTWO - dtmp1;

fmuld %f46,%f28,%f28 ! (7_1) dtmp0 *= dres;

ld [%o4],%f20 ! (2_0) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f46 ! (4_0) x_hi0 = x0 + D2ON36;

lda [%i0]0x82,%o4 ! (6_0) hy0 = *(int*)py;

mov %i0,%o0

faddd %f60,D2ON36,%f50 ! (4_0) y_hi0 = y0 + D2ON36;

and %o1,_0x7fffffff,%o7 ! (6_0) hx0 &= 0x7fffffff;

faddd %f44,%f38,%f22 ! (3_0) dres = res0_hi + res0_lo;

fmuld %f54,%f52,%f16 ! (1_0) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (6_0) hx0 ? 0x7ff00000

stx %g1,[%fp+dtmp12] ! (5_0) *(long long*)&scl0 = ll;

fpsub32 %f20,%f2,%f52 ! (2_0) dd = vis_fpsub32(dtmp0, dexp0);

and %o4,_0x7fffffff,%l7 ! (6_0) hy0 &= 0x7fffffff;

st %f22,[%fp+ftmp0] ! (3_0) iarr = ((int*)&dres)[0];

bge,pn %icc,.update21 ! (6_0) if ( hx0 >= 0x7ff00000 )

fsubd %f46,D2ON36,%f46 ! (4_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (6_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (6_0) hy0 ? 0x7ff00000

bge,pn %icc,.update22 ! (6_0) if ( hy0 >= 0x7ff00000 )

fsubd %f50,D2ON36,%f54 ! (4_0) y_hi0 -= D2ON36;

fmuld %f52,%f14,%f50 ! (2_0) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (6_0) j0 = diff0 >> 31;

faddd %f48,%f28,%f48 ! (7_1) res0 += dtmp0;

and %o1,%o3,%o1 ! (6_0) j0 &= diff0;

fand %f26,DA0,%f28 ! (0_0) res0 = vis_fand(dres,DA0);

fmuld %f46,%f46,%f0 ! (4_0) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (6_0) j0 = hy0 - j0;

cmp %o7,_0x00100000 ! (6_0) hx0 ? 0x00100000

fsubd %f10,%f46,%f2 ! (4_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f20 ! (4_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (6_0) j0 &= 0x7ff00000;

bl,pn %icc,.update23 ! (6_0) if ( hx0 < 0x00100000 )

faddd %f10,%f46,%f62 ! (4_0) res0_lo = x0 + x_hi0;

.cont23a:

fmuld %f16,%f18,%f18 ! (1_0) dtmp2 = dd * dres;

sub %l0,%o4,%g1 ! (6_0) j0 = 0x7ff00000 - j0;

fsubd DTWO,%f50,%f10 ! (2_0) dtmp0 = DTWO - dtmp0;

.cont23b:

fmuld %f32,%f28,%f50 ! (0_0) dtmp0 = res0_hi * res0;

sllx %g1,32,%g1 ! (6_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp13] ! (6_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f46 ! (4_0) dtmp1 = y0 + y_hi0;

fmuld %f36,%f28,%f36 ! (0_0) dtmp1 = res0_lo * res0;

fsubd %f60,%f54,%f60 ! (4_0) y_lo0 = y0 - y_hi0;

.cont24:

fmuld %f62,%f2,%f2 ! (4_0) res0_lo *= x_lo0;

ldd [%fp+dtmp11],%f62 ! (5_0) *(long long*)&scl0 = ll;

faddd %f0,%f20,%f32 ! (4_0) res0_hi += dtmp0;

lda [%i4]%asi,%f0 ! (5_0) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f52,%f10,%f10 ! (2_0) dd *= dtmp0;

lda [%i4+4]%asi,%f1 ! (5_0) ((float*)&x0)[1] = ((float*)px)[1];

fsubd DONE,%f50,%f52 ! (0_0) dtmp0 = DONE - dtmp0;

fmuld %f46,%f60,%f46 ! (4_0) dtmp1 *= y_lo0;

lda [%i3]%asi,%f12 ! (5_0) ((float*)&y0)[0] = ((float*)py)[0];

fsubd DTWO,%f18,%f18 ! (1_0) dtmp2 = DTWO - dtmp2;

add %i1,stridex,%i4 ! px += stridex

lda [%i3+4]%asi,%f13 ! (5_0) ((float*)&y0)[1] = ((float*)py)[1];

fmuld %f10,%f14,%f50 ! (2_0) dtmp1 = dd * dres;

add %i4,stridex,%i1 ! px += stridex

ld [%fp+ftmp0],%o2 ! (3_0) iarr = ((int*)&dres)[0];

fand %f22,DA1,%f54 ! (3_0) dexp0 = vis_fand(dres,DA1);

fmuld %f0,%f62,%f60 ! (5_0) x0 *= scl0;

fsubd %f52,%f36,%f20 ! (0_0) dtmp0 -= dtmp1;

fmuld %f12,%f62,%f52 ! (5_0) y0 *= scl0;

sra %o2,11,%i3 ! (3_0) iarr >>= 11;

faddd %f2,%f46,%f36 ! (4_0) res0_lo += dtmp1;

and %i3,0x1fc,%i3 ! (3_0) iarr &= 0x1fc;

fmuld %f16,%f18,%f16 ! (1_0) dres = dd * dtmp2;

fsqrtd %f48,%f18 ! (7_1) res0 = sqrt ( res0 );

add %i3,TBL,%o4 ! (3_0) (char*)dll1 + iarr

lda [%i4]0x82,%o1 ! (7_0) hx0 = *(int*)px;

fsubd DTWO,%f50,%f46 ! (2_0) dtmp1 = DTWO - dtmp1;

fmuld %f20,%f26,%f48 ! (0_0) dtmp0 *= dres;

add %i0,stridey,%i3 ! py += stridey

ld [%o4],%f20 ! (3_0) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f60,D2ON36,%f50 ! (5_0) x_hi0 = x0 + D2ON36;

lda [%i3]0x82,%o4 ! (7_0) hy0 = *(int*)py;

add %i3,stridey,%i0 ! py += stridey

faddd %f52,D2ON36,%f12 ! (5_0) y_hi0 = y0 + D2ON36;

and %o1,_0x7fffffff,%o7 ! (7_0) hx0 &= 0x7fffffff;

faddd %f32,%f36,%f24 ! (4_0) dres = res0_hi + res0_lo;

fmuld %f10,%f46,%f26 ! (2_0) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (7_0) hx0 ? 0x7ff00000

stx %g1,[%fp+dtmp14] ! (6_0) *(long long*)&scl0 = ll;

fpsub32 %f20,%f54,%f10 ! (3_0) dd = vis_fpsub32(dtmp0, dexp0);

and %o4,_0x7fffffff,%l7 ! (7_0) hy0 &= 0x7fffffff;

st %f24,[%fp+ftmp0] ! (4_0) iarr = ((int*)&dres)[0];

bge,pn %icc,.update25 ! (7_0) if ( hx0 >= 0x7ff00000 )

fsubd %f50,D2ON36,%f20 ! (5_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (7_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (7_0) hy0 ? 0x7ff00000

bge,pn %icc,.update26 ! (7_0) if ( hy0 >= 0x7ff00000 )

fsubd %f12,D2ON36,%f54 ! (5_0) y_hi0 -= D2ON36;

fmuld %f10,%f22,%f50 ! (3_0) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (7_0) j0 = diff0 >> 31;

faddd %f28,%f48,%f48 ! (0_0) res0 += dtmp0;

and %o1,%o3,%o1 ! (7_0) j0 &= diff0;

fand %f16,DA0,%f28 ! (1_0) res0 = vis_fand(dres,DA0);

fmuld %f20,%f20,%f0 ! (5_0) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (7_0) j0 = hy0 - j0;

cmp %o7,_0x00100000 ! (7_0) hx0 ? 0x00100000

fsubd %f60,%f20,%f2 ! (5_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (5_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (7_0) j0 &= 0x7ff00000;

bl,pn %icc,.update27 ! (7_0) if ( hx0 < 0x00100000 )

faddd %f60,%f20,%f62 ! (5_0) res0_lo = x0 + x_hi0;

.cont27a:

fmuld %f26,%f14,%f14 ! (2_0) dtmp2 = dd * dres;

sub %l0,%o4,%g1 ! (7_0) j0 = 0x7ff00000 - j0;

fsubd DTWO,%f50,%f20 ! (3_0) dtmp0 = DTWO - dtmp0;

.cont27b:

fmuld %f42,%f28,%f60 ! (1_0) dtmp0 = res0_hi * res0;

sllx %g1,32,%g1 ! (7_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp15] ! (7_0) *(long long*)&scl0 = ll;

faddd %f52,%f54,%f50 ! (5_0) dtmp1 = y0 + y_hi0;

fmuld %f34,%f28,%f34 ! (1_0) dtmp1 = res0_lo * res0;

fsubd %f52,%f54,%f54 ! (5_0) y_lo0 = y0 - y_hi0;

.cont28:

fmuld %f62,%f2,%f2 ! (5_0) res0_lo *= x_lo0;

ldd [%fp+dtmp13],%f62 ! (6_0) *(long long*)&scl0 = ll;

faddd %f0,%f46,%f42 ! (5_0) res0_hi += dtmp0;

fmuld %f10,%f20,%f52 ! (3_0) dd *= dtmp0;

lda [%i2]%asi,%f10 ! (6_0) ((float*)&x0)[0] = ((float*)px)[0];

lda [%i2+4]%asi,%f11 ! (6_0) ((float*)&x0)[1] = ((float*)px)[1];

fsubd DONE,%f60,%f60 ! (1_0) dtmp0 = DONE - dtmp0;

fmuld %f50,%f54,%f46 ! (5_0) dtmp1 *= y_lo0;

lda [%o0]%asi,%f12 ! (6_0) ((float*)&y0)[0] = ((float*)py)[0];

fsubd DTWO,%f14,%f14 ! (2_0) dtmp2 = DTWO - dtmp2;

lda [%o0+4]%asi,%f13 ! (6_0) ((float*)&y0)[1] = ((float*)py)[1];

fmuld %f52,%f22,%f50 ! (3_0) dtmp1 = dd * dres;

ld [%fp+ftmp0],%o2 ! (4_0) iarr = ((int*)&dres)[0];

fand %f24,DA1,%f54 ! (4_0) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (6_0) x0 *= scl0;

ldd [%fp+dtmp0],%f0 ! (7_1) *(long long*)&scl0 = ll;

fsubd %f60,%f34,%f20 ! (1_0) dtmp0 -= dtmp1;

fmuld %f12,%f62,%f60 ! (6_0) y0 *= scl0;

sra %o2,11,%o4 ! (4_0) iarr >>= 11;

faddd %f2,%f46,%f34 ! (5_0) res0_lo += dtmp1;

and %o4,0x1fc,%o4 ! (4_0) iarr &= 0x1fc;

fmuld %f26,%f14,%f26 ! (2_0) dres = dd * dtmp2;

cmp counter,8

bl,pn %icc,.tail

nop

ba .main_loop

sub counter,8,counter

.align 16

.main_loop:

fsqrtd %f48,%f14 ! (0_1) res0 = sqrt ( res0 );

add %o4,TBL,%o4 ! (4_1) (char*)dll1 + iarr

lda [%i1]0x82,%o1 ! (0_0) hx0 = *(int*)px;

fsubd DTWO,%f50,%f46 ! (3_1) dtmp1 = DTWO - dtmp1;

fmuld %f20,%f16,%f48 ! (1_1) dtmp0 *= dres;

mov %i1,%i2

ld [%o4],%f20 ! (4_1) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f50 ! (6_1) x_hi0 = x0 + D2ON36;

nop

mov %i0,%o0

lda [%i0]0x82,%o4 ! (0_0) hy0 = *(int*)py;

faddd %f60,D2ON36,%f2 ! (6_1) y_hi0 = y0 + D2ON36;

faddd %f42,%f34,%f16 ! (5_1) dres = res0_hi + res0_lo;

and %o1,_0x7fffffff,%o7 ! (0_0) hx0 &= 0x7fffffff;

st %f16,[%fp+ftmp0] ! (5_1) iarr = ((int*)&dres)[0];

fmuld %f0,%f18,%f0 ! (7_2) res0 = scl0 * res0;

fmuld %f52,%f46,%f18 ! (3_1) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (0_0) hx0 ? 0x7ff00000

st %f0,[%i5] ! (7_2) ((float*)pz)[0] = ((float*)&res0)[0];

fpsub32 %f20,%f54,%f54 ! (4_1) dd = vis_fpsub32(dtmp0, dexp0);

and %o4,_0x7fffffff,%l7 ! (0_0) hy0 &= 0x7fffffff;

st %f1,[%i5+4] ! (7_2) ((float*)pz)[1] = ((float*)&res0)[1];

bge,pn %icc,.update29 ! (0_0) if ( hx0 >= 0x7ff00000 )

fsubd %f50,D2ON36,%f20 ! (6_1) x_hi0 -= D2ON36;

cmp %l7,_0x7ff00000 ! (0_0) hy0 ? 0x7ff00000

sub %l7,%o7,%o1 ! (0_0) diff0 = hy0 - hx0;

bge,pn %icc,.update30 ! (0_0) if ( hy0 >= 0x7ff00000 )

fsubd %f2,D2ON36,%f2 ! (6_1) y_hi0 -= D2ON36;

fmuld %f54,%f24,%f50 ! (4_1) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (0_0) j0 = diff0 >> 31;

stx %g1,[%fp+dtmp0] ! (7_1) *(long long*)&scl0 = ll;

faddd %f28,%f48,%f52 ! (1_1) res0 += dtmp0;

and %o1,%o3,%o1 ! (0_0) j0 &= diff0;

cmp %o7,_0x00100000 ! (0_0) hx0 ? 0x00100000

bl,pn %icc,.update31 ! (0_0) if ( hx0 < 0x00100000 )

fand %f26,DA0,%f48 ! (2_1) res0 = vis_fand(dres,DA0);

.cont31:

fmuld %f20,%f20,%f0 ! (6_1) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (0_0) j0 = hy0 - j0;

nop

fsubd %f10,%f20,%f28 ! (6_1) x_lo0 = x0 - x_hi0;

fmuld %f2,%f2,%f46 ! (6_1) dtmp0 = y_hi0 * y_hi0;

add %i5,stridez,%i5 ! pz += stridez

and %o4,%l0,%o4 ! (0_0) j0 &= 0x7ff00000;

faddd %f10,%f20,%f62 ! (6_1) res0_lo = x0 + x_hi0;

fmuld %f18,%f22,%f22 ! (3_1) dtmp2 = dd * dres;

sub %l0,%o4,%o4 ! (0_0) j0 = 0x7ff00000 - j0;

nop

fsubd DTWO,%f50,%f20 ! (4_1) dtmp0 = DTWO - dtmp0;

.cont32:

fmuld %f30,%f48,%f12 ! (2_1) dtmp0 = res0_hi * res0;

sllx %o4,32,%o4 ! (0_0) ll = (long long)j0 << 32;

stx %o4,[%fp+dtmp1] ! (0_0) *(long long*)&scl0 = ll;

faddd %f60,%f2,%f50 ! (6_1) dtmp1 = y0 + y_hi0;

fmuld %f40,%f48,%f40 ! (2_1) dtmp1 = res0_lo * res0;

nop

bn,pn %icc,.exit

fsubd %f60,%f2,%f2 ! (6_1) y_lo0 = y0 - y_hi0;

fmuld %f62,%f28,%f28 ! (6_1) res0_lo *= x_lo0;

nop

ldd [%fp+dtmp15],%f62 ! (7_1) *(long long*)&scl0 = ll;

faddd %f0,%f46,%f30 ! (6_1) res0_hi += dtmp0;

nop

nop

lda [%i4]%asi,%f10 ! (7_1) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f54,%f20,%f54 ! (4_1) dd *= dtmp0;

nop

nop

lda [%i4+4]%asi,%f11 ! (7_1) ((float*)&x0)[1] = ((float*)px)[1];

fsubd DONE,%f12,%f60 ! (2_1) dtmp0 = DONE - dtmp0;

fmuld %f50,%f2,%f46 ! (6_1) dtmp1 *= y_lo0;

nop

lda [%i3]%asi,%f12 ! (7_1) ((float*)&y0)[0] = ((float*)py)[0];

fsubd DTWO,%f22,%f22 ! (3_1) dtmp2 = DTWO - dtmp2;

add %i1,stridex,%i4 ! px += stridex

nop

lda [%i3+4]%asi,%f13 ! (7_1) ((float*)&y0)[1] = ((float*)py)[1];

bn,pn %icc,.exit

fmuld %f54,%f24,%f50 ! (4_1) dtmp1 = dd * dres;

add %i4,stridex,%i1 ! px += stridex

ld [%fp+ftmp0],%o2 ! (5_1) iarr = ((int*)&dres)[0];

fand %f16,DA1,%f2 ! (5_1) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (7_1) x0 *= scl0;

nop

ldd [%fp+dtmp2],%f0 ! (0_1) *(long long*)&scl0 = ll;

fsubd %f60,%f40,%f20 ! (2_1) dtmp0 -= dtmp1;

fmuld %f12,%f62,%f60 ! (7_1) y0 *= scl0;

sra %o2,11,%i3 ! (5_1) iarr >>= 11;

nop

faddd %f28,%f46,%f40 ! (6_1) res0_lo += dtmp1;

and %i3,0x1fc,%i3 ! (5_1) iarr &= 0x1fc;

nop

bn,pn %icc,.exit

fmuld %f18,%f22,%f28 ! (3_1) dres = dd * dtmp2;

fsqrtd %f52,%f22 ! (1_1) res0 = sqrt ( res0 );

lda [%i4]0x82,%o1 ! (1_0) hx0 = *(int*)px;

add %i3,TBL,%g1 ! (5_1) (char*)dll1 + iarr

fsubd DTWO,%f50,%f62 ! (4_1) dtmp1 = DTWO - dtmp1;

fmuld %f20,%f26,%f52 ! (2_1) dtmp0 *= dres;

add %i0,stridey,%i3 ! py += stridey

ld [%g1],%f26 ! (5_1) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f46 ! (7_1) x_hi0 = x0 + D2ON36;

nop

add %i3,stridey,%i0 ! py += stridey

lda [%i3]0x82,%g1 ! (1_0) hy0 = *(int*)py;

faddd %f60,D2ON36,%f50 ! (7_1) y_hi0 = y0 + D2ON36;

faddd %f30,%f40,%f18 ! (6_1) dres = res0_hi + res0_lo;

and %o1,_0x7fffffff,%o7 ! (1_0) hx0 &= 0x7fffffff;

st %f18,[%fp+ftmp0] ! (6_1) iarr = ((int*)&dres)[0];

fmuld %f0,%f14,%f0 ! (0_1) res0 = scl0 * res0;

fmuld %f54,%f62,%f14 ! (4_1) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (1_0) hx0 ? 0x7ff00000

st %f0,[%i5] ! (0_1) ((float*)pz)[0] = ((float*)&res0)[0];

fpsub32 %f26,%f2,%f26 ! (5_1) dd = vis_fpsub32(dtmp0, dexp0);

and %g1,_0x7fffffff,%l7 ! (1_0) hy0 &= 0x7fffffff;

nop

bge,pn %icc,.update33 ! (1_0) if ( hx0 >= 0x7ff00000 )

fsubd %f46,D2ON36,%f20 ! (7_1) x_hi0 -= D2ON36;

cmp %l7,_0x7ff00000 ! (1_0) hy0 ? 0x7ff00000

sub %l7,%o7,%o1 ! (1_0) diff0 = hy0 - hx0;

st %f1,[%i5+4] ! (0_1) ((float*)pz)[1] = ((float*)&res0)[1];

fsubd %f50,D2ON36,%f54 ! (7_1) y_hi0 -= D2ON36;

fmuld %f26,%f16,%f50 ! (5_1) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (1_0) j0 = diff0 >> 31;

bge,pn %icc,.update34 ! (1_0) if ( hy0 >= 0x7ff00000 )

faddd %f48,%f52,%f52 ! (2_1) res0 += dtmp0;

and %o1,%o3,%o1 ! (1_0) j0 &= diff0;

add %i5,stridez,%i5 ! pz += stridez

stx %o4,[%fp+dtmp2] ! (0_0) *(long long*)&scl0 = ll;

fand %f28,DA0,%f48 ! (3_1) res0 = vis_fand(dres,DA0);

fmuld %f20,%f20,%f2 ! (7_1) res0_hi = x_hi0 * x_hi0;

sub %l7,%o1,%o4 ! (1_0) j0 = hy0 - j0;

cmp %o7,_0x00100000 ! (1_0) hx0 ? 0x00100000

fsubd %f10,%f20,%f0 ! (7_1) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (7_1) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (1_0) j0 &= 0x7ff00000;

bl,pn %icc,.update35 ! (1_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (7_1) res0_lo = x0 + x_hi0;

.cont35a:

fmuld %f44,%f48,%f10 ! (3_1) dtmp0 = res0_hi * res0;

nop

sub %l0,%o4,%o4 ! (1_0) j0 = 0x7ff00000 - j0;

fsubd DTWO,%f50,%f20 ! (5_1) dtmp0 = DTWO - dtmp0;

.cont35b:

fmuld %f14,%f24,%f24 ! (4_1) dtmp2 = dd * dres;

sllx %o4,32,%o4 ! (1_0) ll = (long long)j0 << 32;

stx %o4,[%fp+dtmp3] ! (1_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (7_1) dtmp1 = y0 + y_hi0;

fmuld %f38,%f48,%f38 ! (3_1) dtmp1 = res0_lo * res0;

nop

nop

fsubd %f60,%f54,%f12 ! (7_1) y_lo0 = y0 - y_hi0;

.cont36:

fmuld %f62,%f0,%f0 ! (7_1) res0_lo *= x_lo0;

nop

ldd [%fp+dtmp1],%f62 ! (0_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f44 ! (7_1) res0_hi += dtmp0;

fsubd DONE,%f10,%f60 ! (3_1) dtmp0 = DONE - dtmp0;

nop

lda [%i2]%asi,%f10 ! (0_0) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f26,%f20,%f54 ! (5_1) dd *= dtmp0;

nop

nop

lda [%i2+4]%asi,%f11 ! (0_0) ((float*)&x0)[1] = ((float*)px)[1];

bn,pn %icc,.exit

fmuld %f50,%f12,%f26 ! (7_1) dtmp1 *= y_lo0;

nop

lda [%o0]%asi,%f12 ! (0_0) ((float*)&y0)[0] = ((float*)py)[0];

fsubd DTWO,%f24,%f24 ! (4_1) dtmp2 = DTWO - dtmp2;

nop

nop

lda [%o0+4]%asi,%f13 ! (0_0) ((float*)&y0)[1] = ((float*)py)[1];

bn,pn %icc,.exit

fmuld %f54,%f16,%f46 ! (5_1) dtmp1 = dd * dres;

nop

ld [%fp+ftmp0],%o2 ! (6_1) iarr = ((int*)&dres)[0];

fand %f18,DA1,%f2 ! (6_1) dexp0 = vis_fand(dres,DA1);

fmuld %f10,%f62,%f10 ! (0_0) x0 *= scl0;

nop

ldd [%fp+dtmp4],%f50 ! (1_1) *(long long*)&scl0 = ll;

fsubd %f60,%f38,%f20 ! (3_1) dtmp0 -= dtmp1;

fmuld %f12,%f62,%f60 ! (0_0) y0 *= scl0;

sra %o2,11,%g1 ! (6_1) iarr >>= 11;

nop

faddd %f0,%f26,%f38 ! (7_1) res0_lo += dtmp1;

nop

and %g1,0x1fc,%g1 ! (6_1) iarr &= 0x1fc;

bn,pn %icc,.exit

fmuld %f14,%f24,%f26 ! (4_1) dres = dd * dtmp2;

fsqrtd %f52,%f24 ! (2_1) res0 = sqrt ( res0 );

lda [%i1]0x82,%o1 ! (2_0) hx0 = *(int*)px;

add %g1,TBL,%g1 ! (6_1) (char*)dll1 + iarr

fsubd DTWO,%f46,%f62 ! (5_1) dtmp1 = DTWO - dtmp1;

fmuld %f20,%f28,%f52 ! (3_1) dtmp0 *= dres;

mov %i1,%i2

ld [%g1],%f28 ! (6_1) dtmp0 = ((double*)((char*)dll1 + iarr))[0];

faddd %f10,D2ON36,%f46 ! (0_0) x_hi0 = x0 + D2ON36;

nop

mov %i0,%o0

lda [%i0]0x82,%g1 ! (2_0) hy0 = *(int*)py;

faddd %f60,D2ON36,%f12 ! (0_0) y_hi0 = y0 + D2ON36;

faddd %f44,%f38,%f14 ! (7_1) dres = res0_hi + res0_lo;

and %o1,_0x7fffffff,%o7 ! (2_0) hx0 &= 0x7fffffff;

st %f14,[%fp+ftmp0] ! (7_1) iarr = ((int*)&dres)[0];

fmuld %f50,%f22,%f0 ! (1_1) res0 = scl0 * res0;

fmuld %f54,%f62,%f22 ! (5_1) dd *= dtmp1;

cmp %o7,_0x7ff00000 ! (2_0) hx0 ? 0x7ff00000

st %f0,[%i5] ! (1_1) ((float*)pz)[0] = ((float*)&res0)[0];

fpsub32 %f28,%f2,%f28 ! (6_1) dd = vis_fpsub32(dtmp0, dexp0);

and %g1,_0x7fffffff,%l7 ! (2_0) hx0 &= 0x7fffffff;

nop

bge,pn %icc,.update37 ! (2_0) if ( hx0 >= 0x7ff00000 )

fsubd %f46,D2ON36,%f20 ! (0_0) x_hi0 -= D2ON36;

sub %l7,%o7,%o1 ! (2_0) diff0 = hy0 - hx0;

cmp %l7,_0x7ff00000 ! (2_0) hy0 ? 0x7ff00000

st %f1,[%i5+4] ! (1_1) ((float*)pz)[1] = ((float*)&res0)[1];

fsubd %f12,D2ON36,%f54 ! (0_0) y_hi0 -= D2ON36;

fmuld %f28,%f18,%f50 ! (6_1) dtmp0 = dd * dres;

sra %o1,31,%o3 ! (2_0) j0 = diff0 >> 31;

bge,pn %icc,.update38 ! (2_0) if ( hy0 >= 0x7ff00000 )

faddd %f48,%f52,%f52 ! (3_1) res0 += dtmp0;

and %o1,%o3,%o1 ! (2_0) j0 &= diff0;

add %i5,stridez,%i5 ! pz += stridez

stx %o4,[%fp+dtmp4] ! (1_0) *(long long*)&scl0 = ll;

fand %f26,DA0,%f48 ! (4_1) res0 = vis_fand(dres,DA0);

fmuld %f20,%f20,%f2 ! (0_0) res0_hi = x_hi0 * x_hi0;

cmp %o7,_0x00100000 ! (2_0) hx0 ? 0x00100000

sub %l7,%o1,%o4 ! (2_0) j0 = hy0 - j0;

fsubd %f10,%f20,%f0 ! (0_0) x_lo0 = x0 - x_hi0;

fmuld %f54,%f54,%f46 ! (0_0) dtmp0 = y_hi0 * y_hi0;

and %o4,%l0,%o4 ! (2_0) j0 &= 0x7ff00000;

bl,pn %icc,.update39 ! (2_0) if ( hx0 < 0x00100000 )

faddd %f10,%f20,%f62 ! (0_0) res0_lo = x0 + x_hi0;

.cont39a:

fmuld %f32,%f48,%f10 ! (4_1) dtmp0 = res0_hi * res0;

sub %l0,%o4,%g1 ! (2_0) j0 = 0x7ff00000 - j0;

nop

fsubd DTWO,%f50,%f20 ! (6_1) dtmp0 = DTWO - dtmp0;

.cont39b:

fmuld %f22,%f16,%f16 ! (5_1) dtmp2 = dd * dres;

sllx %g1,32,%g1 ! (2_0) ll = (long long)j0 << 32;

stx %g1,[%fp+dtmp5] ! (2_0) *(long long*)&scl0 = ll;

faddd %f60,%f54,%f50 ! (0_0) dtmp1 = y0 + y_hi0;

fmuld %f36,%f48,%f36 ! (4_1) dtmp1 = res0_lo * res0;

nop

nop

fsubd %f60,%f54,%f12 ! (0_0) y_lo0 = y0 - y_hi0;

.cont40:

fmuld %f62,%f0,%f0 ! (0_0) res0_lo *= x_lo0;

nop

ldd [%fp+dtmp3],%f62 ! (1_0) *(long long*)&scl0 = ll;

faddd %f2,%f46,%f32 ! (0_0) res0_hi += dtmp0;

fsubd DONE,%f10,%f60 ! (4_1) dtmp0 = DONE - dtmp0;

nop

lda [%i4]%asi,%f10 ! (1_0) ((float*)&x0)[0] = ((float*)px)[0];

fmuld %f28,%f20,%f54 ! (6_1) dd *= dtmp0;