Logo Search packages:      
Sourcecode: gretl version File versions

ssheet.c

/* 
 *  gretl -- Gnu Regression, Econometrics and Time-series Library
 *  Copyright (C) 2001 Allin Cottrell and Riccardo "Jack" Lucchetti
 * 
 *  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 3 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, see <http://www.gnu.org/licenses/>.
 * 
 */

#include "gretl.h"
#include "treeutils.h"
#include "ssheet.h"
#include "dlgutils.h"
#include "menustate.h"
#include "selector.h"
#include "session.h"
#include "winstack.h"
#include "toolbar.h"
#include "usermat.h"
#include "matrix_extra.h"
#include "gretl_scalar.h"

#include <errno.h>
#include <ctype.h>
#include <float.h>

#define SSDEBUG 0
#define CELLDEBUG 0

typedef enum {
    SHEET_SUBSAMPLED    = 1 << 0,
    SHEET_SHORT_VARLIST = 1 << 1,
    SHEET_INSERT_OBS_OK = 1 << 2,
    SHEET_ADD_OBS_OK    = 1 << 3,
    SHEET_USE_COMMA     = 1 << 4,
    SHEET_MODIFIED      = 1 << 5,
    SHEET_COLNAME_MOD   = 1 << 6
} SheetFlags;

enum {
    SHEET_AT_END,
    SHEET_AT_POINT
};

enum {
    NEXT_DOWN,
    NEXT_RIGHT,
    NEXT_UP,
    NEXT_LEFT,
    NEXT_SAME
};

typedef struct {
    GtkWidget *view;
    GtkWidget *win;
    GtkWidget *locator;
    GtkWidget *entry;
    GtkWidget *popup;
    GtkWidget *save;
    GtkWidget *apply;
    GtkUIManager *ui;
    GtkCellRenderer *dumbcell;
    GtkCellRenderer *datacell;
    GdkPixbuf *pbuf;
    gchar location[64];
    gretl_matrix *matrix;
    gretl_matrix *oldmat;
    char mname[VNAMELEN];
    const char **colnames;
    int *varlist;
    int datacols, datarows;
    int totcols;
    int orig_nobs;
    int added_vars;
    int orig_main_v;
    int next;
    SheetCmd cmd;
    SheetFlags flags;
    guint cid;
    guint point;
} Spreadsheet;

#define MATRIX_DIGITS DBL_DIG

static void set_up_sheet_column (GtkTreeViewColumn *column, gint width, 
                         gboolean expand);
static gint get_data_col_width (void);
static int add_data_column (Spreadsheet *sheet);
static void create_sheet_cell_renderers (Spreadsheet *sheet);

static void matrix_fill_callback (GtkAction *action, gpointer data);
static void matrix_props_callback (GtkAction *action, gpointer data);
static void matrix_edit_callback (GtkAction *action, gpointer data);
static int update_sheet_from_matrix (Spreadsheet *sheet);
static void size_matrix_window (Spreadsheet *sheet);
static void set_ok_transforms (Spreadsheet *sheet);

static void sheet_show_popup (GtkWidget *w, Spreadsheet *sheet);
static void get_data_from_sheet (GtkWidget *w, Spreadsheet *sheet);
static gint maybe_exit_sheet (GtkWidget *w, Spreadsheet *sheet);

static void 
spreadsheet_scroll_to_foot (Spreadsheet *sheet, int row, int col);

enum {
    SERIES_ADD_BTN,
    SERIES_APPLY_BTN
};

static GretlToolItem series_items[] = {
    { N_("Add..."), GTK_STOCK_ADD,    G_CALLBACK(sheet_show_popup),    SERIES_ADD_BTN },
    { N_("Apply"),  GTK_STOCK_APPLY,  G_CALLBACK(get_data_from_sheet), SERIES_APPLY_BTN },
    { N_("Close"),  GTK_STOCK_CLOSE,  G_CALLBACK(maybe_exit_sheet),    0 }
};

static int n_series_items = G_N_ELEMENTS(series_items);

const gchar *matrix_ui = 
    "<ui>"
    "  <menubar>"
    "    <menu action='Fill'>"
    "      <menuitem action='FillIdentity'/>"
    "      <menuitem action='FillUniform'/>"
    "      <menuitem action='FillNormal'/>"
    "    </menu>"
    "    <menu action='Properties'>"
    "      <menuitem action='PropsView'/>"
    "    </menu>"
    "    <menu action='Transform'>"
    "      <menuitem action='XTX'/>"
    "      <menuitem action='Transpose'/>"
    "      <menuitem action='Cholesky'/>"
    "      <menuitem action='Invert'/>"
    "      <menuitem action='ScalarMult'/>"
    "      <menuitem action='ScalarDiv'/>"
    "    </menu>"
    "  </menubar>"
    "</ui>";

static GtkActionEntry matrix_items[] = {
    { "Fill", NULL, N_("_Fill"), NULL, NULL, NULL },
    { "FillIdentity", NULL, N_("_Identity matrix"), NULL, NULL, G_CALLBACK(matrix_fill_callback) },
    { "FillUniform", NULL, N_("_Uniform random"), NULL, NULL, G_CALLBACK(matrix_fill_callback) },
    { "FillNormal", NULL, N_("_Normal random"), NULL, NULL, G_CALLBACK(matrix_fill_callback) },
    { "Properties", NULL, N_("_Properties"), NULL, NULL, NULL },
    { "PropsView", NULL, N_("_View"), NULL, NULL, G_CALLBACK(matrix_props_callback) },
    { "Transform", NULL, N_("_Transform"), NULL, NULL, NULL },
    { "XTX", NULL, N_("_X'X"), NULL, NULL, G_CALLBACK(matrix_edit_callback) },
    { "Transpose", NULL, N_("_Transpose"), NULL, NULL, G_CALLBACK(matrix_edit_callback) },
    { "Cholesky", NULL, N_("_Cholesky"), NULL, NULL, G_CALLBACK(matrix_edit_callback) },
    { "Invert", NULL, N_("_Invert"), NULL, NULL, G_CALLBACK(matrix_edit_callback) },
    { "ScalarMult", NULL, N_("_Multiply by scalar"), NULL, NULL, G_CALLBACK(matrix_edit_callback) },
    { "ScalarDiv", NULL, N_("_Divide by scalar"), NULL, NULL, G_CALLBACK(matrix_edit_callback) }
};

#define sheet_is_modified(s) (s->flags & SHEET_MODIFIED)

static void sheet_set_modified (Spreadsheet *sheet, gboolean s)
{
    if (s) {
      sheet->flags |= SHEET_MODIFIED;
    } else {
      sheet->flags &= ~SHEET_MODIFIED;
    }

    if (sheet->save != NULL) {
      gtk_widget_set_sensitive(sheet->save, s);
    } else if (sheet->apply != NULL) {
      gtk_widget_set_sensitive(sheet->apply, s);
    }
}

static char *single_underscores (char *targ, const char *src)
{
    char *p = targ;

    while (*src) {
      if (*src == '_' && *(src + 1) == '_') {
          src++;
          *p++ = '_';
      } else {
          *p++ = *src;
      }
      src++;
    }

    *p = '\0';

    return targ;
}

static void set_locator_label (Spreadsheet *sheet, GtkTreePath *path,
                         GtkTreeViewColumn *column)
{
    GtkTreeModel *model = gtk_tree_view_get_model(GTK_TREE_VIEW(sheet->view));
    GtkTreeIter iter;
    gchar *rstr;
    const gchar *cstr;
    char tmp[VNAMELEN];

    gtk_tree_model_get_iter(model, &iter, path);
    gtk_tree_model_get(model, &iter, 0, &rstr, -1);
    cstr = gtk_tree_view_column_get_title(column);
    if (sheet->matrix != NULL) {
      sprintf(sheet->location, "%s, %s", rstr, cstr);
    } else {
      single_underscores(tmp, cstr);
      sprintf(sheet->location, "%s, %s", tmp, rstr);
    }
    gtk_statusbar_pop(GTK_STATUSBAR(sheet->locator), sheet->cid);
    gtk_statusbar_push(GTK_STATUSBAR(sheet->locator), 
                   sheet->cid, sheet->location);
    g_free(rstr);
}

static void move_to_next_column (Spreadsheet *sheet, GtkTreePath *path,
                         GtkTreeViewColumn *col)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    gpointer p = g_object_get_data(G_OBJECT(col), "colnum");
    GtkTreeViewColumn *newcol;
    int colnum;

    if (p == NULL) {
      return;
    }

    colnum = GPOINTER_TO_INT(p);

    if (sheet->next == NEXT_RIGHT) {
      colnum++;
    } else {
      colnum--;
    }  

#if CELLDEBUG
    fprintf(stderr, "move_to_next_column: colnum = %d\n", colnum);
#endif

    if (colnum < 1 || colnum > sheet->datacols) {
      return;
    }
    
    newcol = gtk_tree_view_get_column(view, colnum);

    if (newcol != NULL) {
      gtk_tree_view_set_cursor(view, path, newcol, FALSE);
    }    
}

static void move_to_next_cell (Spreadsheet *sheet, GtkTreePath *path,
                         GtkTreeViewColumn *column)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    gint newrow;

    if (sheet->next == NEXT_SAME) {
      sheet->next = NEXT_DOWN;
      return;
    }

    if (sheet->next == NEXT_RIGHT) {
#if CELLDEBUG
      fprintf(stderr, "move_to_next: got NEXT_RIGHT\n");
#endif
      move_to_next_column(sheet, path, column);
      return;
    }

    if (sheet->next == NEXT_UP) {
      newrow = gtk_tree_path_get_indices(path)[0] - 1;
    } else {
      newrow = gtk_tree_path_get_indices(path)[0] + 1;
    }

#if CELLDEBUG
    fprintf(stderr, "move_to_next: newrow = %d\n", newrow);
#endif

    if (newrow >= 0 && newrow < sheet->datarows) {
      GtkTreePath *newpath;
      gchar pstr[8];

      sprintf(pstr, "%d", newrow);
      newpath = gtk_tree_path_new_from_string(pstr);
      if (newpath != NULL) {
          gtk_tree_view_set_cursor(view, newpath, column, FALSE);
          gtk_tree_path_free(newpath);
      }
    } else {
      move_to_next_column(sheet, path, column);
    }

    sheet->next = NEXT_DOWN;
}

static void update_sheet_matrix (Spreadsheet *sheet)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeIter iter;
    GtkTreeModel *model;
    gchar *numstr;
    double x;
    int i, j;

    model = gtk_tree_view_get_model(view);
    gtk_tree_model_get_iter_first(model, &iter);

    for (i=0; i<sheet->datarows; i++) {
      for (j=0; j<sheet->datacols; j++) {
          gtk_tree_model_get(model, &iter, j+1, &numstr, -1);
          x = atof(numstr); 
          gretl_matrix_set(sheet->matrix, i, j, x);
          g_free(numstr);
      }
      gtk_tree_model_iter_next(model, &iter);
    }
}

static void maybe_update_column_names (Spreadsheet *sheet)
{
    if (sheet->flags & SHEET_COLNAME_MOD) {
      gretl_matrix *M = sheet->oldmat;
      GtkTreeViewColumn *col;
      const char *title;
      char *cnames, tmp[13];
      int j, err = 0;

      cnames = malloc(M->cols * 13 + 1);

      if (cnames == NULL) {
          err = 1;
      } else {
          *cnames = '\0';
          for (j=0; j<M->cols && !err; j++) {
            title = NULL;
            col = gtk_tree_view_get_column(GTK_TREE_VIEW(sheet->view), j + 1);
            if (col != NULL) {
                title = gtk_tree_view_column_get_title(col);
            }
            if (title != NULL) {
                *tmp = '\0';
                single_underscores(tmp, title);
                strcat(cnames, tmp);
                strcat(cnames, " ");
            } else {
                err = 1;
            }
          }
          if (!err) {
            umatrix_set_colnames_from_string(M, cnames);
          }
          free(cnames);
      }

      sheet->flags &= ~SHEET_COLNAME_MOD;
    }
}

/* in case we're editing a pre-existing matrix, carry the
   modifications back */

static void update_saved_matrix (Spreadsheet *sheet)
{
    if (sheet->oldmat != NULL) {
      if (sheet->oldmat->rows == sheet->matrix->rows &&
          sheet->oldmat->cols == sheet->matrix->cols) {
          gretl_matrix_copy_values(sheet->oldmat, sheet->matrix);
      } else if (sheet->oldmat->rows * sheet->oldmat->cols ==
               sheet->matrix->rows * sheet->matrix->cols) {
          sheet->oldmat->rows = sheet->matrix->rows;
          sheet->oldmat->cols = sheet->matrix->cols;
          gretl_matrix_copy_values(sheet->oldmat, sheet->matrix);
      } else {
          gretl_matrix *m = gretl_matrix_copy(sheet->matrix);

          if (m == NULL) {
            nomem();
          } else {
            user_matrix_replace_matrix_by_name(sheet->mname, m);
            sheet->oldmat = m;
          }
      }

      maybe_update_column_names(sheet);

      /* record the fact that a matrix has been changed */
      mark_session_changed();
    }
}

static void update_sheet_matrix_element (Spreadsheet *sheet,
                               const gchar *new_text,
                               const gchar *path_string,
                               int colnum)
{
    int i = atoi(path_string);
    int j = colnum - 1;

    gretl_matrix_set(sheet->matrix, i, j, atof(new_text));
}

