Move some fields from the HID* structure to HID_DRAW* and HID_DRAW_CLASS*

[geda-pcb/pcjc2.git] / src / hid / nelma / nelma.c

/*
 *                            COPYRIGHT
 *
 *  PCB, interactive printed circuit board design
 *
 *  NELMA (Numerical capacitance calculator) export HID
 *  Copyright (C) 2006 Tomaz Solc (tomaz.solc@tablix.org)
 *
 *  PNG export code is based on the PNG export HID
 *  Copyright (C) 2006 Dan McMahill
 *
 *  This program is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program; if not, write to the Free Software
 *  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 *
 */

/*
 * This HID exports a PCB layout into: o One layer mask file (PNG format) per
 * copper layer. o Nelma configuration file that contains netlist and pin
 * information.
 */

/*
 * FIXME:
 * 
 * If you have a section of a net that does not contain any pins then that
 * section will be missing from the Nelma's copper geometry.
 * 
 * For example:
 * 
 * this section will be ignored by Nelma |       |
 * 
 * ||             ||=======||            ||     component layer ||
 * ||       ||            || ||=============||       ||============||
 * solder layer
 * 
 * pin1           via      via           pin2
 * 
 * Single layer layouts are always exported correctly.
 * 
 */

#ifdef HAVE_CONFIG_H
#include "config.h"
#endif

#include <stdio.h>
#include <stdarg.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>

#include <time.h>

#include "global.h"
#include "error.h" /* Message() */
#include "data.h"
#include "misc.h"
#include "rats.h"

#include "hid.h"
#include "hid_draw.h"
#include "../hidint.h"
#include "hid/common/hidnogui.h"
#include "hid/common/draw_helpers.h"

#include <gd.h>

#include "hid/common/hidinit.h"

#ifdef HAVE_LIBDMALLOC
#include <dmalloc.h>
#endif

#define CRASH fprintf(stderr, "HID error: pcb called unimplemented PNG function %s.\n", __FUNCTION__); abort()

/* Needed for PNG export */

struct color_struct {
        /* the descriptor used by the gd library */
        int             c;

        /* so I can figure out what rgb value c refers to */
        unsigned int    r, g, b;
};

typedef struct nelma_gc_struct {
        struct hid_gc_struct hid_gc; /* Parent */

        EndCapStyle     cap;
        Coord           width;
        unsigned char   r, g, b;
        int             erase;
        struct color_struct *color;
        gdImagePtr      brush;
} *nelmaGC;

static HID nelma_hid;
static HID_DRAW nelma_graphics;
static HID_DRAW_CLASS nelma_graphics_class;

static struct color_struct *black = NULL, *white = NULL;
static Coord    linewidth = -1;
static gdImagePtr lastbrush = (gdImagePtr)((void *) -1);

/* gd image and file for PNG export */
static gdImagePtr nelma_im = NULL;
static FILE    *nelma_f = NULL;

static int      is_mask;
static int      is_drill;

/*
 * Which groups of layers to export into PNG layer masks. 1 means export, 0
 * means do not export.
 */
static int      nelma_export_group[MAX_GROUP];

/* Group that is currently exported. */
static int      nelma_cur_group;

/* Filename prefix that will be used when saving files. */
static const char *nelma_basename = NULL;

/* Horizontal DPI (grid points per inch) */
static int      nelma_dpi = -1;

/* Height of the copper layers in micrometers. */

/*
 * The height of the copper layer is currently taken as the vertical grid
 * step, since this is the smallest vertical feature in the layout.
 */
static int      nelma_copperh = -1;
/* Height of the substrate layers in micrometers. */
static int      nelma_substrateh = -1;
/* Relative permittivity of the substrate. */
static double   nelma_substratee = -1;

/* Permittivity of empty space (As/Vm) */
static const double nelma_air_epsilon = 8.85e-12;

