/* Copyright 2003-2008 Multigig Ltd
 *
 * Copied and written by Stefan Jones (stefan.jones@multigig.com) at Multigig Ltd
 *
 * Code based on and copied from ScriptEDA ( http://www-cad.eecs.berkeley.edu/~pinhong/scriptEDA )
 *
 * Under LGPLv2 licence since 2008, December 1st
 */

/*******************/
/*   Defines       */
/*******************/

#define TCLSPICE_name    "spice"
#define TCLSPICE_prefix  "spice::"
#define TCLSPICE_namespace "spice"
#ifdef _MSC_VER
#define TCLSPICE_version "25.1"
#define STDIN_FILENO    0
#define STDOUT_FILENO   1
#define STDERR_FILENO   2
#endif

/**********************************************************************/
/*              Header files for C functions                          */
/**********************************************************************/

#include "ngspice/ngspice.h"
#include "ngspice/randnumb.h"
#include "misc/misc_time.h"
#include <tcl.h>

/*Use Tcl threads if on W32 without pthreads*/
#ifndef HAVE_LIBPTHREAD

#if defined(__MINGW32__) || defined(_MSC_VER)

#define mutex_lock(a) Tcl_MutexLock(a)
#define mutex_unlock(a) Tcl_MutexUnlock(a)
#define thread_self() Tcl_GetCurrentThread()
typedef Tcl_Mutex mutexType;
typedef Tcl_ThreadId threadId_t;
#define TCL_THREADS
#define THREADS

#endif

#else

#include <pthread.h>
#define mutex_lock(a) pthread_mutex_lock(a)
#define mutex_unlock(a) pthread_mutex_unlock(a)
#define thread_self() pthread_self()
typedef pthread_mutex_t mutexType;
typedef pthread_t threadId_t;
#define THREADS

#endif


/* Copied from main.c in ngspice*/
#include <stdio.h>
#if defined(__MINGW32__)
#include <stdarg.h>
/* remove type incompatibility with winnt.h*/
#undef BOOLEAN
#include <windef.h>
#include <winbase.h>  /* Sleep */
#elif defined(_MSC_VER)
#include <stdarg.h>
/* remove type incompatibility with winnt.h*/
#undef BOOLEAN
#include <windows.h> /* Sleep */
#include <process.h> /* _getpid */
#define dup _dup
#define dup2 _dup2
#define open _open
#define close _close
#else
#include <unistd.h> /* usleep */
#endif /* __MINGW32__ */

#include "ngspice/iferrmsg.h"
#include "ngspice/ftedefs.h"
#include "ngspice/devdefs.h"
#include <spicelib/devices/dev.h>
#include <spicelib/analysis/analysis.h>
#include <misc/ivars.h>
#include <frontend/resource.h>
#include <frontend/com_measure2.h>
#ifdef _MSC_VER
#include <stdio.h>
#define snprintf _snprintf
#endif
#include <frontend/outitf.h>
#include "ngspice/memory.h"
#include <frontend/com_measure2.h>

#ifndef HAVE_GETRUSAGE
#ifdef HAVE_FTIME
#include <sys/timeb.h>
#endif
#endif

/* To interupt a spice run */
#include <signal.h>
typedef void (*sighandler)(int);

#include <setjmp.h>
#include "frontend/signal_handler.h"

/*Included for the module to access data*/
#include "ngspice/dvec.h"
#include "ngspice/plot.h"

#ifdef __CYGWIN__
#undef WIN32
#endif
#include <blt.h>
#include  "ngspice/sim.h"

/* defines for Tcl support
 * Tcl 8.3 and Tcl 8.4 support,
 * suggested by http://mini.net/tcl/3669, 07.03.03 */
#ifndef CONST84
#define CONST84
#endif
/* Arguments of Tcl_CmpProc for Tcl/Tk 8.4.x */
#define TCL_CMDPROCARGS(clientData, interp, argc, argv)                 \
    (ClientData clientData, Tcl_Interp *interp, int argc, CONST84 char *argv[])

/*For get_output*/
#include <fcntl.h>
#include <sys/stat.h>

#ifdef _MSC_VER
#define S_IRWXU _S_IWRITE
#endif

extern IFfrontEnd nutmeginfo;

extern struct comm spcp_coms[ ];
extern void DevInit(void);
extern int SIMinit(IFfrontEnd *frontEnd, IFsimulator **simulator);
extern wordlist *cp_varwl(struct variable *var);

/*For blt spice to use*/
typedef struct {
    char *name;
#ifdef THREADS
    mutexType mutex;            /*lock for this vector*/
#endif
    double *data;               /* vector data*/
    int size;                   /*data it can store*/
    int length;                 /*actual amount of data*/
} vector;

/*The current run (to get variable names, etc)*/
static runDesc *cur_run;

static vector *vectors;

static int ownVectors;

/* save this each time called */
static Tcl_Interp *spice_interp;
#define save_interp()                           \
    do {                                        \
        spice_interp = interp;                  \
    } while(0)


void tcl_stdflush(FILE *f);
int  tcl_vfprintf(FILE *f, const char *fmt, va_list args);
#if defined(__MINGW32__) || defined(_MSC_VER)
__declspec(dllexport)
#endif
int  Spice_Init(Tcl_Interp *interp);
int  Tcl_ExecutePerLoop(void);
void triggerEventCheck(ClientData clientData, int flags);
void triggerEventSetup(ClientData clientData, int flags);
int  triggerEventHandler(Tcl_Event *evPtr, int flags);
void stepEventCheck(ClientData clientData, int flags);
int  stepEventHandler(Tcl_Event *evPtr, int flags);
void stepEventSetup(ClientData clientData, int flags);
int  sp_Tk_Update(void);
int  sp_Tk_SetColor(int colorid);
int  sp_Tk_SetLinestyle(int linestyleid);
int  sp_Tk_DefineLinestyle(int linestyleid, int mask);
int  sp_Tk_DefineColor(int colorid, double red, double green, double blue);
int  sp_Tk_Text(char *text, int x, int y, int angle);
int  sp_Tk_Arc(int x0, int y0, int radius, double theta, double delta_theta);
int  sp_Tk_DrawLine(int x1, int y1, int x2, int y2);
int  sp_Tk_Clear(void);
int  sp_Tk_Close(void);
int  sp_Tk_NewViewport(GRAPH *graph);
int  sp_Tk_Init(void);
int  get_mod_param TCL_CMDPROCARGS(clientData, interp, argc, argv);
void sighandler_tclspice(int num);
void blt_relink(int index, void *tmp);
void blt_lockvec(int index);
void blt_add(int index, double value);
void blt_init(void *run);
int  blt_plot(struct dvec *y, struct dvec *x, int new);


/****************************************************************************/
/*                          BLT and data routines                           */
/****************************************************************************/

/*helper function*/
/*inline*/
static struct plot *
get_plot_by_index(int plot)
{
    struct plot *pl;
    pl = plot_list;
    for (; 0 < plot; plot--) {
        pl = pl->pl_next;
        if (!pl)
            return NULL;
    }
    return pl;
}


/*this holds the number of time points done (altered by spice)*/
int steps_completed;

/* number of bltvectors*/
static int blt_vnum;


/*Native Tcl functions */

static int
spice_header TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    char buf[256];
    char *date, *name, *title;
    NG_IGNORE(clientData);
    NG_IGNORE(argv);
    if (argc != 1) {
        Tcl_SetResult(interp, "Wrong # args. spice::spice_header", TCL_STATIC);
        return TCL_ERROR;
    }
    if (cur_run) {
        Tcl_ResetResult(interp);
        date = datestring();
        title = cur_run->name;
        name = cur_run->type;
        sprintf(buf, "{title \"%s\"} {name \"%s\"} {date \"%s\"} {variables %u}", title, name, date, cur_run->numData);
        Tcl_AppendResult(interp, buf, TCL_STATIC);
        return TCL_OK;
    } else {
        return TCL_ERROR;
    }
}


static int
spice_data TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    char buf[256];
    int i, type;
    char *name;
    NG_IGNORE(clientData);
    if (argc > 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::spice_data ?plot?",
                      TCL_STATIC);
        return TCL_ERROR;
    }
    if (argc == 1) {
        if (blt_vnum) {
            Tcl_ResetResult(interp);
            for (i = 0; i < blt_vnum; i++) {
                name = vectors[i].name;
                if (substring("#branch", name))
                    type = SV_CURRENT;
                else if (cieq(name, "time"))
                    type = SV_TIME;
                else if (cieq(name, "frequency"))
                    type = SV_FREQUENCY;
                else
                    type = SV_VOLTAGE;
                sprintf(buf, "{%s %s} ", name, ft_typenames(type));
                Tcl_AppendResult(interp, buf, TCL_STATIC);
            }
            return TCL_OK;
        } else {
            return TCL_ERROR;
        }
    } else {
        struct plot *pl;
        struct dvec *v;
        if (!(pl = get_plot_by_index(atoi(argv[1])))) {
            Tcl_SetResult(interp, "Bad plot number", TCL_STATIC);
            return TCL_ERROR;
        }
        for (v = pl->pl_dvecs; v; v = v->v_next) {
            name = v->v_name;
            if (substring("#branch", name))
                type = SV_CURRENT;
            else if (cieq(name, "time"))
                type = SV_TIME;
            else if (cieq(name, "frequency"))
                type = SV_FREQUENCY;
            else
                type = SV_VOLTAGE;
            sprintf(buf, "{%s %s} ", name, ft_typenames(type));
            Tcl_AppendResult(interp, buf, TCL_STATIC);
        }
        return TCL_OK;
    }
}