static void 
maybe_update_store (Spreadsheet *sheet, const gchar *new_text,
                const gchar *path_string)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeModel *model = gtk_tree_view_get_model(view);
    GtkTreeViewColumn *col = NULL;
    GtkTreePath *path;
    GtkTreeIter iter;
    gchar *old_text;
    gint colnum;

    gtk_tree_view_get_cursor(view, NULL, &col);
    if (col == NULL) {
      return;
    }

    colnum = GPOINTER_TO_INT(g_object_get_data(G_OBJECT(col), "colnum"));
    path = gtk_tree_path_new_from_string(path_string);
    gtk_tree_model_get_iter(model, &iter, path);
    gtk_tree_model_get(model, &iter, colnum, &old_text, -1);

#if CELLDEBUG
    fprintf(stderr, "maybe_update_store: old='%s', new='%s'\n",
          old_text, new_text);
#endif    

    if (old_text != NULL && strcmp(old_text, new_text)) {
      gtk_list_store_set(GTK_LIST_STORE(model), &iter, 
                     colnum, new_text, -1);
      if (sheet->matrix != NULL) {
          update_sheet_matrix_element(sheet, new_text, path_string, colnum);
      }
      sheet_set_modified(sheet, TRUE);
    }

    move_to_next_cell(sheet, path, col);
    gtk_tree_path_free(path);
    g_free(old_text);
}

static int scalar_try_genr (gchar **ps)
{
    double x;
    int err = 0;

    x = generate_scalar(*ps + 1, &Z, datainfo, &err);

    if (!err) {
      g_free(*ps);
      *ps = g_strdup_printf("%.*g", DBL_DIG, x);
    }

    return err;
}

static void sheet_cell_edited (GtkCellRendererText *cell,
                         const gchar *path_string,
                         const gchar *user_text,
                         Spreadsheet *sheet)
{
    gchar *new_text = NULL;
    int err = 0;

#if CELLDEBUG
    fprintf(stderr, "*** sheet_cell_edited\n");
#endif

    if (sheet->matrix == NULL && 
      (!strcmp(user_text, "na") || !strcmp(user_text, "NA"))) {
      /* allow conversion to "missing" */
      new_text = g_strdup("");
    } else {
      new_text = g_strdup(user_text);
      if (*new_text == '=' && sheet->cmd == SHEET_EDIT_SCALARS) {
          err = scalar_try_genr(&new_text);
      } else {
          if (sheet->flags & SHEET_USE_COMMA) {
            /* accept point also: convert to locale */
            charsub(new_text, '.', ',');
          }
          err = check_atof(new_text);
      }
      if (err) {
          errbox(gretl_errmsg_get());
          g_free(new_text);
      } 
    }

    if (!err && new_text != NULL) {
      maybe_update_store(sheet, new_text, path_string);
      g_free(new_text);
    }
}

static void 
spreadsheet_scroll_to_new_col (Spreadsheet *sheet, GtkTreeViewColumn *column)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreePath *path;
    gchar *pstr;
    GtkAdjustment *adj;
    GtkWidget *sw;

    pstr = g_strdup("0");
    path = gtk_tree_path_new_from_string(pstr);
    gtk_tree_view_set_cursor(view, path, column, TRUE); /* ?? */
    sw = gtk_widget_get_ancestor(GTK_WIDGET(view), GTK_TYPE_BIN);
    if (sw != NULL) {
      adj = gtk_scrolled_window_get_hadjustment(GTK_SCROLLED_WINDOW(sw));
      gtk_adjustment_set_value(adj, adj->upper);
    }
    gtk_tree_path_free(path);
    g_free(pstr);
}

static GtkTreeViewColumn *
add_treeview_column_with_title (Spreadsheet *sheet, const char *name)
{
    GtkTreeViewColumn *column;
    gint cols, colnum;

    column = gtk_tree_view_column_new();
    gtk_tree_view_column_set_title(column, name);
    set_up_sheet_column(column, get_data_col_width(), TRUE);

    cols = gtk_tree_view_insert_column(GTK_TREE_VIEW(sheet->view), column, -1);
    colnum = cols - 1;

    gtk_tree_view_column_pack_start(column, sheet->datacell, TRUE);
    gtk_tree_view_column_set_attributes(column, 
                              sheet->datacell,
                              "text", colnum,
                              NULL);
    g_object_set_data(G_OBJECT(column), "colnum", GINT_TO_POINTER(colnum));

    return column;
}

static int real_add_new_series (Spreadsheet *sheet, const char *varname)
{
    GtkTreeViewColumn *column;
    char tmp[32];

    if (add_data_column(sheet)) {
      return 1;
    }

#if SSDEBUG
    fprintf(stderr, "real_add_new_series, after add_data_column: sheet->totcols=%d\n", 
          sheet->totcols);
#endif

    double_underscores(tmp, varname);
    column = add_treeview_column_with_title(sheet, tmp);

    /* scroll to editing position if need be */
    spreadsheet_scroll_to_new_col(sheet, column);

    sheet_set_modified(sheet, TRUE);

    return 0;
}

static int real_add_new_scalar (Spreadsheet *sheet, const char *vname,
                        double val)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkListStore *store;
    GtkTreeIter iter;

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));
    gtk_list_store_append(store, &iter);

    if (na(val)) {
      gtk_list_store_set(store, &iter, 0, vname, 1, "", -1);
    } else {
      gchar *numstr;

      numstr = g_strdup_printf("%.*g", DBL_DIG, val);
      gtk_list_store_set(store, &iter, 0, vname, 1, numstr, 
                     2, sheet->pbuf, -1);
      g_free(numstr);
    }

    sheet->datarows += 1;
    spreadsheet_scroll_to_foot(sheet, sheet->datarows - 1, 1);

    sheet_set_modified(sheet, TRUE);

    return 0;
}

static void 
spreadsheet_scroll_to_foot (Spreadsheet *sheet, int row, int col)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreePath *path;
    gchar *pstr;
    GtkTreeViewColumn *column;
    GtkAdjustment *adj;
    GtkWidget *sw;

    pstr = g_strdup_printf("%d", row);
    path = gtk_tree_path_new_from_string(pstr);
    column = gtk_tree_view_get_column(view, col);
    gtk_tree_view_set_cursor(view, path, column, FALSE);
    sw = gtk_widget_get_ancestor(GTK_WIDGET(view), GTK_TYPE_BIN);
    if (sw != NULL) {
      adj = gtk_scrolled_window_get_vadjustment(GTK_SCROLLED_WINDOW(sw));
      gtk_adjustment_set_value(adj, adj->upper);
    }
    gtk_tree_path_free(path);
    g_free(pstr);
}

static void 
real_add_new_obs (Spreadsheet *sheet, const char *obsname, int n)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    gint rownum = 0;
    gint oldrows = sheet->datarows;
    GtkListStore *store;
    gchar *pstr = NULL;
    GtkTreeIter iter;
    gchar rowlabel[10];
    gint i, j;

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));

    if (sheet->point == SHEET_AT_END) {
      rownum = sheet->datarows - 1;
      pstr = g_strdup_printf("%d", rownum);
      for (i=0; i<n; i++) {
          gtk_list_store_append(store, &iter);
      }
    } else if (sheet->point == SHEET_AT_POINT) {
      GtkTreePath *path;

      gtk_tree_view_get_cursor(view, &path, NULL);
      if (path != NULL) {
          rownum = gtk_tree_path_get_indices(path)[0];
          gtk_tree_model_get_iter(GTK_TREE_MODEL(store), &iter, path);
          gtk_list_store_insert(store, &iter, rownum);
          gtk_tree_path_free(path);
      }
    } else {
      return;
    }

    if (datainfo->markers && obsname != NULL) {
      gtk_list_store_set(store, &iter, 0, obsname, -1);
    } else if (sheet->point == SHEET_AT_END) {
      gtk_tree_model_get_iter_from_string(GTK_TREE_MODEL(store), &iter,
                                  pstr);
      for (j=0; j<n; j++) {
          gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
          get_obs_string(rowlabel, sheet->datarows + j, datainfo);
          gtk_list_store_set(store, &iter, 0, rowlabel, -1);
      }
    }

    sheet->datarows += n;

    if (pstr != NULL) {
      gtk_tree_model_get_iter_from_string(GTK_TREE_MODEL(store), &iter,
                                  pstr);
      for (j=0; j<n; j++) {
          gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
          for (i=1; i<=sheet->datacols; i++) {
            gtk_list_store_set(store, &iter, i, "", -1);
          } 
      }
    } else {
      for (i=1; i<=sheet->datacols; i++) {
          gtk_list_store_set(store, &iter, i, "", -1);
      }
    } 

    if (sheet->point == SHEET_AT_POINT && !datainfo->markers) {
      for (i=rownum; i<sheet->datarows; i++) {
          get_obs_string(rowlabel, i, datainfo);
          gtk_list_store_set(store, &iter, 0, rowlabel, -1);
          gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
      }
    } 

    if (sheet->point == SHEET_AT_END) {
      spreadsheet_scroll_to_foot(sheet, oldrows, 1);
    } else {
      GtkTreePath *path;
      GtkTreeIter insiter;

      pstr = g_strdup_printf("%d", rownum);
      path = gtk_tree_path_new_from_string(pstr);
      gtk_tree_model_get_iter(GTK_TREE_MODEL(store), &insiter, path);
      gtk_list_store_set(store, &insiter, 1, "", -1);
      gtk_tree_path_free(path);
    }

    if (pstr != NULL) {
      g_free(pstr);
    }

    sheet_set_modified(sheet, TRUE);
}

static void add_new_scalar (GtkWidget *widget, dialog_t *dlg) 
{
    Spreadsheet *sheet = (Spreadsheet *) edit_dialog_get_data(dlg);
    const gchar *buf;
    gchar *vname = NULL;
    gchar *expr = NULL;
    char varname[VNAMELEN];
    double val = NADBL;
    int err = 0;

    buf = edit_dialog_get_text(dlg);

    if (buf == NULL || *buf == '\0') {
      return;
    }

    if (strchr(buf, '=')) {
      gchar *p, *tmp = g_strdup(buf);

      p = strchr(tmp, '=');
      *p = '\0';
      tailstrip(tmp);
      vname = g_strdup(tmp);
      expr = g_strdup(p + 1);
      g_free(tmp);

      if (gui_validate_varname(vname, GRETL_TYPE_DOUBLE)) {
          err = 1;
      } else {
          *varname = 0;
          strncat(varname, vname, VNAMELEN - 1);
          val = generate_scalar(expr, &Z, datainfo, &err);
      }
    } else if (gui_validate_varname(buf, GRETL_TYPE_DOUBLE)) {
      return;
    } else {
      *varname = 0;
      strncat(varname, buf, VNAMELEN - 1);
    }

    g_free(vname);
    g_free(expr);

    if (!err) {
      close_dialog(dlg);
      if (real_add_new_scalar(sheet, varname, val)) {
          nomem();
      }
    }
}

static void name_new_var (GtkWidget *widget, dialog_t *dlg) 
{
    Spreadsheet *sheet = (Spreadsheet *) edit_dialog_get_data(dlg);
    const gchar *buf;
    char varname[VNAMELEN];

    buf = edit_dialog_get_text(dlg);

    if (buf == NULL || gui_validate_varname(buf, GRETL_TYPE_SERIES)) {
      return;
    }

    *varname = 0;
    strncat(varname, buf, VNAMELEN - 1);

    close_dialog(dlg);

    if (real_add_new_series(sheet, varname)) {
      nomem();
    }
}

static void name_new_obs (GtkWidget *widget, dialog_t *dlg) 
{
    Spreadsheet *sheet = (Spreadsheet *) edit_dialog_get_data(dlg);
    const gchar *buf;
    char obsmarker[OBSLEN];

    buf = edit_dialog_get_text(dlg);
    if (buf == NULL) return;

    *obsmarker = 0;
    strncat(obsmarker, buf, OBSLEN - 1);

    close_dialog(dlg);
    real_add_new_obs(sheet, obsmarker, 1);
}

static void name_var_dialog (Spreadsheet *sheet) 
{
    int cancel = 0;
    
    edit_dialog(_("gretl: name variable"), 
            _("Enter name for new variable\n"
              "(max. 15 characters)"),
            NULL, name_new_var, sheet, 
            0, VARCLICK_NONE, &cancel);
}

static void new_scalar_dialog (Spreadsheet *sheet) 
{
    int cancel = 0;
    
    edit_dialog(_("gretl: scalar"), 
            _("Enter name, or name=value,\n"
              "for new scalar"),
            NULL, add_new_scalar, sheet, 
            0, VARCLICK_NONE, &cancel);
}

static void new_case_dialog (Spreadsheet *sheet) 
{
    int cancel = 0;

    edit_dialog(_("gretl: case marker"), 
            _("Enter case marker for new obs\n"
              "(max. 8 characters)"),
            NULL, name_new_obs, sheet, 
            0, VARCLICK_NONE, &cancel);
}

static void name_matrix_col (GtkWidget *widget, dialog_t *dlg) 
{
    GtkTreeViewColumn *col = (GtkTreeViewColumn *) edit_dialog_get_data(dlg);
    const gchar *buf, *old;
    char tmp[13], colname[24];

    buf = edit_dialog_get_text(dlg);
    if (buf == NULL || gui_validate_varname(buf, GRETL_TYPE_NONE)) {
      return;
    }

    *tmp = 0;
    strncat(tmp, buf, 12);
    double_underscores(colname, tmp);

    close_dialog(dlg);

    old = gtk_tree_view_column_get_title(col);

    if (strcmp(old, colname)) {
      Spreadsheet *sheet = g_object_get_data(G_OBJECT(col), "sheet");

      gtk_tree_view_column_set_title(col, colname);
      sheet->flags |= SHEET_COLNAME_MOD;
      sheet_set_modified(sheet, TRUE);
    }
}

