NetBSD/lib/libc/arch/ns32k/gen/modf.S

73 lines
2.3 KiB
ArmAsm

/* $NetBSD: modf.S,v 1.2 1997/05/08 13:38:40 matthias 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.2 1997/05/08 13:38:40 matthias 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 it's 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