gem5 / arm / linux / a579e949378eab7b992d2bf942b655ca407f2948 / . / arch / mips / math-emu / dp_flong.c

/* IEEE754 floating point arithmetic | |

* double precision: common utilities | |

*/ | |

/* | |

* MIPS floating point support | |

* Copyright (C) 1994-2000 Algorithmics Ltd. | |

* | |

* This program is free software; you can distribute it and/or modify it | |

* under the terms of the GNU General Public License (Version 2) as | |

* published by the Free Software Foundation. | |

* | |

* This program is distributed in the hope it will be useful, but WITHOUT | |

* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |

* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |

* for more details. | |

* | |

* You should have received a copy of the GNU General Public License along | |

* with this program; if not, write to the Free Software Foundation, Inc., | |

* 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. | |

*/ | |

#include "ieee754dp.h" | |

union ieee754dp ieee754dp_flong(s64 x) | |

{ | |

u64 xm; | |

int xe; | |

int xs; | |

ieee754_clearcx(); | |

if (x == 0) | |

return ieee754dp_zero(0); | |

if (x == 1 || x == -1) | |

return ieee754dp_one(x < 0); | |

if (x == 10 || x == -10) | |

return ieee754dp_ten(x < 0); | |

xs = (x < 0); | |

if (xs) { | |

if (x == (1ULL << 63)) | |

xm = (1ULL << 63); /* max neg can't be safely negated */ | |

else | |

xm = -x; | |

} else { | |

xm = x; | |

} | |

/* normalize */ | |

xe = DP_FBITS + 3; | |

if (xm >> (DP_FBITS + 1 + 3)) { | |

/* shunt out overflow bits */ | |

while (xm >> (DP_FBITS + 1 + 3)) { | |

XDPSRSX1(); | |

} | |

} else { | |

/* normalize in grs extended double precision */ | |

while ((xm >> (DP_FBITS + 3)) == 0) { | |

xm <<= 1; | |

xe--; | |

} | |

} | |

return ieee754dp_format(xs, xe, xm); | |

} |