linforward.cc

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.
//
//=============================================================
// Author: H. Lydia Deng
// Implementation of a objective function class
// Implementation of linear operators as objective function for optimization
//=============================================================

#include <Forward.hh>
#include <LinForward.hh>

static const char*  myNameIs =  "LinearForward";

const char* LinearForward::className() const{ return (myNameIs);}
    
LinearForward::LinearForward(Matrix<double>* pm)
: Forward()
{
	int m = pm->numOfRows();
//	int n = pm->numOfCols();

	ndata = m;
	jacob = pm;
}

LinearForward::LinearForward(Matrix<double>* pm, int verbose)
: Forward(verbose)
{
	int m = pm->numOfRows();
//	int n = pm->numOfCols();

	ndata = m;
	jacob = pm;
}

Vector<double>	LinearForward::adjointOp(const Vector<double>& v)
{
    return jacob->atdotx(v);
}

Vector<double>	LinearForward::adjointOp(const Vector<long>& v)
{
   Vector<double> u(v);
   return jacob->atdotx(u);
}

double	LinearForward::forwardOp(int i, const Vector<double>& v)
{
   Vector<double> a(jacob->rowVector(i));
   return a*v;
}

double	LinearForward::forwardOp(int i, const Vector<long>& v)
{
   Vector<double> a(jacob->rowVector(i)), u(v);
   return a*u;
}

Vector<double>	LinearForward::forwardOp(const Vector<double>& v)
{
    return jacob->adotx(v);
}

Vector<double>	LinearForward::forwardOp(const Vector<long>& v)
{
   Vector<double> u(v);
   return jacob->adotx(u);
}

List<double> LinearForward::dataList(const Model<double>& p)
{
    Vector<double> v(forwardOp(p.modParam()));
    List<double> l(ToList(v));
    return l;
}

List<double> LinearForward::dataList(const Model<long>& p)
{
   Vector<double> v(forwardOp(p.modParam()));
   List<double> l(ToList(v));
   return l;
}