static GlpkFileVersion
gnm_glpk_detect_version (GnmGlpk *lp,
GsfInputTextline *tl)
{
GnmSubSolver *subsol = lp->parent;
gsf_off_t cur = gsf_input_tell (GSF_INPUT (tl));
GlpkFileVersion ver = GLPK_UNKNOWN;
const char *line;
unsigned cols, rows;
if ((line = gsf_input_textline_utf8_gets (tl)) == NULL)
goto out;
if (sscanf (line, "%u %u", &rows, &cols) == 2 &&
cols == g_hash_table_size (subsol->cell_from_name)) {
ver = GLPK_457;
if (gnm_solver_debug ())
g_printerr ("Detected version 4.57 file format\n");
goto out;
}
if ((line[0] == 'c' || line[0] == 's') && line[1] == ' ') {
ver = GLPK_458;
if (gnm_solver_debug ())
g_printerr ("Detected version 4.58 file format\n");
goto out;
}
out:
// Extra seek due to gsf bug
gsf_input_seek (GSF_INPUT (tl), cur + 1, G_SEEK_SET);
gsf_input_seek (GSF_INPUT (tl), cur, G_SEEK_SET);
return ver;
}
GnmSolver *
nlsolve_solver_factory (GnmSolverFactory *factory, GnmSolverParameters *params)
{
GnmSolver *res = g_object_new (GNM_SOLVER_TYPE,
"params", params,
NULL);
GnmNlsolve *nl = g_new0 (GnmNlsolve, 1);
GSList *input_cells, *l;
int n;
GnmValue const *vinput = gnm_solver_param_get_input (params);
GnmEvalPos ep;
GnmCellRef origin;
nl->parent = GNM_SOLVER (res);
nl->maximize = (params->problem_type == GNM_SOLVER_MAXIMIZE);
eval_pos_init_sheet (&ep, params->sheet);
if (vinput) {
gnm_cellref_make_abs (&origin, &vinput->v_range.cell.a, &ep);
nl->origin.col = origin.col;
nl->origin.row = origin.row;
nl->input_width = value_area_get_width (vinput, &ep);
nl->input_height = value_area_get_height (vinput, &ep);
}
nl->debug = gnm_solver_debug ();
nl->max_iter = params->options.max_iter;
nl->min_factor = 1e-10;
nl->target = gnm_solver_param_get_target_cell (params);
nl->vars = g_ptr_array_new ();
input_cells = gnm_solver_param_get_input_cells (params);
for (l = input_cells; l; l = l->next)
g_ptr_array_add (nl->vars, l->data);
g_slist_free (input_cells);
n = nl->vars->len;
nl->x0 = g_new (gnm_float, n);
nl->xk = g_new (gnm_float, n);
g_signal_connect (res, "prepare", G_CALLBACK (gnm_nlsolve_prepare), nl);
g_signal_connect (res, "start", G_CALLBACK (gnm_nlsolve_start), nl);
g_signal_connect (res, "stop", G_CALLBACK (gnm_nlsolve_stop), nl);
g_object_set_data_full (G_OBJECT (res), PRIVATE_KEY, nl,
(GDestroyNotify)gnm_nlsolve_final);
return res;
}
GnmSolver *
nlsolve_solver_factory (GnmSolverFactory *factory, GnmSolverParameters *params)
{
GnmIterSolver *isol = g_object_new
(GNM_ITER_SOLVER_TYPE,
"params", params,
"flip-sign", (params->problem_type == GNM_SOLVER_MAXIMIZE),
NULL);
GnmSolver *sol = GNM_SOLVER (isol);
GnmNlsolve *nl = g_new0 (GnmNlsolve, 1);
GnmSolverIteratorCompound *citer;
GnmSolverIterator *iter;
citer = g_object_new (GNM_SOLVER_ITERATOR_COMPOUND_TYPE, NULL);
iter = gnm_solver_iterator_new_func (G_CALLBACK (gnm_nlsolve_iterate), nl);
gnm_solver_iterator_compound_add (citer, iter, 1);
gnm_solver_iterator_compound_add (citer, gnm_solver_iterator_new_polish (isol), 0);
gnm_iter_solver_set_iterator (isol, GNM_SOLVER_ITERATOR (citer));
g_object_unref (citer);
nl->sol = sol;
nl->isol = isol;
nl->debug = gnm_solver_debug ();
nl->min_factor = 1e-10;
nl->n = nl->sol->input_cells->len;
g_signal_connect (isol, "prepare", G_CALLBACK (gnm_nlsolve_prepare), nl);
g_object_set_data_full (G_OBJECT (isol), PRIVATE_KEY, nl,
(GDestroyNotify)gnm_nlsolve_final);
return sol;
}
