ViewVC Help

View File | Revision Log | Show Annotations | Download File | Root Listing

root/Development/ray/src/gen/genrev.c

Comparing ray/src/gen/genrev.c (file contents):
Revision 2.1 by greg, Tue Nov 12 17:05:00 1991 UTC vs.
Revision 2.16 by greg, Fri Apr 18 23:02:53 2025 UTC

#	Line 1 | Line 1
1		#ifndef lint
2	<	static char SCCSid[] = "$SunId$ LBL";
2	>	static const char       RCSid[] = "$Id$";
3		#endif
4	–
5	–	/* Copyright (c) 1989 Regents of the University of California */
6	–
4		/*
5		*  genrev.c - program to generate functions of rotation about z
6		*
#	Line 14 | Line 11 | static char SCCSid[] = "$SunId$ LBL";
11		*      8/6/86
12		*/
13
14	<	#include  <stdio.h>
14	>	#include  <stdlib.h>
15	>	#include  <math.h>
16
17	+	#include  "rtio.h"
18	+	#include  "rterror.h"
19	+	#include  "resolu.h"
20	+	#include  "calcomp.h"
21	+
22		#define  ZNAME          "Z`SYS`"                /* z function name */
23		#define  RNAME          "R`SYS`"                /* r function name */
24
#	Line 28 | Line 31 | static char SCCSid[] = "$SunId$ LBL";
31		#define  UP             04
32		#define  DOWN           010
33
31	–	double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
34
35	+	void
36	+	computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
37	+	double  *nzp, *nrp, z0, r0, z1, r1;
38	+	{
39	+	double  dr, dz, len;
40
41	+	dz = r0 - r1;                           /* right angle vector */
42	+	dr = z1 - z0;
43	+	len = sqrt(dr*dr + dz*dz);
44	+	*nzp = dz/len;
45	+	*nrp = dr/len;
46	+	}
47	+
48	+
49	+	double
50	+	l_hermite(char *nm)
51	+	{
52	+	double  t;
53	+
54	+	t = argument(5);
55	+	return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
56	+	argument(2)*(-2.0*t+3.0)*t*t +
57	+	argument(3)*((t-2.0)*t+1.0)*t +
58	+	argument(4)*(t-1.0)*t*t );
59	+	}
60	+
61	+
62	+	double
63	+	l_bezier(char *nm)
64	+	{
65	+	double  t;
66	+
67	+	t = argument(5);
68	+	return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
69	+	argument(2) * 3.*t*(1.+t*(-2.+t)) +
70	+	argument(3) * 3.*t*t*(1.-t) +
71	+	argument(4) * t*t*t );
72	+	}
73	+
74	+
75	+	double
76	+	l_bspline(char *nm)
77	+	{
78	+	double  t;
79	+
80	+	t = argument(5);
81	+	return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
82	+	argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
83	+	argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
84	+	argument(4) * (1./6.*t*t*t) );
85	+	}
86	+
87	+
88	+	int
89		main(argc, argv)
90		int  argc;
91		char  *argv[];
92		{
38	–	extern long     eclock;
93		char  stmp[256];
94		char  *modname;
95		int  smooth = 0;
#	Line 44 | Line 98 | char  *argv[];
98		int  i, nseg;
99		int  orient;
100
101	+	esupport |= E_VARIABLE|E_FUNCTION|E_RCONST;
102	+	esupport &= ~(E_OUTCHAN|E_INCHAN);
103		varset("PI", ':', PI);
104		funset("hermite", 5, ':', l_hermite);
105		funset("bezier", 5, ':', l_bezier);
#	Line 55 | Line 111 | char  *argv[];
111		for (i = 6; i < argc; i++)
112		if (!strcmp(argv[i], "-e"))
113		scompile(argv[++i], NULL, 0);
114	<	else if (!strcmp(argv[i], "-f"))
115	<	fcompile(argv[++i]);
116	<	else if (!strcmp(argv[i], "-s"))
114	>	else if (!strcmp(argv[i], "-f")) {
115	>	char  *fpath = getpath(argv[++i], getrlibpath(), 0);
116	>	if (fpath == NULL) {
117	>	fprintf(stderr, "%s: cannot find file '%s'\n",
118	>	argv[0], argv[i]);
119	>	quit(1);
120	>	}
121	>	fcompile(fpath);
122	>	} else if (!strcmp(argv[i], "-s"))
123		smooth = 1;
124		else
125		goto userror;
#	Line 66 | Line 128 | char  *argv[];
128		scompile(stmp, NULL, 0);
129		sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
130		scompile(stmp, NULL, 0);
131	<	nseg = atoi(argv[5]);
131	>	nseg = eval(argv[5]) + .5;
132		if (nseg <= 0)
133		goto userror;
134		modname = smooth ? "Phong" : argv[1];
135
136	<	printhead(argc, argv);
137	<	eclock = 0;
136	>	fputs("# ", stdout);
137	>	printargs(argc, argv, stdout);
138	>	doptimize(1);
139	>	eclock++;
140
141		lastnz = lastnr = 0.0;
142		t = 0.0;
#	Line 107 | Line 171 | char  *argv[];
171		printf("0\n4\n");
172		if (orient&(UP|DOWN)) {
173		t = (nextnz - lastnz)/(z - lastz);
174	<	printf("\t%15.9g\t%15.9g\n",
174	>	printf("\t%18.12g\t%18.12g\n",
175		t, lastnz - t*lastz);
176		} else
177		printf("\t0\t%d\n", orient&IN ? 1 : -1);
178		if (orient&(OUT|IN))  {
179		t = (nextnr - lastnr)/(r - lastr);
180	<	printf("\t%15.9g\t%15.9g\n",
180	>	printf("\t%18.12g\t%18.12g\n",
181		t, lastnr - t*lastr);
182		} else
183		printf("\t0\t%d\n", orient&UP ? 1 : -1);
#	Line 123 | Line 187 | char  *argv[];
187		orient&DOWN ? "tube" : "cylinder",
188		argv[2], i);
189		printf("0\n0\n7\n");
190	<	printf("\t0\t0\t%15.9g\n", lastz);
191	<	printf("\t0\t0\t%15.9g\n", z);
192	<	printf("\t%15.9g\n", r);
190	>	printf("\t0\t0\t%18.12g\n", lastz);
191	>	printf("\t0\t0\t%18.12g\n", z);
192	>	printf("\t%18.12g\n", r);
193		} else if (!(orient&(UP|DOWN))) {
194		printf("\n%s ring %s.%d\n", modname, argv[2], i);
195		printf("0\n0\n8\n");
196	<	printf("\t0\t0\t%15.9g\n", z);
197	<	printf("\t0\t0\t%15.9g\n", orient&IN ? 1.0 : -1.0);
198	<	printf("\t%15.9g\t%15.9g\n", lastr, r);
196	>	printf("\t0\t0\t%18.12g\n", z);
197	>	printf("\t0\t0\t%18.12g\n", orient&IN ? 1.0 : -1.0);
198	>	printf("\t%18.12g\t%18.12g\n", lastr, r);
199		} else {
200		printf("\n%s %s %s.%d\n", modname,
201		orient&DOWN ? "cup" : "cone",
202		argv[2], i);
203		printf("0\n0\n8\n");
204	<	printf("\t0\t0\t%15.9g\n", lastz);
205	<	printf("\t0\t0\t%15.9g\n", z);
206	<	printf("\t%15.9g\t%15.9g\n", lastr, r);
204	>	printf("\t0\t0\t%18.12g\n", lastz);
205	>	printf("\t0\t0\t%18.12g\n", z);
206	>	printf("\t%18.12g\t%18.12g\n", lastr, r);
207		}
208		endfor:
209		lastz = z; lastr = r;
#	Line 147 | Line 211 | char  *argv[];
211		lastnz = nz; lastnr = nr;
212		nz = nextnz; nr = nextnr;
213		}
214	<	quit(0);
214	>	return 0;
215
216		userror:
217		fprintf(stderr,
218		"Usage: %s material name z(t) r(t) nseg [-e expr] [-f file] [-s]\n",
219		argv[0]);
220	<	quit(1);
220	>	return 1;
221		}
222
223
160	–	computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
161	–	double  *nzp, *nrp, z0, r0, z1, r1;
162	–	{
163	–	extern double  sqrt();
164	–	double  dr, dz, len;
165	–
166	–	dz = r0 - r1;                           /* right angle vector */
167	–	dr = z1 - z0;
168	–	len = sqrt(dr*dr + dz*dz);
169	–	*nzp = dz/len;
170	–	*nrp = dr/len;
171	–	}
172	–
173	–
174	–	eputs(msg)
175	–	char  *msg;
176	–	{
177	–	fputs(msg, stderr);
178	–	}
179	–
180	–
181	–	wputs(msg)
182	–	char  *msg;
183	–	{
184	–	eputs(msg);
185	–	}
186	–
187	–
188	–	quit(code)
189	–	{
190	–	exit(code);
191	–	}
192	–
193	–
194	–	printhead(ac, av)               /* print command header */
195	–	register int  ac;
196	–	register char  **av;
197	–	{
198	–	putchar('#');
199	–	while (ac--) {
200	–	putchar(' ');
201	–	fputs(*av++, stdout);
202	–	}
203	–	putchar('\n');
204	–	}
205	–
206	–
207	–	double
208	–	l_hermite()
209	–	{
210	–	double  t;
211	–
212	–	t = argument(5);
213	–	return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
214	–	argument(2)*(-2.0*t+3.0)*t*t +
215	–	argument(3)*((t-2.0)*t+1.0)*t +
216	–	argument(4)*(t-1.0)*t*t );
217	–	}
218	–
219	–
220	–	double
221	–	l_bezier()
222	–	{
223	–	double  t;
224	–
225	–	t = argument(5);
226	–	return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
227	–	argument(2) * 3.*t*(1.+t*(-2.+t)) +
228	–	argument(3) * 3.*t*t*(1.-t) +
229	–	argument(4) * t*t*t );
230	–	}
231	–
232	–
233	–	double
234	–	l_bspline()
235	–	{
236	–	double  t;
237	–
238	–	t = argument(5);
239	–	return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
240	–	argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
241	–	argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
242	–	argument(4) * (1./6.*t*t*t) );
243	–	}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines (old)
>	Changed lines (new)