xref: /linux/arch/powerpc/math-emu/fsqrts.c (revision 44f57d78)
1 // SPDX-License-Identifier: GPL-2.0
2 #include <linux/types.h>
3 #include <linux/errno.h>
4 #include <linux/uaccess.h>
5 
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>
10 
11 int
12 fsqrts(void *frD, void *frB)
13 {
14 	FP_DECL_D(B);
15 	FP_DECL_D(R);
16 	FP_DECL_EX;
17 
18 #ifdef DEBUG
19 	printk("%s: %p %p %p %p\n", __func__, frD, frB);
20 #endif
21 
22 	FP_UNPACK_DP(B, frB);
23 
24 #ifdef DEBUG
25 	printk("B: %ld %lu %lu %ld (%ld)\n", B_s, B_f1, B_f0, B_e, B_c);
26 #endif
27 
28 	if (B_s && B_c != FP_CLS_ZERO)
29 		FP_SET_EXCEPTION(EFLAG_VXSQRT);
30 	if (B_c == FP_CLS_NAN)
31 		FP_SET_EXCEPTION(EFLAG_VXSNAN);
32 
33 	FP_SQRT_D(R, B);
34 
35 #ifdef DEBUG
36 	printk("R: %ld %lu %lu %ld (%ld)\n", R_s, R_f1, R_f0, R_e, R_c);
37 #endif
38 
39 	__FP_PACK_DS(frD, R);
40 
41 	return FP_CUR_EXCEPTIONS;
42 }
43