1*0Sstevel@tonic-gate /* 2*0Sstevel@tonic-gate * CDDL HEADER START 3*0Sstevel@tonic-gate * 4*0Sstevel@tonic-gate * The contents of this file are subject to the terms of the 5*0Sstevel@tonic-gate * Common Development and Distribution License, Version 1.0 only 6*0Sstevel@tonic-gate * (the "License"). You may not use this file except in compliance 7*0Sstevel@tonic-gate * with the License. 8*0Sstevel@tonic-gate * 9*0Sstevel@tonic-gate * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 10*0Sstevel@tonic-gate * or http://www.opensolaris.org/os/licensing. 11*0Sstevel@tonic-gate * See the License for the specific language governing permissions 12*0Sstevel@tonic-gate * and limitations under the License. 13*0Sstevel@tonic-gate * 14*0Sstevel@tonic-gate * When distributing Covered Code, include this CDDL HEADER in each 15*0Sstevel@tonic-gate * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 16*0Sstevel@tonic-gate * If applicable, add the following below this CDDL HEADER, with the 17*0Sstevel@tonic-gate * fields enclosed by brackets "[]" replaced with your own identifying 18*0Sstevel@tonic-gate * information: Portions Copyright [yyyy] [name of copyright owner] 19*0Sstevel@tonic-gate * 20*0Sstevel@tonic-gate * CDDL HEADER END 21*0Sstevel@tonic-gate */ 22*0Sstevel@tonic-gate /* Copyright (c) 1984, 1986, 1987, 1988, 1989 AT&T */ 23*0Sstevel@tonic-gate /* All Rights Reserved */ 24*0Sstevel@tonic-gate 25*0Sstevel@tonic-gate 26*0Sstevel@tonic-gate /* 27*0Sstevel@tonic-gate * Copyright 2005 Sun Microsystems, Inc. All rights reserved. 28*0Sstevel@tonic-gate * Use is subject to license terms. 29*0Sstevel@tonic-gate */ 30*0Sstevel@tonic-gate 31*0Sstevel@tonic-gate #pragma ident "%Z%%M% %I% %E% SMI" /* SVr4.0 1.4 */ 32*0Sstevel@tonic-gate /*LINTLIBRARY*/ 33*0Sstevel@tonic-gate 34*0Sstevel@tonic-gate #include <stdlib.h> 35*0Sstevel@tonic-gate #include <math.h> 36*0Sstevel@tonic-gate #include <plot.h> 37*0Sstevel@tonic-gate #include "con.h" 38*0Sstevel@tonic-gate 39*0Sstevel@tonic-gate static int del = 20; 40*0Sstevel@tonic-gate 41*0Sstevel@tonic-gate static void 42*0Sstevel@tonic-gate step(int d) 43*0Sstevel@tonic-gate { 44*0Sstevel@tonic-gate del = d; 45*0Sstevel@tonic-gate } 46*0Sstevel@tonic-gate 47*0Sstevel@tonic-gate /* local definition of quad */ 48*0Sstevel@tonic-gate static int 49*0Sstevel@tonic-gate quad_l(short x, short y, short xp, short yp) 50*0Sstevel@tonic-gate { 51*0Sstevel@tonic-gate if (x < xp) 52*0Sstevel@tonic-gate if (y <= yp) 53*0Sstevel@tonic-gate return (1); 54*0Sstevel@tonic-gate else return (4); 55*0Sstevel@tonic-gate else if (x > xp) 56*0Sstevel@tonic-gate if (y < yp) 57*0Sstevel@tonic-gate return (2); 58*0Sstevel@tonic-gate else return (3); 59*0Sstevel@tonic-gate else if (y < yp) 60*0Sstevel@tonic-gate return (2); 61*0Sstevel@tonic-gate else return (4); 62*0Sstevel@tonic-gate } 63*0Sstevel@tonic-gate 64*0Sstevel@tonic-gate void 65*0Sstevel@tonic-gate arc(short x, short y, short x0, short y0, short x1, short y1) 66*0Sstevel@tonic-gate { 67*0Sstevel@tonic-gate double pc; 68*0Sstevel@tonic-gate int flg, m, xc, yc, xs, ys, qs, qf, qt, qtctr = 0; 69*0Sstevel@tonic-gate int m0, m1; 70*0Sstevel@tonic-gate float dx, dy, r; 71*0Sstevel@tonic-gate char use; 72*0Sstevel@tonic-gate dx = x - x0; 73*0Sstevel@tonic-gate dy = y - y0; 74*0Sstevel@tonic-gate r = dx * dx + dy * dy; 75*0Sstevel@tonic-gate pc = r; 76*0Sstevel@tonic-gate pc = sqrt(pc); 77*0Sstevel@tonic-gate flg = (short)(pc / 4); 78*0Sstevel@tonic-gate if (flg == 0) 79*0Sstevel@tonic-gate step(1); 80*0Sstevel@tonic-gate else if (flg < del) 81*0Sstevel@tonic-gate step(flg); 82*0Sstevel@tonic-gate xc = xs = x0; 83*0Sstevel@tonic-gate yc = ys = y0; 84*0Sstevel@tonic-gate move((short)xs, (short)ys); 85*0Sstevel@tonic-gate if ((x0 == x1) && (y0 == y1)) 86*0Sstevel@tonic-gate flg = 0; 87*0Sstevel@tonic-gate else flg = 1; 88*0Sstevel@tonic-gate qs = quad_l(x, y, x0, y0); 89*0Sstevel@tonic-gate qf = quad_l(x, y, x1, y1); 90*0Sstevel@tonic-gate if (abs(x - x1) < abs(y - y1)) { 91*0Sstevel@tonic-gate use = 'x'; 92*0Sstevel@tonic-gate if ((qs == 2) || (qs == 3)) 93*0Sstevel@tonic-gate m = -1; 94*0Sstevel@tonic-gate else m = 1; 95*0Sstevel@tonic-gate } else { 96*0Sstevel@tonic-gate use = 'y'; 97*0Sstevel@tonic-gate if (qs > 2) 98*0Sstevel@tonic-gate m = -1; 99*0Sstevel@tonic-gate else m = 1; 100*0Sstevel@tonic-gate } 101*0Sstevel@tonic-gate if (qs == qf) { 102*0Sstevel@tonic-gate m0 = (y0 - y) / (x0 - x); 103*0Sstevel@tonic-gate m1 = (y1 - y) / (x1 - x); 104*0Sstevel@tonic-gate if (m0 >= m1) 105*0Sstevel@tonic-gate qt = 4; 106*0Sstevel@tonic-gate else qt = 0; 107*0Sstevel@tonic-gate } else if ((qt = qf - qs) < 0) 108*0Sstevel@tonic-gate qt += 4; 109*0Sstevel@tonic-gate /*LINTED*/ 110*0Sstevel@tonic-gate while (1) { 111*0Sstevel@tonic-gate switch (use) { 112*0Sstevel@tonic-gate case 'x': 113*0Sstevel@tonic-gate if ((qs == 2) || (qs == 3)) 114*0Sstevel@tonic-gate yc -= del; 115*0Sstevel@tonic-gate else yc += del; 116*0Sstevel@tonic-gate dy = yc - y; 117*0Sstevel@tonic-gate pc = r - dy * dy; 118*0Sstevel@tonic-gate xc = (int)(m * sqrt(pc) + x); 119*0Sstevel@tonic-gate if (((x < xs) && (x >= xc)) || 120*0Sstevel@tonic-gate ((x > xs) && (x <= xc)) || 121*0Sstevel@tonic-gate ((y < ys) && (y >= yc)) || 122*0Sstevel@tonic-gate ((y > ys) && (y <= yc))) { 123*0Sstevel@tonic-gate if (++qtctr > qt) 124*0Sstevel@tonic-gate return; 125*0Sstevel@tonic-gate if (++qs > 4) 126*0Sstevel@tonic-gate qs = 1; 127*0Sstevel@tonic-gate if ((qs == 2) || (qs == 3)) 128*0Sstevel@tonic-gate m = -1; 129*0Sstevel@tonic-gate else m = 1; 130*0Sstevel@tonic-gate flg = 1; 131*0Sstevel@tonic-gate } 132*0Sstevel@tonic-gate cont((short)xc, (short)yc); 133*0Sstevel@tonic-gate xs = xc; 134*0Sstevel@tonic-gate ys = yc; 135*0Sstevel@tonic-gate if ((qs == qf) && (flg == 1)) 136*0Sstevel@tonic-gate switch (qf) { 137*0Sstevel@tonic-gate case 3: 138*0Sstevel@tonic-gate case 4: 139*0Sstevel@tonic-gate if (xs >= x1) 140*0Sstevel@tonic-gate return; 141*0Sstevel@tonic-gate continue; 142*0Sstevel@tonic-gate case 1: 143*0Sstevel@tonic-gate case 2: 144*0Sstevel@tonic-gate if (xs <= x1) 145*0Sstevel@tonic-gate return; 146*0Sstevel@tonic-gate } 147*0Sstevel@tonic-gate continue; 148*0Sstevel@tonic-gate case 'y': 149*0Sstevel@tonic-gate if (qs > 2) 150*0Sstevel@tonic-gate xc += del; 151*0Sstevel@tonic-gate else xc -= del; 152*0Sstevel@tonic-gate dx = xc - x; 153*0Sstevel@tonic-gate pc = r - dx * dx; 154*0Sstevel@tonic-gate yc = (int)(m * sqrt(pc) + y); 155*0Sstevel@tonic-gate if (((x < xs) && (x >= xc)) || 156*0Sstevel@tonic-gate ((x > xs) && (x <= xc)) || 157*0Sstevel@tonic-gate ((y < ys) && (y >= yc)) || 158*0Sstevel@tonic-gate ((y > ys) && (y <= yc))) { 159*0Sstevel@tonic-gate if (++qtctr > qt) 160*0Sstevel@tonic-gate return; 161*0Sstevel@tonic-gate if (++qs > 4) 162*0Sstevel@tonic-gate qs = 1; 163*0Sstevel@tonic-gate if (qs > 2) 164*0Sstevel@tonic-gate m = -1; 165*0Sstevel@tonic-gate else m = 1; 166*0Sstevel@tonic-gate flg = 1; 167*0Sstevel@tonic-gate } 168*0Sstevel@tonic-gate cont((short)xc, (short)yc); 169*0Sstevel@tonic-gate xs = xc; 170*0Sstevel@tonic-gate ys = yc; 171*0Sstevel@tonic-gate if ((qs == qf) && (flg == 1)) 172*0Sstevel@tonic-gate switch (qs) { 173*0Sstevel@tonic-gate case 1: 174*0Sstevel@tonic-gate case 4: 175*0Sstevel@tonic-gate if (ys >= y1) 176*0Sstevel@tonic-gate return; 177*0Sstevel@tonic-gate continue; 178*0Sstevel@tonic-gate case 2: 179*0Sstevel@tonic-gate case 3: 180*0Sstevel@tonic-gate if (ys <= y1) 181*0Sstevel@tonic-gate return; 182*0Sstevel@tonic-gate } 183*0Sstevel@tonic-gate } 184*0Sstevel@tonic-gate } 185*0Sstevel@tonic-gate } 186