fp20_write_lp.c

Go to the documentation of this file.

/****************************************************************************/
/*                                                                          */
/*  This file is part of QSopt_ex.                                          */
/*                                                                          */
/*  (c) Copyright 2006 by David Applegate, William Cook, Sanjeeb Dash,      */
/*  and Daniel Espinoza                                                     */
/*                                                                          */
/*  Sanjeeb Dash ownership of copyright in QSopt_ex is derived from his     */
/*  copyright in QSopt.                                                     */
/*                                                                          */
/*  This code may be used under the terms of the GNU General Public License */
/*  (Version 2.1 or later) as published by the Free Software Foundation.    */
/*                                                                          */
/*  Alternatively, use is granted for research purposes only.               */
/*                                                                          */
/*  It is your choice of which of these two licenses you are operating      */
/*  under.                                                                  */
/*                                                                          */
/*  We make no guarantees about the correctness or usefulness of this code. */
/*                                                                          */
/****************************************************************************/

/* RCS_INFO = "$RCSfile: wr_lp.c,v $ $Revision: 1.2 $ $Date: 2003/11/05 16:49:52 $"; */

/****************************************************************************/
/*                                                                          */
/*               Routines to support writing of LP files                    */
/*                                                                          */
/****************************************************************************/

/* 
 * -) anything after '\' is comment 
 * -) variables consist of a-z A-Z 0-9!"#$%(),;.?@_`'{}|~ 
 *    don't start with a digit or '.'
 */

#include "config.h"
#include "fp20_iqsutil.h"
#include "fp20_lpdefs.h"
#include "fp20_write_lp.h"
#ifdef USEDMALLOC
#include "dmalloc.h"
#endif

void fp20_ILLwrite_lp_state_init (
  fp20_ILLwrite_lp_state * line,
  const char *str)
{
  line->total = 0;
  line->p = line->buf;
  *line->p = '\0';
  if (str != NULL)
  {
    fp20_ILLwrite_lp_state_append (line, str);
  }
}

void fp20_ILLwrite_lp_state_append (
  fp20_ILLwrite_lp_state * line,
  const char *str)
{
  int len, rval = 0;

  ILL_FAILfalse (str, "Must have non NULL string");
  sprintf (line->p, str);
  len = strlen (line->p);
  line->total += len;
  line->p += len;
CLEANUP:
  return;
}

void fp20_ILLwrite_lp_state_append_coef (
  fp20_ILLwrite_lp_state * line,
  EGfp20_t v,
  int cnt)
{
  EGfp20_t ntmp;
  int len = 0;

  fp20_EGlpNumInitVar (ntmp);
  fp20_EGlpNumCopy (ntmp, v);
  if (fp20_EGlpNumIsLessZero (ntmp))
  {
    sprintf (line->p, " - ");
    len = 3;
    fp20_EGlpNumSign (ntmp);
  }
  else
  {
    if (cnt > 0)
    {
      sprintf (line->p, " + ");
      len = 3;
    }
    else
    {
      sprintf (line->p, " ");
      len = 1;
    }
  }
  line->p += len;
  line->total += len;
  if (fp20_EGlpNumIsNeqq (ntmp, fp20_oneLpNum))
  {
    fp20_ILLwrite_lp_state_append_number (line, ntmp);
  }
  fp20_EGlpNumClearVar (ntmp);
}

/* so that diff will not stumble over too many number format differences 
 * between c and java generated lp files we make here sure that doubles 
 * which are printed without a ".xxx" part get ".0" from us
 */
static void fp20_append_number (
  fp20_ILLwrite_lp_state * line,
  EGfp20_t v);
void fp20_ILLwrite_lp_state_append_number (
  fp20_ILLwrite_lp_state * line,
  EGfp20_t v)
{
  /* write a blank after 'inf' in case it is used as a coefficient and 
   * a variable follows */
  if (fp20_EGlpNumIsEqqual (v, fp20_ILL_MAXDOUBLE))
  {
    fp20_ILLwrite_lp_state_append (line, "inf ");
  }
  else
  {
    if (fp20_EGlpNumIsEqqual (v, fp20_ILL_MINDOUBLE))
    {
      fp20_ILLwrite_lp_state_append (line, "-inf ");
    }
    else
      fp20_append_number (line, v);
  }
}

static void fp20_append_number (
  fp20_ILLwrite_lp_state * line,
  EGfp20_t v)
{
  int len = 0;
  char *numstr = fp20_EGlpNumGetStr (v);

  sprintf (line->p, "%s%n", numstr, &len);
  EGfree (numstr);
  line->p += len;
  line->total += len;
}

#if 0
#define D_SCALE (1e9)
static void fp20_append_number (
  fp20_ILLwrite_lp_state * line,
  double x)
{
  /* Better code for writing rational problems */
  int i, k;
  int got = 0;
  double w, t;
  int nnum = 0;
  int nden = 0;

  if (x != 0.0)
  {
    if (x < 0.0)
      w = -x;
    else
      w = x;

    for (i = -9, t = 0.000000001; i <= 7; i++, t *= 10.0)
    {
      if (w >= t && w <= t * 10)
      {
        got = 1;
        break;
      }
    }
    if (got == 0)
    {
      fprintf (stderr, "Out-of-range number: %f\n", x);
      exit (1);
    }
  }

  if (x < 0.0)
  {
    sprintf (line->p, "-");
    line->p++;
    line->total++;
    x = -x;
  }

  while (x >= 10.0 * D_SCALE)
  {
    x /= 10.0;
    nnum++;
  }

  /* (x != (double) (int) x) is a hack to let small integers print nicely */

  while (x < D_SCALE && x != (double) (int) x)
  {
    x *= 10.0;
    nden++;
  }

  sprintf (line->p, "%.0f%n", x, &k);
  line->p += k;
  line->total += k;

  for (i = 0; i < nnum; i++)
  {
    sprintf (line->p, "0");
    line->p++;
    line->total++;
  }

  if (nden)
  {
    sprintf (line->p, "/1%n", &k);
    line->p += k;
    line->total += k;
    for (i = 0; i < nden; i++)
    {
      sprintf (line->p, "0");
      line->p++;
      line->total++;
    }
  }
}
#endif

void fp20_ILLwrite_lp_state_save_start (
  fp20_ILLwrite_lp_state * line)
{
  line->startlen = line->total;
}

void fp20_ILLwrite_lp_state_start (
  fp20_ILLwrite_lp_state * line)
{
  int j;

  for (j = 0; j < line->startlen; j++)
  {
    line->buf[j] = ' ';
  }
  line->buf[j] = '\0';
  line->p = line->buf + j;
  line->total = j;
}

---