HID_Attribute   nelma_attribute_list[] = {
        /* other HIDs expect this to be first.  */

/* %start-doc options "nelma Options"
@ftable @code
@item -- basename <string>
File name prefix.
@end ftable
%end-doc
*/
        {"basename", "File name prefix",
        HID_String, 0, 0, {0, 0, 0}, 0, 0},
#define HA_basename 0

/* %start-doc options "nelma Options"
@ftable @code
@item --dpi <num>
Horizontal scale factor (grid points/inch).
@end ftable
%end-doc
*/
        {"dpi", "Horizontal scale factor (grid points/inch)",
        HID_Integer, 0, 1000, {100, 0, 0}, 0, 0},
#define HA_dpi 1

/* %start-doc options "nelma Options"
@ftable @code
@item --copper-height <num>
Copper layer height (um).
@end ftable
%end-doc
*/
        {"copper-height", "Copper layer height (um)",
        HID_Integer, 0, 200, {100, 0, 0}, 0, 0},
#define HA_copperh 2

/* %start-doc options "nelma Options"
@ftable @code
@item --substrate-height <num>
Substrate layer height (um).
@end ftable
%end-doc
*/
        {"substrate-height", "Substrate layer height (um)",
        HID_Integer, 0, 10000, {2000, 0, 0}, 0, 0},
#define HA_substrateh 3

/* %start-doc options "nelma Options"
@ftable @code
@item --substrate-epsilon <num>
Substrate relative epsilon.
@end ftable
%end-doc
*/
        {"substrate-epsilon", "Substrate relative epsilon",
        HID_Real, 0, 100, {0, 0, 4.0}, 0, 0},
#define HA_substratee 4
};

#define NUM_OPTIONS (sizeof(nelma_attribute_list)/sizeof(nelma_attribute_list[0]))

REGISTER_ATTRIBUTES(nelma_attribute_list)
        static HID_Attr_Val nelma_values[NUM_OPTIONS];

/* *** Utility funcions **************************************************** */

/* convert from default PCB units to nelma units */
static int pcb_to_nelma (Coord pcb)
{
  return COORD_TO_INCH(pcb) * nelma_dpi;
}

static char    *
nelma_get_png_name(const char *basename, const char *suffix)
{
        char           *buf;
        int             len;

        len = strlen(basename) + strlen(suffix) + 6;
        buf = (char *)malloc(sizeof(*buf) * len);

        sprintf(buf, "%s.%s.png", basename, suffix);

        return buf;
}

/* Retrieves coordinates (in default PCB units) of a pin or pad. */
/* Copied from netlist.c */
static int 
pin_name_to_xy (LibraryEntryType * pin, Coord *x, Coord *y)
{
        ConnectionType  conn;
        if (!SeekPad(pin, &conn, false))
                return 1;
        switch (conn.type) {
        case PIN_TYPE:
                *x = ((PinType *) (conn.ptr2))->X;
                *y = ((PinType *) (conn.ptr2))->Y;
                return 0;
        case PAD_TYPE:
                *x = ((PadType *) (conn.ptr2))->Point1.X;
                *y = ((PadType *) (conn.ptr2))->Point1.Y;
                return 0;
        }
        return 1;
}

/* *** Exporting netlist data and geometry to the nelma config file ******** */

static void 
nelma_write_space(FILE * out)
{
        double          xh, zh;

        int             z;
        int             i, idx;
        const char     *ext;

        xh = 2.54e-2 / ((double) nelma_dpi);
        zh = nelma_copperh * 1e-6;

        fprintf(out, "\n/* **** Space **** */\n\n");

        fprintf(out, "space pcb {\n");
        fprintf(out, "\tstep = { %e, %e, %e }\n", xh, xh, zh);
        fprintf(out, "\tlayers = {\n");

        fprintf(out, "\t\t\"air-top\",\n");
        fprintf(out, "\t\t\"air-bottom\"");

        z = 10;
        for (i = 0; i < MAX_GROUP; i++)
                if (nelma_export_group[i]) {
                        idx = (i >= 0 && i < max_group) ?
                                PCB->LayerGroups.Entries[i][0] : i;
                        ext = layer_type_to_file_name(idx, FNS_fixed);

                        if (z != 10) {
                                fprintf(out, ",\n");
                                fprintf(out, "\t\t\"substrate-%d\"", z);
                                z++;
                        }
                        fprintf(out, ",\n");
                        fprintf(out, "\t\t\"%s\"", ext);
                        z++;
                }
        fprintf(out, "\n\t}\n");
        fprintf(out, "}\n");
}


