---

Functional.hh

Go to the documentation of this file.

#ifndef FUNCTIONAL
#define FUNCTIONAL

#include "IotrRefCount.hh"
#include "IotrHandle.hh"
#include "PrimalDualSystem.hh"
#include <iostream>

class EvaluationContext;
class ConcreteEvaluationContext;
class IotrVector;
class GenMatrix;
class SymMatrix;
class IotrMask;

class Functional : public IotrRefCount {
public:
  virtual EvaluationContext * evaluate(IotrVector * x, IotrVector * Lm ) = 0;

  virtual IotrVector * newPrimalVector()   const = 0;
  virtual IotrVector * newDualVector()     const = 0;
  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;
};
typedef Handle<Functional> FunctionalHandle;


class ConcreteFunctional : public Functional {
friend class ConcreteEvaluationContext;
protected:
  virtual void concreteValues( double& obj, IotrVector & c,
                               IotrVector & x )  = 0;
  virtual void concreteGrads( IotrVector & g, GenMatrix & J,
                              IotrVector & x ) = 0;
  virtual void concreteHess( SymMatrix & H,
                             IotrVector & x, IotrVector & Lm) = 0;
public:
  // specializes Functional::values
  //virtual void values( double& obj, IotrVector & c, IotrVector & x );
  virtual EvaluationContext * evaluate(IotrVector * x, IotrVector * Lm );
};
typedef Handle<ConcreteFunctional> ConcreteFunctionalHandle;


class EvaluationContext : public IotrRefCount {
public:
  virtual void values( double& obj, IotrVector & c) = 0;
  virtual void grads( IotrVector & g, GenMatrix & J ) = 0;
  virtual void hessian( SymMatrix & H) = 0;
};

typedef Handle<EvaluationContext> EvaluationContextHandle;

class ConcreteEvaluationContext : public EvaluationContext {
protected:
  IotrVectorHandle mX;
  IotrVectorHandle mLm;
  ConcreteFunctionalHandle mFunc;
public:
  ConcreteEvaluationContext(IotrVector * x, IotrVector * Lm,
                           ConcreteFunctional * func );
  //specializes EvaluationContext::values
  virtual void values( double & obj, IotrVector & c);
  //specializes EvaluationContext::grads
  virtual void grads( IotrVector & g, GenMatrix & J );
  //specializes EvalutationContext::hessian
  virtual void hessian( SymMatrix & H);
};
typedef Handle<ConcreteEvaluationContext> ConcreteEvaluationContextHandle;

#endif

---