src/meta/gcalc.c

#ifndef lint
static const char       RCSid[] = "$Id: gcalc.c,v 1.3 2003/11/15 02:13:37 schorsch Exp $";
#endif
/*
 *  gcalc.c - routines to do calculations on graph files.
 *
 *     7/7/86
 *     Greg Ward Larson
 */

#include  <stdio.h>
#include  <string.h>

#include  "mgvars.h"


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

static void gcheader(char  *types);
static void calcpoint(int  c, double  x, double  y);
static void gcvalue(int  c, char  *types);

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

        if (strchr(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);
        }
}


void
gcheader(                       /* 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");
}


void
gcvalue(                /* 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(" %10.6f %11.4f %7.3f", d1, d2, d3);
                        break;
                default:
                        break;
                }
        printf("\n");
}


void
calcpoint(              /* add a point to our calculation */
        int  c,
        double  x,
        double  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.4
Committed:	Thu Dec 1 18:47:35 2011 UTC (13 years, 10 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad5R0, rad5R1, rad4R2, rad4R2P1, rad4R2P2
Changes since 1.3:	+3 -3 lines
Log Message:	Increased significant digits for intercept
#	Content
1	#ifndef lint
2	static const char RCSid[] = "$Id: gcalc.c,v 1.3 2003/11/15 02:13:37 schorsch Exp $";
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	#include <string.h>
13
14	#include "mgvars.h"
15
16
17	static double xsum, xxsum, ysum, yysum, xysum;
18	static double xmin, xmax, ymin, ymax;
19	static double lastx, lasty, rsum;
20	static int npts;
21
22	static void gcheader(char *types);
23	static void calcpoint(int c, double x, double y);
24	static void gcvalue(int c, char *types);
25
26	void
27	gcalc( /* do a calculation */
28	char *types
29	)
30	{
31	int i;
32
33	if (strchr(types, 'h') == NULL)
34	gcheader(types);
35
36	xmin = gparam[XMIN].flags & DEFINED ?
37	varvalue(gparam[XMIN].name) :
38	-1e10;
39	xmax = gparam[XMAX].flags & DEFINED ?
40	varvalue(gparam[XMAX].name) :
41	1e10;
42
43	for (i = 0; i < MAXCUR; i++) {
44	xsum = xxsum = ysum = yysum = xysum = 0.0;
45	rsum = 0.0;
46	npts = 0;
47	mgcurve(i, calcpoint);
48	gcvalue(i, types);
49	}
50	}
51
52
53	void
54	gcheader( /* print header */
55	register char *types
56	)
57	{
58	printf("__");
59	while (*types)
60	switch (*types++) {
61	case 'n':
62	printf("\|_Label___________");
63	break;
64	case 'a':
65	printf("\|____Mean______S.D._");
66	break;
67	case 'm':
68	printf("\|_____Min_______Max_");
69	break;
70	case 'i':
71	printf("\|___Integ_");
72	break;
73	case 'l':
74	printf("\|____Slope_______Intcp____Corr_");
75	break;
76	default:
77	break;
78	}
79	printf("\n");
80	}
81
82
83	void
84	gcvalue( /* print the values for the curve */
85	int c,
86	register char *types
87	)
88	{
89	double d1, d2, d3, sqrt();
90
91	if (npts < 1)
92	return;
93
94	printf("%c:", c+'A');
95	while (*types)
96	switch (*types++) {
97	case 'n':
98	printf(" %-16s", cparam[c][CLABEL].flags & DEFINED ?
99	cparam[c][CLABEL].v.s : "");
100	break;
101	case 'a':
102	d1 = sqrt((yysum - ysum*ysum/npts)/(npts-1));
103	printf(" %9.2f %9.3f", ysum/npts, d1);
104	break;
105	case 'm':
106	printf(" %9.2f %9.2f", ymin, ymax);
107	break;
108	case 'i':
109	printf(" %9.2f", rsum);
110	break;
111	case 'l':
112	d3 = xxsum - xsum*xsum/npts;
113	d1 = (xysum - xsum*ysum/npts)/d3;
114	d2 = (ysum - d1*xsum)/npts;
115	d3 = d1sqrt(d3/(yysum - ysumysum/npts));
116	printf(" %10.6f %11.4f %7.3f", d1, d2, d3);
117	break;
118	default:
119	break;
120	}
121	printf("\n");
122	}
123
124
125	void
126	calcpoint( /* add a point to our calculation */
127	int c,
128	double x,
129	double y
130	)
131	{
132	if (x < xmin \|\| x > xmax)
133	return;
134
135	xsum += x;
136	xxsum += x*x;
137	ysum += y;
138	yysum += y*y;
139	xysum += x*y;
140
141	if (npts) {
142	if (y < ymin)
143	ymin = y;
144	else if (y > ymax)
145	ymax = y;
146	} else
147	ymin = ymax = y;
148
149	if (npts)
150	rsum += ( y + lasty )*( x - lastx )/2.0;
151	lastx = x;
152	lasty = y;
153
154	npts++;
155	}