xref: /titanic_51/usr/src/cmd/eqn/sqrt.c (revision 779fc935796a940997d18b31d64a9fec9c6b40f6)
1*779fc935Sceastha /*
2*779fc935Sceastha  * Copyright 2005 Sun Microsystems, Inc.  All rights reserved.
3*779fc935Sceastha  * Use is subject to license terms.
4*779fc935Sceastha  */
5*779fc935Sceastha 
67c478bd9Sstevel@tonic-gate /*	Copyright (c) 1984, 1986, 1987, 1988, 1989 AT&T	*/
77c478bd9Sstevel@tonic-gate /*	  All Rights Reserved  	*/
87c478bd9Sstevel@tonic-gate 
97c478bd9Sstevel@tonic-gate /*
107c478bd9Sstevel@tonic-gate  * Copyright (c) 1980 Regents of the University of California.
117c478bd9Sstevel@tonic-gate  * All rights reserved. The Berkeley software License Agreement
127c478bd9Sstevel@tonic-gate  * specifies the terms and conditions for redistribution.
137c478bd9Sstevel@tonic-gate  */
147c478bd9Sstevel@tonic-gate 
15*779fc935Sceastha #pragma ident	"%Z%%M%	%I%	%E% SMI"
167c478bd9Sstevel@tonic-gate 
177c478bd9Sstevel@tonic-gate #include "e.h"
187c478bd9Sstevel@tonic-gate 
19*779fc935Sceastha void
20*779fc935Sceastha sqrt(int p2)
21*779fc935Sceastha {
227c478bd9Sstevel@tonic-gate #ifndef NEQN
237c478bd9Sstevel@tonic-gate 	int nps;
247c478bd9Sstevel@tonic-gate 
257c478bd9Sstevel@tonic-gate 	nps = EFFPS(((eht[p2]*9)/10+(resolution/POINT-1))/(resolution/POINT));
26*779fc935Sceastha #endif	/* NEQN */
277c478bd9Sstevel@tonic-gate 	yyval = p2;
287c478bd9Sstevel@tonic-gate #ifndef NEQN
297c478bd9Sstevel@tonic-gate 	eht[yyval] = VERT(EM(1.2, nps));
30*779fc935Sceastha 	if (dbg)
31*779fc935Sceastha 		printf(".\tsqrt: S%d <- S%d;b=%d, h=%d\n",
327c478bd9Sstevel@tonic-gate 		    yyval, p2, ebase[yyval], eht[yyval]);
337c478bd9Sstevel@tonic-gate 	if (rfont[yyval] == ITAL)
347c478bd9Sstevel@tonic-gate 		printf(".as %d \\|\n", yyval);
35*779fc935Sceastha #endif	/* NEQN */
367c478bd9Sstevel@tonic-gate 	nrwid(p2, ps, p2);
377c478bd9Sstevel@tonic-gate #ifndef NEQN
38*779fc935Sceastha 	printf(".ds %d \\v'%du'\\s%d\\v'-.2m'\\(sr\\l'\\n(%du\\"
39*779fc935Sceastha 	    "(rn'\\v'.2m'\\s%d", yyval, ebase[p2], nps, p2, ps);
407c478bd9Sstevel@tonic-gate 	printf("\\v'%du'\\h'-\\n(%du'\\*(%d\n", -ebase[p2], p2, p2);
417c478bd9Sstevel@tonic-gate 	lfont[yyval] = ROM;
42*779fc935Sceastha #else	/* NEQN */
437c478bd9Sstevel@tonic-gate 	printf(".ds %d \\v'%du'\\e\\L'%du'\\l'\\n(%du'",
447c478bd9Sstevel@tonic-gate 	    p2, ebase[p2], -eht[p2], p2);
457c478bd9Sstevel@tonic-gate 	printf("\\v'%du'\\h'-\\n(%du'\\*(%d\n", eht[p2]-ebase[p2], p2, p2);
467c478bd9Sstevel@tonic-gate 	eht[p2] += VERT(1);
47*779fc935Sceastha 	if (dbg)
48*779fc935Sceastha 		printf(".\tsqrt: S%d <- S%d;b=%d, h=%d\n",
497c478bd9Sstevel@tonic-gate 		    p2, p2, ebase[p2], eht[p2]);
50*779fc935Sceastha #endif	/* NEQN */
517c478bd9Sstevel@tonic-gate }
52