static void name_column_dialog (GtkTreeViewColumn *col, gpointer p) 
{
    int cancel = 0;

    edit_dialog(_("gretl: name column"), 
            _("Enter name for column\n"
              "(max. 12 characters)"),
            gtk_tree_view_column_get_title(col),
            name_matrix_col, col, 
            0, VARCLICK_NONE, &cancel);
}

static GtkListStore *make_sheet_liststore (Spreadsheet *sheet)
{
    GtkListStore *store;
    GType *types;
    int ncols;
    int i;

    /* obs col, data cols */
    ncols = sheet->totcols = sheet->datacols + 1;

    if (sheet->cmd == SHEET_EDIT_SCALARS) {
      ncols++;
    }

    types = mymalloc(ncols * sizeof *types);
    if (types == NULL) {
      return NULL;
    }

    for (i=0; i<sheet->totcols; i++) {
      types[i] = G_TYPE_STRING;
    }

    if (sheet->cmd == SHEET_EDIT_SCALARS) {
      types[i] = GDK_TYPE_PIXBUF;
    }

    store = gtk_list_store_newv(ncols, types);
    free(types);

    return store;
}

/* This is relatively complex because, so far as I can tell, you can't
   append or insert additional columns in a GtkListStore: we have to
   create a whole new liststore and copy the old info across.
*/

static int add_data_column (Spreadsheet *sheet)
{
    GtkListStore *old_store, *new_store;
    GtkTreeIter old_iter, new_iter;
    gint row;

    sheet->datacols += 1;

    /* get pointers to original and new stores */
    new_store = make_sheet_liststore(sheet);
    if (new_store == NULL) {
      sheet->datacols -= 1;
      sheet->totcols -= 1;
      return E_ALLOC;
    }

    old_store = GTK_LIST_STORE(gtk_tree_view_get_model(GTK_TREE_VIEW(sheet->view)));

    /* go to start of old and new lists */
    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(old_store), &old_iter);
    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(new_store), &new_iter);

    /* construct the new table */
    for (row=0; row<sheet->datarows; row++) {
      gchar *str;
      int col;

      gtk_list_store_append(new_store, &new_iter);

      for (col=0; col<sheet->totcols; col++) {
          if (col < sheet->datacols) {
            /* copy labels and original data */
            gtk_tree_model_get(GTK_TREE_MODEL(old_store), &old_iter, col, &str, -1);
            gtk_list_store_set(new_store, &new_iter, col, str, -1);
            g_free(str);
          } else { 
            /* new data values: set blank */
            gtk_list_store_set(new_store, &new_iter, col, "", -1);
          }
      }
      gtk_tree_model_iter_next(GTK_TREE_MODEL(old_store), &old_iter);
    }    

    gtk_tree_view_set_model(GTK_TREE_VIEW(sheet->view), GTK_TREE_MODEL(new_store));
    g_object_unref(G_OBJECT(new_store));

    sheet->added_vars += 1;

#if SSDEBUG
    fprintf(stderr, "add_data_column: sheet->added_vars now = %d\n",
          sheet->added_vars);
#endif

    return 0;
}

static void sheet_get_scalar (GtkWidget *w, dialog_t *dlg)
{
    double x, *px = (double *) edit_dialog_get_data(dlg);
    const gchar *buf;
    int err = 0;

    buf = edit_dialog_get_text(dlg);
    if (buf == NULL) return;

    x = gui_double_from_string(buf, &err);
    if (!err) {
      *px = x;
      close_dialog(dlg);
    }
}

static void matrix_edit_callback (GtkAction *action, gpointer data)
{
    const gchar *s = gtk_action_get_name(action);
    Spreadsheet *sheet = (Spreadsheet *) data;
    double x = NADBL;
    int err = 0;

    if (!strcmp(s, "ScalarMult") || !strcmp(s, "ScalarDiv")) {
      int cancel = 0;

      edit_dialog(_("gretl: specify scalar"), 
                _("Enter a numerical value"),
                NULL, sheet_get_scalar, &x, 
                0, VARCLICK_NONE, &cancel);
      if (cancel || na(x)) {
          return;
      }
    }

    if (!strcmp(s, "XTX")) {
      err = matrix_XTX_in_place(sheet->matrix);
    } else if (!strcmp(s, "Transpose")) {
      err = matrix_transpose_in_place(sheet->matrix);
    } else if (!strcmp(s, "Cholesky")) {
      err = matrix_cholesky_in_place(sheet->matrix);
    } else if (!strcmp(s, "Invert")) {
      err = matrix_invert_in_place(sheet->matrix);
    } else if (!strcmp(s, "ScalarMult")) {
      gretl_matrix_multiply_by_scalar(sheet->matrix, x);
    } else if (!strcmp(s, "ScalarDiv")) {
      gretl_matrix_divide_by_scalar(sheet->matrix, x);
    }

    if (err) {
      gui_errmsg(err);
    } else {
      update_sheet_from_matrix(sheet);
    }
}

static void matrix_props_callback (GtkAction *action, gpointer data)
{
    Spreadsheet *sheet = (Spreadsheet *) data;

    view_matrix_properties(sheet->matrix, sheet->mname);
}

static void matrix_fill_callback (GtkAction *action, gpointer data)
{
    const gchar *s = gtk_action_get_name(action);
    Spreadsheet *sheet = (Spreadsheet *) data;
    gretl_matrix *A = sheet->matrix;
    int mindim;

    if (!strcmp(s, "FillIdentity")) {
      mindim = (A->rows < A->cols)? A->rows : A->cols;
      gretl_matrix_inscribe_I(A, 0, 0, mindim);
    } else if (!strcmp(s, "FillUniform")) {
      gretl_matrix_random_fill(A, D_UNIFORM);
    } else if (!strcmp(s, "FillNormal")) {
      gretl_matrix_random_fill(A, D_NORMAL);
    }

    update_sheet_from_matrix(sheet);
}

static void sheet_add_obs_direct (Spreadsheet *sheet)
{
    if (datainfo->markers) {
      new_case_dialog(sheet);
    } else if (sheet->point == SHEET_AT_END) {
      int n = add_obs_dialog(NULL, 1);

      if (n > 0) {
          real_add_new_obs(sheet, NULL, n);
      }
    } else {
      real_add_new_obs(sheet, NULL, 1);
    }
}

static void popup_sheet_add_obs (GtkWidget *w, Spreadsheet *sheet)
{
    sheet->point = SHEET_AT_END;
    sheet_add_obs_direct(sheet);
}

static void popup_sheet_insert_obs (GtkWidget *w, Spreadsheet *sheet)
{
    sheet->point = SHEET_AT_POINT;
    sheet_add_obs_direct(sheet);
}

static void popup_sheet_add_var (GtkWidget *w, Spreadsheet *sheet)
{
    name_var_dialog(sheet);
}

static void popup_sheet_add_scalar (GtkWidget *w, Spreadsheet *sheet)
{
    new_scalar_dialog(sheet);
}

static void sheet_delete_scalar (Spreadsheet *sheet, GtkTreePath *path)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeIter iter;
    GtkListStore *store;

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));
    gtk_tree_model_get_iter(GTK_TREE_MODEL(store), &iter, path);
    gtk_list_store_remove(store, &iter);
    sheet->datarows--;
    sheet_set_modified(sheet, TRUE);
}

static void build_sheet_popup (Spreadsheet *sheet)
{
    sheet->popup = gtk_menu_new();

    if (sheet->cmd == SHEET_EDIT_SCALARS) {
      add_popup_item(_("Add..."), 
                   sheet->popup, 
                   G_CALLBACK(popup_sheet_add_scalar),
                   sheet);
    } else {
      add_popup_item(_("Add Variable"), sheet->popup, 
                   G_CALLBACK(popup_sheet_add_var),
                   sheet);
    }

    if (sheet->flags & SHEET_ADD_OBS_OK) {
      add_popup_item(_("Add Observation"), sheet->popup,
                   G_CALLBACK(popup_sheet_add_obs),
                   sheet);
    } 

    if (sheet->flags & SHEET_INSERT_OBS_OK) {
      add_popup_item(_("Insert Observation"), sheet->popup,
                   G_CALLBACK(popup_sheet_insert_obs),
                   sheet);
    }
}

/* this is connected to the "cursor-changed" signal
   on the treeview */

static void update_cell_position (GtkTreeView *view, 
                          Spreadsheet *sheet)
{
    GtkTreePath *path = NULL;
    GtkTreeViewColumn *col = NULL;
    static int i0, j0;

#if CELLDEBUG > 1
    fprintf(stderr, "*** cursor-changed\n");
#endif

    gtk_tree_view_get_cursor(view, &path, &col);

    if (path != NULL && col != NULL) {
      int i = gtk_tree_path_get_indices(path)[0];
      int j = 
          GPOINTER_TO_INT(g_object_get_data(G_OBJECT(col), "colnum"));

      if (j > 0 && (i != i0 || j != j0)) {
#if CELLDEBUG > 1
          fprintf(stderr, " now in cell(%d, %d)\n", i, j);
#endif
          set_locator_label(sheet, path, col);
          i0 = i;
          j0 = j;
      } else {
#if CELLDEBUG > 1
         fprintf(stderr, " still in cell(%d, %d)\n", i0, j0); 
#endif
      }
    }

    if (path != NULL) {
      gtk_tree_path_free(path);
    }
}

/* put modified values from the spreadsheet into the attached
   matrix */

static void update_matrix_from_sheet_full (Spreadsheet *sheet)
{
    update_sheet_matrix(sheet);

    /* in case we're editing a pre-existing matrix, carry the
       modifications back */
    if (sheet->oldmat != NULL) {
      update_saved_matrix(sheet);
    }

    sheet_set_modified(sheet, FALSE);
    set_ok_transforms(sheet);
}

/* put modified values from the spreadsheet into the array
   of saved scalars */

static void update_scalars_from_sheet (Spreadsheet *sheet)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeModel *model;
    GtkTreeIter iter;
    gchar *vname, *val;
    const char *sname;
    double x;
    int i, j, ns;
    int err = 0;

    model = gtk_tree_view_get_model(view);
    gtk_tree_model_get_iter_first(model, &iter);

    /* deal with added or modified scalars fist */
    for (i=0; i<sheet->datarows && !err; i++) {
      gtk_tree_model_get(model, &iter, 0, &vname, 1, &val, -1);
      x = (*val == '\0')? NADBL : atof(val);
      if (gretl_is_scalar(vname)) {
          gretl_scalar_set_value(vname, x);
      } else {
          err = gretl_scalar_add(vname, x);
      }
      g_free(vname);
      g_free(val);
      gtk_tree_model_iter_next(model, &iter);
    }

    ns = n_saved_scalars();

    if (ns > sheet->datarows) {
      /* must have deleted some scalar(s) via GUI */
      int match, n = 0, ndel = ns - sheet->datarows;

      i = 0;
      while ((sname = gretl_scalar_get_name(i)) && n < ndel) {
          gtk_tree_model_get_iter_first(model, &iter);
          match = 0;
          for (j=0; j<sheet->datarows && !match; j++) {
            gtk_tree_model_get(model, &iter, 0, &vname, -1);
            match = !strcmp(sname, vname);
            g_free(vname);
          }
          if (!match) {
            gretl_scalar_delete(sname, NULL);
            n++;
          } else {
            i++;
          }
      }
    }

    if (err) {
      gui_errmsg(err);
    }

    sheet_set_modified(sheet, FALSE);
}

/* pull modified values from the data-editing spreadsheet
   into the main dataset */

static void update_dataset_from_sheet (Spreadsheet *sheet)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeIter iter;
    GtkTreeViewColumn *column;
    GtkTreeModel *model;
    int oldv = datainfo->v;
    int newvars = sheet->added_vars;
    int newobs = sheet->datarows - sheet->orig_nobs;
    int missobs = 0;
    int i, colnum, s, t;

#if SSDEBUG
    fprintf(stderr, "get_data_from_sheet: oldv = %d, newvars = %d, newobs = %d\n",
          oldv, newvars, newobs);
