73 lines
2.3 KiB
ArmAsm
73 lines
2.3 KiB
ArmAsm
/* $NetBSD: modf.S,v 1.3 2006/04/24 19:00:30 snj Exp $ */
|
|
|
|
/*
|
|
* Copyright (c) 1992 Helsinki University of Technology
|
|
* All Rights Reserved.
|
|
*
|
|
* Permission to use, copy, modify and distribute this software and its
|
|
* documentation is hereby granted, provided that both the copyright
|
|
* notice and this permission notice appear in all copies of the
|
|
* software, derivative works or modified versions, and any portions
|
|
* thereof, and that both notices appear in supporting documentation.
|
|
*
|
|
* HELSINKI UNIVERSITY OF TECHNOLOGY ALLOWS FREE USE OF THIS SOFTWARE IN
|
|
* ITS "AS IS" CONDITION. HELSINKI UNIVERSITY OF TECHNOLOGY DISCLAIMS ANY
|
|
* LIABILITY OF ANY KIND FOR ANY DAMAGES WHATSOEVER RESULTING FROM THE
|
|
* USE OF THIS SOFTWARE.
|
|
*/
|
|
/*
|
|
* HISTORY
|
|
* 29-Apr-92 Tero Kivinen (kivinen) at Helsinki University of Technology
|
|
* Created.
|
|
*/
|
|
|
|
/*
|
|
* double modf (value, iptr)
|
|
* double value, *iptr;
|
|
*
|
|
* Modf returns the fractional part of "value",
|
|
* and stores the integer part indirectly through "iptr".
|
|
*/
|
|
|
|
#include <machine/asm.h>
|
|
|
|
#if defined(LIBC_SCCS)
|
|
RCSID("$NetBSD: modf.S,v 1.3 2006/04/24 19:00:30 snj Exp $")
|
|
#endif
|
|
|
|
ENTRY(modf)
|
|
enter [],0
|
|
movl 8(fp),f0 /* value */
|
|
movd 12(fp),r0 /* higher 32 bit of value */
|
|
lshd -20,r0 /* extract exponent */
|
|
andd 0x7ff,r0 /* 11 lower bits */
|
|
cmpd r0,1023+30 /* compare if its int part can fit in int */
|
|
bgt 0f /* nope else it's ok to truncld it to int */
|
|
truncld f0,r0 /* get integer part */
|
|
movdl r0,f2 /* convert back to float */
|
|
movl f2,0(16(fp)) /* move integer part to *iptr */
|
|
subl f2,f0 /* return fract. part = value - *iptr */
|
|
exit []
|
|
ret 0
|
|
|
|
0: subd 1023+20,r0 /* bias 1023, and upper part of
|
|
floating point mantissa part */
|
|
movqd -1,r2 /* generate mask to get fract. part */
|
|
cmpd r0,32 /* if value > 2^52 (20+32) then all int part */
|
|
bhi 1f
|
|
negd r0,r0 /* shift right */
|
|
lshd r0,r2 /* here */
|
|
comd r2,r2 /* get fractional part complement mask */
|
|
movd 8(fp),r1 /* get lower 32 bit of value */
|
|
andd r2,r1 /* mask fractional part off leave ip part */
|
|
movd r1,0(16(fp)) /* store ip part to *iptr */
|
|
movd 12(fp),4(16(fp)) /* store other half to *iptr */
|
|
subl 0(16(fp)),f0 /* return fract. part = value - *iptr */
|
|
exit []
|
|
ret 0
|
|
|
|
1: movl 8(fp),0(16(fp)) /* copy value to *iptr */
|
|
movdl 0,f0 /* return 0 for fract. part */
|
|
exit []
|
|
ret 0
|