[mirror_ubuntu-focal-kernel.git] / arch / powerpc / math-emu / fnmsubs.c

// SPDX-License-Identifier: GPL-2.0
#include <linux/types.h>
#include <linux/errno.h>
#include <linux/uaccess.h>

#include <asm/sfp-machine.h>
#include <math-emu/soft-fp.h>
#include <math-emu/double.h>
#include <math-emu/single.h>

int
fnmsubs(void *frD, void *frA, void *frB, void *frC)
{
	FP_DECL_D(R);
	FP_DECL_D(A);
	FP_DECL_D(B);
	FP_DECL_D(C);
	FP_DECL_D(T);
	FP_DECL_EX;

#ifdef DEBUG
	printk("%s: %p %p %p %p\n", __func__, frD, frA, frB, frC);
#endif

	FP_UNPACK_DP(A, frA);
	FP_UNPACK_DP(B, frB);
	FP_UNPACK_DP(C, frC);

#ifdef DEBUG
	printk("A: %ld %lu %lu %ld (%ld)\n", A_s, A_f1, A_f0, A_e, A_c);
	printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c);
	printk("C: %ld %lu %lu %ld (%ld)\n", C_s, C_f1, C_f0, C_e, C_c);
#endif

	if ((A_c == FP_CLS_INF && C_c == FP_CLS_ZERO) ||
	    (A_c == FP_CLS_ZERO && C_c == FP_CLS_INF))
		FP_SET_EXCEPTION(EFLAG_VXIMZ);

	FP_MUL_D(T, A, C);

	if (B_c != FP_CLS_NAN)
		B_s ^= 1;

	if (T_s != B_s && T_c == FP_CLS_INF && B_c == FP_CLS_INF)
		FP_SET_EXCEPTION(EFLAG_VXISI);

	FP_ADD_D(R, T, B);

	if (R_c != FP_CLS_NAN)
		R_s ^= 1;

#ifdef DEBUG
	printk("D: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c);
#endif

	__FP_PACK_DS(frD, R);

	return FP_CUR_EXCEPTIONS;
}
Commit	Line	Data
b2441318	1	// SPDX-License-Identifier: GPL-2.0
1da177e4 LT	2	#include <linux/types.h>
1da177e4 LT	3	#include <linux/errno.h>
7c0f6ba6	4	#include <linux/uaccess.h>
1da177e4	5
d2b194ed KG	6	#include <asm/sfp-machine.h>
	7	#include <math-emu/soft-fp.h>
	8	#include <math-emu/double.h>
	9	#include <math-emu/single.h>
1da177e4 LT	10
	11	int
	12	fnmsubs(void frD, void frA, void frB, void frC)
	13	{
	14	FP_DECL_D(R);
	15	FP_DECL_D(A);
	16	FP_DECL_D(B);
	17	FP_DECL_D(C);
	18	FP_DECL_D(T);
d2b194ed	19	FP_DECL_EX;
1da177e4 LT	20
1da177e4 LT	21	#ifdef DEBUG
e48b1b45	22	printk("%s: %p %p %p %p\n", __func__, frD, frA, frB, frC);
1da177e4 LT	23	#endif
1da177e4 LT	24
d2b194ed KG	25	FP_UNPACK_DP(A, frA);
	26	FP_UNPACK_DP(B, frB);
	27	FP_UNPACK_DP(C, frC);
1da177e4 LT	28
	29	#ifdef DEBUG
	30	printk("A: %ld %lu %lu %ld (%ld)\n", A_s, A_f1, A_f0, A_e, A_c);
	31	printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c);
	32	printk("C: %ld %lu %lu %ld (%ld)\n", C_s, C_f1, C_f0, C_e, C_c);
	33	#endif
	34
	35	if ((A_c == FP_CLS_INF && C_c == FP_CLS_ZERO) \|\|
	36	(A_c == FP_CLS_ZERO && C_c == FP_CLS_INF))
033b8a33	37	FP_SET_EXCEPTION(EFLAG_VXIMZ);
1da177e4 LT	38
	39	FP_MUL_D(T, A, C);
	40
	41	if (B_c != FP_CLS_NAN)
	42	B_s ^= 1;
	43
	44	if (T_s != B_s && T_c == FP_CLS_INF && B_c == FP_CLS_INF)
033b8a33	45	FP_SET_EXCEPTION(EFLAG_VXISI);
1da177e4 LT	46
	47	FP_ADD_D(R, T, B);
	48
	49	if (R_c != FP_CLS_NAN)
	50	R_s ^= 1;
	51
	52	#ifdef DEBUG
	53	printk("D: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c);
	54	#endif
	55
d2b194ed KG	56	__FP_PACK_DS(frD, R);
	57
	58	return FP_CUR_EXCEPTIONS;
1da177e4	59	}