22 /*
23 * Copyright 2008 Sun Microsystems, Inc. All rights reserved.
24 * Use is subject to license terms.
25 */
26
27 #pragma ident "%Z%%M% %I% %E% SMI"
28
29 #include "lint.h"
30 #include "base_conversion.h"
31 #include <sys/isa_defs.h>
32
33 #define CSR_DEFAULT 0x1f80
34
35 /* The following should be coded as inline expansion templates. */
36
37 /*
38 * Multiplies two normal or subnormal doubles, returns result and exceptions.
39 */
40 double
41 __mul_set(double x, double y, int *pe) {
42 extern void _putmxcsr(), _getmxcsr();
43 int csr;
44 double z;
45
46 _putmxcsr(CSR_DEFAULT);
47 z = x * y;
48 _getmxcsr(&csr);
49 if ((csr & 0x3f) == 0) {
50 *pe = 0;
51 } else {
52 /* Result may not be exact. */
53 *pe = 1;
54 }
55 return (z);
56 }
57
58 /*
59 * Divides two normal or subnormal doubles x/y, returns result and exceptions.
60 */
61 double
62 __div_set(double x, double y, int *pe) {
63 extern void _putmxcsr(), _getmxcsr();
64 int csr;
65 double z;
66
67 _putmxcsr(CSR_DEFAULT);
68 z = x / y;
|
22 /*
23 * Copyright 2008 Sun Microsystems, Inc. All rights reserved.
24 * Use is subject to license terms.
25 */
26
27 #pragma ident "%Z%%M% %I% %E% SMI"
28
29 #include "lint.h"
30 #include "base_conversion.h"
31 #include <sys/isa_defs.h>
32
33 #define CSR_DEFAULT 0x1f80
34
35 /* The following should be coded as inline expansion templates. */
36
37 /*
38 * Multiplies two normal or subnormal doubles, returns result and exceptions.
39 */
40 double
41 __mul_set(double x, double y, int *pe) {
42 int csr = CSR_DEFAULT;
43 double z;
44
45 __asm__ __volatile__(
46 "ldmxcsr %4\n"
47 "mulsd %2,%3\n"
48 "stmxcsr %0\n"
49 : "=m" (csr), "=C" (z)
50 : "C" (x), "C" (y), "m" (csr));
51
52 if ((csr & 0x3f) == 0) {
53 *pe = 0;
54 } else {
55 /* Result may not be exact. */
56 *pe = 1;
57 }
58 return (z);
59 }
60
61 /*
62 * Divides two normal or subnormal doubles x/y, returns result and exceptions.
63 */
64 double
65 __div_set(double x, double y, int *pe) {
66 extern void _putmxcsr(), _getmxcsr();
67 int csr;
68 double z;
69
70 _putmxcsr(CSR_DEFAULT);
71 z = x / y;
|