blob: ca61513734eacba639a6d442ec2354be0d154673 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
|
/* $OpenBSD: e_sqrt.c,v 1.3 2012/12/05 23:20:03 deraadt Exp $ */
/*
* Written by Martynas Venckus. Public domain
*/
/* LINTLIBRARY */
#include <sys/types.h>
#include <math.h>
#define FPSCR_PR (1 << 19)
#define FPSCR_SZ (1 << 20)
double
sqrt(double d)
{
register_t fpscr, nfpscr;
__asm__ __volatile__ ("sts fpscr, %0" : "=r" (fpscr));
/* Set floating-point mode to double-precision. */
nfpscr = fpscr | FPSCR_PR;
/* Do not set SZ and PR to 1 simultaneously. */
nfpscr &= ~FPSCR_SZ;
__asm__ __volatile__ ("lds %0, fpscr" : : "r" (nfpscr));
__asm__ __volatile__ ("fsqrt %0" : "+f" (d));
/* Restore fp status/control register. */
__asm__ __volatile__ ("lds %0, fpscr" : : "r" (fpscr));
return (d);
}
/* No extended-precision is present. */
#ifdef lint
/* PROTOLIB1 */
long double sqrtl(long double);
#else /* lint */
__weak_alias(sqrtl,sqrt);
#endif /* lint */
|
