---

AmplFunctional.hh

Go to the documentation of this file.

#ifndef AMPLFUNCTIONAL
#define AMPLFUNCTIONAL

#include "Functional.hh"
#include "IotrConfiguration.hh"

#include<string>

struct ASL;
struct Option_Info;
struct keyword;

class ElementGatherSeq;

#include "SimpleVector.hh"

class AmplAdaptor : public IotrRefCount {
protected:
  ASL * asl;
  Option_Info * oi;
  std::string mSname;
  std::string mBsname;

  SimpleVectorHandle mX;
  SimpleVectorHandle mYZ;

  IotrConfigurationHandle mCommandLineOptions;

public:

  IotrVector & x() { return *mX; };
  IotrVector & yz() { return *mYZ; };

  ASL * itsAsl() { return asl; }

  char * processKeyword( keyword * kw, char value[] );
  AmplAdaptor( const char sname[], const char bsname[] );
  ~AmplAdaptor();
  virtual IotrConfiguration * initialize( char ** argv, int &ierr,
                                          keyword * kw, int nkw );
  virtual void writeSolution( char * msg, IotrVector * primal,
                              IotrVector * dual );
  void getBounds( IotrVector & rxlow,
                  IotrVector & rxupp,
                  IotrVector & rclow,
                  IotrVector & rcupp );

  void getInitialPrimal( IotrVector & primal ) const
  {
    primal.copyFrom( *mX );
  };

  void getInitialDual( IotrVector & dual ) {
    dual.copyFrom( *mYZ );
  };

};

typedef Handle<AmplAdaptor> AmplAdaptorHandle;

class AmplFunctional : public ConcreteFunctional {
protected:
  ASL * asl;
  AmplAdaptorHandle adaptor;
public:
  AmplFunctional( AmplAdaptor * adaptor );
  //specialize ConcreteFunctional::concreteValues
  virtual void concreteValues( double & obj, IotrVector & c,
                               IotrVector & x);
  //specialize ConcreteFunctional::concreteGrads
  virtual void concreteGrads( IotrVector & g, GenMatrix & J,
                              IotrVector & x ) = 0;
  //specializes ConcreteFunctional::concreteHess
  virtual void concreteHess( SymMatrix & H,
                             IotrVector & x, IotrVector & Lm) = 0;

  virtual IotrVector * newPrimalVector() const;
  virtual IotrVector * newDualVector() const;
  virtual GenMatrix  * newJacobianMatrix() const = 0;
  virtual SymMatrix  * newHessianMatrix() const = 0;

  virtual PrimalDualSystem *
  newPrimalDualSystem( ElementGather * iclow, ElementGather * icupp,
                       ElementGather * ieq,
                       ElementGather * ixlow, ElementGather * ixupp,
                       IotrMask * aXmask ) = 0;

  virtual void writeConfiguration( std::ostream & out, char indent[] ) = 0;
};

typedef Handle<AmplFunctional> AmplFunctionalHandle;

char * ampladaptor_process_keyword( Option_Info * oi,
                                       keyword * kw, char value[] );

#endif

---