 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_add.c,v 1.4 2005/12/11 12:18:42 christos Exp $ */
+/*	$NetBSD: fpu_add.c,v 1.5 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -37,36 +37,36 @@
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_add.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * Perform an FPU add (return x + y).
+ *
  * To subtract, negate y and call add.
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_add.c,v 1.4 2005/12/11 12:18:42 christos Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_add.c,v 1.5 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #if defined(DIAGNOSTIC)||defined(DEBUG)
 #include <sys/systm.h>
 #endif
 #include <machine/reg.h>
 #include <powerpc/instr.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 #include <powerpc/fpu/fpu_extern.h>
 struct fpn *
 fpu_add(struct fpemu *fe)
+{
 	struct fpn *x = &fe->fe_f1, *y = &fe->fe_f2, *r;
 	u_int r0, r1, r2, r3;
 	int rd;
 	/*

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_compare.c,v 1.4 2005/12/11 12:18:42 christos Exp $ */
+/*	$NetBSD: fpu_compare.c,v 1.5 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -38,32 +38,32 @@
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_compare.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * FCMPU and FCMPO instructions.
+ *
  * These rely on the fact that our internal wide format is achieved by
  * adding zero bits to the end of narrower mantissas.
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_compare.c,v 1.4 2005/12/11 12:18:42 christos Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_compare.c,v 1.5 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 /*
  * Perform a compare instruction (with or without unordered exception).
  * This updates the fcc field in the fsr.
+ *
  * If either operand is NaN, the result is unordered.  For ordered, this
  * causes an NV exception.  Everything else is ordered:
  *	|Inf| > |numbers| > |0|.
  * We already arranged for fp_class(Inf) > fp_class(numbers) > fp_class(0),
  * so we get this directly.  Note, however, that two zeros compare equal

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_div.c,v 1.4 2005/12/11 12:18:42 christos Exp $ */
+/*	$NetBSD: fpu_div.c,v 1.5 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -35,35 +35,35 @@
  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_div.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * Perform an FPU divide (return x / y).
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_div.c,v 1.4 2005/12/11 12:18:42 christos Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_div.c,v 1.5 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #if defined(DIAGNOSTIC)||defined(DEBUG)
 #include <sys/systm.h>
 #endif
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 /*
  * Division of normal numbers is done as follows:
+ *
  * x and y are floating point numbers, i.e., in the form 1.bbbb * 2^e.
  * If X and Y are the mantissas (1.bbbb's), the quotient is then:
+ *
  *	q = (X / Y) * 2^((x exponent) - (y exponent))
+ *
  * Since X and Y are both in [1.0,2.0), the quotient's mantissa (X / Y)

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_mul.c,v 1.4 2005/12/11 12:18:42 christos Exp $ */
+/*	$NetBSD: fpu_mul.c,v 1.5 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -35,35 +35,35 @@
  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_mul.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * Perform an FPU multiply (return x * y).
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_mul.c,v 1.4 2005/12/11 12:18:42 christos Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_mul.c,v 1.5 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #if defined(DIAGNOSTIC)||defined(DEBUG)
 #include <sys/systm.h>
 #endif
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 /*
  * The multiplication algorithm for normal numbers is as follows:
+ *
  * The fraction of the product is built in the usual stepwise fashion.
  * Each step consists of shifting the accumulator right one bit
  * (maintaining any guard bits) and, if the next bit in y is set,
  * adding the multiplicand (x) to the accumulator.  Then, in any case,
  * we advance one bit leftward in y.  Algorithmically:
+ *

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_subr.c,v 1.4 2005/12/11 12:18:42 christos Exp $ */
+/*	$NetBSD: fpu_subr.c,v 1.5 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -35,36 +35,36 @@
  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_subr.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * FPU subroutines.
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_subr.c,v 1.4 2005/12/11 12:18:42 christos Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_subr.c,v 1.5 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #if defined(DIAGNOSTIC)||defined(DEBUG)
 #include <sys/systm.h>
 #endif
 #include <machine/reg.h>
 #include <powerpc/instr.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 #include <powerpc/fpu/fpu_extern.h>
 /*
  * Shift the given number right rsh bits.  Any bits that `fall off' will get
  * shoved into the sticky field; we return the resulting sticky.  Note that
  * shifting NaNs is legal (this will never shift all bits out); a NaN's
  * sticky field is ignored anyway.
  */
 int
 fpu_shr(struct fpn *fp, int rsh)

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_emu.c,v 1.19 2016/12/28 10:52:30 rin Exp $ */
+/*	$NetBSD: fpu_emu.c,v 1.20 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright 2001 Wasabi Systems, Inc.
  * All rights reserved.
+ *
  * Written by Eduardo Horvath and Simon Burge for Wasabi Systems, Inc.
+ *
  * Redistribution and use in source and binary forms, with or without
  * modification, are permitted provided that the following conditions
  * are met:
  * 1. Redistributions of source code must retain the above copyright
  *    notice, this list of conditions and the following disclaimer.
  * 2. Redistributions in binary form must reproduce the above copyright
 @@ -66,42 +66,42 @@
  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu.c	8.1 (Berkeley) 6/11/93
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_emu.c,v 1.19 2016/12/28 10:52:30 rin Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_emu.c,v 1.20 2020/06/27 03:07:57 rin Exp $");
 #include "opt_ddb.h"
 #include <sys/param.h>
 #include <sys/systm.h>
 #include <sys/evcnt.h>
 #include <sys/proc.h>
 #include <sys/siginfo.h>
 #include <sys/signal.h>
 #include <sys/signalvar.h>
 #include <sys/siginfo.h>
 #include <sys/systm.h>
 #include <sys/syslog.h>
 #include <sys/evcnt.h>
 #include <powerpc/instr.h>
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <machine/trap.h>
 #include <powerpc/fpu/fpu_emu.h>
 #include <powerpc/fpu/fpu_extern.h>
 #define	FPU_EMU_EVCNT_DECL(name)					\
 static struct evcnt fpu_emu_ev_##name =					\
     EVCNT_INITIALIZER(EVCNT_TYPE_TRAP, NULL, "fpemu", #name);		\
 EVCNT_ATTACH_STATIC(fpu_emu_ev_##name)
 #define	FPU_EMU_EVCNT_INCR(name)					\
     fpu_emu_ev_##name.ev_count++

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_explode.c,v 1.7 2016/12/28 10:30:04 rin Exp $ */
+/*	$NetBSD: fpu_explode.c,v 1.8 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -36,35 +36,35 @@
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_explode.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * FPU subroutines: `explode' the machine's `packed binary' format numbers
  * into our internal format.
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_explode.c,v 1.7 2016/12/28 10:30:04 rin Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_explode.c,v 1.8 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #include <sys/systm.h>
 #include <machine/ieee.h>
 #include <powerpc/instr.h>
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/ieee.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 #include <powerpc/fpu/fpu_extern.h>
 /*
  * N.B.: in all of the following, we assume the FP format is
+ *
  *	---------------------------
  *	| s | exponent | fraction |
  *	---------------------------
+ *
  * (which represents -1**s * 1.fraction * 2**exponent), so that the

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_implode.c,v 1.6 2005/12/11 12:18:42 christos Exp $ */
+/*	$NetBSD: fpu_implode.c,v 1.7 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -36,35 +36,35 @@
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_implode.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * FPU subroutines: `implode' internal format numbers into the machine's
  * `packed binary' format.
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_implode.c,v 1.6 2005/12/11 12:18:42 christos Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_implode.c,v 1.7 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #include <sys/systm.h>
 #include <machine/ieee.h>
 #include <powerpc/instr.h>
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/ieee.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 #include <powerpc/fpu/fpu_extern.h>
 static int round(struct fpemu *, struct fpn *);
 static int toinf(struct fpemu *, int);
 /*
  * Round a number (algorithm from Motorola MC68882 manual, modified for
  * our internal format).  Set inexact exception if rounding is required.
  * Return true iff we rounded up.
+ *

 @@ -1,14 +1,14 @@
-/*	$NetBSD: fpu_sqrt.c,v 1.5 2014/02/25 14:16:25 martin Exp $ */
+/*	$NetBSD: fpu_sqrt.c,v 1.6 2020/06/27 03:07:57 rin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  All rights reserved.
+ *
  * This software was developed by the Computer Systems Engineering group
  * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
  * contributed to Berkeley.
+ *
  * All advertising materials mentioning features or use of this software
  * must display the following acknowledgement:
  *	This product includes software developed by the University of
  *	California, Lawrence Berkeley Laboratory.
 @@ -35,35 +35,35 @@
  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  * SUCH DAMAGE.
+ *
  *	@(#)fpu_sqrt.c	8.1 (Berkeley) 6/11/93
  */
 /*
  * Perform an FPU square root (return sqrt(x)).
  */
 #include <sys/cdefs.h>
-__KERNEL_RCSID(0, "$NetBSD: fpu_sqrt.c,v 1.5 2014/02/25 14:16:25 martin Exp $");
+__KERNEL_RCSID(0, "$NetBSD: fpu_sqrt.c,v 1.6 2020/06/27 03:07:57 rin Exp $");
 #include <sys/types.h>
 #if defined(DIAGNOSTIC)||defined(DEBUG)
 #include <sys/systm.h>
 #endif
 #include <machine/reg.h>
 #include <machine/fpu.h>
 #include <machine/reg.h>
 #include <powerpc/fpu/fpu_arith.h>
 #include <powerpc/fpu/fpu_emu.h>
 /*
  * Our task is to calculate the square root of a floating point number x0.
  * This number x normally has the form:
+ *
  *		    exp
  *	x = mant * 2		(where 1 <= mant < 2 and exp is an integer)
+ *
  * This can be left as it stands, or the mantissa can be doubled and the
  * exponent decremented: