source/basicsimplex.h

   1 /* This file is part of Shapes.
   2  *
   3  * Shapes is free software: you can redistribute it and/or modify
   4  * it under the terms of the GNU General Public License as published by
   5  * the Free Software Foundation, either version 3 of the License, or
   6  * any later version.
   7  *
   8  * Shapes is distributed in the hope that it will be useful,
   9  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11  * GNU General Public License for more details.
  12  *
  13  * You should have received a copy of the GNU General Public License
  14  * along with Shapes.  If not, see <http://www.gnu.org/licenses/>.
  15  *
  16  * Copyright 2008 Henrik Tidefelt
  17  */
  18
  19 #pragma once
  20
  21 #include "Shapes_Computation_decls.h"
  22
  23 #include <cstddef>
  24
  25 namespace Shapes
  26 {
  27         namespace Computation
  28         {
  29                 class BasicSimplex
  30                 {
  31                         size_t nVars_;
  32                         size_t nEqns_;
  33                         size_t nExt_;
  34
  35                         mutable double * a_;
  36                         mutable double * b_;
  37                         mutable double * c_;
  38
  39                         mutable size_t * varSet_;
  40                         mutable size_t * nonBasicBegin_;
  41                         mutable size_t * nonBasicEnd_;
  42                         mutable size_t * basicBegin_;
  43                         mutable size_t * basicEnd_;
  44
  45                 public:
  46                         BasicSimplex( size_t nVars, size_t nEqns );
  47                         ~BasicSimplex( );
  48
  49                         bool minimize( double * xdst,
  50                                                                                  double * objdst, double objGoal,
  51                                                                                  const double * c, const double * a, const double * b,
  52                                                                                  bool changeSign = false ) const;
  53                         bool maximize( double * xdst,
  54                                                                                  double * objdst, double objGoal,
  55                                                                                  const double * c, const double * a, const double * b ) const;
  56
  57                 protected:
  58                         bool phaseOne( ) const;
  59                         double phaseTwo( double * xdst, double objGoal ) const;
  60                 };
  61         }
  62 }