.file "__vatanf.S"

#include "libm.h"

RO_DATA

.align 64

.CONST_TBL:

.word 0x3fefffff, 0xfffccbbc ! K0 = 9.99999999976686608841e-01

.word 0xbfd55554, 0x51c6b90f ! K1 = -3.33333091601972730504e-01

.word 0x3fc98d6d, 0x926596cc ! K2 = 1.99628540499523379702e-01

.word 0x00020000, 0x00000000 ! DC1

.word 0xfffc0000, 0x00000000 ! DC2

.word 0x7ff00000, 0x00000000 ! DC3

.word 0x3ff00000, 0x00000000 ! DONE = 1.0

.word 0x40000000, 0x00000000 ! DTWO = 2.0

! parr0 = *(int*)&(1.0 / *(double*)&(((long long)i << 45) | 0x3ff0100000000000ULL)) + 0x3ff00000, i = [0, 127]

.word 0x7fdfe01f, 0x7fdfa11c, 0x7fdf6310, 0x7fdf25f6

.word 0x7fdee9c7, 0x7fdeae80, 0x7fde741a, 0x7fde3a91

.word 0x7fde01e0, 0x7fddca01, 0x7fdd92f2, 0x7fdd5cac

.word 0x7fdd272c, 0x7fdcf26e, 0x7fdcbe6d, 0x7fdc8b26

.word 0x7fdc5894, 0x7fdc26b5, 0x7fdbf583, 0x7fdbc4fd

.word 0x7fdb951e, 0x7fdb65e2, 0x7fdb3748, 0x7fdb094b

.word 0x7fdadbe8, 0x7fdaaf1d, 0x7fda82e6, 0x7fda5741

.word 0x7fda2c2a, 0x7fda01a0, 0x7fd9d79f, 0x7fd9ae24

.word 0x7fd9852f, 0x7fd95cbb, 0x7fd934c6, 0x7fd90d4f

.word 0x7fd8e652, 0x7fd8bfce, 0x7fd899c0, 0x7fd87427

.word 0x7fd84f00, 0x7fd82a4a, 0x7fd80601, 0x7fd7e225

.word 0x7fd7beb3, 0x7fd79baa, 0x7fd77908, 0x7fd756ca

.word 0x7fd734f0, 0x7fd71378, 0x7fd6f260, 0x7fd6d1a6

.word 0x7fd6b149, 0x7fd69147, 0x7fd6719f, 0x7fd6524f

.word 0x7fd63356, 0x7fd614b3, 0x7fd5f664, 0x7fd5d867

.word 0x7fd5babc, 0x7fd59d61, 0x7fd58056, 0x7fd56397

.word 0x7fd54725, 0x7fd52aff, 0x7fd50f22, 0x7fd4f38f

.word 0x7fd4d843, 0x7fd4bd3e, 0x7fd4a27f, 0x7fd48805

.word 0x7fd46dce, 0x7fd453d9, 0x7fd43a27, 0x7fd420b5

.word 0x7fd40782, 0x7fd3ee8f, 0x7fd3d5d9, 0x7fd3bd60

.word 0x7fd3a524, 0x7fd38d22, 0x7fd3755b, 0x7fd35dce

.word 0x7fd34679, 0x7fd32f5c, 0x7fd31877, 0x7fd301c8

.word 0x7fd2eb4e, 0x7fd2d50a, 0x7fd2bef9, 0x7fd2a91c

.word 0x7fd29372, 0x7fd27dfa, 0x7fd268b3, 0x7fd2539d

.word 0x7fd23eb7, 0x7fd22a01, 0x7fd21579, 0x7fd20120

.word 0x7fd1ecf4, 0x7fd1d8f5, 0x7fd1c522, 0x7fd1b17c

.word 0x7fd19e01, 0x7fd18ab0, 0x7fd1778a, 0x7fd1648d

.word 0x7fd151b9, 0x7fd13f0e, 0x7fd12c8b, 0x7fd11a30

.word 0x7fd107fb, 0x7fd0f5ed, 0x7fd0e406, 0x7fd0d244

.word 0x7fd0c0a7, 0x7fd0af2f, 0x7fd09ddb, 0x7fd08cab

.word 0x7fd07b9f, 0x7fd06ab5, 0x7fd059ee, 0x7fd04949

.word 0x7fd038c6, 0x7fd02864, 0x7fd01824, 0x7fd00804

.word 0x3ff00000, 0x00000000 ! 1.0

.word 0xbff00000, 0x00000000 ! -1.0

! parr1[i] = atan((double)*(float*)&((i + 460) << 21)), i = [0, 155]

.word 0x3f2fffff, 0xf555555c, 0x3f33ffff, 0xf595555f

.word 0x3f37ffff, 0xee000018, 0x3f3bffff, 0xe36aaadf

.word 0x3f3fffff, 0xd55555bc, 0x3f43ffff, 0xd65555f2

.word 0x3f47ffff, 0xb8000185, 0x3f4bffff, 0x8daaadf3

.word 0x3f4fffff, 0x55555bbc, 0x3f53ffff, 0x59555f19

.word 0x3f57fffe, 0xe000184d, 0x3f5bfffe, 0x36aadf30

.word 0x3f5ffffd, 0x5555bbbc, 0x3f63fffd, 0x6555f195

.word 0x3f67fffb, 0x800184cc, 0x3f6bfff8, 0xdaadf302

.word 0x3f6ffff5, 0x555bbbb7, 0x3f73fff5, 0x955f194a

.word 0x3f77ffee, 0x00184ca6, 0x3f7bffe3, 0x6adf2fd1

.word 0x3f7fffd5, 0x55bbba97, 0x3f83ffd6, 0x55f1929c

.word 0x3f87ffb8, 0x0184c30a, 0x3f8bff8d, 0xadf2e78c

.word 0x3f8fff55, 0x5bbb729b, 0x3f93ff59, 0x5f18a700

.word 0x3f97fee0, 0x184a5c36, 0x3f9bfe36, 0xdf291712

.word 0x3f9ffd55, 0xbba97625, 0x3fa3fd65, 0xf169c9d9

.word 0x3fa7fb81, 0x8430da2a, 0x3fabf8dd, 0xf139c444

.word 0x3faff55b, 0xb72cfdea, 0x3fb3f59f, 0x0e7c559d

.word 0x3fb7ee18, 0x2602f10f, 0x3fbbe39e, 0xbe6f07c4

.word 0x3fbfd5ba, 0x9aac2f6e, 0x3fc3d6ee, 0xe8c6626c

.word 0x3fc7b97b, 0x4bce5b02, 0x3fcb90d7, 0x529260a2

.word 0x3fcf5b75, 0xf92c80dd, 0x3fd36277, 0x3707ebcc

.word 0x3fd6f619, 0x41e4def1, 0x3fda64ee, 0xc3cc23fd

.word 0x3fddac67, 0x0561bb4f, 0x3fe1e00b, 0xabdefeb4

.word 0x3fe4978f, 0xa3269ee1, 0x3fe700a7, 0xc5784634

.word 0x3fe921fb, 0x54442d18, 0x3fecac7c, 0x57846f9e

.word 0x3fef730b, 0xd281f69b, 0x3ff0d38f, 0x2c5ba09f

.word 0x3ff1b6e1, 0x92ebbe44, 0x3ff30b6d, 0x796a4da8

.word 0x3ff3fc17, 0x6b7a8560, 0x3ff4ae10, 0xfc6589a5

.word 0x3ff5368c, 0x951e9cfd, 0x3ff5f973, 0x15254857

.word 0x3ff67d88, 0x63bc99bd, 0x3ff6dcc5, 0x7bb565fd

.word 0x3ff7249f, 0xaa996a21, 0x3ff789bd, 0x2c160054

.word 0x3ff7cd6f, 0x6dc59db4, 0x3ff7fde8, 0x0870c2a0

.word 0x3ff82250, 0x768ac529, 0x3ff8555a, 0x2787981f

.word 0x3ff87769, 0xeb8e956b, 0x3ff88fc2, 0x18ace9dc

.word 0x3ff8a205, 0xfd558740, 0x3ff8bb9a, 0x63718f45

