#line 1 "/user/cvsspst/ees1cg/RAVL/RAVL-0.7/Math/LinearAlgebra/FixedSize/Vector4d.hh"
// This file is part of RAVL, Recognition And Vision Library 
// Copyright (C) 2001, 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_VECTOR4D_HEADER
#define RAVL_VECTOR4D_HEADER 1
/////////////////////////////////////////////////////
//! rcsid="$Id: Vector4d.hh,v 1.6 2002/03/07 13:50:15 craftit Exp $"
//! file="Ravl/Math/LinearAlgebra/FixedSize/Vector4d.hh"
//! lib=RavlMath
//! userlevel=Normal
//! author="Charles Galambos"
//! date="24/01/2001"
//! docentry="Ravl.Math.Geometry.4D"

#include "Ravl/FVector.hh"

namespace RavlN {
  
  //! userlevel=Normal
  //: 4D Vector
  
  class Vector4dC
    : public FVectorC<4>
  {
  public:
    Vector4dC()
      {}
    //: Default constructor.
    
    Vector4dC(const TFVectorC<RealT,4> &base)
      : FVectorC<4>(base)
      {}
    //: Base onstructor.
    
    Vector4dC(RealT v1,RealT v2,RealT v3,RealT v4) { 
      data[0] = v1;
      data[1] = v2;
      data[2] = v3;
      data[3] = v4;
    }
    //: Constructor.
    
    RealT &X()
      { return data[0]; }
    //: First component of vector.

    RealT X() const
      { return data[0]; }
    //: First component of vector.

    RealT &Y()
      { return data[1]; }
    //: Second component of vector.

    RealT Y() const
      { return data[1]; }
    //: Second component of vector.
    
    RealT &Z()
      { return data[2]; }
    //: Third component of vector.
    
    RealT Z() const
      { return data[2]; }
    //: Third component of vector.
    
    RealT &W()
      { return data[3]; }
    //: Fourth component of vector.
    
    RealT W() const
      { return data[3]; }
    //: Fourth component of vector.
  };
  
  inline RealT TFVectorC<RealT,4>::Dot(const TFVectorC<RealT,4> &oth) const 
  { return data[0] * oth.data[0] + data[1] * oth.data[1] + data[2] * oth.data[2] + data[3] * oth.data[3]; }
  //: Calculate the dot product of this and 'oth' vector.

}

#endif