art.hh

COOOL：CWP面向对象最优化库（CWP Object Oriented Optimization Library） COOOL是C++类的一个集合

//============================================================
// COOOL           version 1.1           ---     Nov,  1995
//   Center for Wave Phenomena, Colorado School of Mines
//============================================================
//
//   This code is part of a preliminary release of COOOL (CWP
// Object-Oriented Optimization Library) and associated class 
// libraries. 
//
// The COOOL library is a free software. You can do anything you want
// with it including make a fortune.  However, neither the authors,
// the Center for Wave Phenomena, nor anyone else you can think of
// makes any guarantees about anything in this package or any aspect
// of its functionality.
//
// Since you've got the source code you can also modify the
// library to suit your own purposes. We would appreciate it 
// if the headers that identify the authors are kept in the 
// source code.
//
#ifndef ART_HH
#define ART_HH

#include "QuaOptima.hh"

//=============================
// author:  H. Lydia Deng, 08/08/95
// Definition of the Algebraic Reconstruction Technique (ART)
//=============================
//@Man:
//@Memo: Algebraic Reconstruction Technique
/*@Doc:


*/

class ART : public QuadraticOptima {

  protected:
    int	m, n;
    Model<double> upDating(int, const Model<double>& );
    Model<long> upDating(int, const Model<long>& );

 public:				//Constructors and Destructors
    //@ManMemo: a constructor
    ART(/// dimension of model space
			int n, 
			///pointer to the forward operator (matrix)
			LinearForward* p, 
			/// pointer to the observed data vector
			Vector<double>* data, 
			/// maximum number of iterations in solving linear system
			int it, 
			/// the maxmimum tolerable error
			double tol, 
			///vebose or quiet
			int verb);
    //@ManMemo: a constructor
    ART(int n, LinearForward* p, Vector<double>* data, int it, double tol);
    ~ART();

        //@ManMemo: ART search starting from m0, and returns an optimum model
    Model<double> 	optimizer(Model<double>& m0);
        //@ManMemo: ART search starting from m0, and returns an optimum model
    Model<long> 	optimizer(Model<long>& m0);
};

#endif