largeres1.py

Go to the documentation of this file.

 
import sys
import os
import math
 
import conopt as co
 
sys.path.append('../common/')
import std
 
 
class LargeResModelData(co.ModelData):
   def __init__(self):
      super().__init__()
 
   def buildModel(self):
      """
      adding the variables and constraints to the model
      @ingroup PYTHON1THREAD_LARGERES1
      """
      # adding the variables to the model
      self.varidx = self.addVariable(
         -co.Conopt.Infinity, co.Conopt.Infinity, 35.0
      )
 
      # adding the constraints to the model
      self.considx = self.addConstraint(co.ConstraintType_Eq, 10, [0], [0], [1])
 
      # setting the objective constraint
      self.setObjectiveElement(co.ObjectiveElement_Variable, self.varidx)
 
      # setting the optimisation direction
      self.setOptimizationSense(co.Sense_Maximize)
 
   def evaluateNonlinearTerm(self, x, rowno, ignerr, thread):
      """
      @copydoc conopt.ModelData.evaluateNonlinearTerm
      @ingroup PYTHON1THREAD_LARGERES1
      """
      g = 0
      if rowno == self.considx:
         g = math.exp(x[self.varidx])
 
      return g
 
   def evaluateNonlinearJacobian(self, x, rowno, jacnum, ignerr, thread):
      """
      @copydoc conopt.ModelData.evaluateNonlinearJacobian
      @ingroup PYTHON1THREAD_LARGERES1
      """
      jac = []
      if rowno == self.considx:
         jac.append(math.exp(x[self.varidx]))
 
      return jac
 
 
if __name__ == '__main__':
   name = os.path.basename(__file__)[:-3]
 
   conopt = co.Conopt(name)
   model = LargeResModelData()
   msghdlr = std.TutMessageHandler(name)
 
   model.buildModel()
 
   conopt.loadModel(model)
   conopt.setMessageHandler(msghdlr)
 
   # getting the license variables
   license_int_1 = os.environ.get('CONOPT_LICENSE_INT_1', None)
   license_int_2 = os.environ.get('CONOPT_LICENSE_INT_2', None)
   license_int_3 = os.environ.get('CONOPT_LICENSE_INT_3', None)
   license_text = os.environ.get('CONOPT_LICENSE_TEXT', None)
   if (
      license_int_1 is not None
      and license_int_2 is not None
      and license_int_3 is not None
      and license_text is not None
   ):
      conopt.setLicense(
         int(license_int_1),
         int(license_int_2),
         int(license_int_3),
         license_text,
      )
 
   coi_error = conopt.solve()
 
   if coi_error != 400:
      std.python_log(conopt, 'Expected return code 400, Saw')
 
   std.python_log(conopt, 'Successful Solve')