.word 0x3ff8cca9, 0x27cf0b3d, 0x3ff8d8d8, 0xbf65316f

.word 0x3ff8e1fc, 0xa98cb633, 0x3ff8eec8, 0xcfd00665

.word 0x3ff8f751, 0x0eba96e6, 0x3ff8fd69, 0x4acf36b0

.word 0x3ff901fb, 0x7eee715e, 0x3ff90861, 0xd082d9b5

.word 0x3ff90ca6, 0x0b9322c5, 0x3ff90fb2, 0x37a7ea27

.word 0x3ff911fb, 0x59997f3a, 0x3ff9152e, 0x8a326c38

.word 0x3ff91750, 0xab2e0d12, 0x3ff918d6, 0xc2f9c9e2

.word 0x3ff919fb, 0x54eed7a9, 0x3ff91b94, 0xee352849

.word 0x3ff91ca5, 0xff216922, 0x3ff91d69, 0x0b3f72ff

.word 0x3ff91dfb, 0x5459826d, 0x3ff91ec8, 0x211be619

.word 0x3ff91f50, 0xa99fd49a, 0x3ff91fb2, 0x2fb5defa

.word 0x3ff91ffb, 0x5446d7c3, 0x3ff92061, 0xbaabf105

.word 0x3ff920a5, 0xfeefa208, 0x3ff920d6, 0xc1fb87e7

.word 0x3ff920fb, 0x5444826e, 0x3ff9212e, 0x87778bfc

.word 0x3ff92150, 0xa9999bb6, 0x3ff92169, 0x0b1faabb

.word 0x3ff9217b, 0x544437c3, 0x3ff92194, 0xedddcc28

.word 0x3ff921a5, 0xfeeedaec, 0x3ff921b2, 0x2fb1e5f1

.word 0x3ff921bb, 0x54442e6e, 0x3ff921c8, 0x2110fa94

.word 0x3ff921d0, 0xa99982d3, 0x3ff921d6, 0xc1fb08c6

.word 0x3ff921db, 0x54442d43, 0x3ff921e1, 0xbaaa9395

.word 0x3ff921e5, 0xfeeed7d0, 0x3ff921e9, 0x0b1f9ad7

.word 0x3ff921eb, 0x54442d1e, 0x3ff921ee, 0x8777604e

.word 0x3ff921f0, 0xa999826f, 0x3ff921f2, 0x2fb1e3f5

.word 0x3ff921f3, 0x54442d19, 0x3ff921f4, 0xedddc6b2

.word 0x3ff921f5, 0xfeeed7c3, 0x3ff921f6, 0xc1fb0886

.word 0x3ff921f7, 0x54442d18, 0x3ff921f8, 0x2110f9e5

.word 0x3ff921f8, 0xa999826e, 0x3ff921f9, 0x0b1f9acf

.word 0x3ff921f9, 0x54442d18, 0x3ff921f9, 0xbaaa937f

.word 0x3ff921f9, 0xfeeed7c3, 0x3ff921fa, 0x2fb1e3f4

.word 0x3ff921fa, 0x54442d18, 0x3ff921fa, 0x8777604b

.word 0x3ff921fa, 0xa999826e, 0x3ff921fa, 0xc1fb0886

.word 0x3ff921fa, 0xd4442d18, 0x3ff921fa, 0xedddc6b2

.word 0x3ff921fa, 0xfeeed7c3, 0x3ff921fb, 0x0b1f9acf

.word 0x3ff921fb, 0x14442d18, 0x3ff921fb, 0x2110f9e5

.word 0x3ff921fb, 0x2999826e, 0x3ff921fb, 0x2fb1e3f4

.word 0x3ff921fb, 0x34442d18, 0x3ff921fb, 0x3aaa937f

.word 0x3ff921fb, 0x3eeed7c3, 0x3ff921fb, 0x41fb0886

.word 0x3ff921fb, 0x44442d18, 0x3ff921fb, 0x4777604b

.word 0x3ff921fb, 0x4999826e, 0x3ff921fb, 0x4b1f9acf

.word 0x3ff921fb, 0x4c442d18, 0x3ff921fb, 0x4dddc6b2

.word 0x3ff921fb, 0x4eeed7c3, 0x3ff921fb, 0x4fb1e3f4

.word 0x3ff921fb, 0x50442d18, 0x3ff921fb, 0x5110f9e5

.word 0x3ff921fb, 0x5199826e, 0x3ff921fb, 0x51fb0886

#define DC2 %f2

#define DTWO %f6

#define DONE %f52

#define K0 %f54

#define K1 %f56

#define K2 %f58

#define DC1 %f60

#define DC3 %f62

#define stridex %o2

#define stridey %o3

#define MASK_0x7fffffff %i1

#define MASK_0x100000 %i5

#define tmp_px STACK_BIAS-32

#define tmp_counter STACK_BIAS-24

#define tmp0 STACK_BIAS-16

#define tmp1 STACK_BIAS-8

#define counter %l1

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

#define tmps 0x20

!--------------------------------------------------------------------

! !!!!! vatanf algorithm !!!!!

! ux = ((int*)px)[0];

! ax = ux & 0x7fffffff;

!

! if ( ax < 0x39b89c55 )

! {

! *(int*)py = ux;

! goto next;

! }

!

! if ( ax > 0x4c700518 )

! {

! if ( ax > 0x7f800000 )

! {

! float fpx = fabsf(*px);

! fpx *= fpx;

! *py = fpx;

! goto next;

! }

!

! sign = ux & 0x80000000;

! sign |= pi_2;

! *(int*)py = sign;

! goto next;

! }

!

! ftmp0 = *px;

! x = (double)ftmp0;

! px += stridex;

! y = vis_fpadd32(x,DC1);

! y = vis_fand(y,DC2);

! div = x * y;

! xx = x - y;

! div += DONE;

! i = ((unsigned long long*)&div)[0];

! y0 = vis_fand(div,DC3);

! i >>= 43;

! i &= 508;

! *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

! y0 = vis_fpsub32(dtmp0, y0);

! dtmp0 = div0 * y0;

! dtmp0 = DTWO - dtmp0;

! y0 *= dtmp0;

! dtmp1 = div0 * y0;

! dtmp1 = DTWO - dtmp1;

! y0 *= dtmp1;

! ax = ux & 0x7fffffff;

! ax += 0x00100000;

! ax >>= 18;

! ax &= -8;

! res = *(double*)((char*)parr1 + ax);

! ux >>= 28;

! ux &= -8;

! dtmp0 = *(double*)((char*)sign_arr + ux);

! res *= dtmp0;

! xx *= y0;

! x2 = xx * xx;

! dtmp0 = K2 * x2;

! dtmp0 += K1;

! dtmp0 *= x2;

! dtmp0 += K0;

! dtmp0 *= xx;

! res += dtmp0;

! ftmp0 = (float)res;

! py[0] = ftmp0;

! py += stridey;

!--------------------------------------------------------------------

ENTRY(__vatanf)

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

PIC_SETUP(l7)

PIC_SET(l7,.CONST_TBL,l2)

st %i0,[%fp+tmp_counter]

sllx %i2,2,stridex

sllx %i4,2,stridey

or %g0,%i3,%o1

stx %i1,[%fp+tmp_px]

ldd [%l2],K0

ldd [%l2+8],K1

ldd [%l2+16],K2

ldd [%l2+24],DC1

ldd [%l2+32],DC2

ldd [%l2+40],DC3

ldd [%l2+48],DONE

ldd [%l2+56],DTWO

add %l2,64,%i4

add %l2,64+512,%l0

add %l2,64+512+16-0x1cc*8,%l7

sethi %hi(0x100000),MASK_0x100000

sethi %hi(0x7ffffc00),MASK_0x7fffffff

add MASK_0x7fffffff,1023,MASK_0x7fffffff

sethi %hi(0x39b89c00),%o4

add %o4,0x55,%o4

sethi %hi(0x4c700400),%o5

add %o5,0x118,%o5

.begin:

ld [%fp+tmp_counter],counter

ldx [%fp+tmp_px],%i3

