#line 1 "/user/cvsspst/ees1cg/RAVL/RAVL-0.7/PatternRec/Modeling/Basic/FuncQuadratic.hh"
// This file is part of RAVL, Recognition And Vision Library 
// Copyright (C) 2002, University of Surrey
// This code may be redistributed under the terms of the GNU Lesser
// General Public License (LGPL). See the lgpl.licence file for details or
// see http://www.gnu.org/copyleft/lesser.html
// file-header-ends-here
#ifndef RAVL_FUNCQUADRATIC_HEADER
#define RAVL_FUNCQUADRATIC_HEADER 1
//! rcsid="$Id: FuncQuadratic.hh,v 1.5 2002/08/08 16:03:22 craftit Exp $"
//! lib=RavlPatternRec
//! author="Robert Crida and Charles Galambos"
//! docentry="Ravl.Pattern Recognition.Numerical Modeling"
//! file="Ravl/PatternRec/Modeling/Basic/FuncQuadratic.hh"

#include "Ravl/PatternRec/FuncLinearCoeff.hh"
#include "Ravl/Matrix.hh"

namespace RavlN {
  
  //! userlevel=Develop
  //: Quadratic function.
  
  class FuncQuadraticBodyC 
    : public FuncLinearCoeffBodyC
  {
  public:
    FuncQuadraticBodyC(int inSize,int outSize)
      : FuncLinearCoeffBodyC(inSize,outSize)
    { noCoeffs = NumberCoeffs(inSize); }
    //: Constructor.
    
    FuncQuadraticBodyC(istream &strm);
    //: Load from stream.
    
    FuncQuadraticBodyC(BinIStreamC &strm);
    //: Load from binary stream.
    
    virtual bool Save (ostream &out) const;
    //: Writes object to stream, can be loaded using constructor
    
    virtual bool Save (BinOStreamC &out) const;
    //: Writes object to stream, can be loaded using constructor
    
    virtual VectorC MakeInput (const VectorC &X) const;
    //: Expand vector to linear coefficients.
    
    virtual VectorC MakeJacobianInput (const VectorC &X, IntT i) const;
    //: Expand vector to linear coefficients.
    
    virtual UIntT NumberCoeffs(UIntT inputSize) const;
    //: Calculate the number of coefficents for a given input size.
    
  protected:
    
    UIntT noCoeffs; // Number of coefficents.
  };
  
  //! userlevel=Normal
  //: Quadratic function.
  //
  // It provides a concrete class implementation for functions of the
  // form Y=f(X) where each element of Y is an 2nd polynomial of the elements
  // of X including all the cross parameter terms.
  
  class FuncQuadraticC
    : public FuncLinearCoeffC
  {
  public:    
    FuncQuadraticC(int inSize,int outSize)
      : FuncLinearCoeffC(*new FuncQuadraticBodyC(inSize,outSize))
    {}
    //: Constructor.
    
    FuncQuadraticC(istream &strm);
    //: Load from stream.
    
    FuncQuadraticC(BinIStreamC &strm);
    //: Load from binary stream.
    
  protected:
    FuncQuadraticC(FuncQuadraticBodyC &bod)
      : FuncLinearCoeffC(bod)
    {}
    //: Body constructor.

    FuncQuadraticC(FuncQuadraticBodyC *bod)
      : FuncLinearCoeffC(bod)
    {}
    //: Body ptr constructor.
    
    FuncQuadraticBodyC &Body()
    { return static_cast<FuncQuadraticBodyC &>(FunctionC::Body()); }
    //: Access body.
    
    const FuncQuadraticBodyC &Body() const
    { return static_cast<const FuncQuadraticBodyC &>(FunctionC::Body()); }
    //: Access body.
    
  public:
    
  };
  
  inline istream &operator>>(istream &strm,FuncQuadraticC &obj) {
    obj = FuncQuadraticC(strm);
    return strm;
  }
  //: Load from a stream.
  // Uses virtual constructor.
  
  inline ostream &operator<<(ostream &out,const FuncQuadraticC &obj) {
    obj.Save(out);
    return out;
  }
  //: Save to a stream.
  // Uses virtual constructor.
  
  inline BinIStreamC &operator>>(BinIStreamC &strm,FuncQuadraticC &obj) {
    obj = FuncQuadraticC(strm);
    return strm;
  }
  //: Load from a binary stream.
  // Uses virtual constructor.
  
  inline BinOStreamC &operator<<(BinOStreamC &out,const FuncQuadraticC &obj) {
    obj.Save(out);
    return out;
  }
  //: Save to a stream.
  // Uses virtual constructor.
  
}

#endif