rawfile.c

ngspice又一个电子CAD仿真软件代码.功能更全

/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1986 Wayne A. Christopher, U. C. Berkeley CAD Group
$Id: rawfile.c,v 1.8 2005/05/24 03:03:00 sjborley Exp $
**********/

/*
 * Read and write the ascii and binary rawfile formats.
 */

#include "ngspice.h"
#include "cpdefs.h"
#include "ftedefs.h"
#include "dvec.h"

#include "rawfile.h"
#include "variable.h"

/* static declarations */
static void fixdims(struct dvec *v, char *s);


int raw_prec = -1;        /* How many sigfigs to use, default 15 (max).  */

#define DEFPREC 15

#ifdef HAS_WINDOWS
#undef fscanf             /* redo I/O from WINMAIN.C here
                          otherwise reading ASCII will not work */
#endif


/* Write a raw file.  We write everything in the plot pointed to. */

void
raw_write(char *name, struct plot *pl, bool app, bool binary)
{
    FILE *fp;
    bool realflag = TRUE, writedims;
    bool raw_padding;
    int length, numdims, dims[MAXDIMS];
    int nvars, i, j, prec;
    struct dvec *v, *lv;
    wordlist *wl;
    struct variable *vv;
    double dd;
    char buf[BSIZE_SP];

    if (!cp_getvar("nopadding", VT_BOOL, (char *) &raw_padding))
        raw_padding = FALSE;
    /* Invert since we want to know if we should pad. */
    raw_padding = !raw_padding;

    /* Why bother printing out an empty plot? */
    if (!pl->pl_dvecs) {
	fprintf(cp_err, "Error: plot is empty, nothing written.\n");
	return;
    }

    if (raw_prec != -1)
        prec = raw_prec;
    else
        prec = DEFPREC;

#ifdef __MINGW32__
/* - Binary file binary write -  hvogt 15.03.2000 ---------------------*/
    if (binary) {
        if (!(fp = fopen(name, app ? "ab" : "wb"))) {
            perror(name);
            return;
        }
        fprintf(cp_out,"binary raw file\n");
    }
    else  {
        if (!(fp = fopen(name, app ? "a" : "w"))) {
            perror(name);
            return;
        }
        fprintf(cp_out,"ASCII raw file\n");
    }
/* --------------------------------------------------------------------*/

#else


    if (!(fp = fopen(name, app ? "a" : "w"))) {
        perror(name);
        return;
    }

#endif

    numdims = nvars = length = 0;
    for (v = pl->pl_dvecs; v; v = v->v_next) {
        if (iscomplex(v))
            realflag = FALSE;
        nvars++;
	/* Find the length and dimensions of the longest vector
	 * in the plot.
	 * Be paranoid and assume somewhere we may have
	 * forgotten to set the dimensions of 1-D vectors.
	 */
	if (v->v_numdims <= 1) {
	    v->v_numdims = 1;
	    v->v_dims[0] = v->v_length;
	}
	if (v->v_length > length) {
	    length = v->v_length;
	    numdims = v->v_numdims;
	    for (j = 0; j < numdims; j++) {
		dims[j] = v->v_dims[j];
	    }
	}
    }

    fprintf(fp, "Title: %s\n", pl->pl_title);
    fprintf(fp, "Date: %s\n", pl->pl_date);
    fprintf(fp, "Plotname: %s\n", pl->pl_name);
    fprintf(fp, "Flags: %s%s\n",
	    realflag ? "real" : "complex", raw_padding ? "" : " unpadded" );
    fprintf(fp, "No. Variables: %d\n", nvars);
    fprintf(fp, "No. Points: %d\n", length);
    if (numdims > 1) {
    	dimstring(dims, numdims, buf);
	fprintf(fp, "Dimensions: %s\n", buf);
    }

    for (wl = pl->pl_commands; wl; wl = wl->wl_next)
        fprintf(fp, "Command: %s\n", wl->wl_word);

    for (vv = pl->pl_env; vv; vv = vv->va_next) {
        wl = cp_varwl(vv);
        if (vv->va_type == VT_BOOL) {
            fprintf(fp, "Option: %s\n", vv->va_name);
        } else {
            fprintf(fp, "Option: %s = ", vv->va_name);
            if (vv->va_type == VT_LIST)
                fprintf(fp, "( ");
            wl_print(wl, fp);
            if (vv->va_type == VT_LIST)
                fprintf(fp, " )");
            (void) putc('\n', fp);
        }
    }

    /* Before we write the stuff out, make sure that the scale is the first
     * in the list.
     */
    for (lv = NULL, v = pl->pl_dvecs; v != pl->pl_scale; v = v->v_next)
        lv = v;
    if (lv) {
        lv->v_next = v->v_next;
        v->v_next = pl->pl_dvecs;
        pl->pl_dvecs = v;
    }

    fprintf(fp, "Variables:\n");
    for (i = 0, v = pl->pl_dvecs; v; v = v->v_next) {
        fprintf(fp, "\t%d\t%s\t%s", i++, v->v_name,
                ft_typenames(v->v_type));
        if (v->v_flags & VF_MINGIVEN)
            fprintf(fp, " min=%e", v->v_minsignal);
        if (v->v_flags & VF_MAXGIVEN)
            fprintf(fp, " max=%e", v->v_maxsignal);
        if (v->v_defcolor)
            fprintf(fp, " color=%s", v->v_defcolor);
        if (v->v_gridtype)
            fprintf(fp, " grid=%d", v->v_gridtype);
        if (v->v_plottype)
            fprintf(fp, " plot=%d", v->v_plottype);
	/* Only write dims if they are different from default. */
	writedims = FALSE;
	if (v->v_numdims != numdims) {
	  writedims = TRUE;
	} else {
	    for (j = 0; j < numdims; j++)
		if (dims[j] != v->v_dims[j])
		    writedims = TRUE;
	}
	if (writedims) {
	    dimstring(v->v_dims, v->v_numdims, buf);	
	    fprintf(fp, " dims=%s",buf);
        }
        (void) putc('\n', fp);
    }

    if (binary) {
        fprintf(fp, "Binary:\n");
        for (i = 0; i < length; i++) {
            for (v = pl->pl_dvecs; v; v = v->v_next) {
		/* Don't run off the end of this vector's data. */
		if (i < v->v_length) {
		    if (realflag) {
			dd = (isreal(v) ? v->v_realdata[i] :
			    realpart(&v->v_compdata[i]));
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
		    } else if (isreal(v)) {
			dd = v->v_realdata[i];
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
			dd = 0.0;
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
		    } else {
			dd = realpart(&v->v_compdata[i]);
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
			dd = imagpart(&v->v_compdata[i]);
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
		    }
                } else if (raw_padding) {
		    dd = 0.0;
		    if (realflag) {
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
		    } else {
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
			(void) fwrite((char *) &dd, sizeof
				(double), 1, fp);
		    }
		}
            }
        }
    } else {
        fprintf(fp, "Values:\n");
        for (i = 0; i < length; i++) {
            fprintf(fp, " %d", i);
	    for (v = pl->pl_dvecs; v; v = v->v_next) {
		if (i < v->v_length) {
		    if (realflag)
			fprintf(fp, "\t%.*e\n", prec,
			    isreal(v) ? v->v_realdata[i] :
			    realpart(&v->v_compdata[i]));
		    else if (isreal(v))
			fprintf(fp, "\t%.*e,0.0\n", prec,
			    v->v_realdata[i]);
		    else
			fprintf(fp, "\t%.*e,%.*e\n", prec,
			    realpart(&v->v_compdata[i]),
			    prec,
			    imagpart(&v->v_compdata[i]));
                } else if (raw_padding) {
		    if (realflag) {
			fprintf(fp, "\t%.*e\n", prec, 0.0);
		    } else {
			fprintf(fp, "\t%.*e,%.*e\n",
				prec, 0.0, prec, 0.0);
		    }
		}
	    }
            (void) putc('\n', fp);
        }
    }
    (void) fclose(fp);
    return;
}

