1 /*
2 * Copyright 2005 Sun Microsystems, Inc. All rights reserved.
3 * Use is subject to license terms.
4 */
5
6 /* Copyright (c) 1984, 1986, 1987, 1988, 1989 AT&T */
7 /* All Rights Reserved */
8
9 /*
10 * Copyright (c) 1980 Regents of the University of California.
11 * All rights reserved. The Berkeley software License Agreement
12 * specifies the terms and conditions for redistribution.
13 */
14
15 #include "e.h"
16
17 void
sqrt(int p2)18 sqrt(int p2)
19 {
20 #ifndef NEQN
21 int nps;
22
23 nps = EFFPS(((eht[p2]*9)/10+(resolution/POINT-1))/(resolution/POINT));
24 #endif /* NEQN */
25 yyval = p2;
26 #ifndef NEQN
27 eht[yyval] = VERT(EM(1.2, nps));
28 if (dbg)
29 printf(".\tsqrt: S%d <- S%d;b=%d, h=%d\n",
30 yyval, p2, ebase[yyval], eht[yyval]);
31 if (rfont[yyval] == ITAL)
32 printf(".as %d \\|\n", yyval);
33 #endif /* NEQN */
34 nrwid(p2, ps, p2);
35 #ifndef NEQN
36 printf(".ds %d \\v'%du'\\s%d\\v'-.2m'\\(sr\\l'\\n(%du\\"
37 "(rn'\\v'.2m'\\s%d", yyval, ebase[p2], nps, p2, ps);
38 printf("\\v'%du'\\h'-\\n(%du'\\*(%d\n", -ebase[p2], p2, p2);
39 lfont[yyval] = ROM;
40 #else /* NEQN */
41 printf(".ds %d \\v'%du'\\e\\L'%du'\\l'\\n(%du'",
42 p2, ebase[p2], -eht[p2], p2);
43 printf("\\v'%du'\\h'-\\n(%du'\\*(%d\n", eht[p2]-ebase[p2], p2, p2);
44 eht[p2] += VERT(1);
45 if (dbg)
46 printf(".\tsqrt: S%d <- S%d;b=%d, h=%d\n",
47 p2, p2, ebase[p2], eht[p2]);
48 #endif /* NEQN */
49 }
50