.file "__vhypotf.S"
#include "libm.h"
RO_DATA
.align 64
.CONST_TBL:
.word 0x3fe00001, 0x80007e00 ! K1 = 5.00000715259318464227e-01
.word 0xbfc00003, 0xc0017a01 ! K2 = -1.25000447037521686593e-01
.word 0x000fffff, 0xffffffff ! DC0 = 0x000fffffffffffff
.word 0x3ff00000, 0x00000000 ! DC1 = 0x3ff0000000000000
.word 0x7ffff000, 0x00000000 ! DC2 = 0x7ffff00000000000
.word 0x7fe00000, 0x00000000 ! DA0 = 0x7fe0000000000000
.word 0x47efffff, 0xe0000000 ! DFMAX = 3.402823e+38
.word 0x7f7fffff, 0x80808080 ! FMAX = 3.402823e+38 , SCALE = 0x80808080
.word 0x20000000, 0x00000000 ! DA1 = 0x2000000000000000
#define DC0 %f12
#define DC1 %f10
#define DC2 %f42
#define DA0 %f6
#define DA1 %f4
#define K2 %f26
#define K1 %f28
#define SCALE %f3
#define FMAX %f2
#define DFMAX %f50
#define stridex %l6
#define stridey %i4
#define stridez %l5
#define _0x7fffffff %o1
#define _0x7f3504f3 %o2
#define _0x1ff0 %l2
#define TBL %l1
#define counter %l0
#define tmp_px STACK_BIAS-0x30
#define tmp_py STACK_BIAS-0x28
#define tmp_counter STACK_BIAS-0x20
#define tmp0 STACK_BIAS-0x18
#define tmp1 STACK_BIAS-0x10
#define tmp2 STACK_BIAS-0x0c
#define tmp3 STACK_BIAS-0x08
#define tmp4 STACK_BIAS-0x04
! sizeof temp storage - must be a multiple of 16 for V9
#define tmps 0x30
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !!!!! algorithm !!!!!
! hx0 = *(int*)px;
! x0 = *px;
! px += stridex;
!
! hy0 = *(int*)py;
! y0 = *py;
! py += stridey;
!
! hx0 &= 0x7fffffff;
! hy0 &= 0x7fffffff;
!
! if ( hx >= 0x7f3504f3 || hy >= 0x7f3504f3 )
! {
! if ( hx >= 0x7f800000 || hy >= 0x7f800000 )
! {
! if ( hx == 0x7f800000 || hy == 0x7f800000 )
! *(int*)pz = 0x7f800000;
! else *pz = x * y;
! }
! else
! {
! hyp = sqrt(x * (double)x + y * (double)y);
! if ( hyp <= DMAX ) ftmp0 = (float)hyp;
! else ftmp0 = FMAX * FMAX;
! *pz = ftmp0;
! }
! pz += stridez;
! continue;
! }
! if ( (hx | hy) == 0 )
! {
! *pz = 0;
! pz += stridez;
! continue;
! }
! dx0 = x0 * (double)x0;
! dy0 = y0 * (double)y0;
! db0 = dx0 + dy0;
!
! iexp0 = ((int*)&db0)[0];
!
! h0 = vis_fand(db0,DC0);
! h0 = vis_for(h0,DC1);
! h_hi0 = vis_fand(h0,DC2);
!
! db0 = vis_fand(db0,DA0);
! db0 = vis_fmul8x16(SCALE, db0);
! db0 = vis_fpadd32(db0,DA1);
!
! iexp0 >>= 8;
! di0 = iexp0 & 0x1ff0;
! si0 = (char*)sqrt_arr + di0;
!
! dtmp0 = ((double*)((char*)div_arr + di0))[0];
! xx0 = h0 - h_hi0;
! xx0 *= dmp0;
!
! dtmp0 = ((double*)si0)[1];
! res0 = K2 * xx0;
! res0 += K1;
! res0 *= xx0;
! res0 += DC1;
! res0 = dtmp0 * res0;
! res0 *= db0;
! ftmp0 = (float)res0;
! *pz = ftmp0;
! pz += stridez;
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ENTRY(__vhypotf)
save %sp,-SA(MINFRAME)-tmps,%sp
PIC_SETUP(l7)
PIC_SET(l7,.CONST_TBL,o3)
PIC_SET(l7,__vlibm_TBL_sqrtf,l1)
#ifdef __sparcv9
ldx [%fp+STACK_BIAS+176],stridez
#else
ld [%fp+STACK_BIAS+92],stridez
#endif
st %i0,[%fp+tmp_counter]
stx %i1,[%fp+tmp_px]
stx %i3,[%fp+tmp_py]
ldd [%o3],K1
sethi %hi(0x7ffffc00),%o1
ldd [%o3+8],K2
sethi %hi(0x7f350400),%o2
ldd [%o3+16],DC0
add %o1,1023,_0x7fffffff
add %o2,0xf3,_0x7f3504f3
ldd [%o3+24],DC1
sll %i2,2,stridex
ld [%o3+56],FMAX
ldd [%o3+32],DC2
sll %i4,2,stridey
ldd [%o3+40],DA0
sll stridez,2,stridez
ldd [%o3+48],DFMAX
ld [%o3+60],SCALE
or %g0,0xff8,%l2
ldd [%o3+64],DA1
sll %l2,1,_0x1ff0
or %g0,%i5,%l7
.begin:
ld [%fp+tmp_counter],counter
ldx [%fp+tmp_px],%i1
ldx [%fp+tmp_py],%i2
st %g0,[%fp+tmp_counter]
.begin1:
cmp counter,0
ble,pn %icc,.exit
lda [%i1]0x82,%l3 ! (3_0) hx0 = *(int*)px;
lda [%i2]0x82,%l4 ! (3_0) hy0 = *(int*)py;
lda [%i1]0x82,%f17 ! (3_0) x0 = *px;
and %l3,_0x7fffffff,%l3 ! (3_0) hx0 &= 0x7fffffff;
cmp %l3,_0x7f3504f3 ! (3_0) hx ? 0x7f3504f3
bge,pn %icc,.spec ! (3_0) if ( hx >= 0x7f3504f3 )
and %l4,_0x7fffffff,%l4 ! (3_0) hy0 &= 0x7fffffff;
cmp %l4,_0x7f3504f3 ! (3_0) hy ? 0x7f3504f3
bge,pn %icc,.spec ! (3_0) if ( hy >= 0x7f3504f3 )
or %g0,%i2,%o7
orcc %l3,%l4,%g0
bz,pn %icc,.spec1
add %i1,stridex,%i1 ! px += stridex
fsmuld %f17,%f17,%f44 ! (3_0) dx0 = x0 * (double)x0;
lda [%i2]0x82,%f17 ! (3_0) y0 = *py;
lda [%i1]0x82,%l3 ! (4_0) hx0 = *(int*)px;
lda [stridey+%o7]0x82,%l4 ! (4_0) hy0 = *(int*)py;
and %l3,_0x7fffffff,%l3 ! (4_0) hx0 &= 0x7fffffff;
fsmuld %f17,%f17,%f24 ! (3_0) dy0 = y0 * (double)y0;
cmp %l3,_0x7f3504f3 ! (4_0) hx ? 0x7f3504f3
bge,pn %icc,.update0 ! (4_0) if ( hx >= 0x7f3504f3 )
and %l4,_0x7fffffff,%l4 ! (4_0) hy0 &= 0x7fffffff;
orcc %l3,%l4,%g0
bz,pn %icc,.update0
lda [%i1]0x82,%f17 ! (4_0) x0 = *px;
.cont0:
faddd %f44,%f24,%f24 ! (3_0) db0 = dx0 + dy0;
fsmuld %f17,%f17,%f40 ! (4_1) dy0 = x0 * (double)x0;
cmp %l4,_0x7f3504f3 ! (4_1) hy ? 0x7f3504f3
lda [stridey+%o7]0x82,%f17 ! (4_1) hy0 = *py;
add %o7,stridey,%i5 ! py += stridey
lda [%i1+stridex]0x82,%l3 ! (0_0) hx0 = *(int*)px;
bge,pn %icc,.update1 ! (4_1) if ( hy >= 0x7f3504f3 )
st %f24,[%fp+tmp0] ! (3_1) iexp0 = ((int*)&db0)[0];
.cont1:
and %l3,_0x7fffffff,%l3 ! (0_0) hx0 &= 0x7fffffff;
fsmuld %f17,%f17,%f48 ! (4_1) dy0 = y0 * (double)y0;
lda [%i1+stridex]0x82,%f8 ! (0_0) x0 = *px;
add %i1,stridex,%i1 ! px += stridex
lda [%i5+stridey]0x82,%l4 ! (0_0) hy0 = *(int*)py;
cmp %l3,_0x7f3504f3 ! (0_0) hx ? 0x7f3504f3
bge,pn %icc,.update2 ! (0_0) if ( hx >= 0x7f3504f3 )
add %i5,stridey,%o4 ! py += stridey
.cont2:
faddd %f40,%f48,%f20 ! (4_1) db0 = dx0 + dy0;
fsmuld %f8,%f8,%f40 ! (0_0) dx0 = x0 * (double)x0;
and %l4,_0x7fffffff,%l4 ! (0_0) hy0 &= 0x7fffffff;
lda [%i5+stridey]0x82,%f17 ! (0_0) hy0 = *py;
cmp %l4,_0x7f3504f3 ! (0_0) hy ? 0x7f3504f3
bge,pn %icc,.update3 ! (0_0) if ( hy >= 0x7f3504f3 )
st %f20,[%fp+tmp1] ! (4_1) iexp0 = ((int*)&db0)[0];
orcc %l3,%l4,%g0
bz,pn %icc,.update3
.cont3:
lda [%i1+stridex]0x82,%l3 ! (1_0) hx0 = *(int*)px;
fand %f24,DC0,%f60 ! (3_1) h0 = vis_fand(db0,DC0);
and %l3,_0x7fffffff,%l3 ! (1_0) hx0 &= 0x7fffffff;
fsmuld %f17,%f17,%f34 ! (0_0) dy0 = y0 * (double)y0;
cmp %l3,_0x7f3504f3 ! (1_0) hx ? 0x7f3504f3
lda [%o4+stridey]0x82,%l4 ! (1_0) hy0 = *(int*)py;
add %i1,stridex,%i1 ! px += stridex
lda [%i1]0x82,%f17 ! (1_0) x0 = *px;
bge,pn %icc,.update4 ! (1_0) if ( hx >= 0x7f3504f3 )
add %o4,stridey,%i5 ! py += stridey
.cont4:
and %l4,_0x7fffffff,%l4 ! (1_0) hy0 &= 0x7fffffff;
for %f60,DC1,%f46 ! (3_1) h0 = vis_for(h0,DC1);
cmp %l4,_0x7f3504f3 ! (1_0) hy ? 0x7f3504f3
ld [%fp+tmp0],%o0 ! (3_1) iexp0 = ((int*)&db0)[0];
faddd %f40,%f34,%f0 ! (0_0) db0 = dx0 + dy0;
fsmuld %f17,%f17,%f40 ! (1_0) dx0 = x0 * (double)x0;
add %i1,stridex,%i1 ! px += stridex
lda [%o4+stridey]0x82,%f17 ! (1_0) y0 = *py;
srax %o0,8,%o0 ! (3_1) iexp0 >>= 8;
bge,pn %icc,.update5 ! (1_0) if ( hy >= 0x7f3504f3 )
fand %f46,DC2,%f38 ! (3_1) h_hi0 = vis_fand(h0,DC2);
orcc %l3,%l4,%g0
bz,pn %icc,.update5
.cont5:
lda [%i1]0x82,%l3 ! (2_0) hx0 = *(int*)px;
and %o0,_0x1ff0,%o0 ! (3_1) di0 = iexp0 & 0x1ff0;
st %f0,[%fp+tmp2] ! (0_0) iexp0 = ((int*)&db0)[0];
fand %f20,DC0,%f60 ! (4_1) h0 = vis_fand(db0,DC0);
ldd [TBL+%o0],%f22 ! (3_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
fsubd %f46,%f38,%f38 ! (3_1) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f32 ! (1_0) dy0 = y0 * (double)y0;
add %i5,stridey,%i2 ! py += stridey
lda [stridey+%i5]0x82,%l4 ! (2_0) hy0 = *(int*)py;
and %l3,_0x7fffffff,%l3 ! (2_0) hx0 &= 0x7fffffff;
lda [%i1]0x82,%f17 ! (2_0) x0 = *px;
cmp %l3,_0x7f3504f3 ! (2_0) hx ? 0x7f3504f3
fmuld %f38,%f22,%f38 ! (3_1) xx0 *= dmp0;
and %l4,_0x7fffffff,%l4 ! (2_0) hy0 &= 0x7fffffff;
for %f60,DC1,%f46 ! (4_1) h0 = vis_for(h0,DC1);
bge,pn %icc,.update6 ! (2_0) if ( hx >= 0x7f3504f3 )
ld [%fp+tmp1],%o3 ! (4_1) iexp0 = ((int*)&db0)[0];
.cont6:
faddd %f40,%f32,%f18 ! (1_0) db0 = dx0 + dy0;
fsmuld %f17,%f17,%f44 ! (2_0) dx0 = x0 * (double)x0;
cmp %l4,_0x7f3504f3 ! (2_0) hy ? 0x7f3504f3
lda [stridey+%i5]0x82,%f17 ! (2_0) y0 = *py;
add %i1,stridex,%i1 ! px += stridex
bge,pn %icc,.update7 ! (2_0) if ( hy >= 0x7f3504f3 )
fand %f46,DC2,%f58 ! (4_1) h_hi0 = vis_fand(h0,DC2);
orcc %l3,%l4,%g0
bz,pn %icc,.update7
nop
.cont7:
fmuld K2,%f38,%f56 ! (3_1) res0 = K2 * xx0;
srax %o3,8,%o3 ! (4_1) iexp0 >>= 8;
lda [%i1]0x82,%l3 ! (3_0) hx0 = *(int*)px;
and %o3,_0x1ff0,%o3 ! (4_1) di0 = iexp0 & 0x1ff0;
st %f18,[%fp+tmp3] ! (1_0) iexp0 = ((int*)&db0)[0];
fand %f0,DC0,%f60 ! (0_0) h0 = vis_fand(db0,DC0);
ldd [TBL+%o3],%f22 ! (4_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %i2,stridey,%o7 ! py += stridey
fsubd %f46,%f58,%f58 ! (4_1) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f30 ! (2_0) dy0 = y0 * (double)y0;
lda [stridey+%i2]0x82,%l4 ! (3_0) hy0 = *(int*)py;
and %l3,_0x7fffffff,%l3 ! (3_0) hx0 &= 0x7fffffff;
faddd %f56,K1,%f54 ! (3_1) res0 += K1;
cmp %l3,_0x7f3504f3 ! (3_0) hx ? 0x7f3504f3
lda [%i1]0x82,%f17 ! (3_0) x0 = *px;
add %i1,stridex,%i1 ! px += stridex
bge,pn %icc,.update8 ! (3_0) if ( hx >= 0x7f3504f3 )
fmuld %f58,%f22,%f58 ! (4_1) xx0 *= dmp0;
.cont8:
and %l4,_0x7fffffff,%l4 ! (3_0) hy0 &= 0x7fffffff;
for %f60,DC1,%f46 ! (0_0) h0 = vis_for(h0,DC1);
cmp %l4,_0x7f3504f3 ! (3_0) hy ? 0x7f3504f3
ld [%fp+tmp2],%g1 ! (0_0) iexp0 = ((int*)&db0)[0];
faddd %f44,%f30,%f30 ! (2_0) db0 = dx0 + dy0;
fsmuld %f17,%f17,%f44 ! (3_0) dx0 = x0 * (double)x0;
bge,pn %icc,.update9 ! (3_0) if ( hy >= 0x7f3504f3 )
lda [stridey+%i2]0x82,%f17 ! (3_0) y0 = *py;
orcc %l3,%l4,%g0
bz,pn %icc,.update9
nop
.cont9:
fmuld %f54,%f38,%f40 ! (3_1) res0 *= xx0;
lda [%i1]0x82,%l3 ! (4_0) hx0 = *(int*)px;
fand %f46,DC2,%f38 ! (0_0) h_hi0 = vis_fand(h0,DC2);
fmuld K2,%f58,%f54 ! (4_1) res0 = K2 * xx0;
srax %g1,8,%o5 ! (0_0) iexp0 >>= 8;
lda [stridey+%o7]0x82,%l4 ! (4_0) hy0 = *(int*)py;
fand %f24,DA0,%f56 ! (3_1) db0 = vis_fand(db0,DA0);
and %o5,_0x1ff0,%o5 ! (0_0) di0 = iexp0 & 0x1ff0;
st %f30,[%fp+tmp4] ! (2_0) iexp0 = ((int*)&db0)[0];
fand %f18,DC0,%f60 ! (1_0) h0 = vis_fand(db0,DC0);
ldd [TBL+%o5],%f22 ! (0_0) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %o0,TBL,%g1 ! (3_1) si0 = (char*)sqrt_arr + di0;
and %l3,_0x7fffffff,%l3 ! (4_0) hx0 &= 0x7fffffff;
fsubd %f46,%f38,%f38 ! (0_0) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f24 ! (3_0) dy0 = y0 * (double)y0;
cmp %l3,_0x7f3504f3 ! (4_0) hx ? 0x7f3504f3
bge,pn %icc,.update10 ! (4_0) if ( hx >= 0x7f3504f3 )
faddd %f40,DC1,%f40 ! (3_1) res0 += DC1;
fmul8x16 SCALE,%f56,%f36 ! (3_1) db0 = vis_fmul8x16(SCALE, db0);
and %l4,_0x7fffffff,%l4 ! (4_0) hy0 &= 0x7fffffff;
ldd [%g1+8],%f56 ! (3_1) dtmp0 = ((double*)si0)[1];
faddd %f54,K1,%f54 ! (4_1) res0 += K1;
lda [%i1]0x82,%f17 ! (4_0) x0 = *px;
.cont10:
fmuld %f38,%f22,%f38 ! (0_0) xx0 *= dmp0;
cmp counter,5
for %f60,DC1,%f46 ! (1_0) h0 = vis_for(h0,DC1);
ld [%fp+tmp3],%g1 ! (1_0) iexp0 = ((int*)&db0)[0];
fmuld %f56,%f40,%f62 ! (3_1) res0 = dtmp0 * res0;
faddd %f44,%f24,%f24 ! (3_0) db0 = dx0 + dy0;
bl,pn %icc,.tail
nop
ba .main_loop
sub counter,5,counter
.align 16
.main_loop:
fsmuld %f17,%f17,%f40 ! (4_1) dy0 = x0 * (double)x0;
cmp %l4,_0x7f3504f3 ! (4_1) hy ? 0x7f3504f3
lda [stridey+%o7]0x82,%f17 ! (4_1) hy0 = *py;
fpadd32 %f36,DA1,%f36 ! (3_2) db0 = vis_fpadd32(db0,DA1);
fmuld %f54,%f58,%f58 ! (4_2) res0 *= xx0;
add %o7,stridey,%i5 ! py += stridey
st %f24,[%fp+tmp0] ! (3_1) iexp0 = ((int*)&db0)[0];
fand %f46,DC2,%f44 ! (1_1) h_hi0 = vis_fand(h0,DC2);
fmuld K2,%f38,%f56 ! (0_1) res0 = K2 * xx0;
srax %g1,8,%g5 ! (1_1) iexp0 >>= 8;
bge,pn %icc,.update11 ! (4_1) if ( hy >= 0x7f3504f3 )
fand %f20,DA0,%f54 ! (4_2) db0 = vis_fand(db0,DA0);
orcc %l3,%l4,%g0
nop
bz,pn %icc,.update11
fzero %f52
.cont11:
fmuld %f62,%f36,%f62 ! (3_2) res0 *= db0;
and %g5,_0x1ff0,%g5 ! (1_1) di0 = iexp0 & 0x1ff0;
lda [%i1+stridex]0x82,%l3 ! (0_0) hx0 = *(int*)px;
fand %f30,DC0,%f60 ! (2_1) h0 = vis_fand(db0,DC0);
ldd [%g5+TBL],%f22 ! (1_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %o3,TBL,%g1 ! (4_2) si0 = (char*)sqrt_arr + di0;
add %i1,stridex,%i0 ! px += stridex
fsubd %f46,%f44,%f44 ! (1_1) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f48 ! (4_1) dy0 = y0 * (double)y0;
nop
lda [%i1+stridex]0x82,%f8 ! (0_0) x0 = *px;
faddd %f58,DC1,%f36 ! (4_2) res0 += DC1;
faddd %f56,K1,%f58 ! (0_1) res0 += K1;
and %l3,_0x7fffffff,%l3 ! (0_0) hx0 &= 0x7fffffff;
ldd [%g1+8],%f56 ! (4_2) dtmp0 = ((double*)si0)[1];
fmul8x16 SCALE,%f54,%f54 ! (4_2) db0 = vis_fmul8x16(SCALE, db0);
lda [%i5+stridey]0x82,%l4 ! (0_0) hy0 = *(int*)py;
cmp %l3,_0x7f3504f3 ! (0_0) hx ? 0x7f3504f3
bge,pn %icc,.update12 ! (0_0) if ( hx >= 0x7f3504f3 )
fdtos %f62,%f14 ! (3_2) ftmp0 = (float)res0;
.cont12:
fmuld %f44,%f22,%f44 ! (1_1) xx0 *= dmp0;
add %l7,stridez,%o7 ! pz += stridez
st %f14,[%l7] ! (3_2) *pz = ftmp0;
for %f60,DC1,%f46 ! (2_1) h0 = vis_for(h0,DC1);
fmuld %f56,%f36,%f36 ! (4_2) res0 = dtmp0 * res0;
add %i5,stridey,%o4 ! py += stridey
ld [%fp+tmp4],%g1 ! (2_1) iexp0 = ((int*)&db0)[0];
faddd %f40,%f48,%f20 ! (4_1) db0 = dx0 + dy0;
fsmuld %f8,%f8,%f40 ! (0_0) dx0 = x0 * (double)x0;
and %l4,_0x7fffffff,%l4 ! (0_0) hy0 &= 0x7fffffff;
lda [%i5+stridey]0x82,%f17 ! (0_0) hy0 = *py;
fpadd32 %f54,DA1,%f62 ! (4_2) db0 = vis_fpadd32(db0,DA1);
fmuld %f58,%f38,%f38 ! (0_1) res0 *= xx0;
cmp %l4,_0x7f3504f3 ! (0_0) hy ? 0x7f3504f3
st %f20,[%fp+tmp1] ! (4_1) iexp0 = ((int*)&db0)[0];
fand %f46,DC2,%f58 ! (2_1) h_hi0 = vis_fand(h0,DC2);
fmuld K2,%f44,%f56 ! (1_1) res0 = K2 * xx0;
srax %g1,8,%g1 ! (2_1) iexp0 >>= 8;
bge,pn %icc,.update13 ! (0_0) if ( hy >= 0x7f3504f3 )
fand %f0,DA0,%f54 ! (0_1) db0 = vis_fand(db0,DA0);
orcc %l3,%l4,%g0
nop
bz,pn %icc,.update13
fzero %f52
.cont13:
fmuld %f36,%f62,%f62 ! (4_2) res0 *= db0;
and %g1,_0x1ff0,%g1 ! (2_1) di0 = iexp0 & 0x1ff0;
lda [%i0+stridex]0x82,%l3 ! (1_0) hx0 = *(int*)px;
fand %f24,DC0,%f60 ! (3_1) h0 = vis_fand(db0,DC0);
ldd [TBL+%g1],%f22 ! (2_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %o5,TBL,%o0 ! (0_1) si0 = (char*)sqrt_arr + di0;
add %i0,stridex,%i1 ! px += stridex
fsubd %f46,%f58,%f58 ! (2_1) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f34 ! (0_0) dy0 = y0 * (double)y0;
add %o7,stridez,%i0 ! pz += stridez
lda [%o4+stridey]0x82,%l4 ! (1_0) hy0 = *(int*)py;
faddd %f38,DC1,%f36 ! (0_1) res0 += DC1;
faddd %f56,K1,%f38 ! (1_1) res0 += K1;
and %l3,_0x7fffffff,%l3 ! (1_0) hx0 &= 0x7fffffff;
ldd [%o0+8],%f56 ! (0_1) dtmp0 = ((double*)si0)[1];
fmul8x16 SCALE,%f54,%f54 ! (0_1) db0 = vis_fmul8x16(SCALE, db0);
lda [%i1]0x82,%f17 ! (1_0) x0 = *px;
cmp %l3,_0x7f3504f3 ! (1_0) hx ? 0x7f3504f3
bge,pn %icc,.update14 ! (1_0) if ( hx >= 0x7f3504f3 )
fdtos %f62,%f14 ! (4_2) ftmp0 = (float)res0;
.cont14:
fmuld %f58,%f22,%f58 ! (2_1) xx0 *= dmp0;
and %l4,_0x7fffffff,%l4 ! (1_0) hy0 &= 0x7fffffff;
add %o4,stridey,%i5 ! py += stridey
for %f60,DC1,%f46 ! (3_1) h0 = vis_for(h0,DC1);
fmuld %f56,%f36,%f36 ! (0_1) res0 = dtmp0 * res0;
cmp %l4,_0x7f3504f3 ! (1_0) hy ? 0x7f3504f3
ld [%fp+tmp0],%o0 ! (3_1) iexp0 = ((int*)&db0)[0];
faddd %f40,%f34,%f0 ! (0_0) db0 = dx0 + dy0;
fsmuld %f17,%f17,%f40 ! (1_0) dx0 = x0 * (double)x0;
add %i1,stridex,%i1 ! px += stridex
lda [%o4+stridey]0x82,%f17 ! (1_0) y0 = *py;
fpadd32 %f54,DA1,%f62 ! (0_1) db0 = vis_fpadd32(db0,DA1);
fmuld %f38,%f44,%f44 ! (1_1) res0 *= xx0;
st %f14,[%o7] ! (4_2) *pz = ftmp0;
bge,pn %icc,.update15 ! (1_0) if ( hy >= 0x7f3504f3 )
fand %f46,DC2,%f38 ! (3_1) h_hi0 = vis_fand(h0,DC2);
orcc %l3,%l4,%g0
bz,pn %icc,.update15
nop
.cont15:
fmuld K2,%f58,%f54 ! (2_1) res0 = K2 * xx0;
srax %o0,8,%o0 ! (3_1) iexp0 >>= 8;
st %f0,[%fp+tmp2] ! (0_0) iexp0 = ((int*)&db0)[0];
fand %f18,DA0,%f56 ! (1_1) db0 = vis_fand(db0,DA0);
fmuld %f36,%f62,%f62 ! (0_1) res0 *= db0;
and %o0,_0x1ff0,%o0 ! (3_1) di0 = iexp0 & 0x1ff0;
lda [%i1]0x82,%l3 ! (2_0) hx0 = *(int*)px;
fand %f20,DC0,%f60 ! (4_1) h0 = vis_fand(db0,DC0);
ldd [TBL+%o0],%f22 ! (3_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %g5,TBL,%o3 ! (1_1) si0 = (char*)sqrt_arr + di0;
add %i0,stridez,%i3 ! pz += stridez
fsubd %f46,%f38,%f38 ! (3_1) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f32 ! (1_0) dy0 = y0 * (double)y0;
add %i5,stridey,%i2 ! py += stridey
lda [stridey+%i5]0x82,%l4 ! (2_0) hy0 = *(int*)py;
faddd %f44,DC1,%f44 ! (1_1) res0 += DC1;
fmul8x16 SCALE,%f56,%f36 ! (1_1) db0 = vis_fmul8x16(SCALE, db0);
and %l3,_0x7fffffff,%l3 ! (2_0) hx0 &= 0x7fffffff;
ldd [%o3+8],%f56 ! (1_1) dtmp0 = ((double*)si0)[1];
faddd %f54,K1,%f54 ! (2_1) res0 += K1;
lda [%i1]0x82,%f17 ! (2_0) x0 = *px;
cmp %l3,_0x7f3504f3 ! (2_0) hx ? 0x7f3504f3
add %i3,stridez,%o4 ! pz += stridez
fdtos %f62,%f14 ! (0_1) ftmp0 = (float)res0;
fmuld %f38,%f22,%f38 ! (3_1) xx0 *= dmp0;
and %l4,_0x7fffffff,%l4 ! (2_0) hy0 &= 0x7fffffff;
st %f14,[%i0] ! (0_1) *pz = ftmp0;
for %f60,DC1,%f46 ! (4_1) h0 = vis_for(h0,DC1);
fmuld %f56,%f44,%f62 ! (1_1) res0 = dtmp0 * res0;
bge,pn %icc,.update16 ! (2_0) if ( hx >= 0x7f3504f3 )
ld [%fp+tmp1],%o3 ! (4_1) iexp0 = ((int*)&db0)[0];
faddd %f40,%f32,%f18 ! (1_0) db0 = dx0 + dy0;
.cont16:
fsmuld %f17,%f17,%f44 ! (2_0) dx0 = x0 * (double)x0;
cmp %l4,_0x7f3504f3 ! (2_0) hy ? 0x7f3504f3
lda [stridey+%i5]0x82,%f17 ! (2_0) y0 = *py;
fpadd32 %f36,DA1,%f36 ! (1_1) db0 = vis_fpadd32(db0,DA1);
fmuld %f54,%f58,%f54 ! (2_1) res0 *= xx0;
add %i1,stridex,%l7 ! px += stridex
bge,pn %icc,.update17 ! (2_0) if ( hy >= 0x7f3504f3 )
fand %f46,DC2,%f58 ! (4_1) h_hi0 = vis_fand(h0,DC2);
orcc %l3,%l4,%g0
nop
bz,pn %icc,.update17
fzero %f52
.cont17:
fmuld K2,%f38,%f56 ! (3_1) res0 = K2 * xx0;
srax %o3,8,%o3 ! (4_1) iexp0 >>= 8;
st %f18,[%fp+tmp3] ! (1_0) iexp0 = ((int*)&db0)[0];
fand %f30,DA0,%f40 ! (2_1) db0 = vis_fand(db0,DA0);
fmuld %f62,%f36,%f62 ! (1_1) res0 *= db0;
and %o3,_0x1ff0,%o3 ! (4_1) di0 = iexp0 & 0x1ff0;
lda [%l7]0x82,%l3 ! (3_0) hx0 = *(int*)px;
fand %f0,DC0,%f60 ! (0_0) h0 = vis_fand(db0,DC0);
ldd [TBL+%o3],%f22 ! (4_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %g1,TBL,%g1 ! (2_1) si0 = (char*)sqrt_arr + di0;
add %i2,stridey,%o7 ! py += stridey
fsubd %f46,%f58,%f58 ! (4_1) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f30 ! (2_0) dy0 = y0 * (double)y0;
lda [stridey+%i2]0x82,%l4 ! (3_0) hy0 = *(int*)py;
add %l7,stridex,%i1 ! px += stridex
faddd %f54,DC1,%f36 ! (2_1) res0 += DC1;
faddd %f56,K1,%f54 ! (3_1) res0 += K1;
and %l3,_0x7fffffff,%l3 ! (3_0) hx0 &= 0x7fffffff;
ldd [%g1+8],%f56 ! (2_1) dtmp0 = ((double*)si0)[1];
fmul8x16 SCALE,%f40,%f40 ! (2_1) db0 = vis_fmul8x16(SCALE, db0);
lda [%l7]0x82,%f17 ! (3_0) x0 = *px;
cmp %l3,_0x7f3504f3 ! (3_0) hx ? 0x7f3504f3
bge,pn %icc,.update18 ! (3_0) if ( hx >= 0x7f3504f3 )
fdtos %f62,%f14 ! (1_1) ftmp0 = (float)res0;
.cont18:
fmuld %f58,%f22,%f58 ! (4_1) xx0 *= dmp0;
and %l4,_0x7fffffff,%l4 ! (3_0) hy0 &= 0x7fffffff;
st %f14,[%i3] ! (1_1) *pz = ftmp0;
for %f60,DC1,%f46 ! (0_0) h0 = vis_for(h0,DC1);
fmuld %f56,%f36,%f36 ! (2_1) res0 = dtmp0 * res0;
cmp %l4,_0x7f3504f3 ! (3_0) hy ? 0x7f3504f3
ld [%fp+tmp2],%g1 ! (0_0) iexp0 = ((int*)&db0)[0];
faddd %f44,%f30,%f30 ! (2_0) db0 = dx0 + dy0;
fsmuld %f17,%f17,%f44 ! (3_0) dx0 = x0 * (double)x0;
bge,pn %icc,.update19 ! (3_0) if ( hy >= 0x7f3504f3 )
lda [stridey+%i2]0x82,%f17 ! (3_0) y0 = *py;
fpadd32 %f40,DA1,%f62 ! (2_1) db0 = vis_fpadd32(db0,DA1);
.cont19:
fmuld %f54,%f38,%f40 ! (3_1) res0 *= xx0;
orcc %l3,%l4,%g0
st %f30,[%fp+tmp4] ! (2_0) iexp0 = ((int*)&db0)[0];
fand %f46,DC2,%f38 ! (0_0) h_hi0 = vis_fand(h0,DC2);
fmuld K2,%f58,%f54 ! (4_1) res0 = K2 * xx0;
srax %g1,8,%o5 ! (0_0) iexp0 >>= 8;
lda [%i1]0x82,%l3 ! (4_0) hx0 = *(int*)px;
fand %f24,DA0,%f56 ! (3_1) db0 = vis_fand(db0,DA0);
fmuld %f36,%f62,%f62 ! (2_1) res0 *= db0;
and %o5,_0x1ff0,%o5 ! (0_0) di0 = iexp0 & 0x1ff0;
bz,pn %icc,.update19a
fand %f18,DC0,%f60 ! (1_0) h0 = vis_fand(db0,DC0);
.cont19a:
ldd [TBL+%o5],%f22 ! (0_0) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %o0,TBL,%g1 ! (3_1) si0 = (char*)sqrt_arr + di0;
and %l3,_0x7fffffff,%l3 ! (4_0) hx0 &= 0x7fffffff;
fsubd %f46,%f38,%f38 ! (0_0) xx0 = h0 - h_hi0;
fsmuld %f17,%f17,%f24 ! (3_0) dy0 = y0 * (double)y0;
cmp %l3,_0x7f3504f3 ! (4_0) hx ? 0x7f3504f3
lda [stridey+%o7]0x82,%l4 ! (4_0) hy0 = *(int*)py;
faddd %f40,DC1,%f40 ! (3_1) res0 += DC1;
fmul8x16 SCALE,%f56,%f36 ! (3_1) db0 = vis_fmul8x16(SCALE, db0);
bge,pn %icc,.update20 ! (4_0) if ( hx >= 0x7f3504f3 )
ldd [%g1+8],%f56 ! (3_1) dtmp0 = ((double*)si0)[1];
faddd %f54,K1,%f54 ! (4_1) res0 += K1;
lda [%i1]0x82,%f17 ! (4_0) x0 = *px;
.cont20:
subcc counter,5,counter ! counter -= 5
add %o4,stridez,%l7 ! pz += stridez
fdtos %f62,%f14 ! (2_1) ftmp0 = (float)res0;
fmuld %f38,%f22,%f38 ! (0_0) xx0 *= dmp0;
and %l4,_0x7fffffff,%l4 ! (4_0) hy0 &= 0x7fffffff;
st %f14,[%o4] ! (2_1) *pz = ftmp0;
for %f60,DC1,%f46 ! (1_0) h0 = vis_for(h0,DC1);
ld [%fp+tmp3],%g1 ! (1_0) iexp0 = ((int*)&db0)[0];
fmuld %f56,%f40,%f62 ! (3_1) res0 = dtmp0 * res0;
bpos,pt %icc,.main_loop
faddd %f44,%f24,%f24 ! (3_0) db0 = dx0 + dy0;
add counter,5,counter
.tail:
subcc counter,1,counter
bneg .begin
nop
fpadd32 %f36,DA1,%f36 ! (3_2) db0 = vis_fpadd32(db0,DA1);
fmuld %f54,%f58,%f58 ! (4_2) res0 *= xx0;
fand %f46,DC2,%f44 ! (1_1) h_hi0 = vis_fand(h0,DC2);
fmuld K2,%f38,%f56 ! (0_1) res0 = K2 * xx0;
srax %g1,8,%g5 ! (1_1) iexp0 >>= 8;
fand %f20,DA0,%f54 ! (4_2) db0 = vis_fand(db0,DA0);
fmuld %f62,%f36,%f62 ! (3_2) res0 *= db0;
and %g5,_0x1ff0,%g5 ! (1_1) di0 = iexp0 & 0x1ff0;
ldd [%g5+TBL],%f22 ! (1_1) dtmp0 = ((double*)((char*)div_arr + di0))[0];
add %o3,TBL,%g1 ! (4_2) si0 = (char*)sqrt_arr + di0;
fsubd %f46,%f44,%f44 ! (1_1) xx0 = h0 - h_hi0;
faddd %f58,DC1,%f36 ! (4_2) res0 += DC1;
faddd %f56,K1,%f58 ! (0_1) res0 += K1;
ldd [%g1+8],%f56 ! (4_2) dtmp0 = ((double*)si0)[1];
fmul8x16 SCALE,%f54,%f54 ! (4_2) db0 = vis_fmul8x16(SCALE, db0);
fdtos %f62,%f14 ! (3_2) ftmp0 = (float)res0;
fmuld %f44,%f22,%f44 ! (1_1) xx0 *= dmp0;
add %l7,stridez,%o7 ! pz += stridez
st %f14,[%l7] ! (3_2) *pz = ftmp0;
subcc counter,1,counter
bneg .begin
or %g0,%o7,%l7
fmuld %f56,%f36,%f36 ! (4_2) res0 = dtmp0 * res0;
fpadd32 %f54,DA1,%f62 ! (4_2) db0 = vis_fpadd32(db0,DA1);
fmuld %f58,%f38,%f38 ! (0_1) res0 *= xx0;
fmuld K2,%f44,%f56 ! (1_1) res0 = K2 * xx0;
fand %f0,DA0,%f54 ! (0_1) db0 = vis_fand(db0,DA0);
fmuld %f36,%f62,%f62 ! (4_2) res0 *= db0;
add %o5,TBL,%o0 ! (0_1) si0 = (char*)sqrt_arr + di0;
faddd %f38,DC1,%f36 ! (0_1) res0 += DC1;
faddd %f56,K1,%f38 ! (1_1) res0 += K1;
ldd [%o0+8],%f56 ! (0_1) dtmp0 = ((double*)si0)[1];
fmul8x16 SCALE,%f54,%f54 ! (0_1) db0 = vis_fmul8x16(SCALE, db0);
add %o7,stridez,%i0 ! pz += stridez
fdtos %f62,%f14 ! (4_2) ftmp0 = (float)res0;
fmuld %f56,%f36,%f36 ! (0_1) res0 = dtmp0 * res0;
fpadd32 %f54,DA1,%f62 ! (0_1) db0 = vis_fpadd32(db0,DA1);
fmuld %f38,%f44,%f44 ! (1_1) res0 *= xx0;
add %i0,stridez,%i3 ! pz += stridez
st %f14,[%o7] ! (4_2) *pz = ftmp0;
subcc counter,1,counter
bneg .begin
or %g0,%i0,%l7
fand %f18,DA0,%f56 ! (1_1) db0 = vis_fand(db0,DA0);
fmuld %f36,%f62,%f62 ! (0_1) res0 *= db0;
add %g5,TBL,%o3 ! (1_1) si0 = (char*)sqrt_arr + di0;
faddd %f44,DC1,%f44 ! (1_1) res0 += DC1;
fmul8x16 SCALE,%f56,%f36 ! (1_1) db0 = vis_fmul8x16(SCALE, db0);
ldd [%o3+8],%f56 ! (1_1) dtmp0 = ((double*)si0)[1];
add %i3,stridez,%o4 ! pz += stridez
fdtos %f62,%f14 ! (0_1) ftmp0 = (float)res0;
st %f14,[%i0] ! (0_1) *pz = ftmp0;
subcc counter,1,counter
bneg .begin
or %g0,%i3,%l7
fmuld %f56,%f44,%f62 ! (1_1) res0 = dtmp0 * res0;
fpadd32 %f36,DA1,%f36 ! (1_1) db0 = vis_fpadd32(db0,DA1);
fmuld %f62,%f36,%f62 ! (1_1) res0 *= db0;
fdtos %f62,%f14 ! (1_1) ftmp0 = (float)res0;
st %f14,[%i3] ! (1_1) *pz = ftmp0;
ba .begin
or %g0,%o4,%l7
.align 16
.spec1:
st %g0,[%l7] ! *pz = 0;
add %l7,stridez,%l7 ! pz += stridez
add %i2,stridey,%i2 ! py += stridey
ba .begin1
sub counter,1,counter ! counter--
.align 16
.spec:
sethi %hi(0x7f800000),%i0
cmp %l3,%i0 ! hx ? 0x7f800000
bge,pt %icc,2f ! if ( hx >= 0x7f800000 )
ld [%i2],%f8
cmp %l4,%i0 ! hy ? 0x7f800000
bge,pt %icc,2f ! if ( hy >= 0x7f800000 )
nop
fsmuld %f17,%f17,%f44 ! x * (double)x
fsmuld %f8,%f8,%f24 ! y * (double)y
faddd %f44,%f24,%f24 ! x * (double)x + y * (double)y
fsqrtd %f24,%f24 ! hyp = sqrt(x * (double)x + y * (double)y);
fcmped %f24,DFMAX ! hyp ? DMAX
fbug,a 1f ! if ( hyp > DMAX )
fmuls FMAX,FMAX,%f20 ! ftmp0 = FMAX * FMAX;
fdtos %f24,%f20 ! ftmp0 = (float)hyp;
1:
st %f20,[%l7] ! *pz = ftmp0;
add %l7,stridez,%l7 ! pz += stridez
add %i1,stridex,%i1 ! px += stridex
add %i2,stridey,%i2 ! py += stridey
ba .begin1
sub counter,1,counter ! counter--
2:
fcmps %f17,%f8 ! exceptions
cmp %l3,%i0 ! hx ? 0x7f800000
be,a %icc,1f ! if ( hx == 0x7f800000 )
st %i0,[%l7] ! *(int*)pz = 0x7f800000;
cmp %l4,%i0 ! hy ? 0x7f800000
be,a %icc,1f ! if ( hy == 0x7f800000
st %i0,[%l7] ! *(int*)pz = 0x7f800000;
fmuls %f17,%f8,%f8 ! x * y
st %f8,[%l7] ! *pz = x * y;
1:
add %l7,stridez,%l7 ! pz += stridez
add %i1,stridex,%i1 ! px += stridex
add %i2,stridey,%i2 ! py += stridey
ba .begin1
sub counter,1,counter ! counter--
.align 16
.update0:
cmp counter,1
ble .cont0
fzeros %f17
stx %i1,[%fp+tmp_px]
add %o7,stridey,%i5
stx %i5,[%fp+tmp_py]
sub counter,1,counter
st counter,[%fp+tmp_counter]
ba .cont0
or %g0,1,counter
.align 16
.update1:
cmp counter,1
ble .cont1
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %i5,[%fp+tmp_py]
sub counter,1,counter
st counter,[%fp+tmp_counter]
ba .cont1
or %g0,1,counter
.align 16
.update2:
cmp counter,2
ble .cont2
fzeros %f8
stx %i1,[%fp+tmp_px]
stx %o4,[%fp+tmp_py]
sub counter,2,counter
st counter,[%fp+tmp_counter]
ba .cont2
or %g0,2,counter
.align 16
.update3:
cmp counter,2
ble .cont3
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %o4,[%fp+tmp_py]
sub counter,2,counter
st counter,[%fp+tmp_counter]
ba .cont3
or %g0,2,counter
.align 16
.update4:
cmp counter,3
ble .cont4
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %i5,[%fp+tmp_py]
sub counter,3,counter
st counter,[%fp+tmp_counter]
ba .cont4
or %g0,3,counter
.align 16
.update5:
cmp counter,3
ble .cont5
fzeros %f17
sub %i1,stridex,%i2
stx %i2,[%fp+tmp_px]
stx %i5,[%fp+tmp_py]
sub counter,3,counter
st counter,[%fp+tmp_counter]
ba .cont5
or %g0,3,counter
.align 16
.update6:
cmp counter,4
ble .cont6
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %i2,[%fp+tmp_py]
sub counter,4,counter
st counter,[%fp+tmp_counter]
ba .cont6
or %g0,4,counter
.align 16
.update7:
cmp counter,4
ble .cont7
fzeros %f17
sub %i1,stridex,%o7
stx %o7,[%fp+tmp_px]
stx %i2,[%fp+tmp_py]
sub counter,4,counter
st counter,[%fp+tmp_counter]
ba .cont7
or %g0,4,counter
.align 16
.update8:
cmp counter,5
ble .cont8
fzeros %f17
sub %i1,stridex,%o5
stx %o5,[%fp+tmp_px]
stx %o7,[%fp+tmp_py]
sub counter,5,counter
st counter,[%fp+tmp_counter]
ba .cont8
or %g0,5,counter
.align 16
.update9:
cmp counter,5
ble .cont9
fzeros %f17
sub %i1,stridex,%o5
stx %o5,[%fp+tmp_px]
stx %o7,[%fp+tmp_py]
sub counter,5,counter
st counter,[%fp+tmp_counter]
ba .cont9
or %g0,5,counter
.align 16
.update10:
fmul8x16 SCALE,%f56,%f36 ! (3_1) db0 = vis_fmul8x16(SCALE, db0);
and %l4,_0x7fffffff,%l4 ! (4_0) hy0 &= 0x7fffffff;
ldd [%g1+8],%f56 ! (3_1) dtmp0 = ((double*)si0)[1];
faddd %f54,K1,%f54 ! (4_1) res0 += K1;
cmp counter,6
ble .cont10
fzeros %f17
stx %i1,[%fp+tmp_px]
add %o7,stridey,%i5
stx %i5,[%fp+tmp_py]
sub counter,6,counter
st counter,[%fp+tmp_counter]
ba .cont10
or %g0,6,counter
.align 16
.update11:
cmp counter,1
ble .cont11
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %i5,[%fp+tmp_py]
sub counter,1,counter
st counter,[%fp+tmp_counter]
ba .cont11
or %g0,1,counter
.align 16
.update12:
cmp counter,2
ble .cont12
fzeros %f8
stx %i0,[%fp+tmp_px]
add %i5,stridey,%o4
stx %o4,[%fp+tmp_py]
sub counter,2,counter
st counter,[%fp+tmp_counter]
ba .cont12
or %g0,2,counter
.align 16
.update13:
cmp counter,2
ble .cont13
fzeros %f17
stx %i0,[%fp+tmp_px]
stx %o4,[%fp+tmp_py]
sub counter,2,counter
st counter,[%fp+tmp_counter]
ba .cont13
or %g0,2,counter
.align 16
.update14:
cmp counter,3
ble .cont14
fzeros %f17
stx %i1,[%fp+tmp_px]
add %o4,stridey,%i5
stx %i5,[%fp+tmp_py]
sub counter,3,counter
st counter,[%fp+tmp_counter]
ba .cont14
or %g0,3,counter
.align 16
.update15:
cmp counter,3
ble .cont15
fzeros %f17
sub %i1,stridex,%i2
stx %i2,[%fp+tmp_px]
stx %i5,[%fp+tmp_py]
sub counter,3,counter
st counter,[%fp+tmp_counter]
ba .cont15
or %g0,3,counter
.align 16
.update16:
faddd %f40,%f32,%f18 ! (1_0) db0 = dx0 + dy0;
cmp counter,4
ble .cont16
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %i2,[%fp+tmp_py]
sub counter,4,counter
st counter,[%fp+tmp_counter]
ba .cont16
or %g0,4,counter
.align 16
.update17:
cmp counter,4
ble .cont17
fzeros %f17
stx %i1,[%fp+tmp_px]
stx %i2,[%fp+tmp_py]
sub counter,4,counter
st counter,[%fp+tmp_counter]
ba .cont17
or %g0,4,counter
.align 16
.update18:
cmp counter,5
ble .cont18
fzeros %f17
stx %l7,[%fp+tmp_px]
stx %o7,[%fp+tmp_py]
sub counter,5,counter
st counter,[%fp+tmp_counter]
ba .cont18
or %g0,5,counter
.align 16
.update19:
fpadd32 %f40,DA1,%f62 ! (2_1) db0 = vis_fpadd32(db0,DA1);
cmp counter,5
ble .cont19
fzeros %f17
stx %l7,[%fp+tmp_px]
stx %o7,[%fp+tmp_py]
sub counter,5,counter
st counter,[%fp+tmp_counter]
ba .cont19
or %g0,5,counter
.align 16
.update19a:
cmp counter,5
ble .cont19a
fzeros %f17
stx %l7,[%fp+tmp_px]
stx %o7,[%fp+tmp_py]
sub counter,5,counter
st counter,[%fp+tmp_counter]
ba .cont19a
or %g0,5,counter
.align 16
.update20:
faddd %f54,K1,%f54 ! (4_1) res0 += K1;
cmp counter,6
ble .cont20
fzeros %f17
stx %i1,[%fp+tmp_px]
add %o7,stridey,%g1
stx %g1,[%fp+tmp_py]
sub counter,6,counter
st counter,[%fp+tmp_counter]
ba .cont20
or %g0,6,counter
.exit:
ret
restore
SET_SIZE(__vhypotf)