/* Read a raw file.  Returns a list of plot structures.  This routine should be
 * very flexible about what it expects to see in the rawfile.  Really all we
 * require is that there be one variables and one values section per plot
 * and that the variables precede the values.
 */

#define skip(s) while (*(s) && !isspace(*(s)))(s)++; while (isspace(*(s)))(s)++
#define nonl(s) r = (s); while (*r && (*r != '\n')) r++; *r = '\0'

struct plot *
raw_read(char *name)
{
    char *title = "default title";
    char *date = 0;
    struct plot *plots = NULL, *curpl = NULL;
    char buf[BSIZE_SP], buf2[BSIZE_SP], *s, *t, *r;
    int flags = 0, nvars = 0, npoints = 0, i, j;
    int ndimpoints, numdims=0, dims[MAXDIMS];
    bool raw_padded = TRUE;
    double junk;
    struct dvec *v, *nv;
    struct variable *vv;
    wordlist *wl, *nwl;
    FILE *fp, *lastin, *lastout, *lasterr;


#ifdef __MINGW32__
    bool binary = TRUE;
#endif

    if (!(fp = fopen(name, "r"))) {
        perror(name);
        return (NULL);
    }

#ifdef __MINGW32__
/* Test, whether file really ASCII, otherwise assume binary  hvogt 15.3.2000 */
    while (fgets(buf, BSIZE_SP, fp)) {
        if (ciprefix("values:", buf)) {
        	binary = FALSE;
        	rewind(fp);                /* rewind */
                fprintf(cp_err, "\nASCII raw file\n");
        	break;
        }
    }

    if (binary) {
        (void) fclose(fp);
        if (!(fp = fopen(name, "rb"))) {
            perror(name);
            return (NULL);
            }
            fprintf(cp_err, "\nbinary raw file\n");
        }
/*--------------------------------------------------------*/
#endif

    /* Since we call cp_evloop() from here, we have to do this junk. */
    lastin = cp_curin;
    lastout = cp_curout;
    lasterr = cp_curerr;
    cp_curin = cp_in;
    cp_curout = cp_out;
    cp_curerr = cp_err;

    cp_pushcontrol();

    while (fgets(buf, BSIZE_SP, fp)) {
        /* Figure out what this line is... */
        if (ciprefix("title:", buf)) {
            s = buf;
            skip(s);
            nonl(s);
            title = copy(s);
        } else if (ciprefix("date:", buf)) {
            s = buf;
            skip(s);
            nonl(s);
            date = copy(s);
        } else if (ciprefix("plotname:", buf)) {
            s = buf;
            skip(s);
            nonl(s);
            if (curpl) {    /* reverse commands list */
                for (wl=curpl->pl_commands,
                        curpl->pl_commands=NULL; wl &&
                        wl->wl_next; wl=nwl) {
                    nwl = wl->wl_next;
                    wl->wl_next = curpl->pl_commands;
                    curpl->pl_commands = wl;
                }
            }
            curpl = alloc(struct plot);
            curpl->pl_next = plots;
            plots = curpl;
            curpl->pl_name = copy(s);
	    if (!date)
		date = copy(datestring( ));
            curpl->pl_date = date;
            curpl->pl_title = copy(title);
            flags = VF_PERMANENT;
            nvars = npoints = 0;
        } else if (ciprefix("flags:", buf)) {
            s = buf;
            skip(s);
            while ((t = gettok(&s))) {
                if (cieq(t, "real"))