src/gen/genmarble.c

#ifndef lint
static const char       RCSid[] = "$Id$";
#endif
/*
 *  genmarble.c - generate a marble with bubbles inside.
 *
 *     1/8/86
 */

#include  <stdio.h>

#include <stdlib.h>

#include  <math.h>

#include  "random.h"

#define  PI     3.14159265359

typedef double  FVECT[3];

double  bubble();       /* pretty cute, huh? */


main(argc, argv)
int  argc;
char  **argv;
{
        char  *cmtype, *cname;
        FVECT  cent;
        double  rad;
        int  nbubbles, i;
        double  bubrad;
        FVECT  v;
        double  brad;
        
        if (argc != 9) {
                fprintf(stderr,
                        "Usage: %s material name cent rad #bubbles bubrad\n",
                                         argv[0]);
                exit(1);
        }
        cmtype = argv[1];
        cname = argv[2];
        cent[0] = atof(argv[3]);
        cent[1] = atof(argv[4]);
        cent[2] = atof(argv[5]);
        rad = atof(argv[6]);
        nbubbles = atoi(argv[7]);
        bubrad = atof(argv[8]);

        if (bubrad >= rad) {
                fprintf(stderr, "%s: bubbles too big for marble\n", argv[0]);
                exit(1);
        }

        printf("\n%s sphere %s\n", cmtype, cname);
        printf("0\n0\n4 %f %f %f %f\n", cent[0], cent[1], cent[2], rad);

        for (i = 0; i < nbubbles; i++) {
                brad = bubble(v, cent, rad, bubrad);
                printf("\n%s bubble %s.%d\n", cmtype, cname, i);
                printf("0\n0\n4 %f %f %f %f\n", v[0], v[1], v[2], brad);
        }
        
        return(0);
}


double
bubble(v, cent, rad, bubrad)    /* compute location of random bubble */
FVECT  v, cent;
double  rad, bubrad;
{
        double  r, ro, theta, phi;

        r = frandom()*bubrad;
        ro = sqrt(frandom())*(rad-r);
        theta = frandom()*(2.0*PI);
        phi = frandom()*PI;
        sphere_cart(v, ro, theta, phi);
        v[0] += cent[0]; v[1] += cent[1]; v[2] += cent[2];
        return(r);
}


sphere_cart(v, ro, theta, phi)  /* spherical to cartesian coord. conversion */
FVECT  v;
double  ro, theta, phi;
{
        double  d;
        
        d = sin(phi);
        v[0] = ro*d*cos(theta);
        v[1] = ro*d*sin(theta);
        v[2] = ro*cos(phi);
}
Revision:	2.5
Committed:	Sat Feb 22 02:07:23 2003 UTC (21 years, 2 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R5
Changes since 2.4:	+3 -8 lines
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	greg	2.5	static const char RCSid[] = "$Id$";
3	greg	1.1	#endif
4			/*
5			* genmarble.c - generate a marble with bubbles inside.
6			*
7			* 1/8/86
8			*/
9
10			#include <stdio.h>
11
12	greg	2.5	#include <stdlib.h>
13
14	greg	2.3	#include <math.h>
15
16	greg	1.1	#include "random.h"
17
18			#define PI 3.14159265359
19
20			typedef double FVECT[3];
21
22			double bubble(); /* pretty cute, huh? */
23	greg	2.4
24	greg	1.1
25			main(argc, argv)
26			int argc;
27			char **argv;
28			{
29			char cmtype, cname;
30			FVECT cent;
31			double rad;
32			int nbubbles, i;
33			double bubrad;
34			FVECT v;
35			double brad;
36
37			if (argc != 9) {
38			fprintf(stderr,
39			"Usage: %s material name cent rad #bubbles bubrad\n",
40			argv[0]);
41			exit(1);
42			}
43			cmtype = argv[1];
44			cname = argv[2];
45			cent[0] = atof(argv[3]);
46			cent[1] = atof(argv[4]);
47			cent[2] = atof(argv[5]);
48			rad = atof(argv[6]);
49			nbubbles = atoi(argv[7]);
50			bubrad = atof(argv[8]);
51
52			if (bubrad >= rad) {
53			fprintf(stderr, "%s: bubbles too big for marble\n", argv[0]);
54			exit(1);
55			}
56
57			printf("\n%s sphere %s\n", cmtype, cname);
58			printf("0\n0\n4 %f %f %f %f\n", cent[0], cent[1], cent[2], rad);
59
60			for (i = 0; i < nbubbles; i++) {
61			brad = bubble(v, cent, rad, bubrad);
62			printf("\n%s bubble %s.%d\n", cmtype, cname, i);
63			printf("0\n0\n4 %f %f %f %f\n", v[0], v[1], v[2], brad);
64			}
65
66			return(0);
67			}
68
69
70			double
71			bubble(v, cent, rad, bubrad) /* compute location of random bubble */
72			FVECT v, cent;
73			double rad, bubrad;
74			{
75			double r, ro, theta, phi;
76
77			r = frandom()*bubrad;
78			ro = sqrt(frandom())*(rad-r);
79			theta = frandom()(2.0PI);
80			phi = frandom()*PI;
81			sphere_cart(v, ro, theta, phi);
82			v[0] += cent[0]; v[1] += cent[1]; v[2] += cent[2];
83			return(r);
84			}
85
86
87			sphere_cart(v, ro, theta, phi) /* spherical to cartesian coord. conversion */
88			FVECT v;
89			double ro, theta, phi;
90			{
91			double d;
92
93			d = sin(phi);
94			v[0] = rodcos(theta);
95			v[1] = rodsin(theta);
96			v[2] = ro*cos(phi);
97			}