#endif

    model = gtk_tree_view_get_model(view);

    /* first extend series length, if needed */

    if (newobs > 0) {
      if (dataset_add_observations(newobs, &Z, datainfo, OPT_A) ||
          dataset_destroy_hidden_variables(&Z, datainfo, 0)) {
          nomem();
          return;
      }
    }

    /* then add any new variables to data set */

    if (newvars > 0) {
      int vi, v0 = sheet->varlist[0];
      const gchar *newname;

      if (dataset_add_series(newvars, &Z, datainfo)) {
          nomem();
          return;
      }

      for (i=0; i<newvars; i++) { 
          newname = NULL;
          vi = oldv + i;
          gretl_list_append_term(&sheet->varlist, vi);
          if (sheet->varlist == NULL) {
            nomem();
            return;
          }
          colnum = v0 + 1 + i;
          column = gtk_tree_view_get_column(view, colnum);
          newname = gtk_tree_view_column_get_title(column);
          strcpy(datainfo->varname[vi], newname); 
#if SSDEBUG
          fprintf(stderr, " added var %d (%s) from column %d\n",
                vi, newname, colnum);
#endif
      }
    }

    /* copy data values from spreadsheet */

    colnum = 1;
    for (i=1; i<=sheet->varlist[0]; i++) {
      int vi = sheet->varlist[i];
      gchar *numstr;

      gtk_tree_model_get_iter_first(model, &iter);

#if SSDEBUG
      fprintf(stderr, " updating data for var %d (%s) from column %d\n",
            vi, datainfo->varname[vi], colnum);
#endif
      t = datainfo->t1;
      for (s=0; s<sheet->datarows; s++) {
          gtk_tree_model_get(model, &iter, colnum, &numstr, -1);
          if (*numstr != '\0') {
            Z[vi][t++] = atof(numstr); 
          } else {
            Z[vi][t++] = NADBL;
            missobs = 1;
          }
          g_free(numstr);
          gtk_tree_model_iter_next(model, &iter);
      }
      colnum++;
    }

    /* copy observation markers, if relevant */

    if (datainfo->markers && datainfo->S != NULL) {
      gchar *marker;

      gtk_tree_model_get_iter_first(model, &iter);
      t = datainfo->t1;
      for (s=0; s<sheet->datarows; s++) {
          gtk_tree_model_get(model, &iter, 0, &marker, -1);
          strcpy(datainfo->S[t++], marker);
          g_free(marker);
          gtk_tree_model_iter_next(model, &iter);
      }
    }

    register_data(DATA_APPENDED);

    if (missobs) {
      infobox(_("Warning: there were missing observations"));
    } 

    sheet_set_modified(sheet, FALSE);
    sheet->added_vars -= newvars; /* record that these are handled */
}

static void matrix_new_name (GtkWidget *w, dialog_t *dlg)
{
    char *newname = (char *) edit_dialog_get_data(dlg);
    const gchar *buf = edit_dialog_get_text(dlg);

    if (buf == NULL || gui_validate_varname(buf, GRETL_TYPE_MATRIX)) {
      return;
    }

    *newname = 0;
    strncat(newname, buf, VNAMELEN - 1);

    close_dialog(dlg);
}

static void matrix_save_as (GtkWidget *w, Spreadsheet *sheet)
{
    char newname[VNAMELEN];
    int cancel = 0;

    edit_dialog(_("gretl: save matrix"), 
            _("Enter a name"),
            NULL, matrix_new_name, newname, 
            0, VARCLICK_NONE, &cancel);
    
    if (!cancel) {
      gretl_matrix *m;
      user_matrix *u;
      gchar *tmp;

      m = gretl_matrix_copy(sheet->matrix);
      add_or_replace_user_matrix(m, newname);
      strcpy(sheet->mname, newname);

      tmp = g_strdup_printf("gretl: %s", sheet->mname);
      gtk_window_set_title(GTK_WINDOW(sheet->win), tmp);
      g_free(tmp);

      u = get_user_matrix_by_name(newname);
      g_object_set_data(G_OBJECT(sheet->win), "object", u);

      sheet->oldmat = m;
      maybe_update_column_names(sheet);
      sheet_set_modified(sheet, FALSE);
    }
}

#define new_matrix(s) (s->cmd == SHEET_EDIT_MATRIX && \
                       s->oldmat == NULL)

static void get_data_from_sheet (GtkWidget *w, Spreadsheet *sheet)
{
    if (!sheet_is_modified(sheet) && !new_matrix(sheet)) {
      infobox(_("No changes were made"));
    } else if (sheet->cmd == SHEET_EDIT_MATRIX) {
      update_matrix_from_sheet_full(sheet);
    } else if (sheet->cmd == SHEET_EDIT_SCALARS) {
      update_scalars_from_sheet(sheet);
    } else {
      update_dataset_from_sheet(sheet);
    }
}

static void select_first_editable_cell (Spreadsheet *sheet)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreePath *path;
    GtkTreeViewColumn *column;

    path = gtk_tree_path_new_from_string("0");
    column = gtk_tree_view_get_column(view, 1);
    gtk_tree_view_set_cursor(view, path, column, FALSE);
    if (sheet->locator != NULL) {
      set_locator_label(sheet, path, column);
    }

    gtk_tree_path_free(path);
}

static void set_ok_transforms (Spreadsheet *sheet)
{
    int z = gretl_is_zero_matrix(sheet->matrix);
    int s = gretl_matrix_get_structure(sheet->matrix);

    flip(sheet->ui, "/menubar/Transform/ScalarMult", !z);
    flip(sheet->ui, "/menubar/Transform/ScalarDiv", !z);
    flip(sheet->ui, "/menubar/Transform/XTX", 
       s == 0 || (!z && s != GRETL_MATRIX_IDENTITY));
    flip(sheet->ui, "/menubar/Transform/Cholesky", 
       s == GRETL_MATRIX_SYMMETRIC && !z);
    flip(sheet->ui, "/menubar/Transform/Invert", s > 0 && !z &&
       s != GRETL_MATRIX_IDENTITY);
    flip(sheet->ui, "/menubar/Transform/Transpose", 
       s < GRETL_MATRIX_SYMMETRIC);
}

static int rejig_sheet_cols (Spreadsheet *sheet)
{
    int n = sheet->matrix->cols - sheet->datacols;
    GtkTreeViewColumn *col;
    GtkListStore *store;
    char cstr[16];
    int i;

    sheet->datacols = sheet->matrix->cols;

    store = make_sheet_liststore(sheet);
    if (store == NULL) {
      return E_ALLOC;
    }

    gtk_tree_view_set_model(GTK_TREE_VIEW(sheet->view), GTK_TREE_MODEL(store));
    g_object_unref(G_OBJECT(store));

    if (n > 0) {
      /* add extra columns */
      for (i=0; i<n; i++) {
          sprintf(cstr, "%d", sheet->datacols - n + i + 1);
          add_treeview_column_with_title(sheet, cstr);
      }
    } else {
      /* remove surplus columns */
      n = -n;
      for (i=0; i<n; i++) {
          col = gtk_tree_view_get_column(GTK_TREE_VIEW(sheet->view),
                                 sheet->datacols + n - i);
          gtk_tree_view_remove_column(GTK_TREE_VIEW(sheet->view),
                              col);
      }
    }

    for (i=1; i<=sheet->matrix->cols; i++) {
      sprintf(cstr, "%d", i);
      col = gtk_tree_view_get_column(GTK_TREE_VIEW(sheet->view), i);
      gtk_tree_view_column_set_title(col, cstr);
      gtk_tree_view_column_set_clickable(col, TRUE);
      g_object_set_data(G_OBJECT(col), "sheet", sheet);
      g_signal_connect(G_OBJECT(col), "clicked",
                   G_CALLBACK(name_column_dialog), col);
    }

    return 0;
}

static int update_sheet_from_matrix (Spreadsheet *sheet)
{
    gchar tmpstr[48];
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkListStore *store;
    GtkTreeIter iter;
    int i, j, resized = 0;
    double x;
    int err = 0;

    if (sheet->matrix->cols != sheet->datacols) {
      err = rejig_sheet_cols(sheet);
      if (err) {
          return err;
      }
      resized = 1;
    }

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));
    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(store), &iter);

    for (i=0; i<sheet->matrix->rows; i++) {
      int i1 = i + 1;

      if (resized) {
          gtk_list_store_append(store, &iter);
          sprintf(tmpstr, "%d", i1);
          gtk_list_store_set(store, &iter, 0, tmpstr, -1);
      }
      for (j=0; j<sheet->matrix->cols; j++) {
          x = gretl_matrix_get(sheet->matrix, i, j);
          if (x == -0.0) {
            x = 0.0;
          }
          sprintf(tmpstr, "%.*g", MATRIX_DIGITS, x);
          gtk_list_store_set(store, &iter, j + 1, tmpstr, -1);
      }
      if (i1 < sheet->matrix->rows) {
          /* there are more rows to handle */
          if (resized) {
            /* will be handled above on next iteration */
            continue;
          } else if (i1 >= sheet->datarows) {
            /* need to add another row */
            gtk_list_store_append(store, &iter);
            sprintf(tmpstr, "%d", i1 + 1);
            gtk_list_store_set(store, &iter, 0, tmpstr, -1);
          } else {
            /* row already present */
            gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
          } 
      }
    }

    if (!resized && (sheet->datarows > sheet->matrix->rows)) {
      int rdel = sheet->datarows - sheet->matrix->rows;

      gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
      for (i=0; i<rdel; i++) {
          gtk_list_store_remove(store, &iter);
      }
    }

    if (sheet->datarows != sheet->matrix->rows) {
      sheet->datarows = sheet->matrix->rows;
      resized = 1;
    } 

    sheet_set_modified(sheet, TRUE);
    set_ok_transforms(sheet);

    if (resized) {
      size_matrix_window(sheet);
    }

    return 0;
}

static int add_matrix_data_to_sheet (Spreadsheet *sheet)
{
    gchar tmpstr[32];
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeIter iter;
    GtkListStore *store;
    double x;
    int i, j;

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));

    /* row labels */
    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(store), &iter);
    for (i=0; i<sheet->datarows; i++) {
      sprintf(tmpstr, "%d", i+1);
      gtk_list_store_append(store, &iter);
      gtk_list_store_set(store, &iter, 0, tmpstr, -1);
    }

    /* now insert data values */
    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(store), &iter);

    for (i=0; i<sheet->datarows; i++) {
      for (j=0; j<sheet->datacols; j++) {
          x = gretl_matrix_get(sheet->matrix, i, j);
          sprintf(tmpstr, "%.*g", DBL_DIG, x);
          gtk_list_store_set(store, &iter, j + 1, tmpstr, -1);
      }
      gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
    }

    set_ok_transforms(sheet);

    return 0;
}

static int add_scalars_to_sheet (Spreadsheet *sheet)
{
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeIter iter;
    GtkListStore *store;
    gchar vname[VNAMELEN];
    gchar val[32];
    double x;
    int i;

    if (sheet->pbuf == NULL) {
      sheet->pbuf = gtk_widget_render_icon(sheet->view, GTK_STOCK_DELETE,
                                   GTK_ICON_SIZE_MENU, NULL);
    }

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));
    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(store), &iter);

    if (n_saved_scalars() == 0) {
      /* starting from scratch */
      gtk_list_store_append(store, &iter);
      gtk_list_store_set(store, &iter, 0, sheet->mname, 1, "", 
                     2, sheet->pbuf, -1);
      return 0;
    }

    /* insert variable names and values */

    for (i=0; i<sheet->datarows; i++) {
      strcpy(vname, gretl_scalar_get_name(i)); /* underscores? */
      x = gretl_scalar_get_value_by_index(i);
      if (na(x)) {
          *val = '\0';
      } else {
          sprintf(val, "%.*g", DBL_DIG, x);
      }
      gtk_list_store_append(store, &iter);
      gtk_list_store_set(store, &iter, 0, vname, 1, val, 
                     2, sheet->pbuf, -1);
    }

    return 0;
}

static int add_data_to_sheet (Spreadsheet *sheet, SheetCmd c)
{
    gchar rowlabel[OBSLEN];
    GtkTreeView *view = GTK_TREE_VIEW(sheet->view);
    GtkTreeIter iter;
    GtkListStore *store;
    int i, t;

#if SSDEBUG
    fprintf(stderr, "Doing add_data_to_sheet\n");
#endif

    store = GTK_LIST_STORE(gtk_tree_view_get_model(view));

    /* insert observation markers */

    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(store), &iter);
    for (t=datainfo->t1; t<=datainfo->t2; t++) {
      get_obs_string(rowlabel, t, datainfo);
      gtk_list_store_append(store, &iter);
      gtk_list_store_set(store, &iter, 0, rowlabel, -1);
    }

    sheet->datarows = datainfo->t2 - datainfo->t1 + 1;

    /* insert data values */

    gtk_tree_model_get_iter_first(GTK_TREE_MODEL(store), &iter);

    for (t=datainfo->t1; t<=datainfo->t2; t++) {
      if (c == SHEET_NEW_DATASET) {
          /* no hidden vars to consider; insert NAs for first var */
          gtk_list_store_set(store, &iter, 1, "", -1);
      } else {
          char numstr[32];
          int colnum = 0;
          int vi;

          for (i=1; i<=sheet->varlist[0]; i++) {
            vi = sheet->varlist[i];
            if (na(Z[vi][t])) {
                *numstr = '\0';
            } else {
                sprintf(numstr, "%.*g", DBL_DIG, Z[vi][t]);
            }
            gtk_list_store_set(store, &iter, ++colnum, numstr, -1);
          }
      }
      gtk_tree_model_iter_next(GTK_TREE_MODEL(store), &iter);
    }

    sheet->orig_main_v = datainfo->v;

#if SSDEBUG
    fprintf(stderr, " datarows=%d, orig_vars=%d, orig_main_v=%d\n", 
          sheet->datarows, sheet->varlist[0], sheet->orig_main_v);
#endif

    return 0;
}

gint get_string_width (const gchar *str)
{
    gint width;

    GtkWidget *w;
    PangoLayout *pl;
    PangoContext *pc;

    w = gtk_label_new(NULL);
    pc = gtk_widget_get_pango_context(w);

    pl = pango_layout_new(pc);
    pango_layout_set_text(pl, str, -1);
    pango_layout_get_pixel_size(pl, &width, NULL);

    gtk_widget_destroy(w);
    g_object_unref(G_OBJECT(pl));

    return width;
}

static gint get_row_label_width (void)
{
    static gint width;

    if (width == 0) {
      width = get_string_width("XXXX");
    }
    return width;
}

