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

Comparing ray/src/cv/mgflib/parser.c (file contents):
Revision 1.1 by greg, Tue Jun 21 14:45:46 1994 UTC vs.
Revision 1.7 by greg, Fri Jun 24 17:15:49 1994 UTC

# Line 56 | Line 56 | int    mg_nqcdivs = MG_NQCD;   /* number of divisions per q
56                                  /* alternate handler routines */
57  
58   static int      e_any_toss(),           /* discard unneeded entity */
59 <                e_ies();                /* IES luminaire file */
59 >                e_ies(),                /* IES luminaire file */
60                  e_include(),            /* include file */
61                  e_sph(),                /* sphere */
62 +                e_cmix(),               /* color mixtures */
63 +                e_cspec(),              /* color spectra */
64                  e_cyl(),                /* cylinder */
65                  e_cone(),               /* cone */
66 +                e_prism(),              /* prism */
67                  e_ring(),               /* ring */
68                  e_torus();              /* torus */
69  
# Line 103 | Line 106 | mg_init()                      /* initialize alternate entity handlers */
106                  ineed |= 1<<MG_E_POINT|1<<MG_E_NORMAL|1<<MG_E_VERTEX;
107          } else
108                  uneed |= 1<<MG_E_POINT|1<<MG_E_NORMAL|1<<MG_E_VERTEX|1<<MG_E_XF;
109 +        if (mg_ehand[MG_E_PRISM] == NULL) {
110 +                mg_ehand[MG_E_PRISM] = e_prism;
111 +                ineed |= 1<<MG_E_POINT|1<<MG_E_VERTEX;
112 +        } else
113 +                uneed |= 1<<MG_E_POINT|1<<MG_E_VERTEX|1<<MG_E_XF;
114          if (mg_ehand[MG_E_TORUS] == NULL) {
115                  mg_ehand[MG_E_TORUS] = e_torus;
116                  ineed |= 1<<MG_E_POINT|1<<MG_E_NORMAL|1<<MG_E_VERTEX;
117          } else
118                  uneed |= 1<<MG_E_POINT|1<<MG_E_NORMAL|1<<MG_E_VERTEX|1<<MG_E_XF;
119 +        if (mg_ehand[MG_E_COLOR] != NULL) {
120 +                if (mg_ehand[MG_E_CMIX] == NULL) {
121 +                        mg_ehand[MG_E_CMIX] = e_cmix;
122 +                        ineed |= 1<<MG_E_COLOR|1<<MG_E_CXY|1<<MG_E_CSPEC|1<<MG_E_CMIX;
123 +                }
124 +                if (mg_ehand[MG_E_CSPEC] == NULL) {
125 +                        mg_ehand[MG_E_CSPEC] = e_cspec;
126 +                        ineed |= 1<<MG_E_COLOR|1<<MG_E_CXY|1<<MG_E_CSPEC|1<<MG_E_CMIX;
127 +                }
128 +        }
129                                          /* check for consistency */
130          if (mg_ehand[MG_E_FACE] != NULL)
131                  uneed |= 1<<MG_E_POINT|1<<MG_E_VERTEX|1<<MG_E_XF;
132 <        if (mg_ehand[MG_E_CXY] != NULL)
132 >        if (mg_ehand[MG_E_CXY] != NULL || mg_ehand[MG_E_CSPEC] != NULL ||
133 >                        mg_ehand[MG_E_CMIX] != NULL)
134                  uneed |= 1<<MG_E_COLOR;
135          if (mg_ehand[MG_E_RD] != NULL || mg_ehand[MG_E_TD] != NULL ||
136                          mg_ehand[MG_E_ED] != NULL ||
# Line 131 | Line 150 | mg_init()                      /* initialize alternate entity handlers */
150                  e_supp[MG_E_POINT] = c_hvertex;
151          if (ineed & 1<<MG_E_NORMAL && mg_ehand[MG_E_NORMAL] != c_hvertex)
152                  e_supp[MG_E_NORMAL] = c_hvertex;
153 +        if (ineed & 1<<MG_E_COLOR && mg_ehand[MG_E_COLOR] != c_hcolor)
154 +                e_supp[MG_E_COLOR] = c_hcolor;
155 +        if (ineed & 1<<MG_E_CXY && mg_ehand[MG_E_CXY] != c_hcolor)
156 +                e_supp[MG_E_CXY] = c_hcolor;
157 +        if (ineed & 1<<MG_E_CSPEC && mg_ehand[MG_E_CSPEC] != c_hcolor)
158 +                e_supp[MG_E_CSPEC] = c_hcolor;
159 +        if (ineed & 1<<MG_E_CMIX && mg_ehand[MG_E_CMIX] != c_hcolor)
160 +                e_supp[MG_E_CMIX] = c_hcolor;
161                                          /* discard remaining entities */
162          for (i = 0; i < MG_NENTITIES; i++)
163                  if (mg_ehand[i] == NULL)
# Line 300 | Line 327 | char   *fn;
327          int     rval;
328  
329          if ((rval = mg_open(&cntxt, fn)) != MG_OK) {
330 <                fprintf("%s: %s\n", fn, mg_err[rval]);
330 >                fprintf(stderr, "%s: %s\n", fn, mg_err[rval]);
331                  return(rval);
332          }
333          while (mg_read())               /* parse each line */
# Line 349 | Line 376 | int    (*f)();
376                  if (i+1 >= ac || !isint(av[i+1]))
377                          return(MG_ETYPE);
378                  niter = atoi(av[i+1]);
379 <                argv[0] = "-i";
379 >                argv[0] = mg_ename[MG_E_OBJECT];
380                  argv[1] = cntbuf;
381                  for (j = 2; j+i < ac; j++)
382                          argv[j] = av[j+i];
383                  argv[j] = NULL;
384                  for (j = 0; j < niter; j++) {
385                          sprintf(cntbuf, "%d", j);
386 +                        if ((rval = handle_it(MG_E_OBJECT, 2, argv)) != MG_OK)
387 +                                return(rval);
388 +                        argv[0] = "-i";
389                          if ((rval = mg_iterate(ac-i, argv, f)) != MG_OK)
390                                  return(rval);
391 +                        argv[0] = mg_ename[MG_E_OBJECT];
392 +                        if ((rval = handle_it(MG_E_OBJECT, 1, argv)) != MG_OK)
393 +                                return(rval);
394                  }
395          } else if ((rval = (*f)()) != MG_OK)    /* else do this instance */
396                          return(rval);
# Line 511 | Line 544 | char   **av;
544                  return(MG_EARGC);
545          if ((cv = c_getvert(av[1])) == NULL)
546                  return(MG_EUNDEF);
547 <        if (cv->n[0]==0. && cv->n[1]==0. && cv->n[2]==0.)
547 >        if (is0vect(cv->n))
548                  return(MG_EILL);
549          if (!isflt(av[2]) || !isflt(av[3]))
550                  return(MG_ETYPE);
551          minrad = atof(av[2]);
552 +        round0(minrad);
553          maxrad = atof(av[3]);
554                                          /* check orientation */
555          if (minrad > 0.)
# Line 623 | Line 657 | char   **av;
657                  return(MG_EARGC);
658          if ((cv = c_getvert(av[1])) == NULL)
659                  return(MG_EUNDEF);
660 <        if (cv->n[0]==0. && cv->n[1]==0. && cv->n[2]==0.)
660 >        if (is0vect(cv->n))
661                  return(MG_EILL);
662          if (!isflt(av[2]) || !isflt(av[3]))
663                  return(MG_ETYPE);
664          minrad = atof(av[2]);
665 +        round0(minrad);
666          maxrad = atof(av[3]);
667          if (minrad < 0. || maxrad <= minrad)
668                  return(MG_EILL);
# Line 653 | Line 688 | char   **av;
688                                  sprintf(p3[j], FLTFMT, cv->p[j] +
689                                                  maxrad*u[j]*cos(theta) +
690                                                  maxrad*v[j]*sin(theta));
691 <                        if ((rv = handle_it(MG_E_VERTEX, 3, v3ent)) != MG_OK)
691 >                        if ((rv = handle_it(MG_E_VERTEX, 2, v3ent)) != MG_OK)
692                                  return(rv);
693                          if ((rv = handle_it(MG_E_POINT, 4, p3ent)) != MG_OK)
694                                  return(rv);
# Line 679 | Line 714 | char   **av;
714                                  sprintf(p3[j], FLTFMT, cv->p[j] + maxrad*d);
715                                  sprintf(p4[j], FLTFMT, cv->p[j] + minrad*d);
716                          }
717 <                        if ((rv = handle_it(MG_E_VERTEX, 3, v3ent)) != MG_OK)
717 >                        if ((rv = handle_it(MG_E_VERTEX, 2, v3ent)) != MG_OK)
718                                  return(rv);
719                          if ((rv = handle_it(MG_E_POINT, 4, p3ent)) != MG_OK)
720                                  return(rv);
721 <                        if ((rv = handle_it(MG_E_VERTEX, 3, v4ent)) != MG_OK)
721 >                        if ((rv = handle_it(MG_E_VERTEX, 2, v4ent)) != MG_OK)
722                                  return(rv);
723                          if ((rv = handle_it(MG_E_POINT, 4, p4ent)) != MG_OK)
724                                  return(rv);
# Line 728 | Line 763 | char   **av;
763          if (!isflt(av[2]) || !isflt(av[4]))
764                  return(MG_ETYPE);
765          rad1 = atof(av[2]);
766 +        round0(rad1);
767          rad2 = atof(av[4]);
768 +        round0(rad2);
769          if (rad1 == 0.) {
770                  if (rad2 == 0.)
771                          return(MG_EILL);
# Line 752 | Line 789 | char   **av;
789          if ((d = normalize(w)) == 0.)
790                  return(MG_EILL);
791          n1off = n2off = (rad2 - rad1)/d;
792 <        if (warpconends)                /* hack for e_sph and e_torus */
793 <                n2off = tan(atan(n2off)-(PI/4)/mg_nqcdivs);
794 <        n2off = sgn*n2off;
792 >        if (warpconends) {              /* hack for e_sph and e_torus */
793 >                d = atan(n2off) - (PI/4)/mg_nqcdivs;
794 >                if (d <= -PI/2+FTINY)
795 >                        n2off = -FHUGE;
796 >                else
797 >                        n2off = tan(d);
798 >        }
799          make_axes(u, v, w);
800          for (j = 0; j < 3; j++) {
801                  sprintf(p3[j], FLTFMT, cv2->p[j] + rad2*u[j]);
802 <                sprintf(n3[j], FLTFMT, u[j] + w[j]*n2off);
802 >                if (n2off <= -FHUGE)
803 >                        sprintf(n3[j], FLTFMT, -w[j]);
804 >                else
805 >                        sprintf(n3[j], FLTFMT, u[j] + w[j]*n2off);
806          }
807          if ((rv = handle_it(MG_E_VERTEX, 3, v3ent)) != MG_OK)
808                  return(rv);
# Line 781 | Line 825 | char   **av;
825                          for (j = 0; j < 3; j++) {
826                                  d = u[j]*cos(theta) + v[j]*sin(theta);
827                                  sprintf(p3[j], FLTFMT, cv2->p[j] + rad2*d);
828 <                                sprintf(n3[j], FLTFMT, d + w[j]*n2off);
828 >                                if (n2off > -FHUGE)
829 >                                        sprintf(n3[j], FLTFMT, d + w[j]*n2off);
830                          }
831 <                        if ((rv = handle_it(MG_E_VERTEX, 3, v3ent)) != MG_OK)
831 >                        if ((rv = handle_it(MG_E_VERTEX, 2, v3ent)) != MG_OK)
832                                  return(rv);
833                          if ((rv = handle_it(MG_E_POINT, 4, p3ent)) != MG_OK)
834                                  return(rv);
835 <                        if ((rv = handle_it(MG_E_NORMAL, 4, n3ent)) != MG_OK)
835 >                        if (n2off > -FHUGE &&
836 >                        (rv = handle_it(MG_E_NORMAL, 4, n3ent)) != MG_OK)
837                                  return(rv);
838                          if ((rv = handle_it(MG_E_FACE, 4, fent)) != MG_OK)
839                                  return(rv);
840                  }
841          } else {                        /* quads */
842                  v1ent[3] = "_cv4";
843 <                if (warpconends)                /* hack for e_sph and e_torus */
844 <                        n1off = tan(atan(n1off)+(PI/4)/mg_nqcdivs);
845 <                n1off = sgn*n1off;
843 >                if (warpconends) {              /* hack for e_sph and e_torus */
844 >                        d = atan(n1off) + (PI/4)/mg_nqcdivs;
845 >                        if (d >= PI/2-FTINY)
846 >                                n1off = FHUGE;
847 >                        else
848 >                                n1off = tan(atan(n1off)+(PI/4)/mg_nqcdivs);
849 >                }
850                  for (j = 0; j < 3; j++) {
851                          sprintf(p4[j], FLTFMT, cv1->p[j] + rad1*u[j]);
852 <                        sprintf(n4[j], FLTFMT, u[j] + w[j]*n1off);
852 >                        if (n1off >= FHUGE)
853 >                                sprintf(n4[j], FLTFMT, w[j]);
854 >                        else
855 >                                sprintf(n4[j], FLTFMT, u[j] + w[j]*n1off);
856                  }
857                  if ((rv = handle_it(MG_E_VERTEX, 3, v4ent)) != MG_OK)
858                          return(rv);
# Line 816 | Line 869 | char   **av;
869                          for (j = 0; j < 3; j++) {
870                                  d = u[j]*cos(theta) + v[j]*sin(theta);
871                                  sprintf(p3[j], FLTFMT, cv2->p[j] + rad2*d);
872 <                                sprintf(n3[j], FLTFMT, d + w[j]*n2off);
872 >                                if (n2off > -FHUGE)
873 >                                        sprintf(n3[j], FLTFMT, d + w[j]*n2off);
874                                  sprintf(p4[j], FLTFMT, cv1->p[j] + rad1*d);
875 <                                sprintf(n4[j], FLTFMT, d + w[j]*n1off);
875 >                                if (n1off < FHUGE)
876 >                                        sprintf(n4[j], FLTFMT, d + w[j]*n1off);
877                          }
878 <                        if ((rv = handle_it(MG_E_VERTEX, 3, v3ent)) != MG_OK)
878 >                        if ((rv = handle_it(MG_E_VERTEX, 2, v3ent)) != MG_OK)
879                                  return(rv);
880                          if ((rv = handle_it(MG_E_POINT, 4, p3ent)) != MG_OK)
881                                  return(rv);
882 <                        if ((rv = handle_it(MG_E_NORMAL, 4, n3ent)) != MG_OK)
882 >                        if (n2off > -FHUGE &&
883 >                        (rv = handle_it(MG_E_NORMAL, 4, n3ent)) != MG_OK)
884                                  return(rv);
885 <                        if ((rv = handle_it(MG_E_VERTEX, 3, v4ent)) != MG_OK)
885 >                        if ((rv = handle_it(MG_E_VERTEX, 2, v4ent)) != MG_OK)
886                                  return(rv);
887                          if ((rv = handle_it(MG_E_POINT, 4, p4ent)) != MG_OK)
888                                  return(rv);
889 <                        if ((rv = handle_it(MG_E_NORMAL, 4, n4ent)) != MG_OK)
889 >                        if (n1off < FHUGE &&
890 >                        (rv = handle_it(MG_E_NORMAL, 4, n4ent)) != MG_OK)
891                                  return(rv);
892                          if ((rv = handle_it(MG_E_FACE, 5, fent)) != MG_OK)
893                                  return(rv);
894                  }
895 +        }
896 +        return(MG_OK);
897 + }
898 +
899 +
900 + static int
901 + e_prism(ac, av)                 /* turn a prism into polygons */
902 + int     ac;
903 + char    **av;
904 + {
905 +        static char     p[3][24];
906 +        static char     *vent[4] = {mg_ename[MG_E_VERTEX],NULL,"="};
907 +        static char     *pent[5] = {mg_ename[MG_E_POINT],p[0],p[1],p[2]};
908 +        char    *newav[MG_MAXARGC], nvn[MG_MAXARGC-1][8];
909 +        double  length;
910 +        FVECT   v1, v2, v3, norm;
911 +        register C_VERTEX       *cv;
912 +        C_VERTEX        *cv0;
913 +        int     rv;
914 +        register int    i, j;
915 +
916 +        if (ac < 5)
917 +                return(MG_EARGC);
918 +        if (!isflt(av[ac-1]))
919 +                return(MG_ETYPE);
920 +        length = atof(av[ac-1]);
921 +        if (length <= FTINY && length >= -FTINY)
922 +                return(MG_EILL);
923 +                                        /* do bottom face */
924 +        newav[0] = mg_ename[MG_E_FACE];
925 +        for (i = 1; i < ac-1; i++)
926 +                newav[i] = av[i];
927 +        newav[i] = NULL;
928 +        if ((rv = handle_it(MG_E_FACE, i, newav)) != MG_OK)
929 +                return(rv);
930 +                                        /* compute face normal */
931 +        if ((cv0 = c_getvert(av[1])) == NULL)
932 +                return(MG_EUNDEF);
933 +        norm[0] = norm[1] = norm[2] = 0.;
934 +        v1[0] = v1[1] = v1[2] = 0.;
935 +        for (i = 2; i < ac-1; i++) {
936 +                if ((cv = c_getvert(av[i])) == NULL)
937 +                        return(MG_EUNDEF);
938 +                v2[0] = cv->p[0] - cv0->p[0];
939 +                v2[1] = cv->p[1] - cv0->p[1];
940 +                v2[2] = cv->p[2] - cv0->p[2];
941 +                fcross(v3, v1, v2);
942 +                norm[0] += v3[0];
943 +                norm[1] += v3[1];
944 +                norm[2] += v3[2];
945 +                VCOPY(v1, v2);
946 +        }
947 +        if (normalize(norm) == 0.)
948 +                return(MG_EILL);
949 +                                        /* create moved vertices */
950 +        for (i = 1; i < ac-1; i++) {
951 +                sprintf(nvn[i-1], "_pv%d", i);
952 +                vent[1] = nvn[i-1];
953 +                if ((rv = handle_it(MG_E_VERTEX, 3, vent)) != MG_OK)
954 +                        return(rv);
955 +                cv = c_getvert(av[i]);          /* checked above */
956 +                for (j = 0; j < 3; j++)
957 +                        sprintf(p[j], FLTFMT, cv->p[j] - length*norm[j]);
958 +                if ((rv = handle_it(MG_E_POINT, 4, pent)) != MG_OK)
959 +                        return(rv);
960 +                newav[ac-1-i] = nvn[i-1];       /* reverse */
961 +        }
962 +                                                /* do top face */
963 +        if ((rv = handle_it(MG_E_FACE, ac-1, newav)) != MG_OK)
964 +                return(rv);
965 +                                                /* do the side faces */
966 +        newav[5] = NULL;
967 +        newav[3] = av[ac-2];
968 +        newav[4] = nvn[ac-3];
969 +        for (i = 1; i < ac-1; i++) {
970 +                newav[1] = nvn[i-1];
971 +                newav[2] = av[i];
972 +                if ((rv = handle_it(MG_E_FACE, 5, newav)) != MG_OK)
973 +                        return(rv);
974 +                newav[3] = newav[2];
975 +                newav[4] = newav[1];
976 +        }
977 +        return(MG_OK);
978 + }
979 +
980 +
981 + static int
982 + e_cspec(ac, av)                 /* handle spectral color */
983 + int     ac;
984 + char    **av;
985 + {
986 +        static char     xbuf[24], ybuf[24];
987 +        static char     *ccom[4] = {mg_ename[MG_E_CXY], xbuf, ybuf};
988 +        int     rv;
989 +
990 +        c_ccvt(c_ccolor, C_CSXY);
991 +                                /* if it's really their handler, use it */
992 +        if (mg_ehand[MG_E_CXY] != c_hcolor) {
993 +                sprintf(xbuf, "%.4f", c_ccolor->cx);
994 +                sprintf(ybuf, "%.4f", c_ccolor->cy);
995 +                if ((rv = handle_it(MG_E_CXY, 3, ccom)) != MG_OK)
996 +                        return(rv);
997 +        }
998 +        return(MG_OK);
999 + }
1000 +
1001 +
1002 + static int
1003 + e_cmix(ac, av)                  /* handle mixing of colors */
1004 + int     ac;
1005 + char    **av;
1006 + {
1007 +        char    wl[2][6], vbuf[C_CNSS][24];
1008 +        char    *newav[C_CNSS+4];
1009 +        int     rv;
1010 +        register int    i;
1011 +        /*
1012 +         * Contorted logic works as follows:
1013 +         *      1. the colors are already mixed in c_hcolor() support function
1014 +         *      2. if we would handle a spectral result, make sure it's not
1015 +         *      3. if c_hcolor() would handle a spectral result, don't bother
1016 +         *      4. otherwise, make cspec entity and pass it to their handler
1017 +         *      5. if we have only xy results, handle it as c_spec() would
1018 +         */
1019 +        if (mg_ehand[MG_E_CSPEC] == e_cspec)
1020 +                c_ccvt(c_ccolor, C_CSXY);
1021 +        else if (c_ccolor->flags & C_CDSPEC) {
1022 +                if (mg_ehand[MG_E_CSPEC] != c_hcolor) {
1023 +                        sprintf(wl[0], "%d", C_CMINWL);
1024 +                        sprintf(wl[1], "%d", C_CMAXWL);
1025 +                        newav[0] = mg_ename[MG_E_CSPEC];
1026 +                        newav[1] = wl[0];
1027 +                        newav[2] = wl[1];
1028 +                        for (i = 0; i < C_CNSS; i++) {
1029 +                                sprintf(vbuf[i], "%.6f",
1030 +                                                (double)c_ccolor->ssamp[i] /
1031 +                                                c_ccolor->ssum);
1032 +                                newav[i+3] = vbuf[i];
1033 +                        }
1034 +                        newav[C_CNSS+3] = NULL;
1035 +                        if ((rv = handle_it(MG_E_CSPEC, C_CNSS+3, newav)) != MG_OK)
1036 +                                return(rv);
1037 +                }
1038 +                return(MG_OK);
1039 +        }
1040 +        if (mg_ehand[MG_E_CXY] != c_hcolor) {
1041 +                sprintf(vbuf[0], "%.4f", c_ccolor->cx);
1042 +                sprintf(vbuf[1], "%.4f", c_ccolor->cy);
1043 +                newav[0] = mg_ename[MG_E_CXY];
1044 +                newav[1] = vbuf[0];
1045 +                newav[2] = vbuf[1];
1046 +                newav[3] = NULL;
1047 +                if ((rv = handle_it(MG_E_CXY, 3, newav)) != MG_OK)
1048 +                        return(rv);
1049          }
1050          return(MG_OK);
1051   }

Diff Legend

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