static int resetTriggers(void);


/*Creates and registers the blt vectors, used by spice*/
void
blt_init(void *run)
{
    int i;
    cur_run = NULL;
    /* reset varaibles and free*/
    if (vectors) {
        resetTriggers();
        for (i = blt_vnum-1, blt_vnum = 0 /*stops vector access*/; i >= 0; i--) {
            if (ownVectors)
                FREE(vectors[i].data);
            FREE(vectors[i].name);
#ifdef HAVE_LIBPTHREAD
            pthread_mutex_destroy(&vectors[i].mutex);
#endif
        }
        FREE(vectors);
    }


    /* initilise */
    cur_run = (runDesc *)run;
    vectors = TMALLOC(vector, cur_run->numData);
    for (i = 0; i < cur_run->numData; i++) {
        vectors[i].name = copy((cur_run->data[i]).name);
#ifdef HAVE_LIBPTHREAD
        pthread_mutex_init(&vectors[i].mutex, NULL);
#endif
        vectors[i].data = NULL;
        vectors[i].size = 0;
        vectors[i].length = 0;
    }
    ownVectors = cur_run->writeOut;
    blt_vnum = i;               /*allows access to vectors*/
    return;
}


/*Adds data to the stored vector*/
void
blt_add(int index, double value)
{
    vector *v;
    v = &vectors[index];
#ifdef THREADS
    mutex_lock(&vectors[index].mutex);
#endif
    if (!(v->length < v->size)) {
        v->size += 100;
        v->data = TREALLOC(double, v->data, v->size);
    }
    v->data[v->length] = value;
    v->length ++;
#ifdef THREADS
    mutex_unlock(&vectors[index].mutex);
#endif
    return;
}


/* Locks the vector data to stop conflicts*/
void
blt_lockvec(int index)
{
#ifdef THREADS
    mutex_lock(&vectors[index].mutex);
#else
    NG_IGNORE(index);
#endif
    return;
}


/*links a dvec to a blt vector, used to stop duplication of data when writing to a plot,
  but makes BLT vectors more unsafe */
void
blt_relink(int index, void *tmp)
{
    struct dvec *v = (struct dvec *)tmp;
    vectors[index].data = v->v_realdata;
    vectors[index].length = v->v_length;
    vectors[index].size = v->v_length; /*silly spice doesn't use v_rlength*/
#ifdef THREADS
    mutex_unlock(&vectors[index].mutex);
#endif
    return;
}


/*        Tcl functions to access spice data   */

/* This copys the last Spice state vector to the given blt_vector
 * arg1: blt_vector
 */
static int
lastVector TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    Blt_Vector *vec;
    char *blt;
    int i;
    double *V;
    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::lastVector vecName", TCL_STATIC);
        return TCL_ERROR;
    }
    Tcl_SetResult(interp, "test2", TCL_STATIC);
    return TCL_ERROR;
    blt = (char *)argv[1];
    if (Blt_GetVector(interp, blt, &vec)) {
        Tcl_SetResult(interp, "Bad blt vector ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)blt, TCL_STATIC);
        return TCL_ERROR;
    }
    if (!(V = TMALLOC(double, blt_vnum))) {
        Tcl_SetResult(interp, "Out of Memory", TCL_STATIC);
        return TCL_ERROR;
    }
    Tcl_SetResult(interp, "test1", TCL_STATIC);
    return TCL_ERROR;

    for (i = 0; i < blt_vnum; i++) {
#ifdef THREADS
        mutex_lock(&vectors[i].mutex);
#endif
        V[i] = vectors[i].data[vectors[i].length-1];
#ifdef THREADS
        mutex_unlock(&vectors[i].mutex);
#endif
    }
    Blt_ResetVector(vec, V, blt_vnum, blt_vnum, TCL_VOLATILE);
    txfree(V);
    return TCL_OK;
}


/*agr1: spice variable name
 *arg2: index
 */
static int
get_value TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    char *var;
    int i, vindex, j;
    double val = 0;
    NG_IGNORE(clientData);
    if (argc != 3) {
        Tcl_SetResult(interp,
                      "Wrong # args. spice::get_value spice_variable index",
                      TCL_STATIC);
        return TCL_ERROR;
    }
    var = (char *)argv[1];

    for (i = 0; i < blt_vnum && strcmp(var, vectors[i].name); i++)
        ;

    if (i == blt_vnum) {
        Tcl_SetResult(interp, "Bad spice variable ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        return TCL_ERROR;
    } else {
        vindex = i;
    }

    j = atoi(argv[2]);

#ifdef THREADS
    mutex_lock(&vectors[vindex].mutex);
#endif

    if (j < 0 || j >= vectors[vindex].length) {
        i = 1;
    } else {
        i = 0;
        val = vectors[vindex].data[j];
    }

#ifdef THREADS
    mutex_unlock(&vectors[vindex].mutex);
#endif

    if (i) {
        Tcl_SetResult(interp, "Index out of range", TCL_STATIC);
        return TCL_ERROR;
    } else {
        Tcl_SetObjResult(interp, Tcl_NewDoubleObj(val));
        return TCL_OK;
    }
}


/*
 * arg1: spice vector name
 * arg2: real data blt vector
 * arg3: Optional: imaginary data blt vector*/