static void 
nelma_write_material(FILE * out, char *name, char *type, double e)
{
        fprintf(out, "material %s {\n", name);
        fprintf(out, "\ttype = \"%s\"\n", type);
        fprintf(out, "\tpermittivity = %e\n", e);
        fprintf(out, "\tconductivity = 0.0\n");
        fprintf(out, "\tpermeability = 0.0\n");
        fprintf(out, "}\n");
}

static void 
nelma_write_materials(FILE * out)
{
        fprintf(out, "\n/* **** Materials **** */\n\n");

        nelma_write_material(out, "copper", "metal", nelma_air_epsilon);
        nelma_write_material(out, "air", "dielectric", nelma_air_epsilon);
        nelma_write_material(out, "composite", "dielectric",
                             nelma_air_epsilon * nelma_substratee);
}

static void 
nelma_write_nets(FILE * out)
{
        LibraryType     netlist;
        LibraryMenuType *net;
        LibraryEntryType *pin;

        int             n, m, i, idx;

        const char     *ext;

        netlist = PCB->NetlistLib;

        fprintf(out, "\n/* **** Nets **** */\n\n");

        for (n = 0; n < netlist.MenuN; n++) {
                net = &netlist.Menu[n];

                /* Weird, but correct */
                fprintf(out, "net %s {\n", &net->Name[2]);

                fprintf(out, "\tobjects = {\n");

                for (m = 0; m < net->EntryN; m++) {
                        pin = &net->Entry[m];

                        /* pin_name_to_xy(pin, &x, &y); */

                        for (i = 0; i < MAX_GROUP; i++)
                                if (nelma_export_group[i]) {
                                        idx = (i >= 0 && i < max_group) ?
                                                PCB->LayerGroups.Entries[i][0] : i;
                                        ext = layer_type_to_file_name(idx, FNS_fixed);

                                        if (m != 0 || i != 0)
                                                fprintf(out, ",\n");
                                        fprintf(out, "\t\t\"%s-%s\"", pin->ListEntry,
                                                ext);
                                }
                }

                fprintf(out, "\n");
                fprintf(out, "\t}\n");
                fprintf(out, "}\n");
        }
}

static void 
nelma_write_layer(FILE * out, int z, int h,
                  const char *name, int full,
                  char *mat)
{
        LibraryType     netlist;
        LibraryMenuType *net;
        LibraryEntryType *pin;

        int             n, m;

        fprintf(out, "layer %s {\n", name);
        fprintf(out, "\theight = %d\n", h);
        fprintf(out, "\tz-order = %d\n", z);
        fprintf(out, "\tmaterial = \"%s\"\n", mat);

        if (full) {
                fprintf(out, "\tobjects = {\n");
                netlist = PCB->NetlistLib;

                for (n = 0; n < netlist.MenuN; n++) {
                        net = &netlist.Menu[n];

                        for (m = 0; m < net->EntryN; m++) {
                                pin = &net->Entry[m];

                                if (m != 0 || n != 0)
                                        fprintf(out, ",\n");
                                fprintf(out, "\t\t\"%s-%s\"", pin->ListEntry,
                                        name);
                        }

                }
                fprintf(out, "\n\t}\n");
        }
        fprintf(out, "}\n");
}