st %g0,[%fp+tmp_counter]

.begin1:

cmp counter,0

ble,pn %icc,.exit

nop

lda [%i3]0x82,%l6 ! (0_0) ux = ((int*)px)[0];

and %l6,MASK_0x7fffffff,%l5 ! (0_0) ax = ux & 0x7fffffff;

lda [%i3]0x82,%f0 ! (0_0) ftmp0 = *px;

cmp %l5,%o4 ! (0_0) ax ? 0x39b89c55

bl,pn %icc,.spec0 ! (0_0) if ( ax < 0x39b89c55 )

nop

cmp %l5,%o5 ! (0_0) ax ? 0x4c700518

bg,pn %icc,.spec1 ! (0_0) if ( ax > 0x4c700518 )

nop

add %i3,stridex,%l5 ! px += stridex;

fstod %f0,%f22 ! (0_0) ftmp0 = *px;

mov %l6,%i3

lda [%l5]0x82,%l6 ! (1_0) ux = ((int*)px)[0];

and %l6,MASK_0x7fffffff,%o7 ! (1_0) ax = ux & 0x7fffffff;

lda [%l5]0x82,%f0 ! (1_0) ftmp0 = *px;

add %l5,stridex,%l4 ! px += stridex;

fpadd32 %f22,DC1,%f24 ! (0_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (1_0) ax ? 0x39b89c55

bl,pn %icc,.update0 ! (1_0) if ( ax < 0x39b89c55 )

nop

.cont0:

cmp %o7,%o5 ! (1_0) ax ? 0x4c700518

bg,pn %icc,.update1 ! (1_0) if ( ax > 0x4c700518 )

nop

.cont1:

fstod %f0,%f20 ! (1_0) x = (double)ftmp0;

mov %l6,%l5

fand %f24,DC2,%f26 ! (0_0) y = vis_fand(y,dconst2);

fmuld %f22,%f26,%f32 ! (0_0) div = x * y;

lda [%l4]0x82,%l6 ! (2_0) ux = ((int*)px)[0];

fsubd %f22,%f26,%f22 ! (0_0) xx = x - y;

and %l6,MASK_0x7fffffff,%o7 ! (2_0) ax = ux & 0x7fffffff;

lda [%l4]0x82,%f0 ! (2_0) ftmp0 = *px;

add %l4,stridex,%l3 ! px += stridex;

fpadd32 %f20,DC1,%f24 ! (1_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (2_0) ax ? 0x39b89c55

bl,pn %icc,.update2 ! (2_0) if ( ax < 0x39b89c55 )

faddd DONE,%f32,%f32 ! (0_0) div += done;

.cont2:

cmp %o7,%o5 ! (2_0) ax ? 0x4c700518

bg,pn %icc,.update3 ! (2_0) if ( ax > 0x4c700518 )

nop

.cont3:

std %f32,[%fp+tmp0] ! (0_0) i = ((unsigned long long*)&div)[0];

mov %l6,%l4

fstod %f0,%f18 ! (2_0) x = (double)ftmp0;

fand %f24,DC2,%f26 ! (1_0) y = vis_fand(y,dconst2);

fmuld %f20,%f26,%f30 ! (1_0) div = x * y;

lda [%l3]0x82,%l6 ! (3_0) ux = ((int*)px)[0];

fsubd %f20,%f26,%f20 ! (1_0) xx = x - y;

and %l6,MASK_0x7fffffff,%o7 ! (3_0) ax = ux & 0x7fffffff;

lda [%l3]0x82,%f0 ! (3_0) ftmp0 = *px;

add %l3,stridex,%i0 ! px += stridex;

fpadd32 %f18,DC1,%f24 ! (2_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (3_0) ax ? 0x39b89c55

bl,pn %icc,.update4 ! (3_0) if ( ax < 0x39b89c55 )

faddd DONE,%f30,%f30 ! (1_0) div += done;

.cont4:

cmp %o7,%o5 ! (3_0) ax ? 0x4c700518

bg,pn %icc,.update5 ! (3_0) if ( ax > 0x4c700518 )

nop

.cont5:

std %f30,[%fp+tmp1] ! (1_0) i = ((unsigned long long*)&div)[0];

mov %l6,%l3

fstod %f0,%f16 ! (3_0) x = (double)ftmp0;

ldx [%fp+tmp0],%o0 ! (0_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (2_0) y = vis_fand(y,dconst2);

fand %f32,DC3,%f24 ! (0_0) y0 = vis_fand(div,dconst3);

srlx %o0,43,%o0 ! (0_0) i >>= 43;

and %o0,508,%l6 ! (0_0) i &= 508;

ld [%i4+%l6],%f0 ! (0_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fmuld %f18,%f26,%f28 ! (2_0) div = x * y;

lda [%i0]0x82,%l6 ! (4_0) ux = ((int*)px)[0];

fsubd %f18,%f26,%f18 ! (2_0) xx = x - y;

fpsub32 %f0,%f24,%f40 ! (0_0) y0 = vis_fpsub32(dtmp0, y0);

and %l6,MASK_0x7fffffff,%o7 ! (4_0) ax = ux & 0x7fffffff;

lda [%i0]0x82,%f0 ! (4_0) ftmp0 = *px;

add %i0,stridex,%i2 ! px += stridex;

fpadd32 %f16,DC1,%f24 ! (3_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (4_0) ax ? 0x39b89c55

bl,pn %icc,.update6 ! (4_0) if ( ax < 0x39b89c55 )

faddd DONE,%f28,%f28 ! (2_0) div += done;

.cont6:

fmuld %f32,%f40,%f42 ! (0_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (4_0) ax ? 0x4c700518

bg,pn %icc,.update7 ! (4_0) if ( ax > 0x4c700518 )

nop

.cont7:

std %f28,[%fp+tmp0] ! (2_0) i = ((unsigned long long*)&div)[0];

mov %l6,%i0

fstod %f0,%f14 ! (4_0) x = (double)ftmp0;

ldx [%fp+tmp1],%g1 ! (1_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (3_0) y = vis_fand(y,dconst2);

fand %f30,DC3,%f24 ! (1_0) y0 = vis_fand(div,dconst3);

fsubd DTWO,%f42,%f44 ! (0_0) dtmp0 = dtwo - dtmp0;

srlx %g1,43,%g1 ! (1_0) i >>= 43;

and %g1,508,%l6 ! (1_0) i &= 508;

ld [%i4+%l6],%f0 ! (1_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fmuld %f16,%f26,%f34 ! (3_0) div = x * y;

lda [%i2]0x82,%l6 ! (5_0) ux = ((int*)px)[0];

fsubd %f16,%f26,%f16 ! (3_0) xx = x - y;

fpsub32 %f0,%f24,%f38 ! (1_0) y0 = vis_fpsub32(dtmp0, y0);

add %i2,stridex,%l2 ! px += stridex;

fmuld %f40,%f44,%f40 ! (0_0) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (5_0) ax = ux & 0x7fffffff;

lda [%i2]0x82,%f0 ! (5_0) ftmp0 = *px;

fpadd32 %f14,DC1,%f24 ! (4_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (5_0) ax ? 0x39b89c55

bl,pn %icc,.update8 ! (5_0) if ( ax < 0x39b89c55 )

faddd DONE,%f34,%f34 ! (3_0) div += done;

.cont8:

fmuld %f30,%f38,%f42 ! (1_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (5_0) ax ? 0x4c700518

bg,pn %icc,.update9 ! (5_0) if ( ax > 0x4c700518 )

nop

.cont9:

std %f34,[%fp+tmp1] ! (3_0) i = ((unsigned long long*)&div)[0];

mov %l6,%i2

fstod %f0,%f36 ! (5_0) x = (double)ftmp0;

fmuld %f32,%f40,%f32 ! (0_0) dtmp1 = div0 * y0;

ldx [%fp+tmp0],%o0 ! (2_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (4_0) y = vis_fand(y,dconst2);

fand %f28,DC3,%f24 ! (2_0) y0 = vis_fand(div,dconst3);

fsubd DTWO,%f42,%f44 ! (1_0) dtmp0 = dtwo - dtmp0;

srlx %o0,43,%o0 ! (2_0) i >>= 43;

and %o0,508,%l6 ! (2_0) i &= 508;

fsubd DTWO,%f32,%f46 ! (0_0) dtmp1 = dtwo - dtmp1;

ld [%i4+%l6],%f0 ! (2_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fmuld %f14,%f26,%f32 ! (4_0) div = x * y;

lda [%l2]0x82,%l6 ! (6_0) ux = ((int*)px)[0];

fsubd %f14,%f26,%f14 ! (4_0) xx = x - y;

fmuld %f40,%f46,%f26 ! (0_0) y0 *= dtmp1;

add %l2,stridex,%g5 ! px += stridex;

fpsub32 %f0,%f24,%f40 ! (2_0) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (1_0) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (6_0) ax = ux & 0x7fffffff;

lda [%l2]0x82,%f0 ! (6_0) ftmp0 = *px;

fpadd32 %f36,DC1,%f24 ! (5_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (6_0) ax ? 0x39b89c55

bl,pn %icc,.update10 ! (6_0) if ( ax < 0x39b89c55 )

faddd DONE,%f32,%f32 ! (4_0) div += done;

.cont10:

fmuld %f28,%f40,%f42 ! (2_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (6_0) ax ? 0x4c700518

bg,pn %icc,.update11 ! (6_0) if ( ax > 0x4c700518 )

nop

.cont11:

fmuld %f22,%f26,%f22 ! (0_0) xx *= y0;

mov %l6,%l2

std %f32,[%fp+tmp0] ! (4_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f10 ! (6_0) x = (double)ftmp0;

fmuld %f30,%f38,%f30 ! (1_0) dtmp1 = div0 * y0;

ldx [%fp+tmp1],%g1 ! (3_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (5_0) y = vis_fand(y,dconst2);

fand %f34,DC3,%f24 ! (3_0) y0 = vis_fand(div,dconst3);

fmuld %f22,%f22,%f50 ! (0_0) x2 = xx * xx;

srlx %g1,43,%g1 ! (3_0) i >>= 43;

fsubd DTWO,%f42,%f44 ! (2_0) dtmp0 = dtwo - dtmp0;

and %g1,508,%l6 ! (3_0) i &= 508;

mov %i3,%o7

fsubd DTWO,%f30,%f46 ! (1_0) dtmp1 = dtwo - dtmp1;

ld [%i4+%l6],%f0 ! (3_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fmuld %f36,%f26,%f30 ! (5_0) div = x * y;

srl %o7,28,%g1 ! (0_0) ux >>= 28;

add %g5,stridex,%i3 ! px += stridex;

fmuld K2,%f50,%f4 ! (0_0) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o0 ! (0_0) ax = ux & 0x7fffffff;

lda [%g5]0x82,%l6 ! (7_0) ux = ((int*)px)[0];

fsubd %f36,%f26,%f36 ! (5_0) xx = x - y;

fmuld %f38,%f46,%f26 ! (1_0) y0 *= dtmp1;

add %o0,MASK_0x100000,%o0 ! (0_0) ax += 0x00100000;

and %g1,-8,%g1 ! (0_0) ux &= -8;

fpsub32 %f0,%f24,%f38 ! (3_0) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f40,%f44,%f40 ! (2_0) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (7_0) ax = ux & 0x7fffffff;

lda [%g5]0x82,%f0 ! (7_0) ftmp0 = *px;

fpadd32 %f10,DC1,%f24 ! (6_0) y = vis_fpadd32(x,dconst1);

cmp %o7,%o4 ! (7_0) ax ? 0x39b89c55

bl,pn %icc,.update12 ! (7_0) if ( ax < 0x39b89c55 )

faddd DONE,%f30,%f30 ! (5_0) div += done;

.cont12:

fmuld %f34,%f38,%f42 ! (3_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (7_0) ax ? 0x4c700518

bg,pn %icc,.update13 ! (7_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (0_0) dtmp0 += K1;

.cont13:

fmuld %f20,%f26,%f20 ! (1_0) xx *= y0;

srl %o0,18,%o7 ! (0_0) ax >>= 18;

std %f30,[%fp+tmp1] ! (5_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f8 ! (7_0) x = (double)ftmp0;

fmuld %f28,%f40,%f28 ! (2_0) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (0_0) ux &= -8;

ldx [%fp+tmp0],%o0 ! (4_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (6_0) y = vis_fand(y,dconst2);

add %o7,%l7,%o7 ! (0_0) (char*)parr1 + ax;

mov %l6,%g5

ldd [%l0+%g1],%f48 ! (0_0) dtmp0 = *(double*)((char*)sign_arr + ux);

fmuld %f4,%f50,%f4 ! (0_0) dtmp0 *= x2;

srlx %o0,43,%o0 ! (4_0) i >>= 43;

ldd [%o7],%f0 ! (0_0) res = *(double*)((char*)parr1 + ax);

fand %f32,DC3,%f24 ! (4_0) y0 = vis_fand(div,dconst3);

fmuld %f20,%f20,%f50 ! (1_0) x2 = xx * xx;

and %o0,508,%l6 ! (4_0) i &= 508;

mov %l5,%o7

fsubd DTWO,%f42,%f44 ! (3_0) dtmp0 = dtwo - dtmp0;

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

fmuld %f0,%f48,%f48 ! (0_0) res *= dtmp0;

srl %o7,28,%l5 ! (1_0) ux >>= 28;

ld [%i4+%l6],%f0 ! (4_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fmuld %f10,%f26,%f28 ! (6_0) div = x * y;

faddd %f4,K0,%f42 ! (0_0) dtmp0 += K0;

subcc counter,8,counter

bneg,pn %icc,.tail

or %g0,%o1,%o0

add %fp,tmp0,%g1

lda [%i3]0x82,%l6 ! (0_0) ux = ((int*)px)[0];

ba .main_loop

add %i3,stridex,%l5 ! px += stridex;

.align 16

.main_loop:

fsubd %f10,%f26,%f10 ! (6_1) xx = x - y;

and %o7,MASK_0x7fffffff,%o1 ! (1_1) ax = ux & 0x7fffffff;

st %f12,[%g1] ! (7_1) py[0] = ftmp0;

fmuld K2,%f50,%f4 ! (1_1) dtmp0 = K2 * x2;

fmuld %f40,%f46,%f26 ! (2_1) y0 *= dtmp1;

srl %o7,28,%o7 ! (1_0) ux >>= 28;

add %o1,MASK_0x100000,%g1 ! (1_1) ax += 0x00100000;

fpsub32 %f0,%f24,%f40 ! (4_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (3_1) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o1 ! (0_0) ax = ux & 0x7fffffff;

lda [%i3]0x82,%f0 ! (0_0) ftmp0 = *px;

fpadd32 %f8,DC1,%f24 ! (7_1) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f22,%f44 ! (0_1) dtmp0 *= xx;

cmp %o1,%o4 ! (0_0) ax ? 0x39b89c55

bl,pn %icc,.update14 ! (0_0) if ( ax < 0x39b89c55 )

faddd DONE,%f28,%f28 ! (6_1) div += done;

.cont14:

fmuld %f32,%f40,%f42 ! (4_1) dtmp0 = div0 * y0;

cmp %o1,%o5 ! (0_0) ax ? 0x4c700518

bg,pn %icc,.update15 ! (0_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (1_1) dtmp0 += K1;

.cont15:

fmuld %f18,%f26,%f18 ! (2_1) xx *= y0;

srl %g1,18,%o1 ! (1_1) ax >>= 18;

std %f28,[%fp+tmp0] ! (6_1) i = ((unsigned long long*)&div)[0];

fstod %f0,%f22 ! (0_0) ftmp0 = *px;

fmuld %f34,%f38,%f34 ! (3_1) dtmp1 = div0 * y0;

and %o1,-8,%o1 ! (1_1) ax &= -8;

ldx [%fp+tmp1],%g1 ! (5_1) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (7_1) y = vis_fand(y,dconst2);

ldd [%o1+%l7],%f0 ! (1_1) res = *(double*)((char*)parr1 + ax);

and %o7,-8,%o7 ! (1_1) ux &= -8;

mov %l6,%i3

faddd %f48,%f44,%f12 ! (0_1) res += dtmp0;

fmuld %f4,%f50,%f4 ! (1_1) dtmp0 *= x2;

nop

ldd [%l0+%o7],%f48 ! (1_1) dtmp0 = *(double*)((char*)sign_arr + ux);

fand %f30,DC3,%f24 ! (5_1) y0 = vis_fand(div,dconst3);

fmuld %f18,%f18,%f50 ! (2_1) x2 = xx * xx;

srlx %g1,43,%g1 ! (5_1) i >>= 43;

mov %l4,%o7

fsubd DTWO,%f42,%f44 ! (4_1) dtmp0 = dtwo - dtmp0;

and %g1,508,%l6 ! (5_1) i &= 508;

nop

bn,pn %icc,.exit

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

fmuld %f0,%f48,%f48 ! (1_1) res *= dtmp0;

add %o0,stridey,%g1 ! py += stridey;

ld [%i4+%l6],%f0 ! (5_1) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (0_1) ftmp0 = (float)res;

fmuld %f8,%f26,%f34 ! (7_1) div = x * y;

srl %o7,28,%o1 ! (2_1) ux >>= 28;

lda [%l5]0x82,%l6 ! (1_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (1_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (2_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (2_1) ax = ux & 0x7fffffff;

st %f12,[%o0] ! (0_1) py[0] = ftmp0;

fsubd %f8,%f26,%f8 ! (7_1) xx = x - y;

fmuld %f38,%f46,%f26 ! (3_1) y0 *= dtmp1;

add %l5,stridex,%l4 ! px += stridex;

add %o7,MASK_0x100000,%o0 ! (2_1) ax += 0x00100000;

fpsub32 %f0,%f24,%f38 ! (5_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f40,%f44,%f40 ! (4_1) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (1_0) ax = ux & 0x7fffffff;

lda [%l5]0x82,%f0 ! (1_0) ftmp0 = *px;

fpadd32 %f22,DC1,%f24 ! (0_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f20,%f44 ! (1_1) dtmp0 *= xx;

cmp %o7,%o4 ! (1_0) ax ? 0x39b89c55

bl,pn %icc,.update16 ! (1_0) if ( ax < 0x39b89c55 )

faddd DONE,%f34,%f34 ! (7_1) div += done;

.cont16:

fmuld %f30,%f38,%f42 ! (5_1) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (1_0) ax ? 0x4c700518

bg,pn %icc,.update17 ! (1_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (2_1) dtmp0 += K1;

.cont17:

fmuld %f16,%f26,%f16 ! (3_1) xx *= y0;

srl %o0,18,%o7 ! (2_1) ax >>= 18;

std %f34,[%fp+tmp1] ! (7_1) i = ((unsigned long long*)&div)[0];

fstod %f0,%f20 ! (1_0) x = (double)ftmp0;

fmuld %f32,%f40,%f32 ! (4_1) dtmp1 = div0 * y0;

ldx [%fp+tmp0],%o0 ! (6_1) i = ((unsigned long long*)&div)[0];

and %o1,-8,%o1 ! (2_1) ux &= -8;

fand %f24,DC2,%f26 ! (0_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (1_1) res += dtmp0;

and %o7,-8,%o7 ! (2_1) ax &= -8;

ldd [%l0+%o1],%f48 ! (2_1) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

ldd [%o7+%l7],%f0 ! (2_1) res = *(double*)((char*)parr1 + ax);

mov %l6,%l5

fmuld %f4,%f50,%f4 ! (2_1) dtmp0 *= x2;

fand %f28,DC3,%f24 ! (6_1) y0 = vis_fand(div,dconst3);

fmuld %f16,%f16,%f50 ! (3_1) x2 = xx * xx;

srlx %o0,43,%o0 ! (6_1) i >>= 43;

mov %l3,%o7

fsubd DTWO,%f42,%f44 ! (5_1) dtmp0 = dtwo - dtmp0;

and %o0,508,%l6 ! (6_1) i &= 508;

add %l4,stridex,%l3 ! px += stridex;

bn,pn %icc,.exit

fsubd DTWO,%f32,%f46 ! (4_1) dtmp1 = dtwo - dtmp1;

fmuld %f0,%f48,%f48 ! (2_1) res *= dtmp0;

add %g1,stridey,%o0 ! py += stridey;

ld [%i4+%l6],%f0 ! (6_1) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (1_1) ftmp0 = (float)res;

fmuld %f22,%f26,%f32 ! (0_0) div = x * y;

srl %o7,28,%o1 ! (3_1) ux >>= 28;

lda [%l4]0x82,%l6 ! (2_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (2_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (3_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (3_1) ax = ux & 0x7fffffff;

st %f12,[%g1] ! (1_1) py[0] = ftmp0;

fsubd %f22,%f26,%f22 ! (0_0) xx = x - y;

fmuld %f40,%f46,%f26 ! (4_1) y0 *= dtmp1;

add %o7,MASK_0x100000,%g1 ! (3_1) ax += 0x00100000;

and %o1,-8,%o1 ! (3_1) ux &= -8;

fpsub32 %f0,%f24,%f40 ! (6_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (5_1) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (2_0) ax = ux & 0x7fffffff;

lda [%l4]0x82,%f0 ! (2_0) ftmp0 = *px;

fpadd32 %f20,DC1,%f24 ! (1_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f18,%f44 ! (2_1) dtmp0 *= xx;

cmp %o7,%o4 ! (2_0) ax ? 0x39b89c55

bl,pn %icc,.update18 ! (2_0) if ( ax < 0x39b89c55 )

faddd DONE,%f32,%f32 ! (0_0) div += done;

.cont18:

fmuld %f28,%f40,%f42 ! (6_1) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (2_0) ax ? 0x4c700518

bg,pn %icc,.update19 ! (2_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (3_1) dtmp0 += K1;

.cont19:

fmuld %f14,%f26,%f14 ! (4_1) xx *= y0;

srl %g1,18,%o7 ! (3_1) ax >>= 18;

std %f32,[%fp+tmp0] ! (0_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f18 ! (2_0) x = (double)ftmp0;

fmuld %f30,%f38,%f30 ! (5_1) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (3_1) ax &= -8;

ldx [%fp+tmp1],%g1 ! (7_1) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (1_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (2_1) res += dtmp0;

mov %l6,%l4

ldd [%l0+%o1],%f48 ! (3_1) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

fmuld %f4,%f50,%f4 ! (3_1) dtmp0 *= x2;

ldd [%o7+%l7],%f0 ! (3_1) res = *(double*)((char*)parr1 + ax)

nop

fand %f34,DC3,%f24 ! (7_1) y0 = vis_fand(div,dconst3);

fmuld %f14,%f14,%f50 ! (4_1) x2 = xx * xx;

srlx %g1,43,%g1 ! (7_1) i >>= 43;

mov %i0,%o7

fsubd DTWO,%f42,%f44 ! (6_1) dtmp0 = dtwo - dtmp0;

and %g1,508,%l6 ! (7_1) i &= 508;

add %l3,stridex,%i0 ! px += stridex;

bn,pn %icc,.exit

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

fmuld %f0,%f48,%f48 ! (3_1) res *= dtmp0;

add %o0,stridey,%g1 ! py += stridey;

ld [%i4+%l6],%f0 ! (7_1) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (2_1) ftmp0 = (float)res;

fmuld %f20,%f26,%f30 ! (1_0) div = x * y;

srl %o7,28,%o1 ! (4_1) ux >>= 28;

lda [%l3]0x82,%l6 ! (3_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (3_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (4_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (4_1) ax = ux & 0x7fffffff;

st %f12,[%o0] ! (2_1) py[0] = ftmp0;

fsubd %f20,%f26,%f20 ! (1_0) xx = x - y;

fmuld %f38,%f46,%f26 ! (5_1) y0 *= dtmp1;

add %o7,MASK_0x100000,%o0 ! (4_1) ax += 0x00100000;

and %o1,-8,%o1 ! (4_1) ux &= -8;

fpsub32 %f0,%f24,%f38 ! (7_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f40,%f44,%f40 ! (6_1) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (3_0) ax = ux & 0x7fffffff;

lda [%l3]0x82,%f0 ! (3_0) ftmp0 = *px;

fpadd32 %f18,DC1,%f24 ! (2_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f16,%f44 ! (3_1) dtmp0 *= xx;

cmp %o7,%o4 ! (3_0) ax ? 0x39b89c55

bl,pn %icc,.update20 ! (3_0) if ( ax < 0x39b89c55 )

faddd DONE,%f30,%f30 ! (1_0) div += done;

.cont20:

fmuld %f34,%f38,%f42 ! (7_1) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (3_0) ax ? 0x4c700518

bg,pn %icc,.update21 ! (3_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (4_1) dtmp0 += K1;

.cont21:

fmuld %f36,%f26,%f36 ! (5_1) xx *= y0;

srl %o0,18,%o7 ! (4_1) ax >>= 18;

std %f30,[%fp+tmp1] ! (1_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f16 ! (3_0) x = (double)ftmp0;

fmuld %f28,%f40,%f28 ! (6_1) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (4_1) ax &= -8;

ldx [%fp+tmp0],%o0 ! (0_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (2_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (3_1) res += dtmp0;

nop

ldd [%l0+%o1],%f48 ! (4_1) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

ldd [%o7+%l7],%f0 ! (4_1) res = *(double*)((char*)parr1 + ax);

mov %l6,%l3

fmuld %f4,%f50,%f4 ! (4_1) dtmp0 *= x2;

fand %f32,DC3,%f24 ! (0_0) y0 = vis_fand(div,dconst3);

fmuld %f36,%f36,%f50 ! (5_1) x2 = xx * xx;

srlx %o0,43,%o0 ! (0_0) i >>= 43;

mov %i2,%o7

fsubd DTWO,%f42,%f44 ! (7_1) dtmp0 = dtwo - dtmp0;

and %o0,508,%l6 ! (0_0) i &= 508;

add %i0,stridex,%i2 ! px += stridex;

bn,pn %icc,.exit

fsubd DTWO,%f28,%f46 ! (6_1) dtmp1 = dtwo - dtmp1;

fmuld %f0,%f48,%f48 ! (4_1) res *= dtmp0;

add %g1,stridey,%o0 ! py += stridey;

ld [%i4+%l6],%f0 ! (0_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (3_1) ftmp0 = (float)res;

fmuld %f18,%f26,%f28 ! (2_0) div = x * y;

srl %o7,28,%o1 ! (5_1) ux >>= 28;

lda [%i0]0x82,%l6 ! (4_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (4_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (5_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (5_1) ax = ux & 0x7fffffff;

st %f12,[%g1] ! (3_1) py[0] = ftmp0;

fsubd %f18,%f26,%f18 ! (2_0) xx = x - y;

fmuld %f40,%f46,%f26 ! (6_1) y0 *= dtmp1;

add %o7,MASK_0x100000,%g1 ! (5_1) ax += 0x00100000;

and %o1,-8,%o1 ! (5_1) ux &= -8;

fpsub32 %f0,%f24,%f40 ! (0_0) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (7_1) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (4_0) ax = ux & 0x7fffffff;

lda [%i0]0x82,%f0 ! (4_0) ftmp0 = *px;

fpadd32 %f16,DC1,%f24 ! (3_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f14,%f44 ! (4_1) dtmp0 *= xx;

cmp %o7,%o4 ! (4_0) ax ? 0x39b89c55

bl,pn %icc,.update22 ! (4_0) if ( ax < 0x39b89c55 )

faddd DONE,%f28,%f28 ! (2_0) div += done;

.cont22:

fmuld %f32,%f40,%f42 ! (0_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (4_0) ax ? 0x4c700518

bg,pn %icc,.update23 ! (4_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (5_1) dtmp0 += K1;

.cont23:

fmuld %f10,%f26,%f10 ! (6_1) xx *= y0;

srl %g1,18,%o7 ! (5_1) ax >>= 18;

std %f28,[%fp+tmp0] ! (2_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f14 ! (4_0) x = (double)ftmp0;

fmuld %f34,%f38,%f34 ! (7_1) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (5_1) ax &= -8;

ldx [%fp+tmp1],%g1 ! (1_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (3_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (4_1) res += dtmp0;

mov %l6,%i0

ldd [%l0+%o1],%f48 ! (5_1) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

ldd [%o7+%l7],%f0 ! (5_1) res = *(double*)((char*)parr1 + ax);

nop

fmuld %f4,%f50,%f4 ! (5_1) dtmp0 *= x2;

fand %f30,DC3,%f24 ! (1_0) y0 = vis_fand(div,dconst3);

fmuld %f10,%f10,%f50 ! (6_1) x2 = xx * xx;

srlx %g1,43,%g1 ! (1_0) i >>= 43;

mov %l2,%o7

fsubd DTWO,%f42,%f44 ! (0_0) dtmp0 = dtwo - dtmp0;

and %g1,508,%l6 ! (1_0) i &= 508;

add %i2,stridex,%l2 ! px += stridex;

bn,pn %icc,.exit

fsubd DTWO,%f34,%f46 ! (7_1) dtmp1 = dtwo - dtmp1;

fmuld %f0,%f48,%f48 ! (5_1) res *= dtmp0;

add %o0,stridey,%g1 ! py += stridey;

ld [%i4+%l6],%f0 ! (1_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (4_1) ftmp0 = (float)res;

fmuld %f16,%f26,%f34 ! (3_0) div = x * y;

srl %o7,28,%o1 ! (6_1) ux >>= 28;

lda [%i2]0x82,%l6 ! (5_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (5_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (6_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (6_1) ax = ux & 0x7fffffff;

st %f12,[%o0] ! (4_1) py[0] = ftmp0;

fsubd %f16,%f26,%f16 ! (3_0) xx = x - y;

fmuld %f38,%f46,%f26 ! (7_1) y0 *= dtmp1;

add %o7,MASK_0x100000,%o0 ! (6_1) ax += 0x00100000;

and %o1,-8,%o1 ! (6_1) ux &= -8;

fpsub32 %f0,%f24,%f38 ! (1_0) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f40,%f44,%f40 ! (0_0) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (5_0) ax = ux & 0x7fffffff;

lda [%i2]0x82,%f0 ! (5_0) ftmp0 = *px;

fpadd32 %f14,DC1,%f24 ! (4_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f36,%f44 ! (5_1) dtmp0 *= xx;

cmp %o7,%o4 ! (5_0) ax ? 0x39b89c55

bl,pn %icc,.update24 ! (5_0) if ( ax < 0x39b89c55 )

faddd DONE,%f34,%f34 ! (3_0) div += done;

.cont24:

fmuld %f30,%f38,%f42 ! (1_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (5_0) ax ? 0x4c700518

bg,pn %icc,.update25 ! (5_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (6_1) dtmp0 += K1;

.cont25:

fmuld %f8,%f26,%f8 ! (7_1) xx *= y0;

srl %o0,18,%o7 ! (6_1) ax >>= 18;

std %f34,[%fp+tmp1] ! (3_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f36 ! (5_0) x = (double)ftmp0;

fmuld %f32,%f40,%f32 ! (0_0) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (6_1) ax &= -8;

ldx [%fp+tmp0],%o0 ! (2_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (4_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (5_1) res += dtmp0;

mov %l6,%i2

ldd [%l0+%o1],%f48 ! (6_1) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

ldd [%o7+%l7],%f0 ! (6_1) res = *(double*)((char*)parr1 + ax);

nop

fmuld %f4,%f50,%f4 ! (6_1) dtmp0 *= x2;

fand %f28,DC3,%f24 ! (2_0) y0 = vis_fand(div,dconst3);

fmuld %f8,%f8,%f50 ! (7_1) x2 = xx * xx;

srlx %o0,43,%o0 ! (2_0) i >>= 43;

mov %g5,%o7

fsubd DTWO,%f42,%f44 ! (1_0) dtmp0 = dtwo - dtmp0;

and %o0,508,%l6 ! (2_0) i &= 508;

add %l2,stridex,%g5 ! px += stridex;

bn,pn %icc,.exit

fsubd DTWO,%f32,%f46 ! (0_0) dtmp1 = dtwo - dtmp1;

fmuld %f0,%f48,%f48 ! (6_1) res *= dtmp0;

add %g1,stridey,%o0 ! py += stridey;

ld [%i4+%l6],%f0 ! (2_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (5_1) ftmp0 = (float)res;

fmuld %f14,%f26,%f32 ! (4_0) div = x * y;

srl %o7,28,%o1 ! (7_1) ux >>= 28;

lda [%l2]0x82,%l6 ! (6_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (6_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (7_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (7_1) ax = ux & 0x7fffffff;

st %f12,[%g1] ! (5_1) py[0] = ftmp0;

fsubd %f14,%f26,%f14 ! (4_0) xx = x - y;

fmuld %f40,%f46,%f26 ! (0_0) y0 *= dtmp1;

add %o7,MASK_0x100000,%g1 ! (7_1) ax += 0x00100000;

and %o1,-8,%o1 ! (7_1) ux &= -8;

fpsub32 %f0,%f24,%f40 ! (2_0) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (1_0) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (6_0) ax = ux & 0x7fffffff;

lda [%l2]0x82,%f0 ! (6_0) ftmp0 = *px;

fpadd32 %f36,DC1,%f24 ! (5_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f10,%f44 ! (6_1) dtmp0 *= xx;

cmp %o7,%o4 ! (6_0) ax ? 0x39b89c55

bl,pn %icc,.update26 ! (6_0) if ( ax < 0x39b89c55 )

faddd DONE,%f32,%f32 ! (4_0) div += done;

.cont26:

fmuld %f28,%f40,%f42 ! (2_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (6_0) ax ? 0x4c700518

bg,pn %icc,.update27 ! (6_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (7_1) dtmp0 += K1;

.cont27:

fmuld %f22,%f26,%f22 ! (0_0) xx *= y0;

srl %g1,18,%o7 ! (7_1) ax >>= 18;

std %f32,[%fp+tmp0] ! (4_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f10 ! (6_0) x = (double)ftmp0;

fmuld %f30,%f38,%f30 ! (1_0) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (7_1) ax &= -8;

ldx [%fp+tmp1],%g1 ! (3_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (5_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (6_1) res += dtmp0;

mov %l6,%l2

ldd [%l0+%o1],%f48 ! (7_1) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

ldd [%o7+%l7],%f0 ! (7_1) res = *(double*)((char*)parr1 + ax);

nop

fmuld %f4,%f50,%f4 ! (7_1) dtmp0 *= x2;

fand %f34,DC3,%f24 ! (3_0) y0 = vis_fand(div,dconst3);

fmuld %f22,%f22,%f50 ! (0_0) x2 = xx * xx;

srlx %g1,43,%g1 ! (3_0) i >>= 43;

mov %i3,%o7

fsubd DTWO,%f42,%f44 ! (2_0) dtmp0 = dtwo - dtmp0;

and %g1,508,%l6 ! (3_0) i &= 508;

add %g5,stridex,%i3 ! px += stridex;

bn,pn %icc,.exit

fsubd DTWO,%f30,%f46 ! (1_0) dtmp1 = dtwo - dtmp1;

fmuld %f0,%f48,%f48 ! (7_1) res *= dtmp0;

add %o0,stridey,%g1 ! py += stridey;

ld [%i4+%l6],%f0 ! (3_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (6_1) ftmp0 = (float)res;

fmuld %f36,%f26,%f30 ! (5_0) div = x * y;

srl %o7,28,%o1 ! (0_0) ux >>= 28;

lda [%g5]0x82,%l6 ! (7_0) ux = ((int*)px)[0];

faddd %f4,K0,%f42 ! (7_1) dtmp0 += K0;

fmuld K2,%f50,%f4 ! (0_0) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o7 ! (0_0) ax = ux & 0x7fffffff;

st %f12,[%o0] ! (6_1) py[0] = ftmp0;

fsubd %f36,%f26,%f36 ! (5_0) xx = x - y;

fmuld %f38,%f46,%f26 ! (1_0) y0 *= dtmp1;

add %o7,MASK_0x100000,%o0 ! (0_0) ax += 0x00100000;

and %o1,-8,%o1 ! (0_0) ux &= -8;

fpsub32 %f0,%f24,%f38 ! (3_0) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f40,%f44,%f40 ! (2_0) y0 *= dtmp0;

and %l6,MASK_0x7fffffff,%o7 ! (7_0) ax = ux & 0x7fffffff;

lda [%g5]0x82,%f0 ! (7_0) ftmp0 = *px;

fpadd32 %f10,DC1,%f24 ! (6_0) y = vis_fpadd32(x,dconst1);

fmuld %f42,%f8,%f44 ! (7_1) dtmp0 *= xx;

cmp %o7,%o4 ! (7_0) ax ? 0x39b89c55

bl,pn %icc,.update28 ! (7_0) if ( ax < 0x39b89c55 )

faddd DONE,%f30,%f30 ! (5_0) div += done;

.cont28:

fmuld %f34,%f38,%f42 ! (3_0) dtmp0 = div0 * y0;

cmp %o7,%o5 ! (7_0) ax ? 0x4c700518

bg,pn %icc,.update29 ! (7_0) if ( ax > 0x4c700518 )

faddd %f4,K1,%f4 ! (0_0) dtmp0 += K1;

.cont29:

fmuld %f20,%f26,%f20 ! (1_0) xx *= y0;

srl %o0,18,%o7 ! (0_0) ax >>= 18;

std %f30,[%fp+tmp1] ! (5_0) i = ((unsigned long long*)&div)[0];

fstod %f0,%f8 ! (7_0) x = (double)ftmp0;

fmuld %f28,%f40,%f28 ! (2_0) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (0_0) ux &= -8;

ldx [%fp+tmp0],%o0 ! (4_0) i = ((unsigned long long*)&div)[0];

fand %f24,DC2,%f26 ! (6_0) y = vis_fand(y,dconst2);

faddd %f48,%f44,%f12 ! (7_1) res += dtmp0;

subcc counter,8,counter

ldd [%l0+%o1],%f48 ! (0_0) dtmp0 = *(double*)((char*)sign_arr + ux);

bn,pn %icc,.exit

fmuld %f4,%f50,%f4 ! (0_0) dtmp0 *= x2;

mov %l6,%g5

ldd [%o7+%l7],%f0 ! (0_0) res = *(double*)((char*)parr1 + ax);

fand %f32,DC3,%f24 ! (4_0) y0 = vis_fand(div,dconst3);

fmuld %f20,%f20,%f50 ! (1_0) x2 = xx * xx;

srlx %o0,43,%l6 ! (4_0) i >>= 43;

mov %l5,%o7

fsubd DTWO,%f42,%f44 ! (3_0) dtmp0 = dtwo - dtmp0;

add %g1,stridey,%o0 ! py += stridey;

and %l6,508,%l6 ! (4_0) i &= 508;

bn,pn %icc,.exit

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

fmuld %f0,%f48,%f48 ! (0_0) res *= dtmp0;

ld [%i4+%l6],%f0 ! (4_0) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

add %i3,stridex,%l5 ! px += stridex;

fdtos %f12,%f12 ! (7_1) ftmp0 = (float)res;

lda [%i3]0x82,%l6 ! (0_0) ux = ((int*)px)[0];

fmuld %f10,%f26,%f28 ! (6_0) div = x * y;

bpos,pt %icc,.main_loop

faddd %f4,K0,%f42 ! (0_0) dtmp0 += K0;

srl %o7,28,%l5 ! (1_0) ux >>= 28;

st %f12,[%g1] ! (7_1) py[0] = ftmp0;

.tail:

addcc counter,7,counter

bneg,pn %icc,.begin

or %g0,%o0,%o1

fsubd %f10,%f26,%f10 ! (6_1) xx = x - y;

and %o7,MASK_0x7fffffff,%g1 ! (1_1) ax = ux & 0x7fffffff;

fmuld K2,%f50,%f4 ! (1_1) dtmp0 = K2 * x2;

fmuld %f40,%f46,%f26 ! (2_1) y0 *= dtmp1;

add %g1,MASK_0x100000,%g1 ! (1_1) ax += 0x00100000;

and %l5,-8,%l5 ! (1_1) ux &= -8;

fpsub32 %f0,%f24,%f40 ! (4_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (3_1) y0 *= dtmp0;

fmuld %f42,%f22,%f44 ! (0_1) dtmp0 *= xx;

faddd DONE,%f28,%f28 ! (6_1) div += done;

fmuld %f32,%f40,%f42 ! (4_1) dtmp0 = div0 * y0;

faddd %f4,K1,%f4 ! (1_1) dtmp0 += K1;

fmuld %f18,%f26,%f18 ! (2_1) xx *= y0;

srl %g1,18,%o7 ! (1_1) ax >>= 18;

std %f28,[%fp+tmp0] ! (6_1) i = ((unsigned long long*)&div)[0];

fmuld %f34,%f38,%f34 ! (3_1) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (1_1) ax &= -8;

ldx [%fp+tmp1],%g1 ! (5_1) i = ((unsigned long long*)&div)[0];

faddd %f48,%f44,%f12 ! (0_1) res += dtmp0;

add %o7,%l7,%o7 ! (1_1) (char*)parr1 + ax;

ldd [%l0+%l5],%f48 ! (1_1) dtmp0 = *(double*)((char*)sign_arr + ux);

fmuld %f4,%f50,%f4 ! (1_1) dtmp0 *= x2;

fand %f30,DC3,%f24 ! (5_1) y0 = vis_fand(div,dconst3);

ldd [%o7],%f0 ! (1_1) res = *(double*)((char*)parr1 + ax);

fmuld %f18,%f18,%f50 ! (2_1) x2 = xx * xx;

fsubd DTWO,%f42,%f44 ! (4_1) dtmp0 = dtwo - dtmp0;

srlx %g1,43,%g1 ! (5_1) i >>= 43;

and %g1,508,%l6 ! (5_1) i &= 508;

mov %l4,%o7

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

fmuld %f0,%f48,%f48 ! (1_1) res *= dtmp0;

add %o0,stridey,%g1 ! py += stridey;

ld [%i4+%l6],%f0 ! (5_1) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (0_1) ftmp0 = (float)res;

srl %o7,28,%l4 ! (2_1) ux >>= 28;

st %f12,[%o0] ! (0_1) py[0] = ftmp0;

faddd %f4,K0,%f42 ! (1_1) dtmp0 += K0;

subcc counter,1,counter

bneg,pn %icc,.begin

or %g0,%g1,%o1

fmuld K2,%f50,%f4 ! (2_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o0 ! (2_1) ax = ux & 0x7fffffff;

fmuld %f38,%f46,%f26 ! (3_1) y0 *= dtmp1;

add %o0,MASK_0x100000,%o0 ! (2_1) ax += 0x00100000;

and %l4,-8,%l4 ! (2_1) ux &= -8;

fpsub32 %f0,%f24,%f38 ! (5_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f40,%f44,%f40 ! (4_1) y0 *= dtmp0;

fmuld %f42,%f20,%f44 ! (1_1) dtmp0 *= xx;

fmuld %f30,%f38,%f42 ! (5_1) dtmp0 = div0 * y0;

faddd %f4,K1,%f4 ! (2_1) dtmp0 += K1;

fmuld %f16,%f26,%f16 ! (3_1) xx *= y0;

srl %o0,18,%o7 ! (2_1) ax >>= 18;

fmuld %f32,%f40,%f32 ! (4_1) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (2_1) ax &= -8;

ldx [%fp+tmp0],%o0 ! (6_1) i = ((unsigned long long*)&div)[0];

faddd %f48,%f44,%f12 ! (1_1) res += dtmp0;

add %o7,%l7,%o7 ! (2_1) (char*)parr1 + ax;

ldd [%l0+%l4],%f48 ! (2_1) dtmp0 = *(double*)((char*)sign_arr + ux);

fmuld %f4,%f50,%f4 ! (2_1) dtmp0 *= x2;

fand %f28,DC3,%f24 ! (6_1) y0 = vis_fand(div,dconst3);

ldd [%o7],%f0 ! (2_1) res = *(double*)((char*)parr1 + ax);

fmuld %f16,%f16,%f50 ! (3_1) x2 = xx * xx;

fsubd DTWO,%f42,%f44 ! (5_1) dtmp0 = dtwo - dtmp0;

srlx %o0,43,%o0 ! (6_1) i >>= 43;

and %o0,508,%l6 ! (6_1) i &= 508;

mov %l3,%o7

fsubd DTWO,%f32,%f46 ! (4_1) dtmp1 = dtwo - dtmp1;

fmuld %f0,%f48,%f48 ! (2_1) res *= dtmp0;

add %g1,stridey,%o0 ! py += stridey;

ld [%i4+%l6],%f0 ! (6_1) *(float*)&dtmp0 = *(float*)((char*)parr0 + i);

fdtos %f12,%f12 ! (1_1) ftmp0 = (float)res;

srl %o7,28,%l3 ! (3_1) ux >>= 28;

st %f12,[%g1] ! (1_1) py[0] = ftmp0;

faddd %f4,K0,%f42 ! (2_1) dtmp0 += K0;

subcc counter,1,counter

bneg,pn %icc,.begin

or %g0,%o0,%o1

fmuld K2,%f50,%f4 ! (3_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%g1 ! (3_1) ax = ux & 0x7fffffff;

fmuld %f40,%f46,%f26 ! (4_1) y0 *= dtmp1;

add %g1,MASK_0x100000,%g1 ! (3_1) ax += 0x00100000;

and %l3,-8,%l3 ! (3_1) ux &= -8;

fpsub32 %f0,%f24,%f40 ! (6_1) y0 = vis_fpsub32(dtmp0, y0);

fmuld %f38,%f44,%f38 ! (5_1) y0 *= dtmp0;

fmuld %f42,%f18,%f44 ! (2_1) dtmp0 *= xx;

fmuld %f28,%f40,%f42 ! (6_1) dtmp0 = div0 * y0;

faddd %f4,K1,%f4 ! (3_1) dtmp0 += K1;

fmuld %f14,%f26,%f14 ! (4_1) xx *= y0;

srl %g1,18,%o7 ! (3_1) ax >>= 18;

fmuld %f30,%f38,%f30 ! (5_1) dtmp1 = div0 * y0;

and %o7,-8,%o7 ! (3_1) ax &= -8;

faddd %f48,%f44,%f12 ! (2_1) res += dtmp0;

add %o7,%l7,%o7 ! (3_1) (char*)parr1 + ax;

ldd [%l0+%l3],%f48 ! (3_1) dtmp0 = *(double*)((char*)sign_arr + ux);

fmuld %f4,%f50,%f4 ! (3_1) dtmp0 *= x2;

ldd [%o7],%f0 ! (3_1) res = *(double*)((char*)parr1 + ax)

fmuld %f14,%f14,%f50 ! (4_1) x2 = xx * xx;

fsubd DTWO,%f42,%f44 ! (6_1) dtmp0 = dtwo - dtmp0;

mov %i0,%o7

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

fmuld %f0,%f48,%f48 ! (3_1) res *= dtmp0;

add %o0,stridey,%g1 ! py += stridey;

fdtos %f12,%f12 ! (2_1) ftmp0 = (float)res;

srl %o7,28,%i0 ! (4_1) ux >>= 28;

st %f12,[%o0] ! (2_1) py[0] = ftmp0;

faddd %f4,K0,%f42 ! (3_1) dtmp0 += K0;

subcc counter,1,counter

bneg,pn %icc,.begin

or %g0,%g1,%o1

fmuld K2,%f50,%f4 ! (4_1) dtmp0 = K2 * x2;

and %o7,MASK_0x7fffffff,%o0 ! (4_1) ax = ux & 0x7fffffff;

fmuld %f38,%f46,%f26 ! (5_1) y0 *= dtmp1;