OOPSE/libmdtools/Constraint.hpp

#ifndef _CONSTRAINT_H_
#define _CONSTRAINT_H_

#include <iostream>
#include <vector>

#include "SymMatrix.hpp"

#define ERROR_CONSTRAINT 10

using namespace std;

typedef enum  {simpleBound = 1, linearEqu = 2, linearInequ = 4,
                                     nonlinearEqu = 8, nonlinearInequ =16} ConsType;

typedef enum {btUpper, btLower, btEqu} BoundType;

/**
 * Abstract class of constraint for nonlinear optimization problem
 */ 
class ConstraintBase{
  public:
    ConstraintBase();
    ConstraintBase(int dim);

    virtual void setDim(int dim);
    bool isDimSet() {return init_ndim;}
    
    int getConsType() { return consType;}
    virtual double calcResidual(vector<double>& x) = 0;
    virtual vector<double> calcConsGrad(vector<double>& x) = 0;
    virtual SymMatrix calcConsHessian(vector<double>& x) = 0;

  protected:

    bool init_ndim;
    int ndim;
    
    double bound;
    BoundType boundType;
    ConsType consType;
  
};

#endif
Revision:	1023
Committed:	Wed Feb 4 22:26:00 2004 UTC (21 years, 9 months ago) by tim
File size:	1006 byte(s)
Log Message:	Fix a bunch of bugs :-) Single version of conjugate gradient with golden search linesearch pass a couple of functions test. Brent's algorithm is still broken
#	User	Rev	Content
1	tim	987	#ifndef _CONSTRAINT_H_
2			#define _CONSTRAINT_H_
3
4	tim	1010	#include <iostream>
5	tim	987	#include <vector>
6
7	tim	1010	#include "SymMatrix.hpp"
8	tim	987
9			#define ERROR_CONSTRAINT 10
10
11			using namespace std;
12
13	tim	1010	typedef enum {simpleBound = 1, linearEqu = 2, linearInequ = 4,
14			nonlinearEqu = 8, nonlinearInequ =16} ConsType;
15	tim	987
16	tim	1010	typedef enum {btUpper, btLower, btEqu} BoundType;
17	tim	987
18	tim	995	/**
19			* Abstract class of constraint for nonlinear optimization problem
20			*/
21	tim	987	class ConstraintBase{
22			public:
23			ConstraintBase();
24			ConstraintBase(int dim);
25
26			virtual void setDim(int dim);
27	tim	1023	bool isDimSet() {return init_ndim;}
28	tim	987
29	tim	1010	int getConsType() { return consType;}
30	tim	987	virtual double calcResidual(vector<double>& x) = 0;
31			virtual vector<double> calcConsGrad(vector<double>& x) = 0;
32			virtual SymMatrix calcConsHessian(vector<double>& x) = 0;
33
34			protected:
35
36			bool init_ndim;
37			int ndim;
38
39			double bound;
40			BoundType boundType;
41			ConsType consType;
42
43			};
44
45			#endif