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, 7 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
#	Content
1	#ifndef lint
2	static const char RCSid[] = "$Id$";
3	#endif
4	/*
5	* genmarble.c - generate a marble with bubbles inside.
6	*
7	* 1/8/86
8	*/
9
10	#include <stdio.h>
11
12	#include <stdlib.h>
13
14	#include <math.h>
15
16	#include "random.h"
17
18	#define PI 3.14159265359
19
20	typedef double FVECT[3];
21
22	double bubble(); /* pretty cute, huh? */
23
24
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	}