static int
vectoblt TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    Blt_Vector *real_BltVector, *imag_BltVector;
    char *realBlt, *imagBlt, *var;
    struct dvec *var_dvec;
    double *realData, *compData;
    int compIndex; //index to loop inside the vectors' data
    NG_IGNORE(clientData);

    if (argc < 3 || argc > 4) {
        Tcl_SetResult(interp, "Wrong # args. spice::vectoblt spice_variable real_bltVector [imag_bltVector]", TCL_STATIC);
        return TCL_ERROR;
    }

    real_BltVector = NULL;
    imag_BltVector = NULL;
    var = (char *)argv[1];
    var_dvec = vec_get(var);
    if (var_dvec == NULL) {
        Tcl_SetResult(interp, "Bad spice vector ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        return TCL_ERROR;
    }
    realBlt = (char *)argv[2];
    if (Blt_GetVector(interp, realBlt, &real_BltVector)) {
        Tcl_SetResult(interp, "Bad real blt vector ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)realBlt, TCL_STATIC);
        return TCL_ERROR;
    }
    if (argc == 4) {
        imagBlt = (char *)argv[3];
        if (Blt_GetVector(interp, imagBlt, &imag_BltVector)) {
            Tcl_SetResult(interp, "Bad imag blt vector ", TCL_STATIC);
            Tcl_AppendResult(interp, (char *)imagBlt, TCL_STATIC);
            return TCL_ERROR;
        }
    }
/*If data is complex, it is harder (more complex :) to export...*/
//  int compIndex; //index to loop inside the vectors' data
    if (var_dvec->v_realdata == NULL) {
        if (var_dvec->v_compdata == NULL) {
            Tcl_SetResult(interp, "The vector contains no data", TCL_STATIC);
            Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        }
        else {
            realData = TMALLOC(double, var_dvec->v_length);
            for (compIndex = 0; compIndex < var_dvec->v_length; compIndex++)
                realData[compIndex] = ((var_dvec->v_compdata+compIndex)->cx_real);

            Blt_ResetVector(real_BltVector, realData, var_dvec->v_length, var_dvec->v_length, TCL_VOLATILE);
            if (imag_BltVector != NULL) {
                compData = TMALLOC(double, var_dvec->v_length);
                for (compIndex = 0; compIndex < var_dvec->v_length; compIndex++)
                    compData[compIndex] = ((var_dvec->v_compdata+compIndex)->cx_imag);

                Blt_ResetVector(imag_BltVector, compData, var_dvec->v_length, var_dvec->v_length, TCL_VOLATILE);
            }
        }
    } else {
        Blt_ResetVector(real_BltVector, var_dvec->v_realdata, var_dvec->v_length, var_dvec->v_length, TCL_VOLATILE);
        if (imag_BltVector != NULL) {
            compData = TMALLOC(double, var_dvec->v_length);
            for (compIndex = 0; compIndex < var_dvec->v_length; compIndex++)
                compData[compIndex] = 0;

            Blt_ResetVector(imag_BltVector, compData, var_dvec->v_length, var_dvec->v_length, TCL_VOLATILE);
        }
    }

    Tcl_SetResult(interp, "finished!", TCL_STATIC);
    return TCL_OK;
}


/*agr1: spice variable name
 *arg2: blt_vector
 *arg3: start copy index, optional
 *arg4: end copy index. optional
 */
static int
spicetoblt TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    Blt_Vector *vec;
    int j, i;
    char *blt, *var;
    int start = 0, end = -1, len;

    NG_IGNORE(clientData);
    if (argc < 3 || argc > 5) {
        Tcl_SetResult(interp, "Wrong # args. spice::spicetoblt spice_variable vecName ?start? ?end?", TCL_STATIC);
        return TCL_ERROR;
    }

    var = (char *)argv[1];
    blt = (char *)argv[2];

    for (i = 0; i < blt_vnum && strcmp(var, vectors[i].name); i++)
        ;

    if (i == blt_vnum) {
        Tcl_SetResult(interp, "Bad spice variable ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        return TCL_ERROR;
    } else {
        j = i;
    }

    if (Blt_GetVector(interp, blt, &vec)) {
        Tcl_SetResult(interp, "Bad blt vector ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)blt, TCL_STATIC);
        return TCL_ERROR;
    }


    if (argc >= 4)
        start = atoi(argv[3]);
    if (argc == 5)
        end   = atoi(argv[4]);
    if (vectors[j].length) {
#ifdef THREADS
        mutex_lock(&vectors[j].mutex);
#endif

        len = vectors[j].length;

        if (start) {
            start = start % len;
            if (start < 0)
                start += len;
        }

        end = end % len;
        if (end < 0)
            end += len;

        len = abs(end - start + 1);

        Blt_ResetVector(vec, (vectors[j].data + start), len,
                        len, TCL_VOLATILE);

#ifdef THREADS
        mutex_unlock(&vectors[j].mutex);
#endif
    }
    return TCL_OK;
}


/******************************************************************/
/*     Main spice command executions and thread control           */
/*****************************************************************/

#ifdef THREADS
static threadId_t tid, bgtid = (threadId_t) 0;

static bool fl_running = FALSE;
static bool fl_exited = TRUE;

#if defined(__MINGW32__) || defined(_MSC_VER)
#define EXPORT_FLAVOR WINAPI
#else
#define EXPORT_FLAVOR
#endif

static void * EXPORT_FLAVOR
_thread_run(void *string)
{
    fl_exited = FALSE;
    bgtid = thread_self();
    cp_evloop((char *)string);
    FREE(string);
    bgtid = (threadId_t)0;
    fl_exited = TRUE;
    return NULL;
}


/*Stops a running thread, hopefully */
static int EXPORT_FLAVOR
_thread_stop(void)
{
    int timeout = 0;
    if (fl_running) {
        while (!fl_exited && timeout < 100) {
            ft_intrpt = TRUE;
            timeout++;
#if defined(__MINGW32__) || defined(_MSC_VER)
            Sleep(100); /* va: windows native */
#else
            usleep(10000);
#endif
        }
        if (!fl_exited) {
            fprintf(stderr, "Couldn't stop tclspice\n");
            return TCL_ERROR;
        }
#ifdef HAVE_LIBPTHREAD
        pthread_join(tid, NULL);
#endif
        fl_running = FALSE;
        ft_intrpt = FALSE;
        return TCL_OK;
    } else {
        fprintf(stderr, "Spice not running\n");
    }
    return TCL_OK;
}


void
sighandler_tclspice(int num)
{
    NG_IGNORE(num);
    if (fl_running)
        _thread_stop();
    return;
}

#endif /*THREADS*/


static int
_run(int argc, char **argv)
{
    char buf[1024] = "";
    int i;
    sighandler oldHandler;
#ifdef THREADS
    char *string;
    bool fl_bg = FALSE;
    /* run task in background if preceeded by "bg"*/
    if (!strcmp(argv[0], "bg")) {
        argc--;
        argv = &argv[1];
        fl_bg = TRUE;
    }
#endif


    /* Catch Ctrl-C to break simulations */
#ifndef _MSC_VER_
    oldHandler = signal(SIGINT, (SIGNAL_FUNCTION) ft_sigintr);
    if (SETJMP(jbuf, 1) != 0) {
        ft_sigintr_cleanup();
        signal(SIGINT, oldHandler);
        return TCL_OK;
    }
#else
    oldHandler = SIG_IGN;
#endif

    /*build a char * to pass to cp_evloop */
    for (i = 0; i < argc; i++) {
        strcat(buf, argv[i]);
        strcat(buf, " ");
    }

#ifdef THREADS
    /* run in the background */
    if (fl_bg) {
        if (fl_running)
            _thread_stop();
        fl_running = TRUE;
        string = copy(buf);     /*as buf gets freed fairly quickly*/
#ifdef HAVE_LIBPTHREAD
        pthread_create(&tid, NULL, _thread_run, (void *)string);
#else
        Tcl_CreateThread(&tid, (Tcl_ThreadCreateProc *)_thread_run, string,
                         TCL_THREAD_STACK_DEFAULT, TCL_THREAD_NOFLAGS);
#endif
    } else
        /* halt (pause) a bg run */
        if (!strcmp(argv[0], "halt")) {
            signal(SIGINT, oldHandler);
            return _thread_stop();
        } else
            /* backwards compatability with old command */
            if (!strcmp(argv[0], "stop"))
                if (argc > 1) {
                    cp_evloop(buf);
                } else {
                    _thread_stop();
                    cp_evloop(buf);
                }
            else {
                /* cannot do anything if spice is running in the bg*/
                if (fl_running) {
                    if (fl_exited) {
                        _thread_stop();
                        cp_evloop(buf);
                    } else {
                        fprintf(stderr, "type \"spice stop\" first\n");
                    }
                } else {
                    /*do the command*/
                    cp_evloop(buf);
                }
            }
#else
    cp_evloop(buf);
#endif /*THREADS*/
    signal(SIGINT, oldHandler);
    return TCL_OK;
}


static int
_tcl_dispatch TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    int i;
    NG_IGNORE(clientData);
    save_interp();
    char *prefix = strstr(argv[0], "spice::");
    if (prefix)
        argv[0] = prefix + 7;
    return _run(argc, (char **)argv);
}


/* Runs the spice command given in spice <cmd>*/
static int
_spice_dispatch TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    save_interp();
    if (argc == 1)
        return TCL_OK;
    return _run(argc-1, (char **)&argv[1]);
}


#ifdef THREADS
/*Checks if spice is runnuing in the background */
static int
running TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    NG_IGNORE(argc);
    NG_IGNORE(argv);
    Tcl_SetObjResult(interp, Tcl_NewIntObj((long) (fl_running && !fl_exited)));
    return TCL_OK;
}
#endif


/**************************************/
/*  plot manipulation functions       */
/*  only usefull if plots are saved   */
/**************************************/

/*Outputs the names of all variables in the plot */
static int
plot_variables TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    int plot;
    struct dvec *v;
    NG_IGNORE(clientData);

    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_variables plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot given", TCL_STATIC);
        return TCL_ERROR;
    }

    for (v = pl->pl_dvecs; v; v = v->v_next)
        Tcl_AppendElement(interp, v->v_name);

    return TCL_OK;
}


static int
plot_variablesInfo TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    int plot;
    struct dvec *v;
    char buf[256];
    char *name;
    int length;
    NG_IGNORE(clientData);

    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_variablesInfo plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot given", TCL_STATIC);
        return TCL_ERROR;
    }

    Tcl_ResetResult(interp);
    for (v = pl->pl_dvecs; v; v = v->v_next) {
        name = v->v_name;
        length = v->v_length;

        sprintf(buf, "{%s %s %i} ", name, ft_typenames(v->v_type), length);
        Tcl_AppendResult(interp, (char *)buf, TCL_STATIC);
    }
    return TCL_OK;
}


/*returns the value of a variable */
static int
plot_get_value TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    struct dvec *v;
    char *name;
    int plot, index;
    NG_IGNORE(clientData);

    if (argc != 4) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_get_value name plot index", TCL_STATIC);
        return TCL_ERROR;
    }

    name = (char *)argv[1];
    plot = atoi(argv[2]);
    index = atoi(argv[3]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }
    for (v = pl->pl_dvecs; v; v = v->v_next)
        if (!strcmp(v->v_name, name)) {
            if (index < v->v_length) {
                Tcl_SetObjResult(interp, Tcl_NewDoubleObj((double) v->v_realdata[index]));
                return TCL_OK;
            } else {
                Tcl_SetResult(interp, "Bad index", TCL_STATIC);
                return TCL_ERROR;
            }
        }

    Tcl_SetResult(interp, "variable not found", TCL_STATIC);
    return TCL_ERROR;
}


/*The length of the first vector in a plot*/
static int
plot_datapoints TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    struct dvec *v;
    int plot;
    NG_IGNORE(clientData);

    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_datapoints plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }

    v = pl->pl_dvecs;

    Tcl_SetObjResult(interp, Tcl_NewIntObj((long) v->v_length)); // could be very dangeous
    return TCL_OK;
}