static void 
nelma_write_layers(FILE * out)
{
        int             i, idx;
        int             z;

        const char     *ext;
        char            buf[100];

        int             subh;

        subh = nelma_substrateh / nelma_copperh;

        fprintf(out, "\n/* **** Layers **** */\n\n");

        /* Air layers on top and bottom of the stack */
        /* Their height is double substrate height. */
        nelma_write_layer(out, 1, 2 * subh, "air-top", 0, "air");
        nelma_write_layer(out, 1000, 2 * subh, "air-bottom", 0, "air");

        z = 10;
        for (i = 0; i < MAX_GROUP; i++)
                if (nelma_export_group[i]) {
                        idx = (i >= 0 && i < max_group) ?
                                PCB->LayerGroups.Entries[i][0] : i;
                        ext = layer_type_to_file_name(idx, FNS_fixed);

                        if (z != 10) {
                                sprintf(buf, "substrate-%d", z);
                                nelma_write_layer(out,
                                                  z,
                                                  subh,
                                                  buf,
                                                  0,
                                                  "composite");
                                z++;
                        }
                        /*
                         * FIXME: for layers that are not on top or bottom,
                         * the material should be "composite"
                         */
                        nelma_write_layer(out, z, 1, ext, 1, "air");

                        z++;
                }
}

static void 
nelma_write_object(FILE * out, LibraryEntryType *pin)
{
        int             i, idx;
        Coord           px = 0, py = 0;
        int             x, y;

        char           *f;
        const char     *ext;

        pin_name_to_xy (pin, &px, &py);

        x = pcb_to_nelma (px);
        y = pcb_to_nelma (py);

        for (i = 0; i < MAX_GROUP; i++)
                if (nelma_export_group[i]) {
                        idx = (i >= 0 && i < max_group) ?
                                PCB->LayerGroups.Entries[i][0] : i;
                        ext = layer_type_to_file_name(idx, FNS_fixed);

                        fprintf(out, "object %s-%s {\n", pin->ListEntry, ext);
                        fprintf(out, "\tposition = { 0, 0 }\n");
                        fprintf(out, "\tmaterial = \"copper\"\n");
                        fprintf(out, "\ttype = \"image\"\n");
                        fprintf(out, "\trole = \"net\"\n");

                        f = nelma_get_png_name(nelma_basename, ext);

                        fprintf(out, "\tfile = \"%s\"\n", f);

                        free(f);

                        fprintf(out, "\tfile-pos = { %d, %d }\n", x, y);
                        fprintf(out, "}\n");
                }
}

static void 
nelma_write_objects(FILE * out)
{
        LibraryType     netlist;
        LibraryMenuType *net;
        LibraryEntryType *pin;

        int             n, m;

        netlist = PCB->NetlistLib;

        fprintf(out, "\n/* **** Objects **** */\n\n");

        for (n = 0; n < netlist.MenuN; n++) {
                net = &netlist.Menu[n];

                for (m = 0; m < net->EntryN; m++) {
                        pin = &net->Entry[m];

                        nelma_write_object(out, pin);
                }
        }
}

/* *** Main export callback ************************************************ */

static void 
nelma_parse_arguments(int *argc, char ***argv)
{
        hid_register_attributes(nelma_attribute_list,
                                sizeof(nelma_attribute_list) /
                                sizeof(nelma_attribute_list[0]));
        hid_parse_command_line(argc, argv);
}

static HID_Attribute *
nelma_get_export_options(int *n)
{
        static char    *last_made_filename = 0;

        if (PCB) {
                derive_default_filename(PCB->Filename,
                                        &nelma_attribute_list[HA_basename],
                                        ".nelma",
                                        &last_made_filename);
        }
        if (n) {
                *n = NUM_OPTIONS;
        }
        return nelma_attribute_list;
}

/* Populates nelma_export_group array */
void 
nelma_choose_groups()
{
        int             n, m;
        LayerType      *layer;

        /* Set entire array to 0 (don't export any layer groups by default */
        memset(nelma_export_group, 0, sizeof(nelma_export_group));

        for (n = 0; n < max_copper_layer; n++) {
                layer = &PCB->Data->Layer[n];

                if (layer->LineN || layer->TextN || layer->ArcN ||
                    layer->PolygonN) {
                        /* layer isn't empty */

                        /*
                         * is this check necessary? It seems that special
                         * layers have negative indexes?
                         */

                        if (SL_TYPE(n) == 0) {
                                /* layer is a copper layer */
                                m = GetLayerGroupNumberByNumber(n);

                                /* the export layer */
                                nelma_export_group[m] = 1;
                        }
                }
        }
}

