swh:1:snp:af87cd67498ef4fe47c76ed3e7caffe5b61facaf
Tip revision: c86ba8382b7954f4a5b3840ec65a9726ba981fd7 authored by Unknown Author on 18 October 2006, 08:53:49 UTC
This commit was manufactured by cvs2svn to create tag 'v5-13-04a'.
This commit was manufactured by cvs2svn to create tag 'v5-13-04a'.
Tip revision: c86ba83
TQpDataDens.h
// @(#)root/quadp:$Name: $:$Id: TQpDataDens.h,v 1.6 2006/06/03 05:42:18 brun Exp $
// Author: Eddy Offermann May 2004
/*************************************************************************
* Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
* All rights reserved. *
* *
* For the licensing terms see $ROOTSYS/LICENSE. *
* For the list of contributors see $ROOTSYS/README/CREDITS. *
*************************************************************************/
/*************************************************************************
* Parts of this file are copied from the OOQP distribution and *
* are subject to the following license: *
* *
* COPYRIGHT 2001 UNIVERSITY OF CHICAGO *
* *
* The copyright holder hereby grants you royalty-free rights to use, *
* reproduce, prepare derivative works, and to redistribute this software*
* to others, provided that any changes are clearly documented. This *
* software was authored by: *
* *
* E. MICHAEL GERTZ gertz@mcs.anl.gov *
* Mathematics and Computer Science Division *
* Argonne National Laboratory *
* 9700 S. Cass Avenue *
* Argonne, IL 60439-4844 *
* *
* STEPHEN J. WRIGHT swright@cs.wisc.edu *
* Computer Sciences Department *
* University of Wisconsin *
* 1210 West Dayton Street *
* Madison, WI 53706 FAX: (608)262-9777 *
* *
* Any questions or comments may be directed to one of the authors. *
* *
* ARGONNE NATIONAL LABORATORY (ANL), WITH FACILITIES IN THE STATES OF *
* ILLINOIS AND IDAHO, IS OWNED BY THE UNITED STATES GOVERNMENT, AND *
* OPERATED BY THE UNIVERSITY OF CHICAGO UNDER PROVISION OF A CONTRACT *
* WITH THE DEPARTMENT OF ENERGY. *
*************************************************************************/
#ifndef ROOT_TQpDataDens
#define ROOT_TQpDataDens
#ifndef ROOT_TROOT
#include "TROOT.h"
#endif
#ifndef ROOT_TClass
#include "TClass.h"
#endif
#ifndef ROOT_TError
#include "TError.h"
#endif
#ifndef ROOT_TMath
#include "TMath.h"
#endif
#ifndef ROOT_TQpDataBase
#include "TQpDataBase.h"
#endif
#ifndef ROOT_TQpVar
#include "TQpVar.h"
#endif
#ifndef ROOT_TMatrixD
#include "TMatrixD.h"
#endif
#ifndef ROOT_TMatrixDSym
#include "TMatrixDSym.h"
#endif
//////////////////////////////////////////////////////////////////////////
// //
// TQpDataDens //
// //
// Data for the dense QP formulation //
// //
//////////////////////////////////////////////////////////////////////////
class TQpDataDens : public TQpDataBase
{
protected:
// these variables will be "Used" not copied
TMatrixDSym fQ; // Quadratic part of Objective function
TMatrixD fA; // Equality constraints
TMatrixD fC; // Inequality constraints
public:
TQpDataDens() {}
// data objects of the specified dimensions
TQpDataDens(Int_t nx,Int_t my,Int_t mz);
// sets up pointers to the data objects that are passed as arguments
TQpDataDens(TVectorD &c,TMatrixDSym &Q,TVectorD &xlow,TVectorD &ixlow,TVectorD &xupp,
TVectorD &ixupp,TMatrixD &A,TVectorD &bA,TMatrixD &C,TVectorD &clow,
TVectorD &iclow,TVectorD &cupp,TVectorD &icupp);
TQpDataDens(const TQpDataDens &another);
virtual ~TQpDataDens() {}
virtual void PutQIntoAt(TMatrixDBase &M,Int_t row,Int_t col);
// insert the Hessian Q into the matrix M for the fundamental
// linear system, where M is stored as a TMatrixD
virtual void PutAIntoAt(TMatrixDBase &M,Int_t row,Int_t col);
// insert the constraint matrix A into the matrix M for the
// fundamental linear system, where M is stored as a TMatrixD
virtual void PutCIntoAt(TMatrixDBase &M,Int_t row,Int_t col);
// insert the constraint matrix C into the matrix M for the
// fundamental linear system, where M is stored as a TMatrixD
virtual void Qmult (Double_t beta,TVectorD& y,Double_t alpha,const TVectorD& x);
// y = beta * y + alpha * Q * x
virtual void Amult (Double_t beta,TVectorD& y,Double_t alpha,const TVectorD& x);
// y = beta * y + alpha * A * x
virtual void Cmult (Double_t beta,TVectorD& y,Double_t alpha,const TVectorD& x);
// y = beta * y + alpha * C * x
virtual void ATransmult(Double_t beta,TVectorD& y,Double_t alpha,const TVectorD& x);
// y = beta * y + alpha * A^T * x
virtual void CTransmult(Double_t beta,TVectorD& y,Double_t alpha,const TVectorD& x);
// y = beta * y + alpha * C^T * x
virtual void GetDiagonalOfQ(TVectorD &dQ); // extract the diagonal of Q and put it in the vector dQ
virtual Double_t DataNorm();
virtual void DataRandom(TVectorD &x,TVectorD &y,TVectorD &z,TVectorD &s);
// Create a random problem (x,y,z,s)
// the solution to the random problem
virtual void Print(Option_t *opt="") const;
virtual Double_t ObjectiveValue(TQpVar *vars);
TQpDataDens &operator= (const TQpDataDens &source);
ClassDef(TQpDataDens,1) // Qp Data class for Dens formulation
};
#endif