mlVector4.h

Go to the documentation of this file.

/*************************************************************************************
**
** Copyright 2007, MeVis Medical Solutions AG
**
** The user may use this file in accordance with the license agreement provided with
** the Software or, alternatively, in accordance with the terms contained in a
** written agreement between the user and MeVis Medical Solutions AG.
**
** For further information use the contact form at https://www.mevislab.de/contact
**
**************************************************************************************/
 
#ifndef ML_VECTOR4_H
#define ML_VECTOR4_H

 
// Include system independent file and project settings.
#include "mlLinearAlgebraSystem.h"
#include "mlLinearAlgebraDefs.h"
 
#include "mlFloatingPointVector.h"
#include "mlVector2.h"
#include "mlVector3.h"
#include "mlMatrix3.h"
#include "mlMatrix4.h"
 
#include <mlErrorOutput.h>
 
// All declarations of this header will be in the ML_LA_NAMESPACE namespace.
ML_LA_START_NAMESPACE
 
//--------------------------------------------------------------------
// This is necessary, because we do not know whether vector or
// matrix header is included first and we cannot move template
// code into the C++ file.
//--------------------------------------------------------------------
template <class DT> class Tvec2;
template <class DT> class Tvec3;
template <class DT> class Tmat3;
template <class DT> class Tmat4;
 
//--------------------------------------------------------------------
//--------------------------------------------------------------------
template <class DT>
class Tvec4 : public FloatingPointVector<DT,4>
{
public:

  typedef FloatingPointVector<DT,4> Superclass;

  typedef DT ComponentType;
 
  //--------------------------------------------------------------------
  //--------------------------------------------------------------------
  inline explicit Tvec4(const DT value=0) : Superclass(value)
  {
  }

  inline Tvec4(const Superclass &v): Superclass(v)
  {
  }

  inline Tvec4(const DT x, const DT y, const DT z, const DT w)
  {
    Superclass::_buffer[0] = x;
    Superclass::_buffer[1] = y;
    Superclass::_buffer[2] = z;
    Superclass::_buffer[3] = w;
  }

  inline Tvec4(const Tvec2<DT>& v, const DT z, const DT w)
  {
    Superclass::_buffer[0] = v[0];
    Superclass::_buffer[1] = v[1];
    Superclass::_buffer[2] = z;
    Superclass::_buffer[3] = w;
  }

  inline Tvec4(const Tvec3<DT>& v, const DT w)
  {
    Superclass::_buffer[0] = v[0];
    Superclass::_buffer[1] = v[1];
    Superclass::_buffer[2] = v[2];
    Superclass::_buffer[3] = w;
  }

  inline void assign(const DT x, const DT y, const DT z, const DT c)
  {
    Superclass::_buffer[0] = x;
    Superclass::_buffer[1] = y;
    Superclass::_buffer[2] = z;
    Superclass::_buffer[3] = c;
  }

  inline Tvec4<DT> divideByLastComp() const
  {
    const DT lastComp = Superclass::_buffer[3];
 
    ML_CHECK_FLOAT_THROW(lastComp);
    const DT div = static_cast<DT>(1) / lastComp;
    return Tvec4<DT>(Superclass::_buffer[0] * div,
                     Superclass::_buffer[1] * div,
                     Superclass::_buffer[2] * div, 1);
  }

};

#define _ML_VEC4_RC(i) a[i][0]*v[0] + a[i][1]*v[1] + a[i][2]*v[2] + a[i][3]*v[3]
//-----------------------------------------------------------------------------------
//-----------------------------------------------------------------------------------
template <class DT>
inline Tvec4<DT> operator*(const Tmat4<DT>& a, const Tvec4<DT>& v)
{
  return Tvec4<DT>(_ML_VEC4_RC(0), _ML_VEC4_RC(1), _ML_VEC4_RC(2), _ML_VEC4_RC(3));
}
#undef _ML_VEC4_RC
 
//-----------------------------------------------------------------------------------
//-----------------------------------------------------------------------------------
template <class DT>
inline Tvec4<DT> operator*(const Tvec4<DT>& v, const Tmat4<DT>& a)
{
  return a.transpose() * v;
}
 
 
//-----------------------------------------------------------------------------------
//-----------------------------------------------------------------------------------
typedef Tvec4<MLfloat>   Vector4f;
typedef Tvec4<MLdouble>  Vector4d;
typedef Tvec4<MLldouble> Vector4ld;
typedef Tvec4<MLdouble>  Vector4;
 
ML_LA_END_NAMESPACE
 
#endif //of __mlVector4_H