src/meta/gcalc.c

#ifndef lint
static const char       RCSid[] = "$Id$";
#endif
/*
 *  gcalc.c - routines to do calculations on graph files.
 *
 *     7/7/86
 *     Greg Ward Larson
 */

#include  <stdio.h>

#include  "mgvars.h"

#ifndef BSD
#define index   strchr
#endif

extern char  *index();

static double  xsum, xxsum, ysum, yysum, xysum;
static double  xmin, xmax, ymin, ymax;
static double  lastx, lasty, rsum;
static int  npts;


gcalc(types)                    /* do a calculation */
char  *types;
{
        int  i, calcpoint();

        if (index(types, 'h') == NULL)
                gcheader(types);
        
        xmin = gparam[XMIN].flags & DEFINED ?
                        varvalue(gparam[XMIN].name) :
                        -1e10;
        xmax = gparam[XMAX].flags & DEFINED ?
                        varvalue(gparam[XMAX].name) :
                        1e10;

        for (i = 0; i < MAXCUR; i++) {
                xsum = xxsum = ysum = yysum = xysum = 0.0;
                rsum = 0.0;
                npts = 0;
                mgcurve(i, calcpoint);
                gcvalue(i, types);
        }
}


gcheader(types)                 /* print header */
register char  *types;
{
        printf("__");
        while (*types)
                switch (*types++) {
                case 'n':
                        printf("|_Label___________");
                        break;
                case 'a':
                        printf("|____Mean______S.D._");
                        break;
                case 'm':
                        printf("|_____Min_______Max_");
                        break;
                case 'i':
                        printf("|___Integ_");
                        break;
                case 'l':
                        printf("|___Slope_____Intcp______Corr_");
                        break;
                default:
                        break;
                }
        printf("\n");
}


gcvalue(c, types)               /* print the values for the curve */
int  c;
register char  *types;
{
        double  d1, d2, d3, sqrt();

        if (npts < 1)
                return;

        printf("%c:", c+'A');
        while (*types)
                switch (*types++) {
                case 'n':
                        printf("  %-16s", cparam[c][CLABEL].flags & DEFINED ?
                                        cparam[c][CLABEL].v.s : "");
                        break;
                case 'a':
                        d1 = sqrt((yysum - ysum*ysum/npts)/(npts-1));
                        printf(" %9.2f %9.3f", ysum/npts, d1);
                        break;
                case 'm':
                        printf(" %9.2f %9.2f", ymin, ymax);
                        break;
                case 'i':
                        printf(" %9.2f", rsum);
                        break;
                case 'l':
                        d3 = xxsum - xsum*xsum/npts;
                        d1 = (xysum - xsum*ysum/npts)/d3;
                        d2 = (ysum - d1*xsum)/npts;
                        d3 = d1*sqrt(d3/(yysum - ysum*ysum/npts));
                        printf(" %9.5f %9.2f %9.5f", d1, d2, d3);
                        break;
                default:
                        break;
                }
        printf("\n");
}


calcpoint(c, x, y)              /* add a point to our calculation */
int  c;
double  x, y;
{
        if (x < xmin || x > xmax)
                return;

        xsum += x;
        xxsum += x*x;
        ysum += y;
        yysum += y*y;
        xysum += x*y;

        if (npts) {
                if (y < ymin)
                        ymin = y;
                else if (y > ymax)
                        ymax = y;
        } else
                ymin = ymax = y;

        if (npts)
                rsum += ( y + lasty )*( x - lastx )/2.0;
        lastx = x;
        lasty = y;

        npts++;
}
Revision:	1.1
Committed:	Sat Feb 22 02:07:26 2003 UTC (21 years, 2 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R5
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2			static const char RCSid[] = "$Id$";
3			#endif
4			/*
5			* gcalc.c - routines to do calculations on graph files.
6			*
7			* 7/7/86
8			* Greg Ward Larson
9			*/
10
11			#include <stdio.h>
12
13			#include "mgvars.h"
14
15			#ifndef BSD
16			#define index strchr
17			#endif
18
19			extern char *index();
20
21			static double xsum, xxsum, ysum, yysum, xysum;
22			static double xmin, xmax, ymin, ymax;
23			static double lastx, lasty, rsum;
24			static int npts;
25
26
27			gcalc(types) /* do a calculation */
28			char *types;
29			{
30			int i, calcpoint();
31
32			if (index(types, 'h') == NULL)
33			gcheader(types);
34
35			xmin = gparam[XMIN].flags & DEFINED ?
36			varvalue(gparam[XMIN].name) :
37			-1e10;
38			xmax = gparam[XMAX].flags & DEFINED ?
39			varvalue(gparam[XMAX].name) :
40			1e10;
41
42			for (i = 0; i < MAXCUR; i++) {
43			xsum = xxsum = ysum = yysum = xysum = 0.0;
44			rsum = 0.0;
45			npts = 0;
46			mgcurve(i, calcpoint);
47			gcvalue(i, types);
48			}
49			}
50
51
52			gcheader(types) /* print header */
53			register char *types;
54			{
55			printf("__");
56			while (*types)
57			switch (*types++) {
58			case 'n':
59			printf("\|_Label___________");
60			break;
61			case 'a':
62			printf("\|____Mean______S.D._");
63			break;
64			case 'm':
65			printf("\|_____Min_______Max_");
66			break;
67			case 'i':
68			printf("\|___Integ_");
69			break;
70			case 'l':
71			printf("\|___Slope_____Intcp______Corr_");
72			break;
73			default:
74			break;
75			}
76			printf("\n");
77			}
78
79
80			gcvalue(c, types) /* print the values for the curve */
81			int c;
82			register char *types;
83			{
84			double d1, d2, d3, sqrt();
85
86			if (npts < 1)
87			return;
88
89			printf("%c:", c+'A');
90			while (*types)
91			switch (*types++) {
92			case 'n':
93			printf(" %-16s", cparam[c][CLABEL].flags & DEFINED ?
94			cparam[c][CLABEL].v.s : "");
95			break;
96			case 'a':
97			d1 = sqrt((yysum - ysum*ysum/npts)/(npts-1));
98			printf(" %9.2f %9.3f", ysum/npts, d1);
99			break;
100			case 'm':
101			printf(" %9.2f %9.2f", ymin, ymax);
102			break;
103			case 'i':
104			printf(" %9.2f", rsum);
105			break;
106			case 'l':
107			d3 = xxsum - xsum*xsum/npts;
108			d1 = (xysum - xsum*ysum/npts)/d3;
109			d2 = (ysum - d1*xsum)/npts;
110			d3 = d1sqrt(d3/(yysum - ysumysum/npts));
111			printf(" %9.5f %9.2f %9.5f", d1, d2, d3);
112			break;
113			default:
114			break;
115			}
116			printf("\n");
117			}
118
119
120			calcpoint(c, x, y) /* add a point to our calculation */
121			int c;
122			double x, y;
123			{
124			if (x < xmin \|\| x > xmax)
125			return;
126
127			xsum += x;
128			xxsum += x*x;
129			ysum += y;
130			yysum += y*y;
131			xysum += x*y;
132
133			if (npts) {
134			if (y < ymin)
135			ymin = y;
136			else if (y > ymax)
137			ymax = y;
138			} else
139			ymin = ymax = y;
140
141			if (npts)
142			rsum += ( y + lasty )*( x - lastx )/2.0;
143			lastx = x;
144			lasty = y;
145
146			npts++;
147			}