alpar@9: /* glpapi04.c (problem scaling routines) */
alpar@9: 
alpar@9: /***********************************************************************
alpar@9: *  This code is part of GLPK (GNU Linear Programming Kit).
alpar@9: *
alpar@9: *  Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008,
alpar@9: *  2009, 2010, 2011 Andrew Makhorin, Department for Applied Informatics,
alpar@9: *  Moscow Aviation Institute, Moscow, Russia. All rights reserved.
alpar@9: *  E-mail: <mao@gnu.org>.
alpar@9: *
alpar@9: *  GLPK is free software: you can redistribute it and/or modify it
alpar@9: *  under the terms of the GNU General Public License as published by
alpar@9: *  the Free Software Foundation, either version 3 of the License, or
alpar@9: *  (at your option) any later version.
alpar@9: *
alpar@9: *  GLPK is distributed in the hope that it will be useful, but WITHOUT
alpar@9: *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
alpar@9: *  or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
alpar@9: *  License for more details.
alpar@9: *
alpar@9: *  You should have received a copy of the GNU General Public License
alpar@9: *  along with GLPK. If not, see <http://www.gnu.org/licenses/>.
alpar@9: ***********************************************************************/
alpar@9: 
alpar@9: #include "glpapi.h"
alpar@9: 
alpar@9: /***********************************************************************
alpar@9: *  NAME
alpar@9: *
alpar@9: *  glp_set_rii - set (change) row scale factor
alpar@9: *
alpar@9: *  SYNOPSIS
alpar@9: *
alpar@9: *  void glp_set_rii(glp_prob *lp, int i, double rii);
alpar@9: *
alpar@9: *  DESCRIPTION
alpar@9: *
alpar@9: *  The routine glp_set_rii sets (changes) the scale factor r[i,i] for
alpar@9: *  i-th row of the specified problem object. */
alpar@9: 
alpar@9: void glp_set_rii(glp_prob *lp, int i, double rii)
alpar@9: {     if (!(1 <= i && i <= lp->m))
alpar@9:          xerror("glp_set_rii: i = %d; row number out of range\n", i);
alpar@9:       if (rii <= 0.0)
alpar@9:          xerror("glp_set_rii: i = %d; rii = %g; invalid scale factor\n",
alpar@9:             i, rii);
alpar@9:       if (lp->valid && lp->row[i]->rii != rii)
alpar@9:       {  GLPAIJ *aij;
alpar@9:          for (aij = lp->row[i]->ptr; aij != NULL; aij = aij->r_next)
alpar@9:          {  if (aij->col->stat == GLP_BS)
alpar@9:             {  /* invalidate the basis factorization */
alpar@9:                lp->valid = 0;
alpar@9:                break;
alpar@9:             }
alpar@9:          }
alpar@9:       }
alpar@9:       lp->row[i]->rii = rii;
alpar@9:       return;
alpar@9: }
alpar@9: 
alpar@9: /***********************************************************************
alpar@9: *  NAME
alpar@9: *
alpar@9: *  glp_set sjj - set (change) column scale factor
alpar@9: *
alpar@9: *  SYNOPSIS
alpar@9: *
alpar@9: *  void glp_set_sjj(glp_prob *lp, int j, double sjj);
alpar@9: *
alpar@9: *  DESCRIPTION
alpar@9: *
alpar@9: *  The routine glp_set_sjj sets (changes) the scale factor s[j,j] for
alpar@9: *  j-th column of the specified problem object. */
alpar@9: 
alpar@9: void glp_set_sjj(glp_prob *lp, int j, double sjj)
alpar@9: {     if (!(1 <= j && j <= lp->n))
alpar@9:          xerror("glp_set_sjj: j = %d; column number out of range\n", j);
alpar@9:       if (sjj <= 0.0)
alpar@9:          xerror("glp_set_sjj: j = %d; sjj = %g; invalid scale factor\n",
alpar@9:             j, sjj);
alpar@9:       if (lp->valid && lp->col[j]->sjj != sjj && lp->col[j]->stat ==
alpar@9:          GLP_BS)
alpar@9:       {  /* invalidate the basis factorization */
alpar@9:          lp->valid = 0;
alpar@9:       }
alpar@9:       lp->col[j]->sjj = sjj;
alpar@9:       return;
alpar@9: }
alpar@9: 
alpar@9: /***********************************************************************
alpar@9: *  NAME
alpar@9: *
alpar@9: *  glp_get_rii - retrieve row scale factor
alpar@9: *
alpar@9: *  SYNOPSIS
alpar@9: *
alpar@9: *  double glp_get_rii(glp_prob *lp, int i);
alpar@9: *
alpar@9: *  RETURNS
alpar@9: *
alpar@9: *  The routine glp_get_rii returns current scale factor r[i,i] for i-th
alpar@9: *  row of the specified problem object. */
alpar@9: 
alpar@9: double glp_get_rii(glp_prob *lp, int i)
alpar@9: {     if (!(1 <= i && i <= lp->m))
alpar@9:          xerror("glp_get_rii: i = %d; row number out of range\n", i);
alpar@9:       return lp->row[i]->rii;
alpar@9: }
alpar@9: 
alpar@9: /***********************************************************************
alpar@9: *  NAME
alpar@9: *
alpar@9: *  glp_get_sjj - retrieve column scale factor
alpar@9: *
alpar@9: *  SYNOPSIS
alpar@9: *
alpar@9: *  double glp_get_sjj(glp_prob *lp, int j);
alpar@9: *
alpar@9: *  RETURNS
alpar@9: *
alpar@9: *  The routine glp_get_sjj returns current scale factor s[j,j] for j-th
alpar@9: *  column of the specified problem object. */
alpar@9: 
alpar@9: double glp_get_sjj(glp_prob *lp, int j)
alpar@9: {     if (!(1 <= j && j <= lp->n))
alpar@9:          xerror("glp_get_sjj: j = %d; column number out of range\n", j);
alpar@9:       return lp->col[j]->sjj;
alpar@9: }
alpar@9: 
alpar@9: /***********************************************************************
alpar@9: *  NAME
alpar@9: *
alpar@9: *  glp_unscale_prob - unscale problem data
alpar@9: *
alpar@9: *  SYNOPSIS
alpar@9: *
alpar@9: *  void glp_unscale_prob(glp_prob *lp);
alpar@9: *
alpar@9: *  DESCRIPTION
alpar@9: *
alpar@9: *  The routine glp_unscale_prob performs unscaling of problem data for
alpar@9: *  the specified problem object.
alpar@9: *
alpar@9: *  "Unscaling" means replacing the current scaling matrices R and S by
alpar@9: *  unity matrices that cancels the scaling effect. */
alpar@9: 
alpar@9: void glp_unscale_prob(glp_prob *lp)
alpar@9: {     int m = glp_get_num_rows(lp);
alpar@9:       int n = glp_get_num_cols(lp);
alpar@9:       int i, j;
alpar@9:       for (i = 1; i <= m; i++) glp_set_rii(lp, i, 1.0);
alpar@9:       for (j = 1; j <= n; j++) glp_set_sjj(lp, j, 1.0);
alpar@9:       return;
alpar@9: }
alpar@9: 
alpar@9: /* eof */