---

HuberData.h

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

#ifndef HuberDATA_H
#define HuberDATA_H

#include "Data.h"
#include "DenseGenMatrixHandle.h"
#include "DenseSymMatrixHandle.h"
#include "SimpleVectorHandle.h"

class HuberData : public Data
{
 public:

  int m, n;

  HuberData(int nobservations_in, int npredictors_in,  double cutoff,
            double * X = 0, double * y = 0);

  ~HuberData();

  DenseGenMatrixHandle Xt;

  SimpleVectorHandle Y;

  double cutoff;

  int nobservations;

  int npredictors;

  static HuberData * textInput(char filename[], double cutoff, int &iErr);

  virtual void XtMult( double beta, SimpleVector& y,
                      double alpha, SimpleVector& x );

  virtual void XtTransMult( double beta, SimpleVector& y,
                           double alpha, SimpleVector& x );

  virtual double datanorm();
  virtual void datarandom();
  virtual void print() {};
};


enum { huberinputok = 0,
       huberfileopenerror = -1, // couldn't find the designated input file 
       huberfileinputerror = -2 // not enough input, or error on input
};

#endif

---