static gint get_obs_col_width (void)
{
    static gint width;

    if (width == 0) {
      width = get_string_width("XXXXXXXXX");
    }
    return width;
}

static gint get_name_col_width (void)
{
    static gint width;

    if (width == 0) {
      width = get_string_width("XXXXXXXXXXXXXXX");
    }
    return width;
}

static gint get_data_col_width (void)
{
    static gint width;

    if (width == 0) {
      width = get_string_width("-00.0000000000");
    }
    return width;
}

static gint get_delete_col_width (void)
{
    static gint width;

    if (width == 0) {
      width = get_string_width("XXXXXX");
    }
    return width;
}

static void 
set_up_sheet_column (GtkTreeViewColumn *column, gint width, gboolean expand)
{
    gtk_tree_view_column_set_alignment(column, 0.5); /* header centered */
    gtk_tree_view_column_set_sizing(column, GTK_TREE_VIEW_COLUMN_AUTOSIZE);
    gtk_tree_view_column_set_resizable(column, TRUE);
    gtk_tree_view_column_set_min_width(column, width);
    gtk_tree_view_column_set_expand(column, expand);
}

static void commit_and_move_right (Spreadsheet *sheet,
                           const char *s)
{
    GtkTreePath *path;
    gchar *pathstr;

    gtk_tree_view_get_cursor(GTK_TREE_VIEW(sheet->view), &path, NULL);
    if (path == NULL) {
      return;
    }

    pathstr = gtk_tree_path_to_string(path);

#if CELLDEBUG
    fprintf(stderr, "commit_and_move_right: calling sheet_cell_edited\n");
#endif
    sheet_cell_edited(NULL, pathstr, s, sheet);
    g_free(pathstr);
    gtk_tree_path_free(path);
}

static gboolean
catch_sheet_edit_key (GtkWidget *view, GdkEventKey *key, 
                  Spreadsheet *sheet)
{
    sheet->next = NEXT_DOWN;

    if (key->keyval == GDK_Up) {
#if CELLDEBUG
      fprintf(stderr, "catch_edit_key: GDK_Up\n");
#endif
      sheet->next = NEXT_UP;
    } else if (key->keyval == GDK_Right && sheet->entry != NULL) {
      int n, pos = gtk_editable_get_position(GTK_EDITABLE(sheet->entry));
      const char *s = gtk_entry_get_text(GTK_ENTRY(sheet->entry));

#if CELLDEBUG
      fprintf(stderr, "catch_edit_key: GDK_Right\n");
#endif

      n = (s != NULL)? strlen(s) : 0;
      if (pos == n) {
          sheet->next = NEXT_RIGHT;
          commit_and_move_right(sheet, s);
          return TRUE;
      }
    } else if (key->keyval == GDK_Return) {
#if CELLDEBUG
      fprintf(stderr, "catch_edit_key: GDK_Return\n");
#endif
    } else if (key->keyval == GDK_Down) {
#if CELLDEBUG
      fprintf(stderr, "catch_edit_key: GDK_Down\n");
#endif
    }

    return FALSE;
}

static void nullify_sheet_entry (GtkObject *o, Spreadsheet *sheet)
{
#if CELLDEBUG
    fprintf(stderr, "editing entry destroyed\n");
#endif
    sheet->entry = NULL;
}

static void cell_edit_start (GtkCellRenderer *r,
                       GtkCellEditable *ed,
                       gchar *path,
                       Spreadsheet *sheet)
{
#if CELLDEBUG
    fprintf(stderr, "*** editing-started\n");
#endif
    if (GTK_IS_ENTRY(ed)) {
      sheet->entry = GTK_WIDGET(ed);
      g_signal_connect(G_OBJECT(ed), "key-press-event",
                   G_CALLBACK(catch_sheet_edit_key), sheet);
      g_signal_connect(GTK_OBJECT(ed), "destroy",
                   G_CALLBACK(nullify_sheet_entry), sheet);
    }
}

static void create_sheet_cell_renderers (Spreadsheet *sheet)
{
    GtkCellRenderer *r;

    r = gtk_cell_renderer_text_new();
    g_object_set(r, "ypad", 1, 
             "xalign", 1.0,
             "background", "gray",
             "editable", FALSE, NULL);
    sheet->dumbcell = r;

    r = gtk_cell_renderer_text_new();
    g_object_set(r, "ypad", 1, 
             "xalign", 1.0, 
             "family", "Monospace",
             "editable", TRUE, 
             "mode", GTK_CELL_RENDERER_MODE_EDITABLE,
             NULL);
#if 1
    g_signal_connect(r, "editing-started",
                 G_CALLBACK(cell_edit_start), sheet);
#endif
    g_signal_connect(r, "edited",
                 G_CALLBACK(sheet_cell_edited), sheet);
    sheet->datacell = r;
}

static void manufacture_keystroke (GtkWidget *widget, guint uval)
{
    GdkKeymapKey *keys;
    gint n_keys;

    if (gdk_keymap_get_entries_for_keyval(NULL, uval, &keys, &n_keys)) {
      guint16 hardware_keycode;
      GdkEvent *event;

      hardware_keycode = keys[0].keycode;
      g_free(keys);

      event = gdk_event_new(GDK_KEY_PRESS);
      event->key.window = g_object_ref(widget->window);
      event->key.hardware_keycode = hardware_keycode;

      event->key.keyval = gdk_unicode_to_keyval(uval);
      event->key.length = 1;

      event->key.send_event = FALSE;
      event->key.time = GDK_CURRENT_TIME;   

      gtk_main_do_event(event);
      gdk_event_free(event);
    }
}

static int numeric_key (guint *kval, Spreadsheet *sheet)
{
    if (*kval >= GDK_KP_0 && *kval <= GDK_KP_9) {
      *kval = GDK_0 + (*kval - GDK_KP_0);
    }

    if (*kval >= GDK_0 && *kval <= GDK_9) {
      return 1;
    } else if (*kval == GDK_minus || *kval == GDK_period) {
      return 1;
    } else if ((sheet->flags & SHEET_USE_COMMA) && *kval == GDK_comma) {
      return 1;
    } 

    return 0;
}

static gint catch_spreadsheet_key (GtkWidget *view, GdkEventKey *key, 
                           Spreadsheet *sheet)
{
    guint kval = key->keyval;

#if CELLDEBUG
    fprintf(stderr, "catch_spreadsheet_key: %d\n", kval);
#endif

    if (kval == GDK_Tab) {
      /* FIXME */
      ;
    }

    if (kval == GDK_Left) {
      GtkTreeViewColumn *col = NULL;
      gpointer p;

      gtk_tree_view_get_cursor(GTK_TREE_VIEW(view), NULL, &col);
      if (col != NULL) {
          p = g_object_get_data(G_OBJECT(col), "colnum");
          if (p != NULL && GPOINTER_TO_INT(p) == 1) {
            return TRUE;
          } 
      }
    } else if (kval == GDK_Up || kval == GDK_Down) {
      GtkTreePath *path = NULL;
      GtkTreeViewColumn *col = NULL;
      int i;

      gtk_tree_view_get_cursor(GTK_TREE_VIEW(view), &path, &col);
      i = (gtk_tree_path_get_indices(path))[0];

      if (kval == GDK_Down && i < sheet->datarows - 1) {
          gtk_tree_path_next(path);
      } else if (kval == GDK_Up && i > 0) {
          gtk_tree_path_prev(path);
      }
      gtk_tree_view_set_cursor(GTK_TREE_VIEW(view), path, col, 
                         FALSE);
      gtk_tree_path_free(path);
      return TRUE;
    } else if (numeric_key(&kval, sheet)) {
      GtkTreePath *path = NULL;
      GtkTreeViewColumn *col = NULL;

#if CELLDEBUG
      fprintf(stderr, "numeric key: start editing, k = %d\n", kval);
#endif

      gtk_tree_view_get_cursor(GTK_TREE_VIEW(view), &path, &col);
      if (path != NULL && col != NULL) {
          gtk_tree_view_set_cursor(GTK_TREE_VIEW(view), path, col, 
                             TRUE);
          gtk_tree_path_free(path);
          manufacture_keystroke(view, kval);
      }
    }

    return FALSE;
}

static gint catch_spreadsheet_click (GtkWidget *view, GdkEvent *event,
                             Spreadsheet *sheet)
{   
    GdkModifierType mods; 
    gint ret = FALSE;

#if CELLDEBUG
    fprintf(stderr, "** catch_spreadsheet_click()\n");
#endif

    if (event->type != GDK_BUTTON_PRESS) {
      return FALSE;
    }

    mods = widget_get_pointer_mask(view);

    if (sheet->matrix == NULL && (mods & GDK_BUTTON3_MASK)) {
      GdkEventButton *bevent = (GdkEventButton *) event;

      if (sheet->popup == NULL) {
          build_sheet_popup(sheet);
      }

      gtk_menu_popup(GTK_MENU(sheet->popup), NULL, NULL, NULL, NULL,
                   bevent->button, bevent->time);
      return TRUE;
    }     
      
    if (mods & GDK_BUTTON1_MASK) {
      GdkEventButton *bevent = (GdkEventButton *) event;
      GtkTreePath *path = NULL;
      GtkTreeViewColumn *column = NULL;   
      GtkTreePath *oldpath = NULL;
      GtkTreeViewColumn *oldcol = NULL;

#if CELLDEBUG
      fprintf(stderr, "Got button 1 click\n");
#endif

      /* where's the cursor at present? */
      gtk_tree_view_get_cursor(GTK_TREE_VIEW(sheet->view),
                         &oldpath, &oldcol);

      if (oldpath != NULL) {
          if (oldcol != NULL && sheet->entry != NULL) {
            const gchar *txt = gtk_entry_get_text(GTK_ENTRY(sheet->entry));
            gchar *pathstr = gtk_tree_path_to_string(oldpath);

#if CELLDEBUG
            fprintf(stderr, "click: calling sheet_cell_edited\n");
#endif
            sheet_cell_edited(NULL, pathstr, txt, sheet);
            g_free(pathstr);
          }
          gtk_tree_path_free(oldpath);
      }

      gtk_tree_view_get_path_at_pos(GTK_TREE_VIEW(sheet->view),
                              (gint) bevent->x, 
                              (gint) bevent->y,
                              &path, &column,
                              NULL, NULL);
      if (column != NULL) {
          gpointer p = g_object_get_data(G_OBJECT(column), "colnum");
          gint colnum = GPOINTER_TO_INT(p);

#if CELLDEBUG
          fprintf(stderr, "Clicked column: colnum = %d\n", colnum);
#endif

          if (colnum == 0) {
            /* don't respond to a click in a non-data column */
            ret = TRUE;
          } else if (sheet->cmd == SHEET_EDIT_SCALARS && colnum == 2) {
            /* call to delete scalar */
            sheet_delete_scalar(sheet, path);
          } else {
            /* activate clicked cell for editing */
            gtk_tree_view_set_cursor(GTK_TREE_VIEW(sheet->view), 
                               path, column, TRUE);
            ret = TRUE;
          }
      }
      gtk_tree_path_free(path);
    }

#if CELLDEBUG
    fprintf(stderr, "catch_spreadsheet_click returning %d\n", ret);
#endif

    return ret;
}

static void sheet_show_popup (GtkWidget *w, Spreadsheet *sheet)
{
    if (sheet->popup == NULL) {
      build_sheet_popup(sheet);
    }

    gtk_menu_popup(GTK_MENU(sheet->popup), NULL, NULL, NULL, NULL,
               1, gtk_get_current_event_time());
}

