square2.java

Go to the documentation of this file.

 
import java.util.*;
import java.lang.Math;
import conopt.*;

public class square2 {
   public static void main(String argv[]) {
      System.loadLibrary("conoptjni4");
 
      String name = "square2";
 
      Conopt conopt = new Conopt(name);
      Square2ModelData model = new Square2ModelData();
      MsgHandler msghdlr = new MsgHandler(name);
 
      // adding the message handler to the conopt interface
      conopt.setMessageHandler(msghdlr);
 
      // building the model
      model.buildModel();
 
      // loading the model in the conopt object
      conopt.loadModel(model);
 
      // tell conopt this is a square2 system
      conopt.squareModel(1);
 
      // try to set the license using the environment variables
      try {
         int license_int_1 = Integer.parseInt(System.getenv("CONOPT_LICENSE_INT_1"));
         int license_int_2 = Integer.parseInt(System.getenv("CONOPT_LICENSE_INT_2"));
         int license_int_3 = Integer.parseInt(System.getenv("CONOPT_LICENSE_INT_3"));
         String license_text = System.getenv("CONOPT_LICENSE_TEXT");
 
         conopt.setLicense(license_int_1, license_int_2, license_int_3,
               license_text);
      } catch (Exception e) {
         System.out.println("Unable to set license: " + e.getMessage());
      }
 
      int coi_error = conopt.solve();
      std s = new std();
      int retcode = 0;
 
      if (coi_error != 0) {
         s.java_log(name, "Errors encountered during solution: " + coi_error);
         retcode = -1;
      }
      // checking the status and objective value
      else if (conopt.solutionStatus() != 1 || conopt.modelStatus() < 4 || conopt.modelStatus() > 5) {
         s.java_log(name, "Solver or Model status not as expected (1,4) or (1,5)");
         retcode = -1;
      }
      else {
         // printing the final status of the optimisation
         conopt.printStatus();
         s.java_log(name, "Successful Solve");
      }
 
      msghdlr.close();
      System.exit(retcode);
   }
}
 
class Square2ModelData extends ModelData {
 
   private int x0;
   private int x1;
   private int cons1;
   private int cons2;
   private int cons3;
 
   public Square2ModelData() {
      super();
   }

   public void buildModel() {
      //   Model  x0 + x1 = 10
      //          x0 - x1 = 0
      //          x0 + x1 <= 9
      // adding the variables to the model
      x0 = addVariable(-Conopt.Infinity, Conopt.Infinity, 0.0);
      x1 = addVariable(-Conopt.Infinity, Conopt.Infinity, 0.0);
 
      // adding the constraints to the model
      // Constraint 1
      // Rhs = 10 and type Equality
      {
         int[] index = {x0, x1};
         double[] value = {1, 1};
         int[] nlflag = {0, 0};
 
         cons1 = addConstraint(ConstraintType.Eq, 10.0, index, value, nlflag);
      }
      
      // Constraint 2
      // Rhs = 0 and type Equality
      {
         int[] index = {x0, x1};
         double[] value = {0, 0};
         int[] nlflag = {1, 1}; 
 
 
         cons2 = addConstraint(ConstraintType.Eq, 0.0, index, value, nlflag);
      }
 
      // Constraint 3
      // Rhs = 9 and type Less Than Inequality
      {
         int[] index = {x0, x1};
         double[] value = {0, 0};
         int[] nlflag = {1, 1}; 
 
 
         cons3 = addConstraint(ConstraintType.Eq, 9.0, index, value, nlflag);
      }
   }

   public double evaluateNonlinearTerm(double[] x, int rowno, boolean ignerr, int thread) {
      double g = 0;
      if (rowno == cons1) {}
 
      else if (rowno == cons2) {
         g = x[x0] - x[x1] ;
      }
      else if (rowno == cons3) {
         g = x[x0] + x[x1] ;
      }
 
      return g;
   }

   public void evaluateNonlinearJacobian(double[] x, double[] jac, int rowno, int[] jacnum, boolean ignerr, int thread) {
      assert x.length == jac.length;
 
      if (rowno == cons1) {}
      
      else if (rowno == cons2) {
         jac[x0] = 1 ;
         jac[x1] = -1 ;
      }
      else if (rowno == cons2) {
         jac[x0] = 1 ;
         jac[x1] = 1 ;
      }
   }
}