/*-
 * Copyright (c) 1985 The Regents of the University of California.
 * All rights reserved.
 *
 * This code is derived from software contributed to Berkeley by
 * Computer Consoles Inc.
 *
 * %sccs.include.redist.c%
 *
 *	@(#)Klogf.s	7.1 (Berkeley) 12/06/90
 */

#include "../tahoe/SYS.h"

	.text
ENTRY(Klogf, R5|R4|R3|R2)
	subl3	$88,fp,sp
	clrl	8(fp)
	tstl	4(fp)
	jgtr	L53
	movl	small+4,r1
	movl	small,r0
	ret
L53:	pushl	20(fp)			# hfs
	subl3	$88,fp,-(sp)		# &exp
	pushl	8(fp)			# arg
	pushl	4(fp)
	callf	$20,_Kfrexpf
	ldd	r0
					# movl	r1,-56(fp)
	std	-60(fp)			# movl	r0,-60(fp)
	jbr	L55
L2000001:
	pushl	20(fp)			# hfs
	ldd	two			# 2.0
	pushd
	ldd	-60(fp)			# x
	pushd	
	callf	$24,_Kmuld
	ldd	r0
	std	-60(fp)
	subl2	$1,-88(fp)
L55:	cmpd2	-60(fp),half
	jlss	L2000001
	cmpd2	-60(fp),_sqrto2
	jgeq	L59
	pushl	20(fp)			# hfs
	ldd	-60(fp)			# x
	pushd
	ldd	two			# 2.0
	pushd
	callf	$24,_Kmuld
	ldd	r0
	std	-60(fp)
	subl2	$1,-88(fp)
L59:	pushl	20(fp)			# hfs
	ldd	one			# 1.0
	pushd
	ldd	-60(fp)			# x
	pushd
	callf	$24,_Ksubd	
	ldd	r0
	std	r2
	pushl	20(fp)			# hfs
	ldd	one			# 1.0
	pushd
	ldd	-60(fp)			# x
	pushd
	callf	$24,_Kaddd
	ldd	r0
	std	r4
	pushl	20(fp)			# hfs
	ldd	r4			# temp result of x-1
	pushd
	ldd	r2			# temp result of x+1
	pushd
	callf	$24,_Kdivd
	ldd	r0
	std	-68(fp)			# z
	pushl	20(fp)			# hfs
	pushd				# z
	pushd				# z
	callf	$24,_Kmuld
	ldd	r0
	std	-76(fp)			# zsq
	pushl	20(fp)			# hfs
	ldd	-76(fp)			# zsq
	pushd
	ldd	_p1			# p1
	pushd
	callf	$24,_Kmuld
	ldd	r0
	std	r2
	pushl	20(fp)			# hfs
	ldd 	_p0			# p0
	pushd
	ldd	r2
	pushd
	callf	$24,_Kaddd
	ldd	r0
	std	-84(fp)			# temp
	pushl	20(fp)			# hfs
	ldd	-76(fp)			# zsq
	pushd
	ldd	_q1
	pushd
	callf	$24,_Kmuld
	ldd	r0
	std	r2
	pushl	20(fp)			# hfs
	ldd	_q0			# q0
	pushd
	ldd	r2
	pushd
	callf	$24,_Kaddd
	ldd	r0
	std	r2			# temp result of (q1*zsq+q0)
	pushl	20(fp)			# hfs
	ldd	r2
	pushd
	ldd	-84(fp)			# temp
	pushd
	callf	$24,_Kdivd
	ldd	r0
	std	-84(fp)			# put in temp
	pushl	20(fp)			# hfs
	ldd	-68(fp)			# z
	pushd
	ldd	-84(fp)
	pushd
	callf	$24,_Kmuld
	ldd	r0
	std	r2
	pushl	20(fp)			#hfs
	pushl	$0			# 0 pad
	pushl	-88(fp)			# exp
	pushl	$0			# dummy
	pushl	$0			# dummy
	callf	$24,_Kcvtld
	pushl	20(fp)			# hfs
	ldd	_log2			# log2
	pushd
	pushl	r1			# exp converted
	pushl	r0			
	callf	$24,_Kmuld
	ldd	r0
	std	r4
	pushl	20(fp)			# hfs
	ldd	r2			# result of temp*z
	pushd
	ldd	r4
	pushd
	callf	$24,_Kaddd
	ldd	r0
	cvdf
	stf	r0
	clrl	r1
	ret

	.data
	.align	2
_log2:	.long	0x40317217, 0xf7d1cf7a # .double 0.69314718055994531
_ln10:	.long	0x41135d8d, 0xddaaa8ac # .double 2.3025850929940457
_sqrto2:.long	0x403504f3, 0x33f9de65 # .double 0.70710678118654753
_p0:	.long	0xc154114d, 0xeb0ba468 # .double -3.31355617479
_p1:	.long	0x40654226, 0x56bd0c4c # .double 0.89554061525
_q0:	.long	0xc0d4114c, 0xfdc7df02 # .double -1.65677797691
_q1:	.long	0x40800000, 0x00000000 # .double 1
small:	.long	0xfffffffe, 0xfffffffe # .double -1.7014117331926443e+38
half:	.long	0x40000000, 0x00000000 # .double 0.5
two:	.long	0x41000000, 0x00000000 # .double 2
one:	.long	0x40800000, 0x00000000 # .double 1