/*These functions give you infomation about a plot*/

static int
plot_title TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    int plot;
    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_title plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }
    Tcl_SetObjResult(interp, Tcl_NewStringObj(pl->pl_title, -1));
    return TCL_OK;
}


static int
plot_date TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;

    int plot;
    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_date plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }
    Tcl_SetObjResult(interp, Tcl_NewStringObj(pl->pl_date, -1));
    return TCL_OK;
}


static int
plot_name TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    int plot;
    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_name plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }
    Tcl_SetObjResult(interp, Tcl_NewStringObj(pl->pl_name, -1));
    return TCL_OK;
}


static int
plot_typename TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    int plot;
    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_typename plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }
    Tcl_SetObjResult(interp, Tcl_NewStringObj(pl->pl_typename, -1));
    return TCL_OK;
}


/*number of variables in a plot*/

static int
plot_nvars TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    struct dvec *v;
    int plot;
    int i = 0;

    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_nvars plot", TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }
    for (v = pl->pl_dvecs; v; v = v->v_next)
        i++;
    Tcl_SetObjResult(interp, Tcl_NewIntObj((long) i));
    return TCL_OK;
}


static int
plot_defaultscale TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct plot *pl;
    int plot;

    NG_IGNORE(clientData);
    if (argc != 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::plot_defaultscale plot",
                      TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }

    if (pl->pl_scale)
        Tcl_SetObjResult(interp, Tcl_NewStringObj(pl->pl_scale->v_name, -1));
    return TCL_OK;
}


/*agr1: plot index
 *agr2: spice variable name
 *arg3: blt_vector
 *arg4: start copy index, optional
 *arg5: end copy index. optional
 */
static int
plot_getvector TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    Blt_Vector *vec;
    char *blt, *var;
    int start = 0, end = -1, len;
    int plot;
    struct dvec *v;
    struct plot *pl;

    NG_IGNORE(clientData);
    if (argc < 4 || argc > 6) {
        Tcl_SetResult(interp,
                      "Wrong # args. spice::plot_getvector plot spice_variable vecName ?start? ?end?",
                      TCL_STATIC);
        return TCL_ERROR;
    }

    plot = atoi(argv[1]);

    if (!(pl = get_plot_by_index(plot))) {
        Tcl_SetResult(interp, "Bad plot", TCL_STATIC);
        return TCL_ERROR;
    }

    var = (char *)argv[2];
    blt = (char *)argv[3];

    for (v = pl->pl_dvecs; v; v = v->v_next)
        if (!strcmp(v->v_name, var))
            break;

    if (v == NULL) {
        Tcl_SetResult(interp, "variable not found: ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        return TCL_ERROR;
    }

    if (Blt_GetVector(interp, blt, &vec)) {
        Tcl_SetResult(interp, "Bad blt vector ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)blt, TCL_STATIC);
        return TCL_ERROR;
    }

    if (argc >= 5)
        start = atoi(argv[4]);
    if (argc == 6)
        end   = atoi(argv[5]);
    if (v->v_length) {

        len = v->v_length;

        if (start) {
            start = start % len;
            if (start < 0)
                start += len;
        }

        end = end % len;
        if (end < 0)
            end += len;

        len = abs(end - start + 1);

        Blt_ResetVector(vec, (v->v_realdata + start), len,
                        len, TCL_VOLATILE);
    }
    return TCL_OK;
}


static int
plot_getplot TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    NG_IGNORE(argc);
    NG_IGNORE(argv);

    if (plot_cur)
        Tcl_SetObjResult(interp, Tcl_NewStringObj(plot_cur->pl_typename, -1));

    return TCL_OK;
}


/*******************************************/
/*           Misc functions                */
/*******************************************/

/*Runs a tcl script and returns the output*/
static int
get_output TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    FILE *pipein;
    int tmp_1, tmp_2 = 0;
    char buf[1024];
    int outfd, outfd2 = 0;
    NG_IGNORE(clientData);
    save_interp();
    if ((argc < 2) || (argc > 3)) {
        Tcl_SetResult(interp, "Wrong # args. spice::get_output script ?var_for_stderr?", TCL_STATIC);
        return TCL_ERROR;
    }
    tmp_1 = dup(1);
    outfd = open("/tmp/tclspice.tmp_out", O_WRONLY|O_CREAT|O_TRUNC, S_IRWXU);
    if (argc == 3) {
        tmp_2 = dup(2);
        outfd2 = open("/tmp/tclspice.tmp_err", O_WRONLY|O_CREAT|O_TRUNC, S_IRWXU);
    }
    freopen("/tmp/tclspice.tmp_out", "w", stdout);
    if (argc == 3)
        freopen("/tmp/tclspice.tmp_err", "w", stderr);
    dup2(outfd, 1);
    if (argc == 3)
        dup2(outfd2, 2);

    Tcl_Eval(interp, argv[1]);

    fclose(stdout);
    close(outfd);
    if (argc == 3) {
        fclose(stderr);
        close(outfd2);
    }
    dup2(tmp_1, 1);
    close(tmp_1);
    if (argc == 3) {
        dup2(tmp_2, 2);
        close(tmp_2);
    }
    freopen("/dev/fd/1", "w", stdout);
    if (argc == 3)
        freopen("/dev/fd/2", "w", stderr);
    pipein = fopen("/tmp/tclspice.tmp_out", "r");
    if (pipein == NULL)
        fprintf(stderr, "pipein==NULL\n");

    Tcl_ResetResult(interp);
    while (fgets(buf, 1024, pipein) != NULL)
        Tcl_AppendResult(interp, (char *)buf, TCL_STATIC);

    fclose(pipein);
    if (argc == 3) {
        pipein = fopen("/tmp/tclspice.tmp_err", "r");
        Tcl_SetVar(interp, argv[2], "", 0);
        while (fgets(buf, 1024, pipein) != NULL)
            Tcl_SetVar(interp, argv[2], buf, TCL_APPEND_VALUE);

        fclose(pipein);
    }
    return TCL_OK;
}


/* Returns the current value of a parameter
 * has lots of memory leaks
 */
static int
get_param TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    wordlist *wl = NULL;
    char *device, *param;
    struct variable *v;
    char buf[128];
    NG_IGNORE(clientData);
    if (argc != 3) {
        Tcl_SetResult(interp, "Wrong # args. spice::get_param device param", TCL_STATIC);
        return TCL_ERROR;
    }
    if (!ft_curckt) {
        Tcl_SetResult(interp, "No circuit loaded ", TCL_STATIC);
        return TCL_ERROR;
    }

    device = (char *)argv[1];
    param  = (char *)argv[2];
    /* copied from old_show(wordlist *) */
    v = if_getparam (ft_curckt->ci_ckt, &device, param, 0, 0);
    if (!v)
        v = if_getparam (ft_curckt->ci_ckt, &device, param, 0, 1);
    if (v) {
        wl = cp_varwl(v);
        Tcl_SetResult(interp, wl->wl_word, TCL_VOLATILE);
        wl_free(wl);
        tfree(v);
        return TCL_OK;

    } else {
        sprintf(buf, "%s in %s not found", param, device);
        Tcl_AppendResult(interp, buf, TCL_STATIC);
    }
    return TCL_ERROR;
}


/* va - added
   call:    s. errormessage
   returns: param == all: list of all model parameters of device/model
   param == name: description of given model parameter
*/
int
get_mod_param TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    char *name;
    char *paramname;
    GENinstance *devptr = NULL;
    GENmodel *modptr = NULL;
    IFdevice *device;
    IFparm *opt;
    IFvalue pv;
    int i, err, typecode = -1;
    char buf[128];
    bool found;

    NG_IGNORE(clientData);
    if (argc < 2 || argc >3) {
        Tcl_SetResult(interp,
                      "Wrong # args. spice::get_mod_param device|model [all|param]", TCL_STATIC);
        return TCL_ERROR;
    }
    if (ft_curckt == NULL) {
        Tcl_SetResult(interp, "No circuit loaded ", TCL_STATIC);
        return TCL_ERROR;
    }

    name = (char *)argv[1];
    if (argc > 2)
        paramname = (char *)argv[2];
    else
        paramname = "all";

    if (name == NULL || name[0] == '\0') {
        Tcl_SetResult(interp, "No model or device name provided.", TCL_STATIC);
        return TCL_ERROR;
    }

    /* get the unique IFuid for name (device/model) */
    INPretrieve(&name, ft_curckt->ci_symtab);
    devptr = ft_sim->findInstance (ft_curckt->ci_ckt, name);
    if (devptr) {
        typecode = devptr->GENmodPtr->GENmodType;
    } else  {
        modptr = ft_sim->findModel (ft_curckt->ci_ckt, name);
        if (modptr) {
            typecode = modptr->GENmodType;
        } else {
            sprintf(buf, "No such device or model name %s", name);
            Tcl_SetResult(interp, buf, TCL_VOLATILE);
            return TCL_ERROR;
        }
    }
    device = ft_sim->devices[typecode];
    found = FALSE;
    for (i = 0; i < *(device->numModelParms); i++) {
        opt = &device->modelParms[i];
        if (opt->dataType != (IF_SET|IF_ASK|IF_REAL))
            continue; /* only real IO-parameter */
        if (strcmp(paramname, "all") == 0) {
            Tcl_AppendElement(interp, opt->keyword);
            found = TRUE;
        } else if (strcmp(paramname, opt->keyword) == 0) {
            if (devptr)
                err = ft_sim->askInstanceQuest (ft_curckt->ci_ckt, devptr,
                                                opt->id, &pv, NULL);
            else
                err = ft_sim->askModelQuest (ft_curckt->ci_ckt, modptr,
                                             opt->id, &pv, NULL);
            if (err == OK) {
                sprintf(buf, "%g", pv.rValue); /* dataType is here always real */
                Tcl_SetResult(interp, buf, TCL_VOLATILE);
                return TCL_OK;
            }
        }
    }
    if (found != TRUE) {
        sprintf(buf, "unknown parameter %s", paramname);
        Tcl_SetResult(interp, buf, TCL_VOLATILE);
        return TCL_ERROR;
    }
    return TCL_OK;
}


