fpu-if.s

操作系统SunOS 4.1.3版本的源码

/* static char  frefsccsid[] = "@(#)fpu-if.s 1.1 7/30/92 Copyright Sun Microsystems"; */
!=============================================================================
! 	File: fpu.s					Date:
!=============================================================================
!
!  Objective:
!
!	This module provides an interface between diagnostic application
!	programs and the Floating Point Co-Processor.
!
!	Each routine is callable from C and is responsible for setting
!	up the floating point registers, moving parameters, initiating
!	the co-processor operation and moving the results into the 
!	correct registers.
!
! Open Issues:
!
!	1) Double precision values are passed as 32 bit (single)
!	2) Extended precision operations are Invalid Operations!.
!
! History:
!
!	BDS  04/28/88  	Adding comments to source.
!			Revising code for compatability with TI & Weitek.
!			Restructuring into an interface module.
!			.. Test code is now in fpu-test.s
!
!===========================================================================
	.seg    "data"
	.align	8
temp:	.skip 4			/* storage location for a single precision */
temp1:	.skip 4			/* storage for double precision */
temp2:	.skip 4			/* storage for extended */

	.seg	"text"
	.global	_int_float_s
	.global _int_float_d
	.global _float_int_s
	.global _float_int_d
	.global _convert_sp_dp
	.global _convert_dp_sp
	.global _negate_value
	.global _absolute_value
	.global _square_sp 
	.global _square_dp
	.global _add_sp
	.global _add_dp
	.global _sub_sp
	.global _sub_dp
	.global _mult_sp
	.global _mult_dp
	.global _div_sp
	.global _div_dp
!*	.global _set_psr 
	.global _get_psr 
	.global _set_fsr 
	.global _get_fsr
	.global _cmp_s
	.global _cmp_d
	.global _int_float_e
	.global _cmp_s_ex
	.global _cmp_d_ex

!**************************************************************************
!*  			PSR & FSR Register Interface			 *
!**************************************************************************

!--------------------------------------------------------------------------
! Name:		Get Program Status Register
! Function:	return a copy of the psr to the user
! Calling:	none
! Returns:	i0 = psr contents
! Convention:	psr_value = get_psr() ** 
!--------------------------------------------------------------------------
_get_psr:
	save	%sp, -88, %sp	! save the registers & stack frame
	mov	%psr, %i0	! .. load the PSR
	nop			! .. delay
	nop			! .. delay
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Get the Floating point Status Register
! Function:	return a copy of the FSR to caller
! Calling:	none
! Returns:	i0 = fsr contents
! Convention:	fsr_value = get_fsr() ** 
!--------------------------------------------------------------------------
_get_fsr:
	save	%sp, -88, %sp	! save the registers & stack frame
	set	temp2, %l0	! .. set the address of the result holder
	st	%fsr, [%l0]	! .. set the contents of the FSR register
	ld	[%l0], %i0	! .. return the fsr to caller
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Set Floating point Status Register
! Function:	Set the FSR
! Calling:	i0 = value to write to fsr
! Returns:	none
! Convention:	set_fsr(get_fsr() ** || 0x1001) ** 
!--------------------------------------------------------------------------
_set_fsr:
	save	%sp, -88, %sp	! save the registers & stack frame
	set	temp2, %l0 	! .. set the address of the result holder
	st	%i0, [%l0]	! .. save the value in memory
	ld	[%l0], %fsr	! .. get the contents of the FSR register
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!**************************************************************************
!*			Data Conversion Functions			 *
!**************************************************************************

!--------------------------------------------------------------------------
! Name:		Integer to Float (Single)
! Function:	Convert an integer value to a single precision floating point 
!		value
! Calling:	in0 = value to convert
! Returns:	in0 = converted value
! Convention:	Real = int_float_s(Int) ** 
!--------------------------------------------------------------------------
_int_float_s:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %o5	! .. set the address of the result holder
	set	temp, %o4	! .. set address of temp. mem reg
	st	%i0, [%o4]	! .. put the passed value into memory
	ld	[%o4], %f0	! .. get the value from memory into FPU register
	fitos   %f0, %f2	! .. get the integer into float into fpu r1
	st	%f2, [%o5]	! .. store into the location
	ld	[%o5], %i0	! .. put the value for return
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Integer to Float (double)
! Function:	Convert an integer value to a double precision floating point 
!		value
! Calling:	in0 = value to convert
! Returns:	in0 = converted value
! Convention:	Real = int_float_d(Int) ** 
!--------------------------------------------------------------------------
_int_float_d:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %o5	! .. get the address of temp2
	set	temp, %o4	! .. get the address of temp
	st	%i0, [%o4]	! .. get the user value
	ld	[%o4], %f0	! .. into the float register
	fitod   %f0, %f2	! .... have the fpu perform the operation
	st	%f2, [%o5]	! .. save the result
	ld	[%o5], %i0	! .. and return it to caller
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Integer to Float Exception 
! Function:	This routine is to create the exception for creating an 
! Warning:	** UNIMPLEMENTED INSTRUCTION TRAP (EXTENDED PRECISION) ** **
! Calling:	
! Returns:	
! Convention:	Q) ** Is fitox implemented or not ?
!--------------------------------------------------------------------------
_int_float_e:
	save	%sp, -88, %sp   ! save the registers, stack
	set	temp2, %o5	! .. Get the address of temp2
	set	temp, %o4	! .. Get the address of temp
	st	%i0, [%o4]	! .. Get the users value
	ld	[%o4], %f0	! .. into the float register
	fitox	%f0, %f4	! .... create an unimplemented instruction
	st	%fsr, [%o5]	! .. return the fsr value 
	ld	[%o5], %i0	! to the caller
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		float to integer (single)
! Function:	Convert a real value to an integer
! Calling:	in0 = Value
! Returns:	in0 = Value
! Convention:	Int = float_int_s(real) ** 
!--------------------------------------------------------------------------
_float_int_s:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %o5	! .. get the address of temp2
	set	temp, %o4	! .... and temp
	st	%i0, [%o4]	! .. get the users value
	ld	[%o4], %f0	! .. into the float register
	fstoi   %f0, %f2	! .... have the fpu perform the operation
	st	%f2, [%o5]	! .. save the result
	ld	[%o5], %i0	! .. and return it to the user
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Float to Integer conversion (double)
! Function:	Convert a real value to an integer
! Calling:	in0 = value
! Returns:	in0 = value
! Convention:	Int = float_int_d(real) ** 
!--------------------------------------------------------------------------
_float_int_d:
	save	%sp, -88, %sp	! save the registers, stack 
	set	temp2, %o5 	! .. get the address of temp2
	set	temp, %o4 	! .. and temp
	std     %i0, [%o4] 	! .. get the callers value
	ldd     [%o4], %f0 	! .. into the float register
	fdtoi   %f0, %f2 	! .... have the fpu perform the operation
	std     %f2, [%o5] 	! .. save the result
	ldd     [%o5], %i0 	! .... and return it to caller
	ret 			! Delayed return (get user ret addr)
	restore 		! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Convert Single to double precision