static void 
nelma_alloc_colors()
{
        /*
         * Allocate white and black -- the first color allocated becomes the
         * background color
         */

        white = (struct color_struct *) malloc(sizeof(*white));
        white->r = white->g = white->b = 255;
        white->c = gdImageColorAllocate(nelma_im, white->r, white->g, white->b);

        black = (struct color_struct *) malloc(sizeof(*black));
        black->r = black->g = black->b = 0;
        black->c = gdImageColorAllocate(nelma_im, black->r, black->g, black->b);
}

static void 
nelma_start_png(const char *basename, const char *suffix)
{
        int             h, w;
        char           *buf;

        buf = nelma_get_png_name(basename, suffix);

        h = pcb_to_nelma(PCB->MaxHeight);
        w = pcb_to_nelma(PCB->MaxWidth);

        /* nelma_im = gdImageCreate (w, h); */

        /* Nelma only works with true color images */
        nelma_im = gdImageCreate(w, h);
        nelma_f = fopen(buf, "wb");

        nelma_alloc_colors();

        free(buf);
}

static void 
nelma_finish_png()
{
#ifdef HAVE_GDIMAGEPNG
        gdImagePng(nelma_im, nelma_f);
#else
        Message("NELMA: PNG not supported by gd. Can't write layer mask.\n");
#endif
        gdImageDestroy(nelma_im);
        fclose(nelma_f);

        free(white);
        free(black);

        nelma_im = NULL;
        nelma_f = NULL;
}

void 
nelma_start_png_export()
{
        BoxType         region;

        region.X1 = 0;
        region.Y1 = 0;
        region.X2 = PCB->MaxWidth;
        region.Y2 = PCB->MaxHeight;

        linewidth = -1;
        lastbrush = (gdImagePtr)((void *) -1);

        hid_expose_callback(&nelma_hid, &region, 0);
}

static void 
nelma_do_export(HID_Attr_Val * options)
{
        int             save_ons[MAX_LAYER + 2];
        int             i, idx;
        FILE           *nelma_config;
        char           *buf;
        int             len;

        time_t          t;

        if (!options) {
                nelma_get_export_options(0);
                for (i = 0; i < NUM_OPTIONS; i++) {
                        nelma_values[i] = nelma_attribute_list[i].default_val;
                }
                options = nelma_values;
        }
        nelma_basename = options[HA_basename].str_value;
        if (!nelma_basename) {
                nelma_basename = "pcb-out";
        }
        nelma_dpi = options[HA_dpi].int_value;
        if (nelma_dpi < 0) {
                fprintf(stderr, "ERROR:  dpi may not be < 0\n");
                return;
        }
        nelma_copperh = options[HA_copperh].int_value;
        nelma_substrateh = options[HA_substrateh].int_value;
        nelma_substratee = options[HA_substratee].real_value;

        nelma_choose_groups();

        for (i = 0; i < MAX_GROUP; i++) {
                if (nelma_export_group[i]) {

                        nelma_cur_group = i;

                        /* magic */
                        idx = (i >= 0 && i < max_group) ?
                                PCB->LayerGroups.Entries[i][0] : i;

                        nelma_start_png(nelma_basename,
                                        layer_type_to_file_name(idx, FNS_fixed));

                        hid_save_and_show_layer_ons(save_ons);
                        nelma_start_png_export();
                        hid_restore_layer_ons(save_ons);

                        nelma_finish_png();
                }
        }

        len = strlen(nelma_basename) + 4;
        buf = (char *)malloc(sizeof(*buf) * len);

        sprintf(buf, "%s.em", nelma_basename);
        nelma_config = fopen(buf, "w");

        free(buf);

        fprintf(nelma_config, "/* Made with PCB Nelma export HID */");
        t = time(NULL);
        fprintf(nelma_config, "/* %s */", ctime(&t));

        nelma_write_nets(nelma_config);
        nelma_write_objects(nelma_config);
        nelma_write_layers(nelma_config);
        nelma_write_materials(nelma_config);
        nelma_write_space(nelma_config);

        fclose(nelma_config);
}

/* *** PNG export (slightly modified code from PNG export HID) ************* */