/* Direct control over the step size
 * Spice will still adjust it to keep accuracy wuithin reltol and abstol
 */
static int
delta TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    if (argc < 1 ||argc > 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::delta ?value?", TCL_STATIC);
        return TCL_ERROR;
    }
    if (ft_curckt == NULL) {
        Tcl_SetResult(interp, "No circuit loaded ", TCL_STATIC);
        return TCL_ERROR;
    }

    if (argc == 2)
        (ft_curckt->ci_ckt)->CKTdelta = atof(argv[1]);

    Tcl_SetObjResult(interp, Tcl_NewDoubleObj((ft_curckt->ci_ckt)->CKTdelta));
    return TCL_OK;
}


#include "ngspice/trandefs.h"
/* Direct control over the maximum stepsize
 * Spice will still adjust it to keep accuracy wuithin reltol and abstol
 */
static int
maxstep TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    TRANan *job;
    NG_IGNORE(clientData);
    if (argc < 1 ||argc > 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::maxstep ?value?", TCL_STATIC);
        return TCL_ERROR;
    }
    if (ft_curckt == NULL) {
        Tcl_SetResult(interp, "No circuit loaded ", TCL_STATIC);
        return TCL_ERROR;
    }

    job = (TRANan*)(ft_curckt->ci_ckt)->CKTcurJob;
    if (argc == 2)
        job->TRANmaxStep = atof(argv[1]);

    Tcl_SetObjResult(interp, Tcl_NewDoubleObj(job->TRANmaxStep));
    return TCL_OK;
}


/* obtain the initial time of a transient analysis */
static int
get_initTime TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    TRANan *job;
    NG_IGNORE(argv);
    NG_IGNORE(clientData);

    if (argc < 1 ||argc > 1) {
        Tcl_SetResult(interp, "Wrong # args. spice::get_initTime", TCL_STATIC);
        return TCL_ERROR;
    }
    if (ft_curckt == NULL) {
        Tcl_SetResult(interp, "No circuit loaded ", TCL_STATIC);
        return TCL_ERROR;
    }

    job = (TRANan*)(ft_curckt->ci_ckt)->CKTcurJob;
    Tcl_SetObjResult(interp, Tcl_NewDoubleObj(job->TRANinitTime));
    return TCL_OK;
}


/* obtain the final time of a transient analysis */
static int
get_finalTime TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    TRANan *job;
    NG_IGNORE(argv);
    NG_IGNORE(clientData);

    if (argc < 1 ||argc > 1) {
        Tcl_SetResult(interp, "Wrong # args. spice::get_finalTime", TCL_STATIC);
        return TCL_ERROR;
    }
    if (ft_curckt == NULL) {
        Tcl_SetResult(interp, "No circuit loaded ", TCL_STATIC);
        return TCL_ERROR;
    }

    job = (TRANan*)(ft_curckt->ci_ckt)->CKTcurJob;
    Tcl_SetObjResult(interp, Tcl_NewDoubleObj(job->TRANfinalTime));
    return TCL_OK;
}


/****************************************/
/*          The Tk frontend for plot    */
/****************************************/

/* Use Tcl_GetStringResult to get canvas size etc. from Tcl */
#include "ngspice/ftedev.h"

int
sp_Tk_Init(void)
{
    /* This is hard coded in C at the mo, use X11 values */
    dispdev->numlinestyles = 8;
    dispdev->numcolors = 20;
    dispdev->width = 1280;
    dispdev->height = 1024;

    return 0;
}


#include "ngspice/graph.h"