! Function:	<as the name says>
! Calling:	in0 = value
! Returns:	in0 = result
! Convention:	result = convert_sp_dp(value) ** 
!--------------------------------------------------------------------------
_convert_sp_dp:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %l0	! .. get the address of temp2
	set	temp, %l1	! .. get the address of temp
	st	%i0, [%l0]	! .. get the callers value
	ld	[%l0], %f0	! .. into the float register
	fstod   %f0, %f2	! .... have the fpu perform the operation
	st	%f2, [%l1]	! .. save the result
	ld	[%l1], %i0	! .... and return it to the caller
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Convert Double to Single precision
! Function:	..
! Calling:	in0 = double precision value
! Returns:	in0 = result
! Convention:	result = convert_dp_sp(value) ** 
!--------------------------------------------------------------------------
_convert_dp_sp:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %l0	! .. get the address of temp2
	set	temp, %l1	! .. and temp
	st	%i0, [%l0]	! .. get the users value
	ld	[%l0], %f0	! .. move it to a float register
	fdtos	%f0, %f2	! .... have the fpu perform the operation
	st	%f2, [%l1]	! .. save the result
	ld	[%l1], %i0	! .... and return it to the caller
	ret			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Convert Single Precision to Extended Precision
! Function:	<see above>
! Warning:	** UNIMPLEMENTED INSTRUCTION TRAP (EXTENDED PRECISION) ** 
! Calling:	in0 = number to convert
! Returns:	in0 = extended result (64bits)
! Convention:	
!--------------------------------------------------------------------------
_convert_sp_ext:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %l0	! point to the temp storage
	set	temp, %l1	! .. the result temp
	st	%i0, [%l0]	! .. copy the data to temp memory
	ld	[%l0], %f0	! .... load from memory to FPU
	fstox	%f0, %f4	! .... have the fpu perform the operation
	st	%f4, [%l1]	! .. copy the 32 bit result into memory
	ld	[%l1], %i0	! then copy to the return registers
	ret 			! .. return to caller
	restore			! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Negate a value
! Function:	Compliments the Sign bit
! Calling:	in0 = number to cross her
! Returns:	in0 = result
! Convention:	result = negate_value(value) ** 
!--------------------------------------------------------------------------
_negate_value:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %l0	! .. get the address of temp2
	set	temp, %l1	! .. and of temp
	st	%i0, [%l0]	! .. get the callers value
	ld	[%l0], %f0	! .. into the float register
	fnegs   %f0, %f2	! .... have the fpu perform the operation
	st	%f2, [%l1]	! .. save the result
	ld	[%l1], %i0 	! .... and return it to the caller
	ret 			! Delayed return (get user ret addr)
	restore 		! .. restore the frame window

!--------------------------------------------------------------------------
! Name:		Absolute Value
! Function:	Convert a value to its absolute value (clears sign bit)
! Calling:	in0 = value
! Returns:	in0 = result
! Convention:	result = absolute_value(value) ** 
!--------------------------------------------------------------------------
_absolute_value:
	save	%sp, -88, %sp	! save the registers, stack 
	set	temp2, %l0 	! .. get the address of temp2
	set	temp, %l1 	! .. and temp
	st	%i0, [%l0] 	! .. get the users value
	ld	[%l0], %f0 	! .. into a float register
	fabss	%f0, %f2 	! .... have the fpu perform the operation
	st	%f2, [%l1] 	! .. save the result
	ld	[%l1], %i0  	! .... and return it to caller
	ret  			! Delayed return (get user ret addr)
	restore			! .. restore the frame window

!**************************************************************************
!*				Arithmetic Functions			 *
!**************************************************************************

!--------------------------------------------------------------------------
! Name:		Square Single
! Function:	Calculate the Square of a Single precision value
! Calling:	in0 = value
! Returns:	in0 = result
! Convention:	result = square_sp(value) ** 
!--------------------------------------------------------------------------
_square_sp:
	save	%sp, -88, %sp	! save the registers, stack
	set	temp2, %l0	! .. get the address of temp2
	set	temp, %l1	! .. and temp
	st	%i0, [%l0]	! .. get the callers value