Function main

snopt-interface/cppexamples/catmixa.cpp:107–293

int main(int argc, char **argv) {
  snoptProblemA catmixa;

  int Cold  = 0;

  int nH    = 1000;
  int n     = 3*(nH+1);
  int nCon  = 2*nH;
  int nF    = nCon + 1;

  int    ObjRow;
  double ObjAdd = -1.0;

  double *x      = new double[n];
  double *xlow   = new double[n];
  double *xupp   = new double[n];
  double *xmul   = new double[n];
  int    *xstate = new    int[n];

  double *F      = new double[nF];
  double *Flow   = new double[nF];
  double *Fupp   = new double[nF];
  double *Fmul   = new double[nF];
  int    *Fstate = new int[nF];

  int lenA   = 2;
  int *iAfun = new int[lenA];
  int *jAvar = new int[lenA];
  double *A  = new double[lenA];

  int lenG   = 14*nH + 1;
  int *iGfun = new int[lenG];
  int *jGvar = new int[lenG];

  int nS = 0, nInf = 0, neA = 0, neG = 0;
  int jx1, jx2, ju, ode1, ode2, Obj;
  double sInf;

  double inf = 1.0e20;

  jx1 = 0;
  jx2 = jx1 + nH + 1;
  ju  = jx2 + nH + 1;

  ode1 = 0;
  ode2 = ode1 + nH;
  Obj  = ode2 + nH;

  ObjRow = Obj;

  // Linear terms first
  iAfun[neA] = ObjRow;
  jAvar[neA] = jx1 + nH;
  A[neA]     = 1.0;
  neA        = neA + 1;

  iAfun[neA] = ObjRow;
  jAvar[neA] = jx2 + nH;