Return to arc.c CVS log

Up to [Research Unix] / researchv10no / libplot / lib4014

Annotation of researchv10no/libplot/lib4014/arc.c, revision 1.1

1.1     ! root        1: #include <stdio.h>
        !             2: #include <math.h>
        !             3: #include "tek.h"
        !             4: /*             arc plotting routine            */
        !             5: /*             from x1,y1 to x2,y2             */
        !             6: /*     with center xc,yc and radius rr */
        !             7: /*     integrates difference equation          */
        !             8: /*     negative rr draws counterclockwise      */
        !             9: #define PI4 0.7854
        !            10: arc(x1, y1, x2, y2, xc, yc, rr)
        !            11: double x1, x2, y1, y2, xc, yc, rr;
        !            12: {
        !            13:        register double dx, dy, a, b;
        !            14:        double  ph, dph, rd, xnext;
        !            15:        register int    n;
        !            16:        dx = x1 - xc;
        !            17:        dy = y1 - yc;
        !            18:        rd = sqrt(dx * dx + dy * dy);
        !            19:        if (rd / e1->quantum < 1.0) { 
        !            20:                move(xc, yc); 
        !            21:                vec(xc, yc); 
        !            22:                return(0);
        !            23:        }
        !            24:        dph = acos(1.0 - (e1->quantum / rd));
        !            25:        if (dph > PI4) 
        !            26:                dph = PI4;
        !            27:        ph=atan2((y2-yc),(x2 - xc)) - atan2(dy, dx);
        !            28:        if (ph < 0) 
        !            29:                ph += 6.2832; 
        !            30:        if (rr < 0) 
        !            31:                ph = 6.2832 - ph;
        !            32:        if (ph < dph) 
        !            33:                line(x1, y1, x2, y2);
        !            34:        else {
        !            35:                n = ph / dph; 
        !            36:                a = cos(dph); 
        !            37:                b = sin(dph); 
        !            38:                if (rr < 0) 
        !            39:                        b = -b;
        !            40:                move(x1, y1);
        !            41:                while ((n--) >= 0) {
        !            42:                        xnext = dx * a - dy * b; 
        !            43:                        dy = dx * b + dy * a; 
        !            44:                        dx = xnext;
        !            45:                        vec(dx + xc, dy + yc);
        !            46:                }
        !            47:        }
        !            48: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.