cci/sys/tahoe/Ksqrtf.x - view

File: [Power 6/32 Unix Tahoe 4.2BSD] / cci / sys / tahoe / Ksqrtf.x
Revision 1.1.1.2 (vendor branch): download - view: text, annotated - select for diffs
Sun Jul 28 12:30:01 2019 UTC (6 years, 11 months ago) by root
Branches: bsd, MAIN
CVS tags: v121, HEAD

Power 6/32 Unix version 1.21

/* file: Ksqrtf.x */ .data .comm _errno,4 .text LL0: .align 1 .globl _Ksqrtf .data .align 2 L22: .long 0x40000000, 0x00000000 # .double .5 .text .data .align 2 L23: .long 0x41000000, 0x00000000 # .double 2 .text .data .align 2 L25: .long 0x41000000, 0x00000000 # .double 2 .text .data .align 2 L26: .long 0x40000000, 0x00000000 # .double .5 .align 2 L27: .long 0x40800000, 0x00000000 # .double 1 .text .data .align 2 L30: .long 0x4F800000, 0x00000000 # .double 1073741824 .text .data .align 2 L33: .long 0x4F800000, 0x00000000 # .double 1073741824 .text .data .align 2 L39: .long 0x40000000, 0x00000000 # .double .5 .text .set L14,0x0 .data .text _Ksqrtf: .word L14 subl3 $76,fp,sp tstl 4(fp) # if (arg <= 0.) { jgtr L18 jgeq L19 # if (arg < 0.) movl $33,_errno # errno = EDOM; L19: clrl r1 clrl r0 ret # return(0.);} L18: pushal 20(fp) # hfs subl3 $72,fp,-(sp) pushl 8(fp) pushl 4(fp) callf $20,_Kfrexpf clrl -60+4(fp) movl r0,-60(fp) # x = frexp(arg,&exp); jbr L20 L2000001: pushal 20(fp) # hfs ldd L23 pushd ldd -60(fp) pushd callf $24,_Kmuld # x * 2 ldd r0 std -60(fp) # x *= 2; decl -72(fp) # exp--;} L20: cmpd2 -60(fp),L22 # while(x < 0.5){ jlss L2000001 bitl $1,-72(fp) # if(exp & 1) { jeql L24 pushal 20(fp) # hfs ldd L25 pushd ldd -60(fp) pushd callf $24,_Kmuld # x * 2 ldd r0 std -60(fp) # x *= 2; decl -72(fp) # exp--;} L24: pushal 20(fp) # hfs ldd -60(fp) pushd ldd L27 pushd callf $24,_Kaddd # (1.0+x) pushal 20(fp) # hfs ldd r0 pushd ldd L26 pushd callf $24,_Kmuld ldd r0 std -68(fp) # temp=0.5*(1.0+x); jbr L28 L2000003: pushal 20(fp) # hfs ldd L30 pushd ldd -68(fp) pushd callf $24,_Kmuld # temp*(1L<<30) ldd r0 std -68(fp) # temp *=(1L<<30); subl2 $60,-72(fp) # exp -= 60;} L28: cmpl -72(fp),$60 # while(exp > 60){ jleq L31 jbr L2000003 L2000005: pushal 20(fp) # hfs ldd L33 pushd ldd -68(fp) pushd callf $24,_Kdivd ldd r0 std -68(fp) # temp /= (1L<<30); addl2 $60,-72(fp) # exp +=60;} L31: cmpl -72(fp),$-60 # while(exp < -60) { jlss L2000005 tstl -72(fp) # if(exp >=0) jlss L34 divl3 $2,-72(fp),r0 shll r0,$1,r0 pushal 20(fp) # hfs pushl $0 pushl r0 pushl $0 pushl $0 callf $24,_Kcvtld pushal 20(fp) # hfs ldd r0 pushd ldd -68(fp) pushd callf $24,_Kmuld # temp * (1L <<(exp/2)) jbr L2000006 # temp *= 1L << (exp/2)); L34: mnegl -72(fp),r0 # -exp divl2 $2,r0 # -exp/2 shll r0,$1,r0 pushal 20(fp) # hfs pushl $0 pushl r0 pushl $0 pushl $0 callf $24,_Kcvtld pushal 20(fp) # hfs ldd r0 pushd ldd -68(fp) pushd callf $24,_Kdivd # temp / (1L <<(-exp/2)) L2000006: ldd r0 std -68(fp) # temp (*/)= 1L << (exp/2)); clrl -76(fp) L2000008: pushal 20(fp) # hfs ldd -68(fp) pushd ldd 4(fp) pushd callf $24,_Kdivd pushal 20(fp) # hfs ldd r0 pushd ldd -68(fp) pushd callf $24,_Kaddd pushal 20(fp) # hfs ldd r0 pushd ldd L39 pushd callf $24,_Kmuld ldd r0 std -68(fp) # temp=0.5*(temp+arg/temp); aobleq $3,-76(fp),L2000008 ldd -68(fp) cvdf stf r0 ret

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.