int
sp_Tk_NewViewport(GRAPH *graph)
{
    const char *result;
    int width, height, fontwidth, fontheight;
    graph->devdep = NULL;

    if (Tcl_GlobalEval(spice_interp, "spice_gr_NewViewport") != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }

    result = Tcl_GetStringResult(spice_interp);
    if (sscanf(result, "%i %i %i %i", &width, &height, &fontwidth, &fontheight) != 4) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    graph->absolute.xpos = 0; /* these always seem sensible, let Tcl adjust coods */
    graph->absolute.ypos = 0;
    graph->absolute.width = width;
    graph->absolute.height = height;
    graph->fontwidth = fontwidth;
    graph->fontheight = fontheight;
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_Close(void)
{
    if (Tcl_Eval(spice_interp, "spice_gr_Close") != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_Clear(void)
{
    if (Tcl_Eval(spice_interp, "spice_gr_Clear") != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_DrawLine(int x1, int y1, int x2, int y2)
{
    char buf[1024];
    sprintf(buf, "spice_gr_DrawLine %i %i %i %i", x1, y1, x2, y2);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_Arc(int x0, int y0, int radius, double theta, double delta_theta)
{
    char buf[1024];
    sprintf(buf, "spice_gr_Arc %i %i %i %f %f", x0, y0, radius, theta, delta_theta);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_Text(char *text, int x, int y, int angle)
{
    char buf[1024];
    NG_IGNORE(angle);
    sprintf(buf, "spice_gr_Text \"%s\" %i %i", text, x, y);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_DefineColor(int colorid, double red, double green, double blue)
{
    char buf[1024];
    sprintf(buf, "spice_gr_DefineColor %i %g %g %g", colorid, red, green, blue);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_DefineLinestyle(int linestyleid, int mask)
{
    char buf[1024];
    sprintf(buf, "spice_gr_DefineLinestyle %i %i", linestyleid, mask);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_SetLinestyle(int linestyleid)
{
    char buf[1024];
    sprintf(buf, "spice_gr_SetLinestyle %i", linestyleid);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_SetColor(int colorid)
{
    char buf[1024];
    sprintf(buf, "spice_gr_SetColor %i", colorid);
    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


int
sp_Tk_Update(void)
{
    if (Tcl_Eval(spice_interp, "spice_gr_Update") != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }
    Tcl_ResetResult(spice_interp);
    return 0;
}


/********************************************************/
/*           The Blt method for plotting                */
/********************************************************/

static void
dvecToBlt(Blt_Vector *Data, struct dvec *x)
{
    if (x->v_flags & VF_REAL) {
        Blt_ResetVector (Data, x->v_realdata , x->v_length,
                         x->v_length, TCL_VOLATILE);
    } else {
        double *data;
        int i;

        data = TMALLOC(double, x->v_length);

        for (i = 0; i < x->v_length; i++)
            data[i] = realpart(x->v_compdata[i]);

        Blt_ResetVector (Data, data, x->v_length, x->v_length, TCL_VOLATILE);

        tfree(data);
    }

    return;
}


static void
escape_brackets(char *string)
{
    int printed = strlen(string), i;

    for (i = 0; i < printed; i++) {
        if (string[i] == ']' || string[i] == '[') {
            int j;
            for (j = printed; j >= i; j--)
                string[j+3] = string[j];

            string[i] = '\\';
            string[i+1] = '\\';
            string[i+2] = '\\';
            i += 3;
            printed += 3;
        }
    }
    return;
}


int
blt_plot(struct dvec *y, struct dvec *x, int new)
{
    static int ctr = -1;
    Blt_Vector *X_Data = NULL, *Y_Data = NULL;
    char buf[1024];

    /* A bug in these functions? , crashes if used so make vectors in Tcl
       Blt_CreateVector(spice_interp, "::spice::X_Data", 1, &X_Data);
       Blt_CreateVector(spice_interp, "::spice::Y_Data", 1, &Y_Data);
    */
    Blt_GetVector(spice_interp, "::spice::X_Data", &X_Data);
    Blt_GetVector(spice_interp, "::spice::Y_Data", &Y_Data);

    if (!X_Data || !Y_Data) {
        fprintf(stderr, "Error: Blt vector X_Data or Y_Data not created\n");
        return 1;
    }

    dvecToBlt(X_Data, x);
    dvecToBlt(Y_Data, y);

    if (new)
        ctr++;

    sprintf(buf, "spice_gr_Plot %s %s %s %s %s %s %d",
            x->v_name, ft_typenames(x->v_type), ft_typabbrev(x->v_type),
            y->v_name, ft_typenames(y->v_type), ft_typabbrev(y->v_type), ctr);
    escape_brackets(buf);

    if (Tcl_Eval(spice_interp, buf) != TCL_OK) {
        Tcl_ResetResult(spice_interp);
        return 1;
    }

    Tcl_ResetResult(spice_interp);

    return 0;
}


/********************************************************/
/*             Triggering stuff                         */
/********************************************************/

struct triggerEvent {
    struct triggerEvent *next;
    int vector;
    int type;
    int stepNumber;
    double time;
    double voltage;
    char ident[16];
};


struct triggerEvent *eventQueue;
struct triggerEvent *eventQueueEnd;

#ifdef THREADS
mutexType triggerMutex;
#endif

struct watch {
    struct watch *next;
    char name[16];
    int vector;                 /* index of vector to watch */
    int type;                   /* +ive or -ive trigger */
    int state;                  /* pretriggered or not */
    double Vmin;                /* the boundaries */
    double Vmax;
    /* To get the exact trigger time */
    double Vavg;
    double oT;
    double oV;
};

struct watch *watches;

char *triggerCallback;
unsigned int triggerPollTime = 1;

char *stepCallback;
unsigned int stepPollTime = 1;
unsigned int stepCount = 1;
int stepCallbackPending;


void
stepEventSetup(ClientData clientData, int flags)
{
    Tcl_Time t;
    NG_IGNORE(clientData);
    NG_IGNORE(flags);
    if (stepCallbackPending) {
        t.sec  = 0;
        t.usec = 0;
    } else {
        t.sec = stepPollTime / 1000;
        t.usec = (stepPollTime % 1000) * 1000;
    }
    Tcl_SetMaxBlockTime(&t);
}


int
stepEventHandler(Tcl_Event *evPtr, int flags)
{
    NG_IGNORE(evPtr);
    NG_IGNORE(flags);
    if (stepCallbackPending) {
        stepCallbackPending = 0;
        Tcl_Preserve((ClientData)spice_interp);
        Tcl_Eval(spice_interp, stepCallback);
        Tcl_ResetResult(spice_interp);
        Tcl_Release((ClientData)spice_interp);
    }
    return TCL_OK;
}


void
stepEventCheck(ClientData clientData, int flags)
{
    NG_IGNORE(clientData);
    NG_IGNORE(flags);
    if (stepCallbackPending) {
        Tcl_Event *tclEvent;
        tclEvent = (Tcl_Event *) ckalloc(sizeof(Tcl_Event));
        tclEvent->proc = stepEventHandler;
        Tcl_QueueEvent(tclEvent, TCL_QUEUE_TAIL);
    }
}


int
triggerEventHandler(Tcl_Event *evPtr, int flags)
{
    static char buf[512];
    int rtn = TCL_OK;
    NG_IGNORE(evPtr);
    NG_IGNORE(flags);
    Tcl_Preserve((ClientData)spice_interp);
#ifdef THREADS
    mutex_lock(&triggerMutex);
#endif
    while (eventQueue) {
        struct triggerEvent *event = eventQueue;
        eventQueue = eventQueue->next;


        snprintf(buf, 512, "%s %s %g %d %d %g %s", triggerCallback, vectors[event->vector].name,
                 event->time, event->stepNumber, event->type, event->voltage, event->ident);

        rtn = Tcl_Eval(spice_interp, buf);
        FREE(event);
        if (rtn)
            goto quit;
    }
    eventQueueEnd = NULL;
quit:
#ifdef THREADS
    mutex_unlock(&triggerMutex);
#endif
    Tcl_ResetResult(spice_interp);
    Tcl_Release((ClientData)spice_interp);
    return TCL_OK;
}


void
triggerEventSetup(ClientData clientData, int flags)
{
    Tcl_Time t;
    NG_IGNORE(clientData);
    NG_IGNORE(flags);
    if (eventQueue) {
        t.sec  = 0;
        t.usec = 0;
    } else {
        t.sec = triggerPollTime / 1000;
        t.usec = (triggerPollTime % 1000) * 1000;
    }
    Tcl_SetMaxBlockTime(&t);
}


void
triggerEventCheck(ClientData clientData, int flags)
{
    NG_IGNORE(clientData);
    NG_IGNORE(flags);
#ifdef THREADS
    mutex_lock(&triggerMutex);
#endif
    if (eventQueue) {
        Tcl_Event *tclEvent;
        tclEvent = (Tcl_Event *) ckalloc(sizeof(Tcl_Event));
        tclEvent->proc = triggerEventHandler;
        Tcl_QueueEvent(tclEvent, TCL_QUEUE_TAIL);
    }
#ifdef THREADS
    mutex_unlock(&triggerMutex);
#endif
}


int
Tcl_ExecutePerLoop(void)
{
    struct watch *current;

#ifdef THREADS
    mutex_lock(&vectors[0].mutex);
    mutex_lock(&triggerMutex);
#endif

    for (current = watches; current; current = current->next) {
        vector *v;
        v = &vectors[current->vector];
#ifdef THREADS
        mutex_lock(&v->mutex);
#endif

        if ((current->type > 0 && current->state && v->data[v->length-1] > current->Vmax) ||
            (current->type < 0 && current->state && v->data[v->length-1] < current->Vmin)) {
            struct triggerEvent *tmp = TMALLOC(struct triggerEvent, 1);

            tmp->next = NULL;

            if (eventQueue) {
                eventQueueEnd->next = tmp;
                eventQueueEnd = tmp;
            } else {
                eventQueue = tmp;
            }

            eventQueueEnd = tmp;

            tmp->vector = current->vector;
            tmp->type = current->type;
            tmp->stepNumber = vectors[0].length;

            {
                double T = vectors[0].data[vectors[0].length-1];
                double V = v->data[v->length-1];

                tmp->time = current->oT +
                    (current->Vavg - current->oV) * (T - current->oT) / (V - current->oV);
                tmp->voltage = current->Vavg;
            }

            strcpy(tmp->ident, current->name);
            current->state = 0;

        } else
            if ((current->type > 0 && v->data[v->length-1] < current->Vmin) ||
                (current->type < 0 && v->data[v->length-1] > current->Vmax))
                current->state = 1;

        current->oT = vectors[0].data[vectors[0].length-1];
        current->oV = v->data[v->length-1];

#ifdef THREADS
        mutex_unlock(&v->mutex);
#endif
    }

    if (stepCallback && vectors[0].length % stepCount == 0)
        stepCallbackPending = 1;

#ifdef THREADS
    mutex_unlock(&triggerMutex);

    mutex_unlock(&vectors[0].mutex);

    if (triggerCallback && eventQueue && bgtid != thread_self())
        triggerEventHandler(NULL, 0);

    if (stepCallback && stepCallbackPending && bgtid != thread_self())
        stepEventHandler(NULL, 0);
#else
    if (triggerCallback && eventQueue)
        triggerEventHandler(NULL, 0);

    if (stepCallback && stepCallbackPending)
        triggerEventHandler(NULL, 0);
#endif

    return 0;
}


static int
resetTriggers(void)
{
#ifdef THREADS
    mutex_lock(&triggerMutex);
#endif

    while (watches) {
        struct watch *tmp = watches;
        watches = tmp->next;
        FREE(tmp);
    }

    while (eventQueue) {
        struct triggerEvent *tmp = eventQueue;
        eventQueue = tmp->next;
        FREE(tmp);
    }

    eventQueueEnd = NULL;

#ifdef THREADS
    mutex_unlock(&triggerMutex);
#endif
    return 0;
}


/* Registers a watch for a trigger
 *arg0: Callback function (optional - none removes callback)
 *arg1: Poll interval usec (optional - defaults to 500000 )
 */
static int
registerTriggerCallback TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    if (argc > 3) {
        Tcl_SetResult(interp,
                      "Wrong # args. spice::registerTriggerCallback ?proc? ?ms?",
                      TCL_STATIC);
        return TCL_ERROR;
    }

    if (triggerCallback) {
        Tcl_DeleteEventSource(triggerEventSetup, triggerEventCheck, NULL);
        free(triggerCallback);
        triggerCallback = NULL;
    }

    if (argc == 1)
        return TCL_OK;

    triggerCallback = strdup(argv[1]);
    Tcl_CreateEventSource(triggerEventSetup, triggerEventCheck, NULL);

    if (argc == 3) {
        triggerPollTime = atoi(argv[2]);
        if (triggerPollTime == 0)
            triggerPollTime = 500;
    }

    return TCL_OK;
}


/* Registers step counter callback
 *arg0: Callback function (optional - none removes callback)
 *arg1: Number of steps per Callback
 *arg2: Poll interval usec (optional - defaults to 500000 )
 */
static int
registerStepCallback TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    if (argc > 4) {
        Tcl_SetResult(interp,
                      "Wrong # args. spice::registerStepCallback ?proc? ?steps? ?ms?",
                      TCL_STATIC);
        return TCL_ERROR;
    }

    if (stepCallback) {
        Tcl_DeleteEventSource(stepEventSetup, stepEventCheck, NULL);
        free(stepCallback);
        stepCallback = NULL;
    }

    if (argc == 1)
        return TCL_OK;

    stepCallback = strdup(argv[1]);
    Tcl_CreateEventSource(stepEventSetup, stepEventCheck, NULL);

    if (argc >= 3) {
        stepCount = atoi(argv[2]);
        if (stepCount == 0)
            stepCount = 1;
    }

    if (argc == 4) {
        stepPollTime = atoi(argv[3]);
        if (stepPollTime == 0)
            stepPollTime = 50;
    }

    return TCL_OK;
}


/* Registers a watch for a trigger
 *arg0: Vector Name to watch
 *arg1: Vmin
 *arg2: Vmax
 *arg3: 1 / -1 for +ive(voltage goes +ive) or -ive trigger
 */
static int
registerTrigger TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    int i, index;
    const char *var;
    char ident[16];
    struct watch *tmp;
    int type;
    double Vavg, Vmin, Vmax;

    NG_IGNORE(clientData);
    if (argc < 4 && argc > 6) {
        Tcl_SetResult(interp, "Wrong # args. spice::registerTrigger vecName Vmin Vmax ?type? ?string?", TCL_STATIC);
        return TCL_ERROR;
    }

    var = argv[1];

    for (i = 0; i < blt_vnum && strcmp(var, vectors[i].name); i++)
        ;

    if (i == blt_vnum) {
        Tcl_SetResult(interp, "Bad spice variable ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        return TCL_ERROR;
    } else {
        index = i;
    }

    if (argc >= 5)
        type = atoi(argv[4]);
    else
        type = 1;

    if (argc >= 6) {
        strncpy(ident, argv[5], sizeof(ident));
        ident[sizeof(ident)-1] = '\0';
    } else {
        ident[0] = '\0';
    }

    Vmin = atof(argv[2]);
    Vmax = atof(argv[3]);
    Vavg = (Vmin + Vmax) / 2 ;

#ifdef THREADS
    mutex_lock(&triggerMutex);
#endif

    for (tmp = watches; tmp != NULL; tmp = tmp->next)
        if (ident[0] != '\0') {
            if (strcmp(ident, tmp->name) == 0) {
                watches->vector = index;
                watches->type = type;
                strcpy(watches->name, ident);

                watches->state = 0;
                watches->Vmin = Vmin;
                watches->Vmax = Vmax;
                watches->Vavg = Vavg;

                break;
            }
        } else {
            if (tmp->vector == index && tmp->type == type
                && tmp->Vavg == Vavg) {
                tmp->Vmin = Vmin;
                tmp->Vmax = Vmax;
                break;
            }
        }

    if (tmp == NULL) {

        tmp = TMALLOC(struct watch, 1);
        tmp->next = watches;
        watches = tmp;

        watches->vector = index;
        watches->type = type;
        strcpy(watches->name, ident);

        watches->state = 0;
        watches->Vmin = Vmin;
        watches->Vmax = Vmax;
        watches->Vavg = Vavg;

    }

#ifdef THREADS
    mutex_unlock(&triggerMutex);
#endif

    return TCL_OK;
}


/*unregisters a trigger
 *arg0: Vector name
 *arg1: type
 */
static int
unregisterTrigger TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    int i, index, type;
    char *var;
    struct watch *tmp;
    struct watch **cut;

    NG_IGNORE(clientData);
    if (argc != 2 && argc != 3) {
        Tcl_SetResult(interp, "Wrong # args. spice::unregisterTrigger vecName ?type?", TCL_STATIC);
        return TCL_ERROR;
    }

    var = (char *)argv[1];

    for (i = 0; i < blt_vnum && strcmp(var, vectors[i].name); i++)
        ;

    if (i == blt_vnum)
        index = -1;
    else
        index = i;

    if (argc == 3)
        type = atoi(argv[4]);
    else
        type = 1;

#ifdef THREADS
    mutex_lock(&triggerMutex);
#endif

    cut = &watches;

    tmp = watches;

    while (tmp)
        if ((tmp->vector == index && tmp->type == type) || strcmp(var, tmp->name) == 0) {
            *cut = tmp->next;
            txfree(tmp);
            break;
        } else {
            cut = &tmp->next;
            tmp = tmp->next;
        }

#ifdef THREADS
    mutex_unlock(&triggerMutex);
#endif

    if (tmp == NULL) {
        Tcl_SetResult(interp, "Could not find trigger ", TCL_STATIC);
        Tcl_AppendResult(interp, (char *)var, TCL_STATIC);
        return TCL_ERROR;
    }

    return TCL_OK;
}


/* returns:
   "vecName" "time" "stepNumber" "type"
*/
static int
popTriggerEvent TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    NG_IGNORE(clientData);
    NG_IGNORE(argv);
    if (argc != 1) {
        Tcl_SetResult(interp, "Wrong # args. spice::popTriggerEvent", TCL_STATIC);
        return TCL_ERROR;
    }

    if (eventQueue) {
        struct triggerEvent *popedEvent;
        Tcl_Obj *list;

#ifdef THREADS
        mutex_lock(&triggerMutex);
#endif

        popedEvent = eventQueue;

        eventQueue = popedEvent->next;
        if (eventQueue == NULL)
            eventQueueEnd = NULL;

        list = Tcl_NewListObj(0, NULL);

        Tcl_ListObjAppendElement(interp, list, Tcl_NewStringObj(vectors[popedEvent->vector].name, strlen(vectors[popedEvent->vector].name)));

        Tcl_ListObjAppendElement(interp, list, Tcl_NewDoubleObj(popedEvent->time));

        Tcl_ListObjAppendElement(interp, list, Tcl_NewIntObj(popedEvent->stepNumber));

        Tcl_ListObjAppendElement(interp, list, Tcl_NewIntObj(popedEvent->type));

        Tcl_ListObjAppendElement(interp, list, Tcl_NewDoubleObj(popedEvent->voltage));

        Tcl_ListObjAppendElement(interp, list, Tcl_NewStringObj(popedEvent->ident, strlen(popedEvent->ident)));

        Tcl_SetObjResult(interp, list);

        FREE(popedEvent);

#ifdef THREADS
        mutex_unlock(&triggerMutex);
#endif
    }

    return TCL_OK;
}


static int
listTriggers TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    struct watch *tmp;
    Tcl_Obj *list;

    NG_IGNORE(clientData);
    NG_IGNORE(argv);
    if (argc != 1) {
        Tcl_SetResult(interp, "Wrong # args. spice::listTriggers", TCL_STATIC);
        return TCL_ERROR;
    }

    list = Tcl_NewListObj(0, NULL);

#ifdef THREADS
    mutex_lock(&triggerMutex);
#endif

    for (tmp = watches; tmp; tmp = tmp->next)
        Tcl_ListObjAppendElement(interp, list, Tcl_NewStringObj(vectors[tmp->vector].name, strlen(vectors[tmp->vector].name)));

#ifdef THREADS
    mutex_unlock(&triggerMutex);
#endif

    Tcl_SetObjResult(interp, list);

    return TCL_OK;
}


static int
tmeasure TCL_CMDPROCARGS(clientData, interp, argc, argv)
{
    wordlist *wl = NULL;
    double mvalue;

    NG_IGNORE(clientData);
    if (argc <= 2) {
        Tcl_SetResult(interp, "Wrong # args. spice::listTriggers", TCL_STATIC);
        return TCL_ERROR;
    }

    wl = wl_build((char **)argv);

    get_measure2(wl, &mvalue, NULL, FALSE);

    printf(" %e \n", mvalue);

    Tcl_ResetResult(spice_interp);

    Tcl_SetObjResult(interp, Tcl_NewDoubleObj((double) mvalue));

    return TCL_OK;
}


/*******************************************************/
/*  Initialise spice and setup native methods          */
/*******************************************************/

#if defined(__MINGW32__) || defined(_MSC_VER)
__declspec(dllexport)
#endif
int
Spice_Init(Tcl_Interp *interp)
{
    if (interp == 0)
        return TCL_ERROR;

#ifdef USE_TCL_STUBS
    if (Tcl_InitStubs(interp, (char *)"8.1", 0) == NULL)
        return TCL_ERROR;
#endif

    Tcl_PkgProvide(interp, (char*) TCLSPICE_name, (char*) TCLSPICE_version);

    Tcl_Eval(interp, "namespace eval " TCLSPICE_namespace " { }");

    save_interp();

    {
        int i;
        char *key;
        Tcl_CmdInfo infoPtr;
        char buf[256];
        sighandler old_sigint;

        ft_rawfile = NULL;
        ivars(NULL);

        cp_in = stdin;
        cp_out = stdout;
        cp_err = stderr;

        /*timer*/
        init_time();

        /*IFsimulator struct initilised*/
        SIMinit(&nutmeginfo, &ft_sim);

        /* program name*/
        cp_program = ft_sim->simulator;

        srand((unsigned int) getpid());
        TausSeed();

        /*parameter fetcher, used in show*/
        if_getparam = spif_getparam_special;

        /* Get startup system limits */
        init_rlimits();

        /*Command prompt stuff */
        ft_cpinit();


        /* Read the user config files */
        /* To catch interrupts during .spiceinit... */
        old_sigint = signal(SIGINT, (SIGNAL_FUNCTION) ft_sigintr);
        if (SETJMP(jbuf, 1) == 1) {
            ft_sigintr_cleanup();
            fprintf(cp_err, "Warning: error executing .spiceinit.\n");
            goto bot;
        }

#ifdef HAVE_PWD_H
        /* Try to source either .spiceinit or ~/.spiceinit. */
        if (access(".spiceinit", 0) == 0) {
            inp_source(".spiceinit");
        } else {
            char *s;
            struct passwd *pw;
            pw = getpwuid(getuid());

            s = tprintf("%s" DIR_PATHSEP "%s", pw->pw_dir, INITSTR);

            if (access(s, 0) == 0)
                inp_source(s);
        }
#else /* ~ HAVE_PWD_H */
        {
            FILE *fp;
            /* Try to source the file "spice.rc" in the current directory.  */
            if ((fp = fopen("spice.rc", "r")) != NULL) {
                (void) fclose(fp);
                inp_source("spice.rc");
            }
        }
#endif /* ~ HAVE_PWD_H */
    bot:
        signal(SIGINT, old_sigint);

        /* initilise Tk display */
        DevInit();

        /* init the mutex */
#ifdef HAVE_LIBPTHREAD
        pthread_mutex_init(&triggerMutex, NULL);
#endif
#ifdef THREADS
        signal(SIGINT, sighandler_tclspice);
#endif

        /*register functions*/
        for (i = 0; (key = cp_coms[i].co_comname); i++) {
            sprintf(buf, "%s%s", TCLSPICE_prefix, key);
            if (Tcl_GetCommandInfo(interp, buf, &infoPtr) != 0)
                printf("Command '%s' can not be registered!\n", buf);
            else
                Tcl_CreateCommand(interp, buf, _tcl_dispatch, NULL, NULL);
        }

        Tcl_CreateCommand(interp, TCLSPICE_prefix "spice_header", spice_header, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "spice_data", spice_data, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "spicetoblt", spicetoblt, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "vectoblt", vectoblt, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "lastVector", lastVector, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "get_value", get_value, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "spice", _spice_dispatch, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "get_output", get_output, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "get_param", get_param, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "get_mod_param", get_mod_param, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "delta", delta, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "maxstep", maxstep, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "get_initTime", get_initTime, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "get_finalTime", get_finalTime, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_variables", plot_variables, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_variablesInfo", plot_variablesInfo, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_get_value", plot_get_value, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_datapoints", plot_datapoints, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_title", plot_title, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_date", plot_date, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_name", plot_name, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_typename", plot_typename, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_nvars", plot_nvars, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_defaultscale", plot_defaultscale, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "plot_getvector", plot_getvector, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "getplot", plot_getplot, NULL, NULL);

        Tcl_CreateCommand(interp, TCLSPICE_prefix "registerTrigger", registerTrigger, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "registerTriggerCallback", registerTriggerCallback, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "popTriggerEvent", popTriggerEvent, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "unregisterTrigger", unregisterTrigger, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "listTriggers", listTriggers, NULL, NULL);

        Tcl_CreateCommand(interp, TCLSPICE_prefix "registerStepCallback", registerTriggerCallback, NULL, NULL);
#ifdef THREADS
        Tcl_CreateCommand(interp, TCLSPICE_prefix "bg", _tcl_dispatch, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "halt", _tcl_dispatch, NULL, NULL);
        Tcl_CreateCommand(interp, TCLSPICE_prefix "running", running, NULL, NULL);
#endif

        Tcl_CreateCommand(interp, TCLSPICE_prefix "tmeasure", tmeasure, NULL, NULL);

        Tcl_CreateCommand(interp, TCLSPICE_prefix "registerStepCallback", registerStepCallback, NULL, NULL);

        Tcl_LinkVar(interp, TCLSPICE_prefix "steps_completed", (char *)&steps_completed, TCL_LINK_READ_ONLY|TCL_LINK_INT);
        Tcl_LinkVar(interp, TCLSPICE_prefix "blt_vnum", (char *)&blt_vnum, TCL_LINK_READ_ONLY|TCL_LINK_INT);
    }
    return TCL_OK;
}


/***************************************/
/* printf wrappers to redirect to puts */
/***************************************/

/* Contributed by Tim Edwards (tim@stravinsky.jhuapl.edu), 2003 */


/*------------------------------------------------------*/
/* Redefine the vfprintf() functions for use with tkcon */
/*------------------------------------------------------*/

int
tcl_vfprintf(FILE *f, const char *fmt, va_list args)
{
    char buf[1024];
    char *p, *s;
    int size, nchars, escapes, result;

    const char * const escape_chars = "$[]\"\\";

    const char * const prolog =
        (f == stderr)
        ? "puts -nonewline stderr \""
        : "puts -nonewline stdout \"";

    const char * const epilog = "\"";

    const int prolog_len = strlen(prolog);
    const int epilog_len = strlen(epilog);

    if ((fileno(f) != STDOUT_FILENO && fileno(f) != STDERR_FILENO &&
         f != stderr && f != stdout)
#ifdef THREADS
        || (fl_running && bgtid == thread_self())
#endif
        )
        return vfprintf(f, fmt, args);

    p = buf;

    // size: how much ist left for chars and terminating '\0'
    size = sizeof(buf) - prolog_len - epilog_len;
    // assert(size > 0);

    for (;;) {
        va_list ap;

        va_copy(ap, args);
        nchars = vsnprintf(p + prolog_len, size, fmt, ap);
        va_end(ap);

        if(nchars == -1) {           /* compatibility to old implementations */
            size *= 2;
        } else if (size < nchars + 1) {
            size = nchars + 1;
        } else {
            break;
        }

        if(p == buf)
            p = Tcl_Alloc(prolog_len + size + epilog_len);
        else
            p = Tcl_Realloc(p, prolog_len + size + epilog_len);
    }

    strncpy(p, prolog, prolog_len);

    s = p + prolog_len;
    for (escapes = 0; ; escapes++) {
        s = strpbrk(s, escape_chars);
        if (!s)
            break;
        s++;
    }

    if (escapes) {

        int new_size = prolog_len + nchars + escapes + epilog_len + 1;
        char *src, *dst;

        if (p != buf) {
            p = Tcl_Realloc(p, new_size);
        } else if (new_size > sizeof(buf)) {
            p = Tcl_Alloc(new_size);
            strcpy(p, buf);
        }

        src = p + prolog_len + nchars;
        dst = src + escapes;

        while (dst > src) {
            char c = *--src;
            *--dst = c;
            if (strchr(escape_chars, c))
                *--dst = '\\';
        }
    }

    strcpy(p + prolog_len + nchars + escapes, epilog);

    result = Tcl_Eval(spice_interp, p);

    if (p != buf)
        Tcl_Free(p);

    return nchars;
}


/*----------------------------------------------------------------------*/
/* Reimplement fprintf() as a call to Tcl_Eval().                       */
/*----------------------------------------------------------------------*/

int
tcl_fprintf(FILE *f, const char *format, ...)
{
    va_list args;
    int rtn;

    va_start (args, format);
    rtn = tcl_vfprintf(f, format, args);
    va_end(args);

    return rtn;
}


/*----------------------------------------------------------------------*/
/* Reimplement printf() as a call to Tcl_Eval().                       */
/*----------------------------------------------------------------------*/

int
tcl_printf(const char *format, ...)
{
    va_list args;
    int rtn;

    va_start (args, format);
    rtn = tcl_vfprintf(stdout, format, args);
    va_end(args);

    return rtn;
}


/*------------------------------------------------------*/
/* Console output flushing which goes along with the    */
/* routine tcl_vprintf() above.                         */
/*------------------------------------------------------*/

void
tcl_stdflush(FILE *f)
{
    Tcl_SavedResult state;
    static char stdstr[] = "flush stdxxx";
    char *stdptr = stdstr + 9;

#ifdef THREADS
    if (fl_running && bgtid == thread_self())
        return;
#endif

    Tcl_SaveResult(spice_interp, &state);
    strcpy(stdptr, (f == stderr) ? "err" : "out");
    Tcl_Eval(spice_interp, stdstr);
    Tcl_RestoreResult(spice_interp, &state);
}