static int 
nelma_set_layer(const char *name, int group, int empty)
{
        int             idx = (group >= 0 && group < max_group) ?
        PCB->LayerGroups.Entries[group][0] : group;

        if (name == 0) {
                name = PCB->Data->Layer[idx].Name;
        }
        if (strcmp(name, "invisible") == 0) {
                return 0;
        }
        is_drill = (SL_TYPE(idx) == SL_PDRILL || SL_TYPE(idx) == SL_UDRILL);
        is_mask = (SL_TYPE(idx) == SL_MASK);

        if (is_mask) {
                /* Don't print masks */
                return 0;
        }
        if (is_drill) {
                /*
                 * Print 'holes', so that we can fill gaps in the copper
                 * layer
                 */
                return 1;
        }
        if (group == nelma_cur_group) {
                return 1;
        }
        return 0;
}

static hidGC 
nelma_make_gc(void)
{
        hidGC gc = (hidGC) calloc (1, sizeof(struct hid_gc_struct));
        nelmaGC nelma_gc = (nelmaGC)gc;

        gc->hid = &nelma_hid;
        gc->hid_draw = &nelma_graphics;

        nelma_gc->cap = Trace_Cap;
        nelma_gc->width = 1;
        nelma_gc->color = (struct color_struct *) malloc(sizeof(*nelma_gc->color));
        nelma_gc->color->r = nelma_gc->color->g = nelma_gc->color->b = 0;
        nelma_gc->color->c = 0;

        return gc;
}

static void 
nelma_destroy_gc(hidGC gc)
{
        free(gc);
}

static void 
nelma_use_mask(enum mask_mode mode)
{
        /* does nothing */
}

static void 
nelma_set_color(hidGC gc, const char *name)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        if (nelma_im == NULL) {
                return;
        }
        if (name == NULL) {
                name = "#ff0000";
        }
        if (!strcmp(name, "drill")) {
                nelma_gc->color = black;
                nelma_gc->erase = 0;
                return;
        }
        if (!strcmp(name, "erase")) {
                /* FIXME -- should be background, not white */
                nelma_gc->color = white;
                nelma_gc->erase = 1;
                return;
        }
        nelma_gc->color = black;
        nelma_gc->erase = 0;
        return;
}

static void
nelma_set_line_cap(hidGC gc, EndCapStyle style)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        nelma_gc->cap = style;
}

static void
nelma_set_line_width(hidGC gc, Coord width)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        nelma_gc->width = width;
}

static void
nelma_set_draw_xor(hidGC gc, int xor_)
{
        ;
}

static void
nelma_set_draw_faded(hidGC gc, int faded)
{
}

static void
use_gc(hidGC gc)
{
        nelmaGC nelma_gc = (nelmaGC)gc;
        int             need_brush = 0;

        if (gc->hid != &nelma_hid) {
                fprintf(stderr, "Fatal: GC from another HID passed to nelma HID\n");
                abort();
        }
        if (linewidth != nelma_gc->width) {
                /* Make sure the scaling doesn't erase lines completely */
                /*
                if (SCALE (nelma_gc->width) == 0 && nelma_gc->width > 0)
                  gdImageSetThickness (im, 1);
                else
                */
                gdImageSetThickness(nelma_im, pcb_to_nelma(nelma_gc->width));
                linewidth = nelma_gc->width;
                need_brush = 1;
        }
        if (lastbrush != nelma_gc->brush || need_brush) {
                static void    *bcache = 0;
                hidval          bval;
                char            name[256];
                char            type;
                int             r;

                switch (nelma_gc->cap) {
                case Round_Cap:
                case Trace_Cap:
                        type = 'C';
                        r = pcb_to_nelma(nelma_gc->width / 2);
                        break;
                default:
                case Square_Cap:
                        r = pcb_to_nelma(nelma_gc->width);
                        type = 'S';
                        break;
                }
                sprintf(name, "#%.2x%.2x%.2x_%c_%d", nelma_gc->color->r, nelma_gc->color->g,
                        nelma_gc->color->b, type, r);

                if (hid_cache_color(0, name, &bval, &bcache)) {
                  nelma_gc->brush = (gdImagePtr)bval.ptr;
                } else {
                        int             bg, fg;
                        if (type == 'C')
                                nelma_gc->brush = gdImageCreate(2 * r + 1, 2 * r + 1);
                        else
                                nelma_gc->brush = gdImageCreate(r + 1, r + 1);
                        bg = gdImageColorAllocate(nelma_gc->brush, 255, 255, 255);
                        fg =
                                gdImageColorAllocate(nelma_gc->brush, nelma_gc->color->r, nelma_gc->color->g,
                                                     nelma_gc->color->b);
                        gdImageColorTransparent(nelma_gc->brush, bg);

                        /*
                         * if we shrunk to a radius/box width of zero, then just use
                         * a single pixel to draw with.
                         */
                        if (r == 0)
                                gdImageFilledRectangle(nelma_gc->brush, 0, 0, 0, 0, fg);
                        else {
                                if (type == 'C')
                                        gdImageFilledEllipse(nelma_gc->brush, r, r, 2 * r, 2 * r, fg);
                                else
                                        gdImageFilledRectangle(nelma_gc->brush, 0, 0, r, r, fg);
                        }
                        bval.ptr = nelma_gc->brush;
                        hid_cache_color(1, name, &bval, &bcache);
                }

                gdImageSetBrush(nelma_im, nelma_gc->brush);
                lastbrush = nelma_gc->brush;

        }
}