static int build_sheet_view (Spreadsheet *sheet)
{
    GtkListStore *store; 
    GtkWidget *view;
    GtkTreeViewColumn *column;
    GtkTreeSelection *select;
    gchar *col0str = NULL;
    gchar tmpstr[32];
    gint i, width, colnum;

    if (sheet->varlist != NULL) {
      sheet->datacols = 0;

      for (i=sheet->varlist[0]; i>0; i--) {
          if (var_is_hidden(datainfo, sheet->varlist[i])) {
            gretl_list_delete_at_pos(sheet->varlist, i);
          } else {
            sheet->datacols += 1;
          }
      }  

      if (sheet->varlist[0] < datainfo->v - 1) {
          sheet->flags |= SHEET_SHORT_VARLIST;
      }
    }

    if (get_local_decpoint() == ',') {
      sheet->flags |= SHEET_USE_COMMA;
    }

    store = make_sheet_liststore(sheet);
    view = gtk_tree_view_new_with_model(GTK_TREE_MODEL(store));
    g_object_unref(G_OBJECT(store));

    gtk_tree_view_set_rules_hint(GTK_TREE_VIEW(view), TRUE);

    /* build and attach the (two) cell renderers */
    create_sheet_cell_renderers(sheet);

    /* construct the first (text) column */
    if (sheet->cmd == SHEET_EDIT_SCALARS) {
      col0str = _("Name");
      width = get_name_col_width();
    } else if (sheet->matrix != NULL) {
      width = get_row_label_width();
    } else {
      width = get_obs_col_width();
    }
    column = gtk_tree_view_column_new_with_attributes(col0str,
                                          sheet->dumbcell,
                                          "text", 0, 
                                          NULL);
    gtk_tree_view_append_column(GTK_TREE_VIEW(view), column);
    set_up_sheet_column(column, width, FALSE);
    g_object_set_data(G_OBJECT(column), "colnum", GINT_TO_POINTER(0));

    width = get_data_col_width();

    if (sheet->matrix != NULL) {
      for (i=1; i<=sheet->datacols; i++) {
          if (sheet->colnames != NULL) {
            double_underscores(tmpstr, sheet->colnames[i-1]);
          } else {
            sprintf(tmpstr, "%d", i);
          }
          column = gtk_tree_view_column_new_with_attributes(tmpstr,
                                                sheet->datacell,
                                                "text", 
                                                i, 
                                                NULL);
          gtk_tree_view_append_column(GTK_TREE_VIEW(view), column);
          set_up_sheet_column(column, width, TRUE);
          g_object_set_data(G_OBJECT(column), "colnum", GINT_TO_POINTER(i));
          g_object_set_data(G_OBJECT(column), "sheet", sheet);
          gtk_tree_view_column_set_clickable(column, TRUE);
          g_signal_connect(G_OBJECT(column), "clicked",
                       G_CALLBACK(name_column_dialog), column);
      }
      sheet->colnames = NULL;
    } else if (sheet->cmd == SHEET_EDIT_SCALARS) {
      GtkCellRenderer *pixcell;
      
      /* numerical value */
      column = gtk_tree_view_column_new_with_attributes(_("Value"),
                                            sheet->datacell,
                                            "text", 
                                            1, 
                                            NULL);
      gtk_tree_view_append_column(GTK_TREE_VIEW(view), column);
      set_up_sheet_column(column, width, TRUE);
      g_object_set_data(G_OBJECT(column), "colnum", GINT_TO_POINTER(1));
      g_object_set_data(G_OBJECT(column), "sheet", sheet);

      /* deletion icon */
      pixcell = gtk_cell_renderer_pixbuf_new();
      column = gtk_tree_view_column_new_with_attributes(_("Delete"),
                                            pixcell,
                                            "pixbuf", 
                                            2, 
                                            NULL);
      gtk_tree_view_append_column(GTK_TREE_VIEW(view), column);
      gtk_tree_view_column_set_resizable(column, FALSE);
      g_object_set_data(G_OBJECT(column), "colnum", GINT_TO_POINTER(2));
      g_object_set_data(G_OBJECT(column), "sheet", sheet);
    } else {
      colnum = 0;
      for (i=1; i<=sheet->varlist[0]; i++) {
          int vi = sheet->varlist[i];

          if (var_is_hidden(datainfo, vi)) {
            continue;
          }
          colnum++;
          double_underscores(tmpstr, datainfo->varname[vi]);
          column = gtk_tree_view_column_new_with_attributes(tmpstr,
                                                sheet->datacell,
                                                "text", 
                                                colnum, 
                                                NULL);
          gtk_tree_view_append_column(GTK_TREE_VIEW(view), column);
          set_up_sheet_column(column, width, TRUE);
          g_object_set_data(G_OBJECT(column), "colnum", GINT_TO_POINTER(colnum));
      }
    }

    /* set the selection property on the tree view */
    select = gtk_tree_view_get_selection(GTK_TREE_VIEW(view));
    gtk_tree_selection_set_mode(select, GTK_SELECTION_NONE);

    if (sheet->cmd != SHEET_EDIT_SCALARS) {
      g_signal_connect(G_OBJECT(view), "cursor-changed",
                   G_CALLBACK(update_cell_position), sheet);
    }

    g_signal_connect(G_OBJECT(view), "key-press-event",
                 G_CALLBACK(catch_spreadsheet_key), sheet);

    /* attach to sheet struct */
    sheet->view = view;

    return 0;
}

static void free_spreadsheet (GtkWidget *widget, Spreadsheet **psheet) 
{
    Spreadsheet *sheet = *psheet;

    if (sheet->popup != NULL) {
      gtk_widget_destroy(sheet->popup);
    }

    free(sheet->varlist);

    if (sheet->matrix == NULL) {
      set_dataset_locked(FALSE);
    }

    if (sheet->matrix != NULL && sheet->oldmat != NULL) {
      /* delete the copied matrix */
      gretl_matrix_free(sheet->matrix);
    }

    free(sheet);

    *psheet = NULL;
}

static void free_matrix_sheet (GtkWidget *widget, Spreadsheet *sheet) 
{
    if (sheet->popup != NULL) {
      gtk_widget_destroy(sheet->popup);
    }

    free(sheet->varlist);

    if (sheet->matrix == NULL) {
      set_dataset_locked(FALSE);
    }

    if (sheet->matrix != NULL && sheet->oldmat != NULL) {
      /* delete the copied matrix */
      gretl_matrix_free(sheet->matrix);
    }

    free(sheet);
}

static int sheet_list_empty (Spreadsheet *sheet)
{
    int ret = 0;

    if (sheet->varlist == NULL) {
      free(sheet);
      ret = 1;
    } else if (sheet->varlist[0] == 0) {
      free(sheet->varlist);
      free(sheet);
      ret = 1;
    }

    if (ret) {
      errbox(_("No series to edit"));
    }

    return ret;
}

static Spreadsheet *spreadsheet_new (SheetCmd c)
{
    Spreadsheet *sheet;

    sheet = mymalloc(sizeof *sheet);
    if (sheet == NULL) return NULL;

    sheet->view = NULL;
    sheet->win = NULL;
    sheet->locator = NULL;
    sheet->entry = NULL;
    sheet->popup = NULL;
    sheet->ui = NULL;
    sheet->save = NULL;
    sheet->apply = NULL;
    sheet->dumbcell = NULL;
    sheet->datacell = NULL;
    sheet->pbuf = NULL;
    sheet->datacols = sheet->datarows = 0;
    sheet->totcols = 0;
    sheet->added_vars = 0;
    sheet->orig_main_v = 0;
    sheet->orig_nobs = 0;
    sheet->next = NEXT_DOWN;
    sheet->cid = 0;
    sheet->varlist = NULL;
    sheet->matrix = NULL;
    sheet->oldmat = NULL;
    sheet->colnames = NULL;
    sheet->cmd = c;
    sheet->flags = 0;

    sheet->mname[0] = '\0';

    if (c == SHEET_EDIT_MATRIX || c == SHEET_EDIT_SCALARS) {
      return sheet;
    }

    if (datainfo->t1 != 0 || datainfo->t2 < datainfo->n - 1) {
      sheet->flags |= SHEET_SUBSAMPLED;
    }

    sheet->orig_nobs = datainfo->t2 - datainfo->t1 + 1;

    if (sheet->cmd == SHEET_NEW_DATASET) {
      sheet->varlist = gretl_list_new(1);
    } else {
      if (sheet->cmd == SHEET_EDIT_VARLIST) {
          sheet->varlist = main_window_selection_as_list();
      } else {
          sheet->varlist = full_var_list(datainfo, NULL);
      }
      if (sheet_list_empty(sheet)) {
          return NULL;
      }
    }    

    if (sheet->varlist == NULL) {
      free(sheet);
      sheet = NULL;
    } else if (sheet->cmd == SHEET_NEW_DATASET) {
      sheet->varlist[1] = 1;
    }

    return sheet;
}

static void empty_dataset_guard (void)
{
    int t, empty = 0, miss = 0;

    if (datainfo->v == 2) {
      empty = 1;
      for (t=0; t<datainfo->n; t++) {
          if (na(Z[1][t])) {
            miss = 1;
          } else {
            empty = 0;
          }
      }
    }

    if (empty) {
      infobox(_("Warning: series %s is empty"), datainfo->varname[1]);
    } else if (miss) {
      infobox(_("Warning: there were missing observations"));
    }

    register_data(DATA_APPENDED);
}

static gint maybe_exit_sheet (GtkWidget *w, Spreadsheet *sheet)
{
    if (sheet_is_modified(sheet)) {
      int resp = yes_no_dialog ("gretl", _("Save changes?"), 1);

      if (resp == GRETL_YES) {
          get_data_from_sheet(NULL, sheet);
      } else if (resp == GRETL_CANCEL) {
          return FALSE;
      }
    } 

    if (sheet->cmd == SHEET_NEW_DATASET) {
      empty_dataset_guard();
    }
  
    gtk_widget_destroy(sheet->win);

    return FALSE;
}

static gint sheet_delete_event (GtkWidget *w, GdkEvent *event,
                        Spreadsheet *sheet)
{
    int resp;

    if (sheet_is_modified(sheet)) {
      resp = yes_no_dialog ("gretl", 
                        (sheet->matrix != NULL)? _("Save changes?") :
                        _("Do you want to save changes you have\n"
                        "made to the current data set?"), 1);
      if (resp == GRETL_YES) {
          get_data_from_sheet(NULL, sheet);
      } else if (resp == GRETL_CANCEL) {
          return TRUE;
      }
    }

    if (sheet->cmd == SHEET_NEW_DATASET) {
      empty_dataset_guard();
    }

    return FALSE;
}

static gint simple_exit_sheet (GtkWidget *w, Spreadsheet *sheet)
{
    gtk_widget_destroy(sheet->win);
    return FALSE;
}

static void size_matrix_window (Spreadsheet *sheet)
{
    int nc = sheet->datacols;
    int w, h;

    if (nc < 2) {
      nc = 2;
    }

    w = get_row_label_width() + nc * get_data_col_width() + 30;
    if (w > 640) {
      w = 640;
    }

    h = sheet->datarows * 20 + 160;
    if (h > 480) {
      h = 480;
    }

    gtk_window_resize(GTK_WINDOW(sheet->win), w, h);
}

static void size_scalars_window (Spreadsheet *sheet)
{
    int nc = 2;
    int w, h;

    w = get_row_label_width() + nc * get_data_col_width() + 
      get_delete_col_width() + 30;
    if (w > 640) {
      w = 640;
    }

    h = sheet->datarows * 20 + 160;
    if (h > 480) {
      h = 480;
    }

    gtk_window_resize(GTK_WINDOW(sheet->win), w, h);    
}

static void size_data_window (Spreadsheet *sheet)
{
    int ocw = get_obs_col_width();
    int dcw = get_data_col_width();
    int extra = 40;
    int nc, w, h = 400;

    nc = (sheet->varlist[0] > 4)? 4 : sheet->varlist[0];
    if (nc < 2) {
      extra += 40;
    }

    w = ocw + nc * dcw + extra;

    gtk_window_set_default_size(GTK_WINDOW(sheet->win), w, h);
}

/* hack to avoid losing a not-yet-committed edit to a cell
   in the sheet, on choosing Save, Apply, OK, etc */

gboolean 
button_entered (GtkWidget *w, GdkEventCrossing *e, Spreadsheet *sheet)
{
    if (sheet->entry != NULL) {
      const gchar *s = gtk_entry_get_text(GTK_ENTRY(sheet->entry));
      GtkTreePath *path;
      gchar *pathstr;

      gtk_tree_view_get_cursor(GTK_TREE_VIEW(sheet->view), &path, NULL);

      if (path != NULL) {
          pathstr = gtk_tree_path_to_string(path);
          sheet->next = NEXT_SAME;
          sheet_cell_edited(NULL, pathstr, s, sheet);
          g_free(pathstr);
          gtk_tree_path_free(path);
      }
    }

    return FALSE;
}

static void adjust_add_menu_state (Spreadsheet *sheet)
{
    sheet->flags |= SHEET_ADD_OBS_OK | SHEET_INSERT_OBS_OK;

    if (complex_subsampled() || datainfo->t2 < datainfo->n - 1) {
      sheet->flags &= ~SHEET_ADD_OBS_OK;
      sheet->flags &= ~SHEET_INSERT_OBS_OK;
    } else if ((sheet->flags & SHEET_SHORT_VARLIST) ||
             dataset_is_panel(datainfo)) {
      sheet->flags &= ~SHEET_INSERT_OBS_OK;
    }
}

static void series_sheet_add_locator (Spreadsheet *sheet, 
                              GtkWidget *hbox)
{
    GtkWidget *vbox = gtk_vbox_new(FALSE, 1);
    GtkWidget *status_box = gtk_hbox_new(FALSE, 1);
    gint w = get_obs_col_width();

    gtk_container_set_border_width(GTK_CONTAINER(status_box), 0);
    gtk_box_pack_start(GTK_BOX(vbox), status_box, FALSE, FALSE, 0);

    sheet->locator = gtk_statusbar_new(); 
    gtk_widget_set_size_request(sheet->locator, 2 * w, 20);
    gtk_statusbar_set_has_resize_grip(GTK_STATUSBAR(sheet->locator), FALSE);
    sheet->cid = gtk_statusbar_get_context_id(GTK_STATUSBAR(sheet->locator), 
                                    "current row and column");
    gtk_box_pack_start(GTK_BOX(status_box), sheet->locator, FALSE, FALSE, 0);
    gtk_box_pack_end(GTK_BOX(hbox), vbox, FALSE, FALSE, 5);
}

static void sheet_add_toolbar (Spreadsheet *sheet, GtkWidget *vbox)
{
    GtkWidget *hbox, *tbar;
    GretlToolItem *item;
    GtkToolItem *button;
    int i;

    hbox = gtk_hbox_new(FALSE, 0);
    tbar = gretl_toolbar_new();

    for (i=0; i<n_series_items; i++) {
      item = &series_items[i];
      button = gretl_toolbar_insert(tbar, item, item->func, sheet, -1);
      if (item->flag == SERIES_APPLY_BTN) {
          sheet->apply = GTK_WIDGET(button);
          gtk_widget_set_sensitive(sheet->apply, FALSE);
      } 
      if (item->flag != SERIES_ADD_BTN) {
          g_signal_connect(G_OBJECT(button), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
      }
    }

    gtk_box_pack_start(GTK_BOX(hbox), tbar, FALSE, FALSE, 0);
    series_sheet_add_locator(sheet, hbox);

    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 0);

    adjust_add_menu_state(sheet);
}

