ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/mgf2rad.c
(Generate patch)

Comparing ray/src/cv/mgf2rad.c (file contents):
Revision 2.1 by greg, Wed Jun 22 15:33:03 1994 UTC vs.
Revision 2.14 by greg, Thu Apr 13 14:43:33 1995 UTC

# Line 17 | Line 17 | static char SCCSid[] = "$SunId$ LBL";
17  
18   #define putv(v)         printf("%18.12g %18.12g %18.12g\n",(v)[0],(v)[1],(v)[2])
19  
20 < #define isgrey(cxy)     ((cxy)->cx > .31 && (cxy)->cx < .35 && \
21 <                        (cxy)->cy > .31 && (cxy)->cy < .35)
20 > #define invert          (xf_context != NULL && xf_context->rev)
21  
22 < #define is0vect(v)      ((v)[0] == 0. && (v)[1] == 0. && (v)[2] == 0.)
22 > double  glowdist = FHUGE;               /* glow test distance */
23  
24 < #define BIGFLT          1e8
24 > double  emult = 1.;                     /* emitter multiplier */
25  
26 < double  glowdist = 1.5*BIGFLT;          /* glow test distance */
26 > FILE    *matfp = stdout;                /* material output file */
27  
29 – double  emult = 1.;                     /* emmitter multiplier */
30 –
28   int     r_comment(), r_cone(), r_cyl(), r_face(), r_ies(), r_ring(), r_sph();
29   char    *material(), *object(), *addarg();
30  
# Line 41 | Line 38 | char   *argv[];
38          mg_ehand[MG_E_COMMENT] = r_comment;
39          mg_ehand[MG_E_COLOR] = c_hcolor;
40          mg_ehand[MG_E_CONE] = r_cone;
41 +        mg_ehand[MG_E_CMIX] = c_hcolor;
42 +        mg_ehand[MG_E_CSPEC] = c_hcolor;
43          mg_ehand[MG_E_CXY] = c_hcolor;
44 +        mg_ehand[MG_E_CCT] = c_hcolor;
45          mg_ehand[MG_E_CYL] = r_cyl;
46          mg_ehand[MG_E_ED] = c_hmaterial;
47          mg_ehand[MG_E_FACE] = r_face;
48          mg_ehand[MG_E_IES] = r_ies;
49 +        mg_ehand[MG_E_IR] = c_hmaterial;
50          mg_ehand[MG_E_MATERIAL] = c_hmaterial;
51          mg_ehand[MG_E_NORMAL] = c_hvertex;
52          mg_ehand[MG_E_OBJECT] = obj_handler;
# Line 53 | Line 54 | char   *argv[];
54          mg_ehand[MG_E_RD] = c_hmaterial;
55          mg_ehand[MG_E_RING] = r_ring;
56          mg_ehand[MG_E_RS] = c_hmaterial;
57 +        mg_ehand[MG_E_SIDES] = c_hmaterial;
58          mg_ehand[MG_E_SPH] = r_sph;
59          mg_ehand[MG_E_TD] = c_hmaterial;
60          mg_ehand[MG_E_TS] = c_hmaterial;
# Line 65 | Line 67 | char   *argv[];
67                  printf(" %s", argv[i]);
68                  switch (argv[i][1]) {
69                  case 'g':                       /* glow distance (meters) */
70 <                        if (argv[i][2] || badarg(argc-i, argv+i, "f"))
70 >                        if (argv[i][2] || badarg(argc-i-1, argv+i+1, "f"))
71                                  goto userr;
72                          glowdist = atof(argv[++i]);
73                          printf(" %s", argv[i]);
74                          break;
75                  case 'e':                       /* emitter multiplier */
76 <                        if (argv[i][2] || badarg(argc-i, argv+i, "f"))
76 >                        if (argv[i][2] || badarg(argc-i-1, argv+i+1, "f"))
77                                  goto userr;
78                          emult = atof(argv[++i]);
79                          printf(" %s", argv[i]);
80                          break;
81 +                case 'm':                       /* materials file */
82 +                        matfp = fopen(argv[++i], "a");
83 +                        if (matfp == NULL) {
84 +                                fprintf(stderr, "%s: cannot append\n", argv[i]);
85 +                                exit(1);
86 +                        }
87 +                        printf(" %s", argv[i]);
88 +                        break;
89                  default:
90                          goto userr;
91                  }
# Line 93 | Line 103 | char   *argv[];
103                  }
104          exit(0);
105   userr:
106 <        fprintf(stderr, "Usage: %s [-g dist][-m mult] [file.mgf] ..\n",
106 >        fprintf(stderr, "Usage: %s [-g dist][-e mult][-m matf] [file.mgf] ..\n",
107                          argv[0]);
108          exit(1);
109   }
# Line 104 | Line 114 | r_comment(ac, av)              /* repeat a comment verbatim */
114   register int    ac;
115   register char   **av;
116   {
117 <        fputs("\n#", stdout);   /* use Radiance comment character */
117 >        putchar('#');           /* use Radiance comment character */
118          while (--ac) {
119                  putchar(' ');
120                  fputs(*++av, stdout);
# Line 280 | Line 290 | char   **av;
290                  return(MG_EARGC);
291          if ((mat = material()) == NULL)
292                  return(MG_EBADMAT);
293 <        if (ac < 5) {                           /* check for surface normals */
293 >        if (ac <= 5) {                          /* check for surface normals */
294                  for (i = 1; i < ac; i++) {
295                          if ((cv = c_getvert(av[i])) == NULL)
296                                  return(MG_EUNDEF);
# Line 297 | Line 307 | char   **av;
307          printf("\n%s polygon %sf%d\n", mat, object(), ++nfaces);
308          printf("0\n0\n%d\n", 3*(ac-1));
309          for (i = 1; i < ac; i++) {
310 <                if ((cv = c_getvert(av[i])) == NULL)
310 >                if ((cv = c_getvert(av[invert ? ac-i : i])) == NULL)
311                          return(MG_EUNDEF);
312                  xf_xfmpoint(v, cv->p);
313                  putv(v);
# Line 344 | Line 354 | char   **av;
354                                          /* put out xform command */
355          printf("\n!xform");
356          oname = object();
357 <        if (*oname)
358 <                printf(" -n %s", oname);
357 >        if (*oname) {
358 >                printf(" -n ");
359 >                for (op = oname; op[1]; op++)   /* remove trailing separator */
360 >                        putchar(*op);
361 >        }
362          for (i = xa0; i < ac; i++)
363                  printf(" %s", av[i]);
364          if (ac > xa0 && xf_argc > 0)
# Line 361 | Line 374 | do_tri(mat, vn1, vn2, vn3)             /* put out smoothed triang
374   char    *mat, *vn1, *vn2, *vn3;
375   {
376          static int      ntris;
364 –        char    *mod = mat;
377          BARYCCM bvecs;
378          FLOAT   bcoor[3][3];
379          C_VERTEX        *cv1, *cv2, *cv3;
# Line 369 | Line 381 | char   *mat, *vn1, *vn2, *vn3;
381          FVECT   n1, n2, n3;
382          register int    i;
383                          /* the following is repeat code, so assume it's OK */
372 –        cv1 = c_getvert(vn1);
384          cv2 = c_getvert(vn2);
385 <        cv3 = c_getvert(vn3);
385 >        if (invert) {
386 >                cv3 = c_getvert(vn1);
387 >                cv1 = c_getvert(vn3);
388 >        } else {
389 >                cv1 = c_getvert(vn1);
390 >                cv3 = c_getvert(vn3);
391 >        }
392          xf_xfmpoint(v1, cv1->p);
393          xf_xfmpoint(v2, cv2->p);
394          xf_xfmpoint(v3, cv3->p);
395 <        if (comp_baryc(&bvecs, v1, v2, v3) == 0) {
396 <                printf("\n%s texfunc T-nor\n", mod);
397 <                mod = "T-nor";
398 <                printf("4 dx dy dz %s\n0\n", TCALNAME);
399 <                xf_rotvect(n1, cv1->n);
400 <                xf_rotvect(n2, cv2->n);
401 <                xf_rotvect(n3, cv3->n);
402 <                for (i = 0; i < 3; i++) {
403 <                        bcoor[i][0] = n1[i];
404 <                        bcoor[i][1] = n2[i];
405 <                        bcoor[i][2] = n3[i];
389 <                }
390 <                put_baryc(&bvecs, bcoor, 3);
395 >        if (comp_baryc(&bvecs, v1, v2, v3) < 0)
396 >                return;                         /* degenerate triangle! */
397 >        printf("\n%s texfunc T-nor\n", mat);
398 >        printf("4 dx dy dz %s\n0\n", TCALNAME);
399 >        xf_rotvect(n1, cv1->n);
400 >        xf_rotvect(n2, cv2->n);
401 >        xf_rotvect(n3, cv3->n);
402 >        for (i = 0; i < 3; i++) {
403 >                bcoor[i][0] = n1[i];
404 >                bcoor[i][1] = n2[i];
405 >                bcoor[i][2] = n3[i];
406          }
407 <        printf("\n%s polygon %st%d\n", mod, object(), ++ntris);
407 >        put_baryc(&bvecs, bcoor, 3);
408 >        printf("\nT-nor polygon %st%d\n", object(), ++ntris);
409          printf("0\n0\n9\n");
410          putv(v1);
411          putv(v2);
# Line 405 | Line 421 | material()                     /* get (and print) current material */
421          double  d;
422          register int    i;
423  
424 <        if (c_cmaterial->name != NULL)
425 <                mname = c_cmaterial->name;
424 >        if (c_cmname != NULL)
425 >                mname = c_cmname;
426          if (!c_cmaterial->clock)
427                  return(mname);          /* already current */
428                                  /* else update output */
429          c_cmaterial->clock = 0;
430          if (c_cmaterial->ed > .1) {     /* emitter */
431                  cvtcolor(radrgb, &c_cmaterial->ed_c,
432 <                                emult*c_cmaterial->ed/WHTEFFICACY);
433 <                if (glowdist < BIGFLT) {        /* do a glow */
434 <                        printf("\nvoid glow %s\n0\n0\n", mname);
435 <                        printf("4 %f %f %f %f\n", colval(radrgb,RED),
432 >                                emult*c_cmaterial->ed/(PI*WHTEFFICACY));
433 >                if (glowdist < FHUGE) {         /* do a glow */
434 >                        fprintf(matfp, "\nvoid glow %s\n0\n0\n", mname);
435 >                        fprintf(matfp, "4 %f %f %f %f\n", colval(radrgb,RED),
436                                          colval(radrgb,GRN),
437                                          colval(radrgb,BLU), glowdist);
438                  } else {
439 <                        printf("\nvoid light %s\n0\n0\n", mname);
440 <                        printf("3 %f %f %f\n", colval(radrgb,RED),
439 >                        fprintf(matfp, "\nvoid light %s\n0\n0\n", mname);
440 >                        fprintf(matfp, "3 %f %f %f\n", colval(radrgb,RED),
441                                          colval(radrgb,GRN),
442                                          colval(radrgb,BLU));
443                  }
# Line 429 | Line 445 | material()                     /* get (and print) current material */
445          }
446          d = c_cmaterial->rd + c_cmaterial->td +
447                          c_cmaterial->rs + c_cmaterial->ts;
448 <        if (d <= 0. | d >= 1.)
448 >        if (d < 0. | d > 1.)
449                  return(NULL);
450 <        if (c_cmaterial->td > .01 || c_cmaterial->ts > .01) {   /* trans */
450 >                                        /* check for glass/dielectric */
451 >        if (c_cmaterial->nr > 1.1 &&
452 >                        c_cmaterial->ts > .25 && c_cmaterial->rs <= .125 &&
453 >                        c_cmaterial->td <= .01 && c_cmaterial->rd <= .01 &&
454 >                        c_cmaterial->rs_a <= .01 && c_cmaterial->ts_a <= .01) {
455 >                cvtcolor(radrgb, &c_cmaterial->ts_c,
456 >                                c_cmaterial->ts + c_cmaterial->rs);
457 >                if (c_cmaterial->sided) {               /* dielectric */
458 >                        colval(radrgb,RED) = pow(colval(radrgb,RED),
459 >                                                        1./C_1SIDEDTHICK);
460 >                        colval(radrgb,GRN) = pow(colval(radrgb,GRN),
461 >                                                        1./C_1SIDEDTHICK);
462 >                        colval(radrgb,BLU) = pow(colval(radrgb,BLU),
463 >                                                        1./C_1SIDEDTHICK);
464 >                        fprintf(matfp, "\nvoid dielectric %s\n0\n0\n", mname);
465 >                        fprintf(matfp, "5 %g %g %g %f 0\n", colval(radrgb,RED),
466 >                                        colval(radrgb,GRN), colval(radrgb,BLU),
467 >                                        c_cmaterial->nr);
468 >                        return(mname);
469 >                }
470 >                                                        /* glass */
471 >                fprintf(matfp, "\nvoid glass %s\n0\n0\n", mname);
472 >                fprintf(matfp, "4 %f %f %f %f\n", colval(radrgb,RED),
473 >                                colval(radrgb,GRN), colval(radrgb,BLU),
474 >                                c_cmaterial->nr);
475 >                return(mname);
476 >                }
477 >                                        /* check for trans */
478 >        if (c_cmaterial->td > .01 || c_cmaterial->ts > .01) {
479                  double  ts, a5, a6;
480  
481 <                ts = sqrt(c_cmaterial->ts);     /* because we use 2 sides */
481 >                if (c_cmaterial->sided) {
482 >                        ts = sqrt(c_cmaterial->ts);     /* approximate */
483 >                        a5 = .5;
484 >                } else {
485 >                        ts = c_cmaterial->ts;
486 >                        a5 = 1.;
487 >                }
488                                                  /* average colors */
489                  d = c_cmaterial->rd + c_cmaterial->td + ts;
490                  cvtcolor(radrgb, &c_cmaterial->rd_c, c_cmaterial->rd/d);
# Line 444 | Line 494 | material()                     /* get (and print) current material */
494                  addcolor(radrgb, c2);
495                  if (c_cmaterial->rs + ts > .0001)
496                          a5 = (c_cmaterial->rs*c_cmaterial->rs_a +
497 <                                        ts*.5*c_cmaterial->ts_a) /
497 >                                        ts*a5*c_cmaterial->ts_a) /
498                                          (c_cmaterial->rs + ts);
499                  a6 = (c_cmaterial->td + ts) /
500                                  (c_cmaterial->rd + c_cmaterial->td + ts);
501 <                if (a6 < .999) {
501 >                if (a6 < .999)
502                          d = c_cmaterial->rd/(1. - c_cmaterial->rs)/(1. - a6);
503 <                        scalecolor(radrgb, d);
504 <                }
505 <                printf("\nvoid trans %s\n0\n0\n", mname);
506 <                printf("7 %f %f %f\n", colval(radrgb,RED),
503 >                else
504 >                        d = c_cmaterial->td + ts;
505 >                scalecolor(radrgb, d);
506 >                fprintf(matfp, "\nvoid trans %s\n0\n0\n", mname);
507 >                fprintf(matfp, "7 %f %f %f\n", colval(radrgb,RED),
508                                  colval(radrgb,GRN), colval(radrgb,BLU));
509 <                printf("\t%f %f %f %f\n", c_cmaterial->rs, a5, a6,
509 >                fprintf(matfp, "\t%f %f %f %f\n", c_cmaterial->rs, a5, a6,
510                                  ts/(ts + c_cmaterial->td));
511                  return(mname);
512          }
513 <        if (c_cmaterial->rs < .01 || isgrey(&c_cmaterial->rs_c)) { /* plastic */
514 <                if (c_cmaterial->rs > .999)
515 <                        cvtcolor(radrgb, &c_cmaterial->rd_c, 1.);
516 <                else
466 <                        cvtcolor(radrgb, &c_cmaterial->rd_c,
513 >                                        /* check for plastic */
514 >        if (c_cmaterial->rs < .1 && (c_cmaterial->rs < .01 ||
515 >                                        c_isgrey(&c_cmaterial->rs_c))) {
516 >                cvtcolor(radrgb, &c_cmaterial->rd_c,
517                                          c_cmaterial->rd/(1.-c_cmaterial->rs));
518 <                printf("\nvoid plastic %s\n0\n0\n", mname);
519 <                printf("5 %f %f %f %f %f\n", colval(radrgb,RED),
518 >                fprintf(matfp, "\nvoid plastic %s\n0\n0\n", mname);
519 >                fprintf(matfp, "5 %f %f %f %f %f\n", colval(radrgb,RED),
520                                  colval(radrgb,GRN), colval(radrgb,BLU),
521                                  c_cmaterial->rs, c_cmaterial->rs_a);
522                  return(mname);
523          }
524                                          /* else it's metal */
525 <        d = c_cmaterial->rd + c_cmaterial->rs;  /* average colors */
526 <        cvtcolor(radrgb, &c_cmaterial->rd_c, c_cmaterial->rd/d);
527 <        cvtcolor(c2, &c_cmaterial->rs_c, c_cmaterial->rs/d);
525 >                                                /* average colors */
526 >        cvtcolor(radrgb, &c_cmaterial->rd_c, c_cmaterial->rd);
527 >        cvtcolor(c2, &c_cmaterial->rs_c, c_cmaterial->rs);
528          addcolor(radrgb, c2);
529 <        if (c_cmaterial->rs < .999) {
530 <                d = c_cmaterial->rd/(1. - c_cmaterial->rs);
481 <                scalecolor(radrgb, d);
482 <        }
483 <        printf("\nvoid metal %s\n0\n0\n", mname);
484 <        printf("5 %f %f %f %f %f\n", colval(radrgb,RED),
529 >        fprintf(matfp, "\nvoid metal %s\n0\n0\n", mname);
530 >        fprintf(matfp, "5 %f %f %f %f %f\n", colval(radrgb,RED),
531                          colval(radrgb,GRN), colval(radrgb,BLU),
532 <                        c_cmaterial->rs, c_cmaterial->rs_a);
532 >                        c_cmaterial->rs/(c_cmaterial->rd + c_cmaterial->rs),
533 >                        c_cmaterial->rs_a);
534          return(mname);
535   }
536  
# Line 495 | Line 542 | double intensity;
542   {
543          static COLOR    ciexyz;
544  
545 +        c_ccvt(ciec, C_CSXY);           /* get xy representation */
546          ciexyz[1] = intensity;
547          ciexyz[0] = ciec->cx/ciec->cy*ciexyz[1];
548          ciexyz[2] = ciexyz[1]*(1./ciec->cy - 1.) - ciexyz[0];

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines