ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/genbeads.c
Revision: 2.1
Committed: Tue Nov 12 17:04:35 1991 UTC (32 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 1.1: +0 -0 lines
Log Message:
updated revision number for release 2.0

File Contents

# User Rev Content
1 greg 1.1 /* Copyright (c) 1986 Regents of the University of California */
2    
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8     * genbeads.c - generate a string of spheres using Hermite
9     * curve specification.
10     *
11     * 10/29/85
12     */
13    
14     #include <stdio.h>
15    
16    
17     char *mtype; /* material type */
18    
19     char *name; /* name */
20    
21    
22     main(argc, argv)
23     int argc;
24     char **argv;
25     {
26     double atof();
27     double p0[3], p1[3], r0[3], r1[3];
28     double rad, inc;
29    
30     if (argc != 17) {
31     fprintf(stderr, "Usage: %s material name p0 p1 r0 r1 rad inc\n",
32     argv[0]);
33     exit(1);
34     }
35     mtype = argv[1];
36     name = argv[2];
37     p0[0] = atof(argv[3]);
38     p0[1] = atof(argv[4]);
39     p0[2] = atof(argv[5]);
40     p1[0] = atof(argv[6]);
41     p1[1] = atof(argv[7]);
42     p1[2] = atof(argv[8]);
43     r0[0] = atof(argv[9]);
44     r0[1] = atof(argv[10]);
45     r0[2] = atof(argv[11]);
46     r1[0] = atof(argv[12]);
47     r1[1] = atof(argv[13]);
48     r1[2] = atof(argv[14]);
49     rad = atof(argv[15]);
50     inc = atof(argv[16]);
51    
52     genstring(mtype, name, p0, p1, r0, r1, rad, inc);
53    
54     return(0);
55     }
56    
57    
58     genstring(mtype, name, p0, p1, r0, r1, rad, inc)
59     char *mtype;
60     char *name;
61     double p0[3];
62     double p1[3];
63     double r0[3];
64     double r1[3];
65     double rad;
66     double inc;
67     {
68     double sqrt();
69     register int i;
70     double v[3];
71     double t;
72    
73     t = 0.0;
74     for (i = 0; t <= 1.0; i++) {
75    
76     printf("\n%s sphere %s.%d\n", mtype, name, i);
77    
78     hermite3(v, p0, p1, r0, r1, t);
79    
80     printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
81     v[0], v[1], v[2], rad);
82    
83     htan3(v, p0, p1, r0, r1, t);
84     t += inc / sqrt(v[0]*v[0] + v[1]*v[1] + v[2]*v[2]);
85     }
86     }