static void sheet_add_matrix_menu (Spreadsheet *sheet, GtkWidget *vbox)
{
    GtkActionGroup *actions;
    GtkWidget *mbar;

    sheet->ui = gtk_ui_manager_new();
    actions = gtk_action_group_new("SheetActions");
    gtk_action_group_set_translation_domain(actions, "gretl");

    gtk_action_group_add_actions(actions, matrix_items, 
                         sizeof matrix_items / sizeof matrix_items[0],
                         sheet);
    gtk_ui_manager_add_ui_from_string(sheet->ui, matrix_ui, -1, NULL);

    gtk_ui_manager_insert_action_group(sheet->ui, actions, 0);
    g_object_unref(actions);

    mbar = gtk_ui_manager_get_widget(sheet->ui, "/menubar");
    gtk_box_pack_start(GTK_BOX(vbox), mbar, FALSE, FALSE, 0);
}

static void sheet_add_matrix_locator (Spreadsheet *sheet, GtkWidget *vbox)
{
    GtkWidget *status_box = gtk_hbox_new(FALSE, 1);
    gint w = get_row_label_width();

    gtk_container_set_border_width(GTK_CONTAINER(status_box), 0);
    gtk_box_pack_start(GTK_BOX(vbox), status_box, FALSE, FALSE, 0);

    sheet->locator = gtk_statusbar_new(); 
    gtk_widget_set_size_request(sheet->locator, 2 * w, 20);
    gtk_statusbar_set_has_resize_grip(GTK_STATUSBAR(sheet->locator), FALSE);
    sheet->cid = gtk_statusbar_get_context_id(GTK_STATUSBAR(sheet->locator), 
                                    "current row and column");
    gtk_box_pack_start(GTK_BOX(status_box), sheet->locator, FALSE, FALSE, 0);
}

/* this has to block when the user is defining a matrix in the
   course of responding to the function call dialog; otherwise
   it should not block 
*/

static void real_show_spreadsheet (Spreadsheet **psheet, SheetCmd c,
                           int block)
{
    Spreadsheet *sheet = *psheet;
    GtkWidget *tmp, *scroller, *main_vbox;
    int err = 0;

    sheet->win = gtk_window_new(GTK_WINDOW_TOPLEVEL);

    if (sheet->matrix != NULL) {
      if (sheet->mname[0] != '\0') {
          gchar *tmp = g_strdup_printf("gretl: %s", sheet->mname);

          gtk_window_set_title(GTK_WINDOW(sheet->win), tmp);
          g_free(tmp);
      } else {
          gtk_window_set_title(GTK_WINDOW(sheet->win), _("gretl: edit matrix"));
      }
    } else if (c == SHEET_EDIT_SCALARS) {
      gtk_window_set_title(GTK_WINDOW(sheet->win), _("gretl: edit scalars"));
    } else {
      gtk_window_set_title(GTK_WINDOW(sheet->win), _("gretl: edit data"));
    }

    if (sheet->matrix != NULL) {
      size_matrix_window(sheet);
    } else if (c == SHEET_EDIT_SCALARS) {
      size_scalars_window(sheet);
    } else {
      size_data_window(sheet);
    }

    if (block) {
      g_signal_connect(G_OBJECT(sheet->win), "destroy",
                   G_CALLBACK(gtk_main_quit), NULL);
    }

    g_signal_connect(G_OBJECT(sheet->win), "delete-event",
                 G_CALLBACK(sheet_delete_event), sheet);

    main_vbox = gtk_vbox_new(FALSE, 5);
    gtk_container_set_border_width(GTK_CONTAINER(main_vbox), 5); 
    gtk_container_add(GTK_CONTAINER(sheet->win), main_vbox);

    if (sheet->matrix != NULL) {
      sheet_add_matrix_menu(sheet, main_vbox);
      sheet_add_matrix_locator(sheet, main_vbox);
    } else if (c != SHEET_EDIT_SCALARS) {
      sheet_add_toolbar(sheet, main_vbox);
    }

    gtk_widget_show_all(main_vbox);

    scroller = gtk_scrolled_window_new(NULL, NULL);
    gtk_scrolled_window_set_policy(GTK_SCROLLED_WINDOW(scroller),
                           GTK_POLICY_AUTOMATIC, GTK_POLICY_AUTOMATIC);
    gtk_scrolled_window_set_shadow_type(GTK_SCROLLED_WINDOW(scroller),
                              GTK_SHADOW_IN);    

    build_sheet_view(sheet);
    gtk_container_add(GTK_CONTAINER(scroller), sheet->view);
    gtk_box_pack_start(GTK_BOX(main_vbox), scroller, TRUE, TRUE, TRUE);

    gtk_widget_show(sheet->view);
    gtk_widget_show(scroller);

    if (sheet->matrix != NULL) {
      g_signal_connect(G_OBJECT(sheet->win), "destroy",
                   G_CALLBACK(free_matrix_sheet), sheet);
    } else if (c == SHEET_EDIT_SCALARS) {
      g_signal_connect(G_OBJECT(sheet->win), "destroy",
                   G_CALLBACK(free_spreadsheet), psheet); /* FIXME? */
    } else {
      g_signal_connect(G_OBJECT(sheet->win), "destroy",
                   G_CALLBACK(free_spreadsheet), psheet);
    } 

    if (sheet->matrix != NULL || c == SHEET_EDIT_SCALARS) {
      /* control buttons (we don't use these for editing series) */
      GtkWidget *button_box;

      button_box = gtk_hbutton_box_new();
      gtk_button_box_set_layout(GTK_BUTTON_BOX(button_box), 
                          GTK_BUTTONBOX_END);
      gtk_box_set_spacing(GTK_BOX(button_box), 10);
      gtk_widget_show(button_box);
      gtk_box_pack_start(GTK_BOX(main_vbox), button_box, FALSE, FALSE, 0);
      gtk_container_set_border_width(GTK_CONTAINER(button_box), 0);

      if (sheet->matrix != NULL && sheet->oldmat == NULL) {
          tmp = gtk_button_new_from_stock(GTK_STOCK_OK);
          gtk_container_add(GTK_CONTAINER(button_box), tmp);
          g_signal_connect(G_OBJECT(tmp), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(get_data_from_sheet), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(simple_exit_sheet), sheet);
          gtk_widget_show(tmp);
      } else if (sheet->matrix != NULL) {
          /* editing a new matrix */
          tmp = gtk_button_new_from_stock(GTK_STOCK_SAVE_AS);
          gtk_container_add(GTK_CONTAINER(button_box), tmp);
          g_signal_connect(G_OBJECT(tmp), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(matrix_save_as), sheet);
          gtk_widget_show(tmp);

          tmp = gtk_button_new_from_stock(GTK_STOCK_SAVE);
          gtk_container_add(GTK_CONTAINER(button_box), tmp);
          g_signal_connect(G_OBJECT(tmp), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(get_data_from_sheet), sheet);
          gtk_widget_show(tmp);
          sheet->save = tmp;
          gtk_widget_set_sensitive(sheet->save, FALSE);

          tmp = gtk_button_new_from_stock(GTK_STOCK_CLOSE);
          gtk_container_add(GTK_CONTAINER(button_box), tmp);
          g_signal_connect(G_OBJECT(tmp), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(maybe_exit_sheet), sheet);
          gtk_widget_show(tmp);
      } else {
          tmp = gtk_button_new_from_stock(GTK_STOCK_APPLY);
          gtk_container_add(GTK_CONTAINER(button_box), tmp);
          g_signal_connect(G_OBJECT(tmp), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(get_data_from_sheet), sheet);
          gtk_widget_show(tmp);
          sheet->apply = tmp;
          gtk_widget_set_sensitive(sheet->apply, FALSE);

          tmp = gtk_button_new_from_stock(GTK_STOCK_CLOSE);
          gtk_container_add(GTK_CONTAINER(button_box), tmp);
          g_signal_connect(G_OBJECT(tmp), "enter-notify-event",
                       G_CALLBACK(button_entered), sheet);
          g_signal_connect(G_OBJECT(tmp), "clicked",
                       G_CALLBACK(maybe_exit_sheet), sheet);
          gtk_widget_show(tmp);
      }
    }

    g_signal_connect(G_OBJECT(sheet->view), "button-press-event",
                 G_CALLBACK(catch_spreadsheet_click),
                 sheet);

    if (c == SHEET_EDIT_MATRIX) {
      err = add_matrix_data_to_sheet(sheet);
    } else if (c == SHEET_EDIT_SCALARS) {
      err = add_scalars_to_sheet(sheet);
    } else {
      err = add_data_to_sheet(sheet, c);
    }

    if (err) {
      gtk_widget_destroy(sheet->win);
      return;
    }

    select_first_editable_cell(sheet);

    gtk_widget_show(sheet->win);

    if (c != SHEET_EDIT_MATRIX && c != SHEET_EDIT_SCALARS) {
      /* we can't have the user making confounding changes elsewhere,
         while editing the dataset here */
      set_dataset_locked(TRUE);
    }

    if (block) {
      gtk_main();
    }
}

void show_spreadsheet (SheetCmd c) 
{
    static Spreadsheet *sheet;    

    if (datainfo->v == 1) {
      warnbox(_("Please add a variable to the dataset first"));
      return;
    }

    if (sheet != NULL) {
      gtk_window_present(GTK_WINDOW(sheet->win));
      return;
    }

    sheet = spreadsheet_new(c);
    if (sheet == NULL) {
      return;
    }

    real_show_spreadsheet(&sheet, c, 0);
}

static void real_edit_scalars (const char *vname);

static void start_scalars (GtkWidget *w, dialog_t *dlg) 
{
    const gchar *buf;
    char varname[VNAMELEN];
 
    buf = edit_dialog_get_text(dlg);

    if (buf == NULL || gui_validate_varname(buf, GRETL_TYPE_DOUBLE)) {
      return;
    }

    *varname = 0;
    strncat(varname, buf, VNAMELEN - 1);

    close_dialog(dlg);
    real_edit_scalars(varname);
}

static void maybe_add_scalar (void)
{
    int resp;

    resp = yes_no_dialog(_("gretl: scalars"),
                   _("No scalar variable are currently defined.\n"
                     "Add one now?"),
                   0);

    if (resp == GRETL_YES) {
      int cancel = 0;

      edit_dialog(_("gretl: name scalar"), 
                _("Enter name for new scalar\n"
                  "(max. 15 characters)"),
                NULL, start_scalars, NULL, 
                0, VARCLICK_NONE, &cancel);
    }
}

static void real_edit_scalars (const char *vname) 
{
    static Spreadsheet *sheet; 
    int ns;

    if (sheet != NULL) {
      gtk_window_present(GTK_WINDOW(sheet->win));
      return;
    }

    if (vname != NULL) {
      ns = 1;
    } else {
      ns = n_saved_scalars();
    }

    if (ns == 0) {
      maybe_add_scalar();
      return;
    }

    sheet = spreadsheet_new(SHEET_EDIT_SCALARS);
    if (sheet == NULL) {
      return;
    }

    sheet->datarows = ns;
    sheet->datacols = 1;

    if (vname != NULL) {
      strcpy(sheet->mname, vname);
    }

    real_show_spreadsheet(&sheet, SHEET_EDIT_SCALARS, 0);
}

void edit_scalars (void) 
{
    real_edit_scalars(NULL);
}

struct gui_matrix_spec {
    char name[VNAMELEN];
    int rows;
    int cols;
    double fill;
    char *formula;
    int uselist;
    gretl_matrix *m;
};

struct mdialog {
    GtkWidget *dlg;
    GtkWidget *nentry;
    GtkWidget *ventry;
    GtkWidget *numerics;
    GtkWidget *formula;
    struct gui_matrix_spec *spec;
};

static void spin_call (GtkWidget *s, int *n)
{
    *n = (int) gtk_spin_button_get_value(GTK_SPIN_BUTTON(s));
}

static void matrix_dialog_ok (GtkWidget *w, struct mdialog *mdlg)
{
    const char *etxt;
    int err;

    etxt = gtk_entry_get_text(GTK_ENTRY(mdlg->nentry));

    if (etxt == NULL || *etxt == '\0') {
      infobox(_("You must give the matrix a name"));
      return;
    }

    err = gui_validate_varname(etxt, GRETL_TYPE_MATRIX);
    if (err) {
      return;
    }

    strcpy(mdlg->spec->name, etxt);
    mdlg->spec->uselist = 0;

    if (GTK_WIDGET_SENSITIVE(mdlg->numerics)) {
      double x;

      etxt = gtk_entry_get_text(GTK_ENTRY(mdlg->ventry));
      x = gui_double_from_string(etxt, &err);
      if (err) {
          return;
      }
      mdlg->spec->fill = x;
    } else if (GTK_WIDGET_SENSITIVE(mdlg->formula)) {
      etxt = gtk_entry_get_text(GTK_ENTRY(mdlg->formula));
      if (etxt == NULL || *etxt == '\0') {
          errbox(_("The matrix formula is empty"));
          return;
      }
      mdlg->spec->formula = g_strdup(etxt);
    } else {
      /* will select data series */
      mdlg->spec->uselist = 1;
    }

    gtk_widget_destroy(mdlg->dlg);
}

