Stars of Shareware: Raytrace & Morphing/SOS-RAYTRACE.ISO/programm/source/radsrc22/src/cv/thf2rad2.c

C/C++ Source or Header | 1989-02-15 | 4.1 KB | 193 lines

/* Copyright (c) 1988 Regents of the University of California */

/*
 *  thf2rad - convert GDS things file to Radiance scene description.
 *
 *    12/21/88
 */

#include <stdio.h>

#define MAXSTR        128        /* maximum string or id length */
#define MAXPTS        1024        /* maximum points per object */

typedef struct {
    double    x, y, z;
} POINT;

double    rad = 0.0;            /* line radius */


main(argc, argv)
int    argc;
char    *argv[];
{
    extern double    atof();
    char    *name = NULL;
    int    an;

    for (an = 1; an < argc; an++)
        if (argv[an][0] == '-')
            switch (argv[an][1]) {
            case 'n':
                name = argv[++an];
                break;
            case 'r':
                rad = atof(argv[++an]);
                break;
            default:
                fprintf(stderr,
                "Usage: %s [-n name][-r rad] [file ..]\n",
                        argv[0]);
                exit(1);
            }
        else
            break;
    header(argc, argv);
    if (an < argc)
        for ( ; an < argc; an++)
            thf2rad(name, argv[an]);
    else
        thf2rad(name, NULL);
    exit(0);
}


thf2rad(name, file)        /* convert things file */
char    *name, *file;
{
    char    nambuf[MAXSTR];
    register char    *cp;
    FILE    *fp;
    int    on;
                /* open file and hack name */
    if (file == NULL) {
        fp = stdin;
        file = "<stdin>";
        if (name == NULL)
            name = "thing";
    } else {
        if ((fp = fopen(file, "r")) == NULL) {
            fprintf(stderr, "%s: cannot open\n", file);
            exit(1);
        }
        if (name == NULL)
            name = file;
    }
    for (cp = nambuf; *cp = *name++; cp++)
        ;
                /* get objects from file */
    on = 0;
    while (cvobject(nambuf, file, fp))
        sprintf(cp, ".%d", ++on);
}


cvobject(name, file, fp)    /* convert next object in things file */
char    *name, *file;
FILE    *fp;
{
    static POINT    parr[MAXPTS];
    char    sbuf[MAXSTR], nambuf[MAXSTR];
    int    npts, n, nv, p1, p2;
    register int    i;
                        /* get points */
    if (fscanf(fp, "%s", sbuf) != 1)
        return(0);
    if (strcmp(sbuf, "YIN"))
        readerr(file, "not a GDS things file");
    if (fscanf(fp, "%d", &npts) != 1 || npts < 0 || npts > MAXPTS)
        readerr(file, "bad number of points");
    for (i = 0; i < npts; i++)
        if (fscanf(fp, "%lf %lf %lf",
                &parr[i].x, &parr[i].y, &parr[i].z) != 3)
            readerr(file, "bad point");
                        /* get lines */
    if (fscanf(fp, "%d", &n) != 1 || n < 0)
        readerr(file, "bad number of lines");
    for (i = 0; i < n; i++) {
        if (fscanf(fp, "%d %d %s", &p1, &p2, sbuf) != 3)
            readerr(file, "error reading vector");
        if (--p1 < 0 || p1 >= npts || --p2 < 0 || p2 >= npts)
            readerr(file, "bad point in vector");
        sprintf(nambuf, "%s.l%d.c", name, i+1);
        line(sbuf, nambuf, &parr[p1], &parr[p2]);
    }
                        /* get surfaces */
    if (fscanf(fp, "%d", &n) != 1 || n < 0)
        readerr(file, "bad number of surfaces");
    for (i = 0; i < n; i++) {
        if (fscanf(fp, "%d %s %d", &p1, sbuf, &nv) != 3)
            readerr(file, "error reading surface");
        if (p1 != 1 || nv < 0)
            readerr(file, "bad surface type (has hole?)");
        if (nv < 3) {
            while (nv--)
                fscanf(fp, "%*d");
            continue;
        }
        printf("\n%s polygon %s.p%d\n", sbuf, name, i+1);
        printf("0\n0\n%d\n", 3*nv);
        while (nv--) {
            if (fscanf(fp, "%d", &p1) != 1)
                readerr(file, "error reading surface vertex");
            loc(&parr[p1-1]);
        }
    }
    return(1);
}


header(ac, av)            /* print header */
int    ac;
char    **av;
{
    putchar('#');
    while (ac--) {
        putchar(' ');
        fputs(*av++, stdout);
    }
    putchar('\n');
}


line(mat, nam, p1, p2)        /* print line */
char    *mat, *nam;
POINT    *p1, *p2;
{
    if (rad <= 0.0)
        return;
    printf("\n%s ring %s.r0\n", mat, nam);
    printf("0\n0\n8\n");
    loc(p1);
    printf("\t%14.8g\t%14.8g\t%14.8g\n",
            p1->x-p2->x, p1->y-p2->y, p1->z-p2->z);
    printf("\t0\t%14.8g\n", rad);
    printf("\n%s cylinder %s\n", mat, nam);
    printf("0\n0\n7\n");
    loc(p1);
    loc(p2);
    printf("\t%g\n", rad);
    printf("\n%s ring %s.r1\n", mat, nam);
    printf("0\n0\n8\n");
    loc(p2);
    printf("\t%14.8g\t%14.8g\t%14.8g\n",
            p2->x-p1->x, p2->y-p1->y, p2->z-p1->z);
    printf("\t0\t%14.8g\n", rad);
}


loc(p)                /* print location */
register POINT    *p;
{
    printf("\t%14.8g\t%14.8g\t%14.8g\n", p->x, p->y, p->z);
}


readerr(file, err)        /* print read error and exit */
char    *file, *err;
{
    fprintf(stderr, "%s: %s\n", file, err);
    exit(1);
}