---

OoqpVector.h

/* OOQP                                                               *
 * Authors: E. Michael Gertz, Stephen J. Wright                       *
 * (C) 2001 University of Chicago. See Copyright Notification in OOQP */

#ifndef OOQPVECTOR_H
#define OOQPVECTOR_H

#include <string.h>
#include <iostream.h>
#include "IotrRefCount.h"


class OoqpVector : public IotrRefCount {
public:
  int n;
  int length() { return n; }

  OoqpVector( int n_ = 0 );
  virtual ~OoqpVector();

  virtual void setToZero() = 0;
  virtual void setToConstant( double c ) = 0;
  virtual void randomize( double alpha, double beta, double *ix ) = 0;
  virtual void copyFrom( OoqpVector& v ) = 0;
  virtual double infnorm() = 0;
  virtual double onenorm() = 0;

  virtual void componentMult( OoqpVector& v ) = 0;
  virtual void componentDiv ( OoqpVector& v ) = 0;
  virtual void scalarMult( double num) = 0;

  virtual void writeToStream(ostream& out) const = 0;
  virtual void writefToStream( ostream& out,
                               const char format[] ) const = 0;

  virtual void scale( double alpha ) = 0;

  virtual void axpy  ( double alpha, OoqpVector& x ) = 0;
  virtual void axzpy ( double alpha, OoqpVector& x, OoqpVector& z ) = 0;
  virtual void axdzpy( double alpha, OoqpVector& x, OoqpVector& z ) = 0;

  virtual void addConstant( double c ) = 0;

  virtual void gondzioProjection( double rmin, double rmax ) = 0;

  virtual void min( double& m, int& index ) = 0;
  virtual double dotProductWith( OoqpVector& v ) = 0;

  virtual double shiftedDotProductWith( double alpha, OoqpVector& mystep,
                                        OoqpVector& yvec,
                                        double beta,  OoqpVector& ystep ) = 0;
  virtual void negate() = 0;

  virtual void invert() = 0;

  virtual int allPositive() = 0;
  
  virtual int numberOfNonzeros() = 0;

  virtual int matchesNonZeroPattern( OoqpVector& select ) = 0;

  virtual void selectNonZeros( OoqpVector& select ) = 0;
  virtual void addSomeConstants( double c, OoqpVector& select ) = 0;

  virtual void writefSomeToStream( ostream& out,
                                   const char format[],
                                   OoqpVector& select ) const = 0;
  virtual void axdzpy( double alpha, OoqpVector& x,
                       OoqpVector& z, OoqpVector& select ) = 0;

  virtual int somePositive( OoqpVector& select ) = 0;
  virtual void divideSome( OoqpVector& div, OoqpVector& select ) = 0;

  virtual int isKindOf( int kind ) = 0;

  virtual double stepbound(OoqpVector & v, double maxStep ) = 0;

  virtual double findBlocking(OoqpVector & wstep_vec, 
                              OoqpVector & u_vec, 
                              OoqpVector & ustep_vec, 
                              double maxStep,
                              double *w_elt, 
                              double *wstep_elt,
                              double *u_elt, 
                              double *ustep_elt,
                              int& first_or_second) = 0;

  virtual void copyIntoArray( double v[] ) const = 0;
  virtual void copyFromArray( double v[] ) = 0;
  virtual void copyFromArray( char v[] ) = 0;
};


enum { kSimpleVector = 0, kPetscVector };

#endif

---