/* $Id: ClpQuadraticObjective.hpp 1525 2010-02-26 17:27:59Z mjs $ */
// Copyright (C) 2003, International Business Machines
// Corporation and others.  All Rights Reserved.
#ifndef ClpQuadraticObjective_H
#define ClpQuadraticObjective_H

#include "ClpObjective.hpp"
#include "CoinPackedMatrix.hpp"

//#############################################################################

/** Quadratic Objective Class

*/

class ClpQuadraticObjective : public ClpObjective {

public:

     ///@name Stuff
     //@{

     /** Returns gradient.  If Quadratic then solution may be NULL,
         also returns an offset (to be added to current one)
         If refresh is false then uses last solution
         Uses model for scaling
         includeLinear 0 - no, 1 as is, 2 as feasible
     */
     virtual double * gradient(const ClpSimplex * model,
                               const double * solution, double & offset, bool refresh,
                               int includeLinear = 2);
     /// Resize objective
     /** Returns reduced gradient.Returns an offset (to be added to current one).
     */
     virtual double reducedGradient(ClpSimplex * model, double * region,
                                    bool useFeasibleCosts);
     /** Returns step length which gives minimum of objective for
         solution + theta * change vector up to maximum theta.

         arrays are numberColumns+numberRows
         Also sets current objective, predicted and at maximumTheta
     */
     virtual double stepLength(ClpSimplex * model,
                               const double * solution,
                               const double * change,
                               double maximumTheta,
                               double & currentObj,
                               double & predictedObj,
                               double & thetaObj);
     /// Return objective value (without any ClpModel offset) (model may be NULL)
     virtual double objectiveValue(const ClpSimplex * model, const double * solution) const ;
     virtual void resize(int newNumberColumns) ;
     /// Delete columns in  objective
     virtual void deleteSome(int numberToDelete, const int * which) ;
     /// Scale objective
     virtual void reallyScale(const double * columnScale) ;
     /** Given a zeroed array sets nonlinear columns to 1.
         Returns number of nonlinear columns
      */
     virtual int markNonlinear(char * which);

     //@}


     ///@name Constructors and destructors
     //@{
     /// Default Constructor
     ClpQuadraticObjective();

     /// Constructor from objective
     ClpQuadraticObjective(const double * linearObjective, int numberColumns,
                           const CoinBigIndex * start,
                           const int * column, const double * element,
                           int numberExtendedColumns_ = -1);

     /** Copy constructor .
         If type is -1 then make sure half symmetric,
         if +1 then make sure full
     */
     ClpQuadraticObjective(const ClpQuadraticObjective & rhs, int type = 0);
     /** Subset constructor.  Duplicates are allowed
         and order is as given.
     */
     ClpQuadraticObjective (const ClpQuadraticObjective &rhs, int numberColumns,
                            const int * whichColumns) ;

     /// Assignment operator
     ClpQuadraticObjective & operator=(const ClpQuadraticObjective& rhs);

     /// Destructor
     virtual ~ClpQuadraticObjective ();

     /// Clone
     virtual ClpObjective * clone() const;
     /** Subset clone.  Duplicates are allowed
         and order is as given.
     */
     virtual ClpObjective * subsetClone (int numberColumns,
                                         const int * whichColumns) const;

     /** Load up quadratic objective.  This is stored as a CoinPackedMatrix */
     void loadQuadraticObjective(const int numberColumns,
                                 const CoinBigIndex * start,
                                 const int * column, const double * element,
                                 int numberExtendedColumns = -1);
     void loadQuadraticObjective (  const CoinPackedMatrix& matrix);
     /// Get rid of quadratic objective
     void deleteQuadraticObjective();
     //@}
     ///@name Gets and sets
     //@{
     /// Quadratic objective
     inline CoinPackedMatrix * quadraticObjective() const     {
          return quadraticObjective_;
     }
     /// Linear objective
     inline double * linearObjective() const     {
          return objective_;
     }
     /// Length of linear objective which could be bigger
     inline int numberExtendedColumns() const {
          return numberExtendedColumns_;
     }
     /// Number of columns in quadratic objective
     inline int numberColumns() const {
          return numberColumns_;
     }
     /// If a full or half matrix
     inline bool fullMatrix() const {
          return fullMatrix_;
     }
     //@}

     //---------------------------------------------------------------------------

private:
     ///@name Private member data
     /// Quadratic objective
     CoinPackedMatrix * quadraticObjective_;
     /// Objective
     double * objective_;
     /// Gradient
     double * gradient_;
     /// Useful to have number of columns about
     int numberColumns_;
     /// Also length of linear objective which could be bigger
     int numberExtendedColumns_;
     /// True if full symmetric matrix, false if half
     bool fullMatrix_;
     //@}
};

#endif