static void
nelma_draw_rect(hidGC gc, Coord x1, Coord y1, Coord x2, Coord y2)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        use_gc(gc);
        gdImageRectangle(nelma_im,
                         pcb_to_nelma(x1), pcb_to_nelma(y1),
                         pcb_to_nelma(x2), pcb_to_nelma(y2), nelma_gc->color->c);
}

static void
nelma_fill_rect(hidGC gc, Coord x1, Coord y1, Coord x2, Coord y2)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        use_gc(gc);
        gdImageSetThickness(nelma_im, 0);
        linewidth = 0;
        gdImageFilledRectangle(nelma_im, pcb_to_nelma(x1), pcb_to_nelma(y1),
                          pcb_to_nelma(x2), pcb_to_nelma(y2), nelma_gc->color->c);
}

static void
nelma_draw_line(hidGC gc, Coord x1, Coord y1, Coord x2, Coord y2)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        if (x1 == x2 && y1 == y2) {
                Coord             w = nelma_gc->width / 2;
                nelma_fill_rect(gc, x1 - w, y1 - w, x1 + w, y1 + w);
                return;
        }
        use_gc(gc);

        gdImageSetThickness(nelma_im, 0);
        linewidth = 0;
        gdImageLine(nelma_im, pcb_to_nelma(x1), pcb_to_nelma(y1),
                    pcb_to_nelma(x2), pcb_to_nelma(y2), gdBrushed);
}

static void
nelma_draw_arc(hidGC gc, Coord cx, Coord cy, Coord width, Coord height,
               Angle start_angle, Angle delta_angle)
{
        Angle sa, ea;

        /*
         * in gdImageArc, 0 degrees is to the right and +90 degrees is down
         * in pcb, 0 degrees is to the left and +90 degrees is down
         */
        start_angle = 180 - start_angle;
        delta_angle = -delta_angle;
        if (delta_angle > 0) {
                sa = start_angle;
                ea = start_angle + delta_angle;
        } else {
                sa = start_angle + delta_angle;
                ea = start_angle;
        }

        /*
         * make sure we start between 0 and 360 otherwise gd does strange
         * things
         */
        sa = NormalizeAngle (sa);
        ea = NormalizeAngle (ea);

#if 0
        printf("draw_arc %d,%d %dx%d %d..%d %d..%d\n",
               cx, cy, width, height, start_angle, delta_angle, sa, ea);
        printf("gdImageArc (%p, %d, %d, %d, %d, %d, %d, %d)\n",
               im, SCALE_X(cx), SCALE_Y(cy),
               SCALE(width), SCALE(height), sa, ea, nelma_gc->color->c);
#endif
        use_gc(gc);
        gdImageSetThickness(nelma_im, 0);
        linewidth = 0;
        gdImageArc(nelma_im, pcb_to_nelma(cx), pcb_to_nelma(cy),
                   pcb_to_nelma(2 * width), pcb_to_nelma(2 * height), sa, ea, gdBrushed);
}

