Changeset 1559 for trunk/Clp/examples


Ignore:
Timestamp:
Jun 5, 2010 3:42:36 PM (10 years ago)
Author:
stefan
Message:

merge split branch into trunk

Location:
trunk
Files:
32 edited

Legend:

Unmodified
Added
Removed
  • trunk

    • Property svn:externals
      •  

        old new  
        1 BuildTools    https://projects.coin-or.org/svn/BuildTools/stable/0.6
        2 ThirdParty/Blas https://projects.coin-or.org/svn/BuildTools/ThirdParty/Blas/stable/1.0
        3 ThirdParty/Lapack https://projects.coin-or.org/svn/BuildTools/ThirdParty/Lapack/stable/1.0
        4 Data/Netlib   https://projects.coin-or.org/svn/Data/stable/1.0/Netlib
        5 Data/Sample   https://projects.coin-or.org/svn/Data/stable/1.0/Sample
        6 CoinUtils     https://projects.coin-or.org/svn/CoinUtils/stable/2.6/CoinUtils
        7 
         1BuildTools    https://projects.coin-or.org/svn/BuildTools/trunk
         2ThirdParty/Blas https://projects.coin-or.org/svn/BuildTools/ThirdParty/Blas/trunk
         3ThirdParty/Lapack https://projects.coin-or.org/svn/BuildTools/ThirdParty/Lapack/trunk
         4#Data/Netlib   https://projects.coin-or.org/svn/Data/trunk/Netlib
         5Data/Sample   https://projects.coin-or.org/svn/Data/trunk/Sample
         6CoinUtils     https://projects.coin-or.org/svn/CoinUtils/trunk/CoinUtils
    • Property svn:mergeinfo changed
      /branches/split (added)merged: 1522,​1527-1531,​1542-1550,​1558
  • trunk/Clp/examples/Makefile.in

    r1265 r1559  
    4242
    4343# C++ Compiler options
    44 CXXFLAGS = @CXXFLAGS@
     44CXXFLAGS = @CXXFLAGS@ -DSAMPLEDIR=\"@SAMPLE_DATA@\" -DNETLIBDIR=\"@NETLIB_DATA@\"
    4545
    4646# additional C++ Compiler options for linking
     
    5151
    5252# C Compiler options
    53 CFLAGS = @CFLAGS@
     53CFLAGS = @CFLAGS@ -DSAMPLEDIR=\"@SAMPLE_DATA@\" -DNETLIBDIR=\"@NETLIB_DATA@\"
    5454
    5555# Directory with COIN header files
     
    5757
    5858# Directory with COIN libraries
    59 COINLIBDIR = @abs_lib_dir@
     59COINLIBDIR = @abs_lib_dir@/coin
    6060
    6161# Libraries necessary to link with Clp
    62 LIBS = -L$(COINLIBDIR) -lClp -lCoinUtils \
    63         @ADDLIBS@ `cat @COINUTILSDOCDIR@/coinutils_addlibs.txt`
     62LIBS = @COINUTILS_LIBS@ -L`$(CYGPATH_W) $(COINLIBDIR)` -lClp @ADDLIBS@
    6463
    6564# Necessary Include dirs (we use the CYGPATH_W variables to allow
    6665# compilation with Windows compilers)
    67 INCL = -I`$(CYGPATH_W) $(COININCDIR)` $(ADDINCFLAGS)
     66INCL = @COINUTILS_CFLAGS@ -I`$(CYGPATH_W) $(COININCDIR)` $(ADDINCFLAGS)
    6867
    6968# The following is necessary under cygwin, if native compilers are used
  • trunk/Clp/examples/addBits.cpp

    r1552 r1559  
    2525  (you can have as many different strings as you want)
    2626*/
     27#include "ClpConfig.h"
    2728#include "ClpSimplex.hpp"
    2829#include "CoinHelperFunctions.hpp"
     
    3233#include <cassert>
    3334
    34 int main (int argc, const char *argv[])
     35int main(int argc, const char *argv[])
    3536{
    3637     // Empty model
    3738     ClpSimplex  model;
    38      std::string mpsFileName = "../../Data/Netlib/25fv47.mps";
     39     std::string mpsFileName;
    3940     if (argc >= 2) mpsFileName = argv[1];
     41     else {
     42#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     43          mpsFileName = NETLIBDIR "/25fv47.mps";
     44#else
     45          fprintf(stderr, "Do not know where to find netlib MPS files.\n");
     46          exit(1);
     47#endif
     48     }
    4049     int status = model.readMps(mpsFileName.c_str(), true);
    4150
     
    160169     int numberErrors = model2.loadProblem(build);
    161170     // should fail as we never set multiplier
    162      assert (numberErrors);
     171     assert(numberErrors);
    163172     time3 = CoinCpuTime() - time2;
    164173     // subtract out unsuccessful times
     
    167176     build.associateElement("multiplier", 0.0);
    168177     numberErrors = model2.loadProblem(build);
    169      assert (!numberErrors);
     178     assert(!numberErrors);
    170179     time3 = CoinCpuTime();
    171180     printf("Time for build using CoinModel is %g (%g for successful loadproblem)\n", time3 - time1,
     
    176185          build.associateElement("multiplier", multiplier);
    177186          numberErrors = model2.loadProblem(build, true);
    178           assert (!numberErrors);
     187          assert(!numberErrors);
    179188          model2.dual();
    180189     }
  • trunk/Clp/examples/addColumns.cpp

    r1552 r1559  
    1313#include <cassert>
    1414
    15 int main (int argc, const char *argv[])
     15int main(int argc, const char *argv[])
    1616{
    1717     {
     
    6161          double time1 = CoinCpuTime();
    6262          int k;
    63           for ( k = 0; k < 10000; k++) {
     63          for (k = 0; k < 10000; k++) {
    6464               int column2Index[] = {0, 1};
    6565               double column2Value[] = {1.0, -5.0};
     
    7373          CoinBuild buildObject;
    7474          time1 = CoinCpuTime();
    75           for ( k = 0; k < 100000; k++) {
     75          for (k = 0; k < 100000; k++) {
    7676               int column2Index[] = {0, 1};
    7777               double column2Value[] = {1.0, -5.0};
     
    8888          CoinBuild buildObject2;
    8989          time1 = CoinCpuTime();
    90           for ( k = 0; k < 10000; k++) {
     90          for (k = 0; k < 10000; k++) {
    9191               int column2Index[] = {0, 1};
    9292               double column2Value[] = {1.0, 1.0, -1.0};
     
    103103          CoinModel modelObject2;
    104104          time1 = CoinCpuTime();
    105           for ( k = 0; k < 10000; k++) {
     105          for (k = 0; k < 10000; k++) {
    106106               int column2Index[] = {0, 1};
    107107               double column2Value[] = {1.0, 1.0, -1.0};
     
    118118          CoinModel modelObject;
    119119          time1 = CoinCpuTime();
    120           for ( k = 0; k < 100000; k++) {
     120          for (k = 0; k < 100000; k++) {
    121121               int column2Index[] = {0, 1};
    122122               double column2Value[] = {1.0, -5.0};
     
    183183          ClpSimplex  model;
    184184          int status;
    185           if (argc < 2)
    186                status = model.readMps("../../Data/Sample/p0033.mps");
    187           else
     185          if (argc < 2) {
     186#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     187               status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     188#else
     189               fprintf(stderr, "Do not know where to find sample MPS files.\n");
     190               exit(1);
     191#endif
     192          } else
    188193               status = model.readMps(argv[1]);
    189194          if (status) {
     
    227232          ClpSimplex  model;
    228233          int status;
    229           if (argc < 2)
    230                status = model.readMps("../../Data/Sample/p0033.mps");
    231           else
     234          if (argc < 2) {
     235#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     236               status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     237#else
     238               fprintf(stderr, "Do not know where to find sample MPS files.\n");
     239               exit(1);
     240#endif
     241          } else
    232242               status = model.readMps(argv[1]);
    233243          if (status) {
  • trunk/Clp/examples/addRows.cpp

    r1552 r1559  
    1313#include <cassert>
    1414
    15 int main (int argc, const char *argv[])
     15int main(int argc, const char *argv[])
    1616{
    1717     try {
     
    6767          double time1 = CoinCpuTime();
    6868          int k;
    69           for ( k = 0; k < 10000; k++) {
     69          for (k = 0; k < 10000; k++) {
    7070               int row2Index[] = {0, 1, 2};
    7171               double row2Value[] = {1.0, -5.0, 1.0};
     
    7979          CoinBuild buildObject;
    8080          time1 = CoinCpuTime();
    81           for ( k = 0; k < 10000; k++) {
     81          for (k = 0; k < 10000; k++) {
    8282               int row2Index[] = {0, 1, 2};
    8383               double row2Value[] = {1.0, -5.0, 1.0};
     
    9494          CoinBuild buildObject2;
    9595          time1 = CoinCpuTime();
    96           for ( k = 0; k < 10000; k++) {
     96          for (k = 0; k < 10000; k++) {
    9797               int row2Index[] = {0, 1, 2};
    9898               double row2Value[] = {1.0, -1.0, 1.0};
     
    108108          CoinModel modelObject2;
    109109          time1 = CoinCpuTime();
    110           for ( k = 0; k < 10000; k++) {
     110          for (k = 0; k < 10000; k++) {
    111111               int row2Index[] = {0, 1, 2};
    112112               double row2Value[] = {1.0, -1.0, 1.0};
     
    121121          CoinModel modelObject3;
    122122          time1 = CoinCpuTime();
    123           for ( k = 0; k < 10000; k++) {
     123          for (k = 0; k < 10000; k++) {
    124124               int row2Index[] = {0, 1, 2};
    125125               double row2Value[] = {1.0, -1.0, 1.0};
     
    135135          CoinModel modelObject;
    136136          time1 = CoinCpuTime();
    137           for ( k = 0; k < 10000; k++) {
     137          for (k = 0; k < 10000; k++) {
    138138               int row2Index[] = {0, 1, 2};
    139139               double row2Value[] = {1.0, -5.0, 1.0};
     
    266266          CoinBuild buildObject3;
    267267          time1 = CoinCpuTime();
    268           for ( k = 0; k < 10000; k++) {
     268          for (k = 0; k < 10000; k++) {
    269269               int row2Index[] = {0, 1, 2};
    270270               double row2Value[] = {1.0, -1.0, 1.0};
  • trunk/Clp/examples/decomp2.cpp

    r1552 r1559  
    33// Corporation and others.  All Rights Reserved.
    44
     5#include "ClpConfig.h"
    56#include "ClpSimplex.hpp"
    67#include "CoinStructuredModel.hpp"
    78#include <iomanip>
    89
    9 int main (int argc, const char *argv[])
     10int main(int argc, const char *argv[])
    1011{
    1112     /* Create a structured model by reading mps file and trying
     
    1314     */
    1415     // At present D-W rows are hard coded - will move stuff from OSL
    15      CoinStructuredModel model((argc < 2) ? "../../Data/Netlib/czprob.mps"
     16#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     17     CoinStructuredModel model((argc < 2) ? NETLIBDIR "/czprob.mps"
    1618                               : argv[1], 1);
     19#else
     20     if (argc<2) {
     21          fprintf(stderr, "Do not know where to find netlib MPS files.\n");
     22          return 1;
     23     }
     24     CoinStructuredModel model(argv[1], 1);
     25#endif
    1726     if (!model.numberRows())
    1827          exit(10);
  • trunk/Clp/examples/decomp3.cpp

    r1552 r1559  
    33// Corporation and others.  All Rights Reserved.
    44
     5#include "ClpConfig.h"
    56#include "ClpSimplex.hpp"
    67#include "ClpPresolve.hpp"
     
    910#include <iomanip>
    1011
    11 int main (int argc, const char *argv[])
     12int main(int argc, const char *argv[])
    1213{
    1314     /* Create a structured model by reading mps file and trying
     
    3132     //#define PRESOLVE
    3233#ifndef PRESOLVE
    33      CoinStructuredModel model((argc < 2) ? "../../Data/Netlib/czprob.mps"
     34#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     35     CoinStructuredModel model((argc < 2) ? NETLIBDIR "/czprob.mps"
    3436                               : argv[1], decompose, maxBlocks);
     37#else
     38     if (argc<2) {
     39          fprintf(stderr, "Do not know where to find netlib MPS files.\n");
     40          return 1;
     41     }
     42     CoinStructuredModel model(argv[1], 1);
     43#endif
    3544     if (!model.numberRows())
    3645          exit(10);
     
    4958#else
    5059     ClpSimplex  model;
    51      int status = model.readMps((argc < 2) ? "../../Data/Netlib/czprob.mps"
     60#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     61     int status = model.readMps((argc < 2) ? NETLIBDIR "/czprob.mps"
    5262                                : argv[1]);
     63#else
     64     if (argc<2) {
     65          fprintf(stderr, "Do not know where to find netlib MPS files.\n");
     66          return 1;
     67     }
     68     int status = model.readMps(argv[1]);
     69#endif
    5370     if (status) {
    5471          fprintf(stdout, "Bad readMps %s\n", argv[1]);
     
    7693     return 0;
    7794     ClpSimplex solver2;
    78      solver2.readMps((argc < 2) ? "../../Data/Netlib/czprob.mps" : argv[1]);
     95#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     96     solver2.readMps((argc < 2) ? NETLIBDIR "/czprob.mps" : argv[1]);
     97#else
     98     if (argc<2) {
     99          fprintf(stderr, "Do not know where to find netlib MPS files.\n");
     100          return 1;
     101     }
     102     solver2.readMps(argv[1]);
     103#endif
    79104     time1 = CoinCpuTime() ;
    80105     solver2.dual();
  • trunk/Clp/examples/decompose.cpp

    r1552 r1559  
    33// Corporation and others.  All Rights Reserved.
    44
     5#include "ClpConfig.h"
    56#include "ClpSimplex.hpp"
    67#include "CoinMpsIO.hpp"
    78#include <iomanip>
    89
    9 int main (int argc, const char *argv[])
     10int main(int argc, const char *argv[])
    1011{
    1112     ClpSimplex  model;
     
    1314     // Keep names
    1415     if (argc < 2) {
    15           //status=model.readMps("/home/forrest/data/ken-18.mps.gz",true);
    16           status = model.readMps("../../Data/Netlib/czprob.mps", true);
     16#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     17          status = model.readMps(NETLIBDIR "/czprob.mps", true);
     18#else
     19          fprintf(stderr, "Do not know where to find netlib MPS files.\n");
     20          return 1;
     21#endif
    1722     } else {
    1823          status = model.readMps(argv[1], true);
     
    115120                                   }
    116121                              } else {
    117                                    assert (columnBlock[iColumn] == numberBlocks - 1);
     122                                   assert(columnBlock[iColumn] == numberBlocks - 1);
    118123                              }
    119124                         }
     
    296301               dual = master.dualRowSolution();
    297302               memset(master.dualRowSolution(),
    298                       0, (numberMasterRows + numberBlocks)*sizeof(double));
     303                      0, (numberMasterRows + numberBlocks) *sizeof(double));
    299304          } else {
    300305               abort();
     
    442447               if (whichBlock[i-numberMasterColumns] == iBlock)
    443448                    sol[i] = solution[i];
    444           memset(lower, 0, (numberMasterRows + numberBlocks)*sizeof(double));
     449          memset(lower, 0, (numberMasterRows + numberBlocks) *sizeof(double));
    445450          master.times(1.0, sol, lower);
    446451          for (iRow = 0; iRow < numberMasterRows; iRow++) {
  • trunk/Clp/examples/defaults.cpp

    r1552 r1559  
    88#include <iomanip>
    99
    10 int main (int argc, const char *argv[])
     10int main(int argc, const char *argv[])
    1111{
    1212     ClpSimplex  model;
    1313     int status;
    1414     // Keep names
    15      if (argc < 2)
    16           status = model.readMps("../../Data/Sample/p0033.mps", true);
    17      else
     15     if (argc < 2) {
     16#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     17          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     18#else
     19          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     20          exit(1);
     21#endif
     22     } else
    1823          status = model.readMps(argv[1], true);
    1924     /*
     
    7075     model.setPrimalTolerance(value);
    7176
    72      model.setDualTolerance( model.dualTolerance()) ;
     77     model.setDualTolerance(model.dualTolerance()) ;
    7378
    7479     // Other Param stuff
  • trunk/Clp/examples/driver.cpp

    r1552 r1559  
    66#include <iomanip>
    77
    8 int main (int argc, const char *argv[])
     8int main(int argc, const char *argv[])
    99{
    1010     ClpSimplex  model;
    1111     int status;
    1212     // Keep names when reading an mps file
    13      if (argc < 2)
    14           status = model.readMps("../../Data/Sample/p0033.mps", true);
    15      else
     13     if (argc < 2) {
     14#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     15          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     16#else
     17          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     18          exit(1);
     19#endif
     20     } else
    1621          status = model.readMps(argv[1], true);
    1722
  • trunk/Clp/examples/driver2.cpp

    r1552 r1559  
    5555     MyMessageHandler(const CoinMessageHandler&);
    5656
    57      MyMessageHandler& operator=(const MyMessageHandler&);
     57     MyMessageHandler& operator= (const MyMessageHandler&);
    5858     /// Clone
    5959     virtual CoinMessageHandler * clone() const ;
     
    7878// Default Constructor
    7979//-------------------------------------------------------------------
    80 MyMessageHandler::MyMessageHandler ()
     80MyMessageHandler::MyMessageHandler()
    8181     : CoinMessageHandler(),
    8282       model_(NULL)
     
    8787// Copy constructor
    8888//-------------------------------------------------------------------
    89 MyMessageHandler::MyMessageHandler (const MyMessageHandler & rhs)
     89MyMessageHandler::MyMessageHandler(const MyMessageHandler & rhs)
    9090     : CoinMessageHandler(rhs),
    9191       model_(rhs.model_)
     
    9393}
    9494
    95 MyMessageHandler::MyMessageHandler (const CoinMessageHandler & rhs)
     95MyMessageHandler::MyMessageHandler(const CoinMessageHandler & rhs)
    9696     : CoinMessageHandler(),
    9797       model_(NULL)
     
    110110// Destructor
    111111//-------------------------------------------------------------------
    112 MyMessageHandler::~MyMessageHandler ()
     112MyMessageHandler::~MyMessageHandler()
    113113{
    114114}
     
    118118//-------------------------------------------------------------------
    119119MyMessageHandler &
    120 MyMessageHandler::operator=(const MyMessageHandler& rhs)
     120MyMessageHandler::operator= (const MyMessageHandler & rhs)
    121121{
    122122     if (this != &rhs) {
    123           CoinMessageHandler::operator=(rhs);
     123          CoinMessageHandler::operator= (rhs);
    124124          model_ = rhs.model_;
    125125     }
     
    174174}
    175175
    176 int main (int argc, const char *argv[])
     176int main(int argc, const char *argv[])
    177177{
    178178     ClpSimplex  model;
     
    183183     int status;
    184184     // Keep names when reading an mps file
    185      if (argc < 2)
    186           status = model.readMps("../../Data/Sample/p0033.mps", true);
    187      else
     185     if (argc < 2) {
     186#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     187          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     188#else
     189          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     190          exit(1);
     191#endif
     192     } else
    188193          status = model.readMps(argv[1], true);
    189194
  • trunk/Clp/examples/driverC.c

    r1552 r1559  
    55/* This example shows the use of the "C" interface */
    66
     7#include "ClpConfig.h"
    78#include "Clp_C_Interface.h"
    89#include <stdio.h>
     
    1920     if (messageNumber == 1000002) {
    2021          /* Coin0002 */
    21           assert (nString == 1 && nInt == 3);
     22          assert(nString == 1 && nInt == 3);
    2223          printf("Name of problem is %s\n", vString[0]);
    2324          printf("row %d col %d el %d\n", vInt[0], vInt[1], vInt[2]);
     
    2526          /* Clp0005 */
    2627          int i;
    27           assert (nInt == 4 && nDouble == 3); /* they may not all print */
     28          assert(nInt == 4 && nDouble == 3);    /* they may not all print */
    2829          for (i = 0; i < 3; i++)
    2930               printf("%d %g\n", vInt[i], vDouble[i]);
     
    3334
    3435
    35 int main (int argc, const char *argv[])
     36int main(int argc, const char *argv[])
    3637{
    3738     /* Get default model */
     
    4142     Clp_registerCallBack(model, callBack);
    4243     /* Keep names when reading an mps file */
    43      if (argc < 2)
    44           status = Clp_readMps(model, "../../Data/Sample/p0033.mps", 1, 0);
    45      else
     44     if (argc < 2) {
     45#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     46          status = Clp_readMps(model, SAMPLEDIR "/p0033.mps", 1, 0);
     47#else
     48          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     49          exit(1);
     50#endif
     51     } else
    4652          status = Clp_readMps(model, argv[1], 1, 0);
    4753
  • trunk/Clp/examples/dualCuts.cpp

    r1552 r1559  
    2020#include <iomanip>
    2121
    22 int main (int argc, const char *argv[])
     22int main(int argc, const char *argv[])
    2323{
    2424     ClpSimplex  model;
  • trunk/Clp/examples/ekk.cpp

    r1552 r1559  
    4949#include <stdarg.h>
    5050
    51 int main (int argc, const char *argv[])
     51int main(int argc, const char *argv[])
    5252{
    5353     const char * name;
     
    6262     EKKContext * context;
    6363
    64      if ( argc > 1 ) {
     64     if (argc > 1) {
    6565          name = argv[1];
    6666     } else {
    67           name = "../../Data/Sample/p0033.mps";
    68      }
     67#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     68          name = (SAMPLEDIR "/p0033.mps";
     69#else
     70          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     71          exit(1);
     72#endif
     73             }
    6974
    70      /* initialize OSL environment */
    71      context = ekk_initializeContext();
     75             /* initialize OSL environment */
     76             context = ekk_initializeContext();
    7277     model = ekk_newModel(context, "");
    7378
     
    143148          }
    144149          if (useosl == 2)
    145                ekk_allSlackBasis(model); // otherwise it would be easy
     150               ekk_allSlackBasis(model);    // otherwise it would be easy
    146151     }
    147152     if ((useosl & 2) == 0) {
  • trunk/Clp/examples/ekk_interface.cpp

    r1552 r1559  
    168168     int numberRows = ekk_getInumrows(model);
    169169     int numberColumns = ekk_getInumcols(model);
    170      int * rowStatus =  (int *) ekk_rowstat(model);
     170     int * rowStatus = (int *) ekk_rowstat(model);
    171171     double * rowSolution = (double *) ekk_rowacts(model);
    172172     int i;
    173173     int * columnStatus = (int *) ekk_colstat(model);
    174      double * columnSolution =  (double *) ekk_colsol(model);
     174     double * columnSolution = (double *) ekk_colsol(model);
    175175     memcpy(rowSolution, clp->primalRowSolution(), numberRows * sizeof(double));
    176176     memcpy(rowDual, clp->dualRowSolution(), numberRows * sizeof(double));
     
    200200     ekk_setRsumdinf(model, clp->sumDualInfeasibilities());
    201201     delete clp;
    202      if (compressInfo )
     202     if (compressInfo)
    203203          ekk_decompressModel(model, compressInfo);
    204204
  • trunk/Clp/examples/hello.cpp

    r1552 r1559  
    1717#include "ClpSimplex.hpp"
    1818
    19 int main (int argc, const char *argv[])
     19int main(int argc, const char *argv[])
    2020{
    2121     ClpSimplex  model;
  • trunk/Clp/examples/makeDual.cpp

    r1552 r1559  
    66#include "CoinHelperFunctions.hpp"
    77#include "CoinBuild.hpp"
    8 int main (int argc, const char *argv[])
     8int main(int argc, const char *argv[])
    99{
    1010     ClpSimplex  model;
    1111     int status;
    12      if (argc < 2)
    13           status = model.readMps("../../Data/Sample/p0033.mps");
    14      else
     12     if (argc < 2) {
     13#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     14          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     15#else
     16          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     17          exit(1);
     18#endif
     19     } else
    1520          status = model.readMps(argv[1]);
    1621     if (status) {
     
    121126     for (iRow = 0; iRow < numberRows; iRow++) {
    122127          if (rowLower[iRow] <= -1.0e20) {
    123                assert (rowUpper[iRow] < 1.0e20);
     128               assert(rowUpper[iRow] < 1.0e20);
    124129               newObjective[kRow] = -rowUpper[iRow];
    125130               fromRowsLower[kRow] = -COIN_DBL_MAX;
  • trunk/Clp/examples/minimum.cpp

    r1552 r1559  
    44
    55#include "ClpSimplex.hpp"
    6 int main (int argc, const char *argv[])
     6int main(int argc, const char *argv[])
    77{
    88     ClpSimplex  model;
    99     int status;
    10      if (argc < 2)
    11           status = model.readMps("../../Data/Sample/p0033.mps");
    12      else
     10     if (argc < 2) {
     11#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     12          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     13#else
     14          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     15          exit(1);
     16#endif
     17     } else
    1318          status = model.readMps(argv[1]);
    1419     if (!status) {
  • trunk/Clp/examples/modify.cpp

    r1552 r1559  
    88#include "ClpSimplex.hpp"
    99#include "CoinHelperFunctions.hpp"
    10 int main (int argc, const char *argv[])
     10int main(int argc, const char *argv[])
    1111{
    1212     ClpSimplex  model;
  • trunk/Clp/examples/myPdco.cpp

    r1552 r1559  
    2121// Default Constructor
    2222//-------------------------------------------------------------------
    23 myPdco::myPdco ()
     23myPdco::myPdco()
    2424     : ClpPdcoBase(),
    2525       rowIndex_(NULL),
     
    4545// Copy constructor
    4646//-------------------------------------------------------------------
    47 myPdco::myPdco (const myPdco & rhs)
     47myPdco::myPdco(const myPdco & rhs)
    4848     : ClpPdcoBase(rhs),
    4949       numlinks_(rhs.numlinks_),
     
    5656// Destructor
    5757//-------------------------------------------------------------------
    58 myPdco::~myPdco ()
     58myPdco::~myPdco()
    5959{
    6060     delete [] rowIndex_;
     
    6565//-------------------------------------------------------------------
    6666myPdco &
    67 myPdco::operator=(const myPdco& rhs)
     67myPdco::operator= (const myPdco & rhs)
    6868{
    6969     if (this != &rhs) {
    70           ClpPdcoBase::operator=(rhs);
     70          ClpPdcoBase::operator= (rhs);
    7171          numlinks_ = rhs.numlinks_;
    7272          numnodes_ = rhs.numnodes_;
     
    9191               y_sum += y_elts[k];
    9292               int i1 = rowIndex_[2*k];
    93                assert (i1 >= 0);
     93               assert(i1 >= 0);
    9494               x_elts[i1] += y_elts[k];
    9595               int i2 = rowIndex_[2*k+1];
    96                assert (i2 >= 0);
     96               assert(i2 >= 0);
    9797               x_elts[i2] -= y_elts[k];
    9898          }
     
    216216     int nonzpt = 0;
    217217     numlinks_ = 0;
    218      while(fscanf(fpData, "%d %d", ifrom, ito) && kct++ < maxlinks) {
     218     while (fscanf(fpData, "%d %d", ifrom, ito) && kct++ < maxlinks) {
    219219          //  while(fread(ipair, 4,2, fpData) && kct++ < maxlinks){
    220220          if ((*ifrom) < 0) {
     
    240240     fclose(fpData);
    241241     fclose(fpParam);
    242      printf ("imax %d  imin %d\n", imax, imin);
     242     printf("imax %d  imin %d\n", imax, imin);
    243243     // Set model size
    244244     numnodes_ = imax + 1;
  • trunk/Clp/examples/network.cpp

    r1552 r1559  
    1818#include <cmath>
    1919
    20 int main (int argc, const char *argv[])
     20int main(int argc, const char *argv[])
    2121{
    2222     int numberColumns;
     
    2424
    2525     FILE * fp;
    26      if ( argc > 1 ) {
     26     if (argc > 1) {
    2727          fp = fopen(argv[1], "r");
    2828          if (!fp) {
     
    4141     // read and skip
    4242     fscanf(fp, "%s", temp);
    43      assert (!strcmp(temp, "BEGIN"));
     43     assert(!strcmp(temp, "BEGIN"));
    4444     fscanf(fp, "%*s %*s %d %d %*s %*s %d %*s", &problem, &numberRows,
    4545            &numberColumns);
  • trunk/Clp/examples/pdco.cpp

    r1552 r1559  
    1616#include <cmath>
    1717
    18 int main (int argc, const char *argv[])
     18int main(int argc, const char *argv[])
    1919{
    2020
     
    5757     // These are already set?
    5858     model.xsize_ = 50.0 / (model.numberColumns());
    59      model.xsize_ = CoinMin( 1.0, model.xsize_);
     59     model.xsize_ = CoinMin(1.0, model.xsize_);
    6060
    6161     /*
    6262      *     Solve the test problem
    6363      */
    64      model.pdco( &stuff, options, info, outfo);
     64     model.pdco(&stuff, options, info, outfo);
    6565
    6666     /*
  • trunk/Clp/examples/piece.cpp

    r1552 r1559  
    1111#include <iomanip>
    1212
    13 int main (int argc, const char *argv[])
     13int main(int argc, const char *argv[])
    1414{
    1515     int status;
     
    115115               double obj2 = objective1[iColumn+1];
    116116               double sol2 = oldSolution[iColumn+1];
    117                if(fabs(up1 - lo2) > 1.0e-8) {
     117               if (fabs(up1 - lo2) > 1.0e-8) {
    118118                    // try other way
    119119                    double temp;
     
    202202     // Create nonlinear objective
    203203     int returnCode = model.createPiecewiseLinearCosts(segstart, breakpt, slope);
    204      assert (!returnCode);
     204     assert(!returnCode);
    205205
    206206     // delete
  • trunk/Clp/examples/rowColumn.cpp

    r1552 r1559  
    1515#include <cassert>
    1616
    17 int main (int argc, const char *argv[])
     17int main(int argc, const char *argv[])
    1818{
    1919     ClpSimplex  modelByRow, modelByColumn;
     
    6161     modelByColumn.dual();
    6262     // check value of objective
    63      assert (fabs(modelByColumn.objectiveValue() - 76000.0) < 1.0e-7);
     63     assert(fabs(modelByColumn.objectiveValue() - 76000.0) < 1.0e-7);
    6464     // b) As two steps - first creating a CoinPackedMatrix
    6565     // NULL for column lengths indicate they are stored without gaps
     
    7575     modelByColumn.dual();
    7676     // check value of objective
    77      assert (fabs(modelByColumn.objectiveValue() - 76000.0) < 1.0e-7);
     77     assert(fabs(modelByColumn.objectiveValue() - 76000.0) < 1.0e-7);
    7878     // Now do by row
    7979     // The false says row ordered so numberRows and numberColumns swapped - see CoinPackedMatrix.hpp
    80      assert(numberElements == rowStart[numberRows]); // check same number of elements in each copy
     80     assert(numberElements == rowStart[numberRows]);    // check same number of elements in each copy
    8181     CoinPackedMatrix byRow(false, numberColumns, numberRows, numberElements,
    8282                            elementByRow, column, rowStart, NULL);
     
    8888     modelByRow.dual();
    8989     // check value of objective
    90      assert (fabs(modelByRow.objectiveValue() - 76000.0) < 1.0e-7);
     90     assert(fabs(modelByRow.objectiveValue() - 76000.0) < 1.0e-7);
    9191     // write solution
    9292     const double * solution = modelByRow.primalColumnSolution();
  • trunk/Clp/examples/sprint.cpp

    r1552 r1559  
    77#include <iomanip>
    88
    9 int main (int argc, const char *argv[])
     9int main(int argc, const char *argv[])
    1010{
    1111     ClpSimplex  model;
     
    6262          // now see what that does to row solution
    6363          double * rowSolution = model.primalRowSolution();
    64           memset (rowSolution, 0, numberRows * sizeof(double));
     64          memset(rowSolution, 0, numberRows * sizeof(double));
    6565          model.times(1.0, columnSolution, rowSolution);
    6666
     
    135135          // now see what variables left out do to row solution
    136136          double * rowSolution = model.primalRowSolution();
    137           memset (rowSolution, 0, numberRows * sizeof(double));
     137          memset(rowSolution, 0, numberRows * sizeof(double));
    138138          int iRow, iColumn;
    139139          // zero out ones in small problem
  • trunk/Clp/examples/sprint2.cpp

    r1552 r1559  
    88#include <iomanip>
    99
    10 int main (int argc, const char *argv[])
     10int main(int argc, const char *argv[])
    1111{
    1212     ClpSimplex  model;
     
    8484          if (iPass) {
    8585               double ratio = ((double) smallNumberRows) / ((double) model2->numberRows());
    86                smallNumberColumns = (int) (smallNumberColumns * ratio);
     86               smallNumberColumns = (int)(smallNumberColumns * ratio);
    8787          }
    8888          delete model2;
  • trunk/Clp/examples/testBarrier.cpp

    r1552 r1559  
    77#include "ClpCholeskyWssmp.hpp"
    88#include "ClpCholeskyDense.hpp"
    9 int main (int argc, const char *argv[])
     9int main(int argc, const char *argv[])
    1010{
    1111     ClpInterior  model;
    1212     int status;
    13      if (argc < 2)
    14           status = model.readMps("../../Data/Sample/p0033.mps");
    15      else
     13     if (argc < 2) {
     14#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     15          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     16#else
     17          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     18          exit(1);
     19#endif
     20     } else
    1621          status = model.readMps(argv[1]);
    1722     if (status) {
  • trunk/Clp/examples/testBasis.cpp

    r1552 r1559  
    88#include <cassert>
    99
    10 int main (int argc, const char *argv[])
     10int main(int argc, const char *argv[])
    1111{
    1212     ClpSimplex  model;
  • trunk/Clp/examples/testGub.cpp

    r1552 r1559  
    33// Corporation and others.  All Rights Reserved.
    44
     5#include "ClpConfig.h"
    56#include "ClpSimplex.hpp"
    67#include "ClpDynamicExampleMatrix.hpp"
     
    1011#include "CoinTime.hpp"
    1112#include "CoinMpsIO.hpp"
    12 int main (int argc, const char *argv[])
     13int main(int argc, const char *argv[])
    1314{
    1415     ClpSimplex  model;
     
    1617     int maxFactor = 100;
    1718     if (argc < 2) {
    18           status = model.readMps("../../Data/Netlib/czprob.mps");
     19#if defined(COIN_HAS_NETLIB) && defined(NETLIBDIR)
     20          status = model.readMps(NETLIBDIR "/czprob.mps");
    1921          if (status) {
    2022               printf("Unable to read matrix - trying gzipped version\n");
    21                status = model.readMps("../../Data/Netlib/czprob.mps.gz");
    22           }
     23               status = model.readMps(NETLIBDIR "/czprob.mps.gz");
     24          }
     25#else
     26          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     27          exit(1);
     28#endif
    2329     } else {
    2430          status = model.readMps(argv[1]);
     
    245251          printf("** While after adding matrix there are %d rows and %d columns\n",
    246252                 model2.numberRows(), model2.numberColumns());
    247           model2.setSpecialOptions(4); // exactly to bound
     253          model2.setSpecialOptions(4);    // exactly to bound
    248254          // For now scaling off
    249255          model2.scaling(0);
     
    261267               ClpDynamicMatrix * oldMatrix =
    262268                    dynamic_cast< ClpDynamicMatrix*>(model2.clpMatrix());
    263                assert (oldMatrix);
     269               assert(oldMatrix);
    264270               ClpDynamicMatrix * newMatrix = new
    265271               ClpDynamicMatrix(&model3, numberGub,
     
    272278               // and ordinary status (but only NON gub rows)
    273279               memcpy(model3.statusArray(), model2.statusArray(),
    274                       (newMatrix->numberStaticRows() + model3.numberColumns())*sizeof(unsigned char));
     280                      (newMatrix->numberStaticRows() + model3.numberColumns()) *sizeof(unsigned char));
    275281               newMatrix->switchOffCheck();
    276282               newMatrix->setRefreshFrequency(1000);
     
    278284               ClpDynamicExampleMatrix * oldMatrix =
    279285                    dynamic_cast< ClpDynamicExampleMatrix*>(model2.clpMatrix());
    280                assert (oldMatrix);
     286               assert(oldMatrix);
    281287               ClpDynamicExampleMatrix * newMatrix = new
    282288               ClpDynamicExampleMatrix(&model3, numberGub,
     
    290296               // and ordinary status (but only NON gub rows)
    291297               memcpy(model3.statusArray(), model2.statusArray(),
    292                       (newMatrix->numberStaticRows() + model3.numberColumns())*sizeof(unsigned char));
     298                      (newMatrix->numberStaticRows() + model3.numberColumns()) *sizeof(unsigned char));
    293299               newMatrix->switchOffCheck();
    294300               newMatrix->setRefreshFrequency(1000);
    295301          }
    296           model3.setSpecialOptions(4); // exactly to bound
     302          model3.setSpecialOptions(4);    // exactly to bound
    297303          // For now scaling off
    298304          model3.scaling(0);
     
    313319          ClpDynamicMatrix * gubMatrix =
    314320               dynamic_cast< ClpDynamicMatrix*>(model2.clpMatrix());
    315           assert (gubMatrix);
     321          assert(gubMatrix);
    316322          ClpDynamicExampleMatrix * gubMatrix2 =
    317323               dynamic_cast< ClpDynamicExampleMatrix*>(model2.clpMatrix());
     
    330336               const int * id = gubMatrix->id();
    331337               int i;
    332                const float * columnLower = gubMatrix->columnLower();
    333                const float * columnUpper = gubMatrix->columnUpper();
     338               const double * columnLower = gubMatrix->columnLower();
     339               const double * columnUpper = gubMatrix->columnUpper();
    334340               for (i = 0; i < numberGubColumns; i++) {
    335341                    if (gubMatrix->getDynamicStatus(i) == ClpDynamicMatrix::atUpperBound) {
     
    409415               }
    410416               // Coding below may not work if gub rows not at end
    411                FILE * fp = fopen ("xx.sol", "w");
     417               FILE * fp = fopen("xx.sol", "w");
    412418               fwrite(gubSolution, sizeof(double), numberTotalColumns, fp);
    413419               fwrite(status, sizeof(char), numberTotalColumns, fp);
     
    417423               model.times(1.0, gubSolution, rowsol2);
    418424               for (i = 0; i < numberStaticRows; i++)
    419                     assert (fabs(rowsol[i] - rowsol2[i]) < 1.0e-3);
     425                    assert(fabs(rowsol[i] - rowsol2[i]) < 1.0e-3);
    420426               for (; i < originalNumberRows; i++)
    421                     assert (rowsol2[i] > lower[i-numberStaticRows] - 1.0e-3 &&
    422                             rowsol2[i] < upper[i-numberStaticRows] + 1.0e-3);
     427                    assert(rowsol2[i] > lower[i-numberStaticRows] - 1.0e-3 &&
     428                           rowsol2[i] < upper[i-numberStaticRows] + 1.0e-3);
    423429               //for (;i<originalNumberRows;i++)
    424430               //printf("%d %g\n",i,rowsol2[i]);
  • trunk/Clp/examples/testGub2.cpp

    r1552 r1559  
    1010#include "CoinTime.hpp"
    1111#include "CoinMpsIO.hpp"
    12 int main (int argc, const char *argv[])
     12int main(int argc, const char *argv[])
    1313{
    1414     ClpSimplex  model;
     
    1616     int maxIts = 0;
    1717     int maxFactor = 100;
    18      if (argc < 2)
    19           status = model.readMps("../../Data/Sample/p0033.mps");
    20      else
     18     if (argc < 2) {
     19#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     20          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     21#else
     22          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     23          exit(1);
     24#endif
     25     } else
    2126          status = model.readMps(argv[1]);
    2227     if (status) {
     
    259264               ClpGubDynamicMatrix * gubMatrix =
    260265                    dynamic_cast< ClpGubDynamicMatrix*>(model2.clpMatrix());
    261                assert (gubMatrix);
     266               assert(gubMatrix);
    262267               const double * solution = model2.primalColumnSolution();
    263268               int numberGubColumns = gubMatrix->numberGubColumns();
     
    271276               const int * id = gubMatrix->id();
    272277               int i;
    273                const float * lowerColumn = gubMatrix->lowerColumn();
    274                const float * upperColumn = gubMatrix->upperColumn();
     278               const double * lowerColumn = gubMatrix->lowerColumn();
     279               const double * upperColumn = gubMatrix->upperColumn();
    275280               for (i = 0; i < numberGubColumns; i++) {
    276281                    if (gubMatrix->getDynamicStatus(i) == ClpGubDynamicMatrix::atUpperBound) {
     
    350355                         abort();
    351356               }
    352                FILE * fp = fopen ("xx.sol", "w");
     357               FILE * fp = fopen("xx.sol", "w");
    353358               fwrite(gubSolution, sizeof(double), numberTotalColumns, fp);
    354359               fwrite(status, sizeof(char), numberTotalColumns, fp);
     
    359364               model.times(1.0, gubSolution, rowsol2);
    360365               for (i = 0; i < numberRows; i++)
    361                     assert (fabs(rowsol[i] - rowsol2[i]) < 1.0e-3);
     366                    assert(fabs(rowsol[i] - rowsol2[i]) < 1.0e-3);
    362367               //for (;i<originalNumberRows;i++)
    363368               //printf("%d %g\n",i,rowsol2[i]);
  • trunk/Clp/examples/testQP.cpp

    r1552 r1559  
    99#include "ClpQuadraticObjective.hpp"
    1010#include <cassert>
    11 int main (int argc, const char *argv[])
     11int main(int argc, const char *argv[])
    1212{
    1313     /* Read quadratic model in two stages to test loadQuadraticObjective.
     
    2222     if (argc < 2) {
    2323          CoinMpsIO  m;
    24           int status = m.readMps("../../Data/Sample/share2qp", "mps");
     24#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     25          int status = m.readMps(SAMPLEDIR "/share2qp", "mps");
     26#else
     27          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     28          exit(1);
     29#endif
    2530          if (status) {
    2631               printf("errors on input\n");
     
    6166          int numberColumns = model.numberColumns();
    6267          int numberRows = model.numberRows();
    63           for ( i = 0; i < numberColumns; i++) {
     68          for (i = 0; i < numberColumns; i++) {
    6469               if (fabs(primal[i]) > 1.0e-8)
    6570                    printf("%d primal %g\n", i, primal[i]);
    6671          }
    67           for ( i = 0; i < numberRows; i++) {
     72          for (i = 0; i < numberRows; i++) {
    6873               if (fabs(dual[i]) > 1.0e-8)
    6974                    printf("%d dual %g\n", i, dual[i]);
     
    103108               const ClpQuadraticObjective * quadraticObj =
    104109                    (dynamic_cast<const ClpQuadraticObjective*>(model.objectiveAsObject()));
    105                assert (quadraticObj);
     110               assert(quadraticObj);
    106111               CoinPackedMatrix * quad = quadraticObj->quadraticObjective();
    107112               const int * columnQuadratic = quad->getIndices();
     
    115120               memcpy(gradient, quadraticObj->linearObjective(), numberColumns * sizeof(double));
    116121               int iColumn;
    117                for ( iColumn = 0; iColumn < numberColumns; iColumn++) {
     122               for (iColumn = 0; iColumn < numberColumns; iColumn++) {
    118123                    double valueI = primal[iColumn];
    119124                    CoinBigIndex j;
     
    149154               const int * columnLength = matrix->getVectorLengths();
    150155               const double * element = matrix->getElements();
    151                for ( iColumn = 0; iColumn < numberColumns; iColumn++) {
     156               for (iColumn = 0; iColumn < numberColumns; iColumn++) {
    152157                    double dj = gradient[iColumn];
    153158                    CoinBigIndex j;
     
    158163                    }
    159164                    if (model.getColumnStatus(iColumn) == ClpSimplex::basic) {
    160                          assert (fabs(dj) < 1.0e-5);
     165                         assert(fabs(dj) < 1.0e-5);
    161166                    } else {
    162                          assert (dj > -1.0e-5);
     167                         assert(dj > -1.0e-5);
    163168                    }
    164169               }
  • trunk/Clp/examples/useVolume.cpp

    r1552 r1559  
    1212private:
    1313     lpHook(const lpHook&);
    14      lpHook& operator=(const lpHook&);
     14     lpHook& operator= (const lpHook&);
    1515private:
    1616     /// Pointer to dense vector of structural variable upper bounds
     
    157157//#############################################################################
    158158
    159 int main (int argc, const char *argv[])
     159int main(int argc, const char *argv[])
    160160{
    161161     ClpSimplex  model;
    162162     int status;
    163163
    164      if (argc < 2)
    165           status = model.readMps("../../Data/Sample/p0033.mps", true);
    166      else
     164     if (argc < 2) {
     165#if defined(COIN_HAS_SAMPLE) && defined(SAMPLEDIR)
     166          status = model.readMps(SAMPLEDIR "/p0033.mps", true);
     167#else
     168          fprintf(stderr, "Do not know where to find sample MPS files.\n");
     169          exit(1);
     170#endif
     171     } else
    167172          status = model.readMps(argv[1], true);
    168173     /*
Note: See TracChangeset for help on using the changeset viewer.