largeres2.cpp

Go to the documentation of this file.

 
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <iostream>
#include "conopt.hpp"
 
class LR2_ModelData : public ConoptModelData
{
public:
   int varidx;
   int considx;
 
   LR2_ModelData() {}

   void buildModel()
   {
      /*                                                               */
      /*   Information about Variables:                                */
      /*   Default: Lower = -Inf, Curr = 0, and Upper = +inf.          */
      /*   Default: the status information in Vsta is not used.        */
      /*                                                               */
      /*   Initial value = 100.0                                       */
      /*                                                               */
      varidx = addVariable(-Conopt::Infinity, Conopt::Infinity, 100.0);
 
      /*                                                               */
      /*   Information about Constraints:                              */
      /*   Default: Rhs = 0                                            */
      /*   Default: the status information in Esta and the function    */
      /*            value in FV are not used.                          */
      /*   Default: Type: There is no default.                         */
      /*            0 = Equality,                                      */
      /*            1 = Greater than or equal,                         */
      /*            2 = Less than or equal,                            */
      /*            3 = Non binding.                                   */
      /*                                                               */
      /*   Constraint 0                                                */
      /*   Rhs = 10.0 and type Equality                                */
      /*                                                               */
      considx = addConstraint(ConoptConstraintType::Eq, 10, {0}, {0}, {1});
      /*                                                               */
      /*   Nonlinearity Structure: L = 0 are linear and NL = 1 are nonlinear */
      /*        X[0]                                                   */
      /*   0:    NL                                                    */
      /*                                                               */
 
      /* setting the objective variable */
      setObjectiveElement(ConoptObjectiveElement::Variable, varidx);
 
      /* setting the optimisation direction */
      setOptimizationSense(ConoptSense::Maximize);
   }

   int FDEval(const double x[], double *g, double jac[], int rowno, const int jacnum[], int mode,
         int ignerr, int *errcnt, int numvar, int numjac, int thread) override
   {
      /*                                                                   */
      /*   Row 0: exp(x) = 10                                              */
      /*                                                                   */
 
      if (rowno == considx)
      {
         /*                                                                   */
         /*   Mode = 1 or 3. Function value: g = P * Out                      */
         /*                                                                   */
         if (mode == 1 || mode == 3)
            *g = exp(x[varidx]);
         /*                                                                   */
         /*   Mode = 2 or 3: Derivative values:                               */
         /*                                                                   */
         if (mode == 2 || mode == 3)
         {
            jac[varidx] = exp(x[varidx]);
         }
      }
      return 0;
   }
};
 
#include "std.cpp"

int main(int argc, char **argv)
{
   int COI_Error = 0;
 
   // getting the program name from the executable path
   std::string pname = getProgramName(argv[0]);
 
   // initialising the Conopt Object
   Conopt conopt(pname);
   LR2_ModelData modeldata;
   Tut_MessageHandler msghandler(pname);
 
   // adding the message handler to the conopt interface
   conopt.setMessageHandler(msghandler);
 
   // building the model
   modeldata.buildModel();
 
   // loading the model in the conopt object
   conopt.loadModel(modeldata);
 
   // Disable debug output for function values
   conopt.debugFV(0);
 
#if defined(CONOPT_LICENSE_INT_1) && defined(CONOPT_LICENSE_INT_2) && defined(CONOPT_LICENSE_INT_3) && defined(CONOPT_LICENSE_TEXT)
   std::string license = CONOPT_LICENSE_TEXT;
   COI_Error += conopt.setLicense(CONOPT_LICENSE_INT_1, CONOPT_LICENSE_INT_2, CONOPT_LICENSE_INT_3, license);
#endif
 
   if (COI_Error)
      cpp_log(
            conopt, "Skipping COI_Solve due to license error. COI_Error = " + std::to_string(COI_Error), COI_Error);
 
   COI_Error = conopt.solve(); /* Optimize                      */
 
   printf("After solving. COI_Error =%d\n", COI_Error);
   if (COI_Error != 400)
   {
      cpp_log(conopt, "Expected return code 400, Saw", COI_Error);
   }
 
   cpp_log(conopt, "Successful Solve", 0);
}