static void
nelma_fill_circle(hidGC gc, Coord cx, Coord cy, Coord radius)
{
        nelmaGC nelma_gc = (nelmaGC)gc;

        use_gc(gc);

        gdImageSetThickness(nelma_im, 0);
        linewidth = 0;
        gdImageFilledEllipse(nelma_im, pcb_to_nelma(cx), pcb_to_nelma(cy),
          pcb_to_nelma(2 * radius), pcb_to_nelma(2 * radius), nelma_gc->color->c);

}

static void
nelma_fill_polygon(hidGC gc, int n_coords, Coord *x, Coord *y)
{
        nelmaGC nelma_gc = (nelmaGC)gc;
        int             i;
        gdPoint        *points;

        points = (gdPoint *) malloc(n_coords * sizeof(gdPoint));
        if (points == NULL) {
                fprintf(stderr, "ERROR:  nelma_fill_polygon():  malloc failed\n");
                exit(1);
        }
        use_gc(gc);
        for (i = 0; i < n_coords; i++) {
                points[i].x = pcb_to_nelma(x[i]);
                points[i].y = pcb_to_nelma(y[i]);
        }
        gdImageSetThickness(nelma_im, 0);
        linewidth = 0;
        gdImageFilledPolygon(nelma_im, points, n_coords, nelma_gc->color->c);
        free(points);
}

static void
nelma_calibrate(double xval, double yval)
{
        CRASH;
}

static void
nelma_set_crosshair(int x, int y, int a)
{
}

/* *** Miscellaneous ******************************************************* */

#include "dolists.h"

void
hid_nelma_init()
{
  memset (&nelma_hid, 0, sizeof (HID));
  memset (&nelma_graphics, 0, sizeof (HID_DRAW));
  memset (&nelma_graphics_class, 0, sizeof (HID_DRAW_CLASS));

  common_nogui_init (&nelma_hid);

  nelma_hid.struct_size         = sizeof (HID);
  nelma_hid.name                = "nelma";
  nelma_hid.description         = "Numerical analysis package export";
  nelma_hid.exporter            = 1;

  nelma_hid.get_export_options  = nelma_get_export_options;
  nelma_hid.do_export           = nelma_do_export;
  nelma_hid.parse_arguments     = nelma_parse_arguments;
  nelma_hid.calibrate           = nelma_calibrate;
  nelma_hid.set_crosshair       = nelma_set_crosshair;

  nelma_hid.graphics            = &nelma_graphics;

  common_draw_helpers_class_init (&nelma_graphics_class);

  nelma_graphics_class.set_layer      = nelma_set_layer;
  nelma_graphics_class.make_gc        = nelma_make_gc;
  nelma_graphics_class.destroy_gc     = nelma_destroy_gc;
  nelma_graphics_class.use_mask       = nelma_use_mask;
  nelma_graphics_class.set_color      = nelma_set_color;
  nelma_graphics_class.set_line_cap   = nelma_set_line_cap;
  nelma_graphics_class.set_line_width = nelma_set_line_width;
  nelma_graphics_class.set_draw_xor   = nelma_set_draw_xor;
  nelma_graphics_class.set_draw_faded = nelma_set_draw_faded;
  nelma_graphics_class.draw_line      = nelma_draw_line;
  nelma_graphics_class.draw_arc       = nelma_draw_arc;
  nelma_graphics_class.draw_rect      = nelma_draw_rect;
  nelma_graphics_class.fill_circle    = nelma_fill_circle;
  nelma_graphics_class.fill_polygon   = nelma_fill_polygon;
  nelma_graphics_class.fill_rect      = nelma_fill_rect;

  nelma_graphics.klass = &nelma_graphics_class;
  nelma_graphics.poly_before = true;
  common_draw_helpers_init (&nelma_graphics);

  hid_register_hid (&nelma_hid);

#include "nelma_lists.h"
}