static gint choose_series (GtkWidget *w, struct mdialog *mdlg)
{
    if (gtk_toggle_button_get_active(GTK_TOGGLE_BUTTON(w))) {
      if (mdlg->numerics != NULL) {
          gtk_widget_set_sensitive(mdlg->numerics, FALSE);
      } 
      if (mdlg->formula != NULL) {
          gtk_widget_set_sensitive(mdlg->formula, FALSE);
      } 
    }

    return FALSE;
}

static gint choose_numeric (GtkWidget *w, struct mdialog *mdlg)
{
    if (gtk_toggle_button_get_active(GTK_TOGGLE_BUTTON(w))) {
      if (mdlg->numerics != NULL) {
          gtk_widget_set_sensitive(mdlg->numerics, TRUE);
      } 
      if (mdlg->formula != NULL) {
          gtk_widget_set_sensitive(mdlg->formula, FALSE);
      } 
    } 

    return FALSE;
}

static gint choose_formula (GtkWidget *w, struct mdialog *mdlg)
{
    if (gtk_toggle_button_get_active(GTK_TOGGLE_BUTTON(w))) {
      if (mdlg->numerics != NULL) {
          gtk_widget_set_sensitive(mdlg->numerics, FALSE);
      } 
      if (mdlg->formula != NULL) {
          gtk_widget_set_sensitive(mdlg->formula, TRUE);
      }
    }

    return FALSE;
}

static int new_matrix_dialog (struct gui_matrix_spec *spec)
{
    struct mdialog mdlg;
    GSList *group;
    GtkWidget *dlg;
    GtkWidget *hbox;
    GtkWidget *vbox;
    GtkWidget *tab;
    GtkWidget *rb;
    GtkWidget *w;
    int maxdim = 1000;
    int canceled = 0;

    dlg = gretl_dialog_new("Matrix", mdata->main, 
                     GRETL_DLG_BLOCK | GRETL_DLG_MODAL);
    vbox = gtk_dialog_get_content_area(GTK_DIALOG(dlg));
    mdlg.dlg = dlg;
    mdlg.spec = spec;
    mdlg.numerics = NULL;
    mdlg.formula = NULL;

    /* top label */
    hbox = gtk_hbox_new(FALSE, 5);
    w = gtk_label_new(_("New matrix"));
    gtk_box_pack_start(GTK_BOX(hbox), w, TRUE, TRUE, 0);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);

    /* matrix name entry */
    hbox = gtk_hbox_new(FALSE, 5);
    w = gtk_label_new(_("Name:"));
    gtk_box_pack_start(GTK_BOX(hbox), w, FALSE, FALSE, 5);
    w = gtk_entry_new();
    gtk_entry_set_max_length(GTK_ENTRY(w), VNAMELEN-1);
    gtk_entry_set_width_chars(GTK_ENTRY(w), VNAMELEN+3);
    if (*spec->name != '\0') {
      gtk_entry_set_text(GTK_ENTRY(w), spec->name);
      gtk_widget_set_sensitive(w, FALSE);
    } else {
      gtk_entry_set_activates_default(GTK_ENTRY(w), TRUE);
    }
    gtk_box_pack_start(GTK_BOX(hbox), w, FALSE, FALSE, 5);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);
    mdlg.nentry = w;

    /* matrix construction options */

    /* option: build from series */

    hbox = gtk_hbox_new(FALSE, 5);
    rb = gtk_radio_button_new_with_label(NULL, _("Build from series"));
    g_signal_connect(G_OBJECT(rb), "clicked", G_CALLBACK(choose_series), &mdlg);
    gtk_box_pack_start(GTK_BOX(hbox), rb, TRUE, TRUE, 0);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);

    /* option: build numerically */

    hbox = gtk_hbox_new(FALSE, 5);
    group = gtk_radio_button_get_group(GTK_RADIO_BUTTON(rb));
    rb = gtk_radio_button_new_with_label(group, _("Build numerically"));
    g_signal_connect(G_OBJECT(rb), "clicked", G_CALLBACK(choose_numeric), &mdlg);
    gtk_box_pack_start(GTK_BOX(hbox), rb, TRUE, TRUE, 0);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);

    hbox = gtk_hbox_new(FALSE, 5);
    tab = gtk_table_new(3, 2, FALSE);
    gtk_table_set_col_spacing(GTK_TABLE(tab), 0, 5);
    
    w = gtk_label_new(_("Number of rows:"));
    gtk_misc_set_alignment(GTK_MISC(w), 0, 1);
    gtk_table_attach_defaults(GTK_TABLE(tab), w, 0, 1, 0, 1);
    w = gtk_spin_button_new_with_range(1, maxdim, 1);
    gtk_spin_button_set_value(GTK_SPIN_BUTTON(w), (gdouble) spec->rows);
    g_signal_connect(G_OBJECT(w), "value-changed",
                 G_CALLBACK(spin_call), &spec->rows);
    gtk_table_attach_defaults(GTK_TABLE(tab), w, 1, 2, 0, 1);

    w = gtk_label_new(_("Number of columns:"));
    gtk_misc_set_alignment(GTK_MISC(w), 0, 1);
    gtk_table_attach_defaults(GTK_TABLE(tab), w, 0, 1, 1, 2);
    w = gtk_spin_button_new_with_range(1, maxdim, 1);
    gtk_spin_button_set_value(GTK_SPIN_BUTTON(w), (gdouble) spec->cols);
    g_signal_connect(G_OBJECT(w), "value-changed",
                 G_CALLBACK(spin_call), &spec->cols);
    gtk_table_attach_defaults(GTK_TABLE(tab), w, 1, 2, 1, 2);

    w = gtk_label_new(_("Initial fill value:"));
    gtk_misc_set_alignment(GTK_MISC(w), 0, 1);
    gtk_table_attach_defaults(GTK_TABLE(tab), w, 0, 1, 2, 3);
    w = gtk_entry_new();
    gtk_entry_set_max_length(GTK_ENTRY(w), VNAMELEN-1);
    gtk_entry_set_width_chars(GTK_ENTRY(w), VNAMELEN+3);
    gtk_entry_set_activates_default(GTK_ENTRY(w), TRUE);
    gtk_entry_set_text(GTK_ENTRY(w), "0");
    gtk_table_attach_defaults(GTK_TABLE(tab), w, 1, 2, 2, 3);
    mdlg.ventry = w;

    w = gtk_label_new(" ");
    gtk_box_pack_start(GTK_BOX(hbox), w, FALSE, FALSE, 5);
    gtk_box_pack_start(GTK_BOX(hbox), tab, FALSE, FALSE, 5);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);

    gtk_widget_set_sensitive(hbox, FALSE);
    mdlg.numerics = hbox;

    /* option: build from formula */

    hbox = gtk_hbox_new(FALSE, 5);
    group = gtk_radio_button_get_group(GTK_RADIO_BUTTON(rb));
    rb = gtk_radio_button_new_with_label(group, _("Build from formula"));
    g_signal_connect(G_OBJECT(rb), "clicked", G_CALLBACK(choose_formula), &mdlg);
    gtk_box_pack_start(GTK_BOX(hbox), rb, TRUE, TRUE, 0);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);
    
    hbox = gtk_hbox_new(FALSE, 5);
    w = gtk_label_new(" ");
    gtk_box_pack_start(GTK_BOX(hbox), w, FALSE, FALSE, 5);
    w = gtk_entry_new();
    gtk_entry_set_max_length(GTK_ENTRY(w), MAXLEN);
    gtk_entry_set_width_chars(GTK_ENTRY(w), 32);
    gtk_entry_set_activates_default(GTK_ENTRY(w), TRUE);
    gtk_box_pack_start(GTK_BOX(hbox), w, FALSE, FALSE, 5);
    gtk_box_pack_start(GTK_BOX(vbox), hbox, FALSE, FALSE, 5);

    gtk_widget_set_sensitive(w, FALSE);
    mdlg.formula = w;

    /* control buttons */
    hbox = gtk_dialog_get_action_area(GTK_DIALOG(dlg));
    w = cancel_delete_button(hbox, dlg, &canceled);
    w = ok_button(hbox);
    g_signal_connect(G_OBJECT(w), "clicked", G_CALLBACK(matrix_dialog_ok), &mdlg);
    gtk_widget_grab_default(w);

    gtk_widget_show_all(dlg);

    return canceled;
}

static int matrix_from_list (selector *sr)
{
    struct gui_matrix_spec *s = selector_get_data(sr);
    const char *buf = selector_list(sr);
    int *list;
    int err = 0;

    if (buf == NULL || *buf == '\0') {
      errbox("No variables are selected");
      return 1;
    } else {
      list = gretl_list_from_string(buf, &err);
    }

    if (err) {
      gui_errmsg(err);
      return err;
    }

    s->m = gretl_matrix_data_subset(list, (const double **) Z, 
                            datainfo->t1, datainfo->t2, 
                            M_MISSING_SKIP, &err);

    if (!err) {
      err = add_or_replace_user_matrix(s->m, s->name);
    }

    if (err) {
      gui_errmsg(err);
    }

    free(list);

    return err;
} 

static int
matrix_from_formula (struct gui_matrix_spec *s)
{
    gchar *genline = g_strdup_printf("matrix %s = %s", s->name, s->formula);
    int err;

    err = generate(genline, &Z, datainfo, OPT_NONE, NULL); 

    if (err) {
      gui_errmsg(err);
    } else {
      s->m = get_matrix_by_name(s->name);
      if (s->m == NULL) {
          err = 1;
      }
    }

    g_free(genline);

    return err;
}

static int
matrix_from_spec (struct gui_matrix_spec *s)
{
    int err = 0;

    s->m = gretl_matrix_alloc(s->rows, s->cols);
    if (s->m == NULL) {
      err = E_ALLOC;
    } else {
      gretl_matrix_fill(s->m, s->fill);
      err = add_or_replace_user_matrix(s->m, s->name);
    }

    if (err) {
      gui_errmsg(err);
    } 

    return err;
}

static void gui_matrix_spec_init (struct gui_matrix_spec *s,
                          gretl_matrix *m,
                          const char *name)
{
    *s->name = '\0';
    s->rows = 2;
    s->cols = 2;
    s->fill = 0;
    s->formula = NULL;
    s->uselist = 1;
    s->m = m;

    if (name != NULL) {
      strncat(s->name, name, VNAMELEN - 1);
    }
}

static void edit_matrix (gretl_matrix *m, const char *name,
                   int block)
{
    Spreadsheet *sheet = NULL;
    int err = 0;

    if (m == NULL || name == NULL) {
      return;
    }

    sheet = spreadsheet_new(SHEET_EDIT_MATRIX);
    if (sheet == NULL) {
      return;
    }

    strcpy(sheet->mname, name);
    sheet->oldmat = m;
    sheet->matrix = gretl_matrix_copy(m);
    if (sheet->matrix == NULL) {
      err = E_ALLOC;
    }

    if (err) {
      nomem();
      free(sheet);
      sheet = NULL;
    } else {
      sheet->colnames = user_matrix_get_column_names(m);
      sheet->datarows = gretl_matrix_rows(sheet->matrix);
      sheet->datacols = gretl_matrix_cols(sheet->matrix);
      real_show_spreadsheet(&sheet, SHEET_EDIT_MATRIX, block);
    }

    if (!block && sheet != NULL) {
      /* protect matrix from deletion while editing */
      user_matrix *u = get_user_matrix_by_name(name);

      if (u != NULL) {
          g_object_set_data(G_OBJECT(sheet->win), "object", u);
          g_signal_connect(G_OBJECT(sheet->win), "destroy", 
                       G_CALLBACK(winstack_remove), sheet->win);
          winstack_add(sheet->win);
      }
    }
}

/* note that both @m and @name may be NULL depending on how
   we are called */

static void real_gui_new_matrix (gretl_matrix *m, const char *name)
{
    struct gui_matrix_spec spec;
    int block = (m == NULL);
    int err = 0;

    gui_matrix_spec_init(&spec, m, name);

    err = new_matrix_dialog(&spec);
    if (err) {
      return;
    }

    if (spec.uselist) {
      /* matrix from listed vars */
      simple_selection(_("Define matrix"), matrix_from_list, 
                   DEFINE_MATRIX, &spec);
      return;
    } else if (spec.formula != NULL) {
      /* matrix from genr-style formula */
      matrix_from_formula(&spec);
      free(spec.formula);
      return;
    } else {
      /* numerical specification: open editor if OK */
      err = matrix_from_spec(&spec);
      if (err) {
          return;
      }
    }

    edit_matrix(spec.m, spec.name, block);
}

/* callback for "Define matrix..." in main window */

void gui_new_matrix (void)
{
    real_gui_new_matrix(NULL, NULL);
}

void edit_user_matrix_by_name (const char *name)
{
    user_matrix *u = get_user_matrix_by_name(name);
    gretl_matrix *m = get_matrix_by_name(name);
    GtkWidget *w;

    /* do we already have a window open, editing this
       matrix? */
    w = match_window_by_data(u);
    if (w != NULL) {
      gtk_window_present(GTK_WINDOW(w));
      return;
    }

    if (m == NULL) {
      errbox(_("Couldn't open '%s'"), name);
    } else if (gretl_is_null_matrix(m)) {
      real_gui_new_matrix(m, name);
    } else {
      edit_matrix(m, name, 0);
    }
}

/* mechanism for locking dataset against changes while editing */

static int locked;

void set_dataset_locked (gboolean s)
{
    locked = s;
    flip(mdata->ui, "/menubar/Sample", !s);
}

int dataset_locked (void)
{
    if (locked) {
      errbox(_("You can't do this while editing the dataset"));
    }

    return locked;
}

Generated by  Doxygen 1.6.0   Back to index