Vector2.h

// Copyright (C) Geom Software e.U, Bernhard Kornberger, Graz/Austria
//
// This file is part of the Fade2D library. The student license is free
// of charge and covers personal non-commercial research. Licensees
// holding a commercial license may use this file in accordance with
// the Commercial License Agreement.
//
// This software is provided AS IS with NO WARRANTY OF ANY KIND,
// INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS
// FOR A PARTICULAR PURPOSE.
//
// Please contact the author if any conditions of this licensing are
// not clear to you.
//
// Author: Bernhard Kornberger, bkorn (at) geom.at
// http://www.geom.at
#pragma once
#include "common.h"
 
#if GEOM_PSEUDO3D==GEOM_TRUE
    namespace GEOM_FADE25D {
#elif GEOM_PSEUDO3D==GEOM_FALSE
    namespace GEOM_FADE2D {
#else
    #error GEOM_PSEUDO3D is not defined
#endif
 
class CLASS_DECLSPEC Vector2;
 
class Vector2
{
public:
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    Vector2(const double x_,const double y_,const double z_);
 
    Vector2();
 
    Vector2(const Vector2& v_);
 
#else
 
    Vector2(const double x_,const double y_);
    Vector2();
    Vector2(const Vector2& v_);
#endif
    Vector2& operator=(const Vector2& other);
 
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
#endif
 
    Vector2 orthogonalVector() const;
 
    bool isDegenerate() const;
 
 
    double x() const;
 
    double y() const;
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    double z() const;
#endif
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    void set(const double x_,const double y_,const double z_);
 
#else
 
    void set(const double x_,const double y_);
#endif
 
    double sqLength() const;
 
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    int getMaxIndex() const;
#endif
 
    double length() const;
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    double operator*(const Vector2& other) const;
#else
 
    double operator*(const Vector2& other) const;
#endif
 
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    Vector2 operator*(double val) const;
#else
 
    Vector2 operator*(double val) const;
#endif
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
    Vector2 operator/(double val) const;
#else
 
    Vector2 operator/(double val) const;
#endif
 
 
 
protected:
    double valX;
    double valY;
#if GEOM_PSEUDO3D==GEOM_TRUE
    double valZ;
#endif
};
 
 
 
 
// Free functions
 
CLASS_DECLSPEC
inline std::ostream &operator<<(std::ostream &stream, const Vector2& vec)
{
#if GEOM_PSEUDO3D==GEOM_TRUE
    stream << "Vector2: "<<vec.x()<<", "<<vec.y()<<", "<<vec.z();
#else
    stream << "Vector2: "<<vec.x()<<", "<<vec.y();
#endif
    return stream;
}
 
#if GEOM_PSEUDO3D==GEOM_TRUE
 
CLASS_DECLSPEC
inline Vector2 crossProduct(const Vector2& vec0,const Vector2& vec1)
{
    double x=vec0.y() * vec1.z() - vec0.z() * vec1.y();
    double y=vec0.z() * vec1.x() - vec0.x() * vec1.z();
    double z=vec0.x() * vec1.y() - vec0.y() * vec1.x();
    return Vector2(x,y,z);
}
#endif
 
 
 
 
CLASS_DECLSPEC
inline Vector2 normalize(const Vector2& other)
{
    double len(other.length());
#if GEOM_PSEUDO3D==GEOM_TRUE
    if(len>0)
    {
        return Vector2(other.x()/len,other.y()/len,other.z()/len);
    }
    else
    {
 
#ifndef NDEBUG
        std::cout<<"\n\n** warning: normalize(const Vector2& other), Null length vector!"<<std::endl;// COUTOK
        //std::cout<<"vec="<<other<<std::endl;
        //std::cout<<"** debug mode -> deliberate segfault, now:"<<std::endl;
        //int* i(NULL);
        //++*i;
#endif
 
        return Vector2(0,0,0);
    }
#else
    if(len>0)
    {
 
        return Vector2(other.x()/len,other.y()/len);
    }
    else
    {
        std::cout<<"warning: normalize(const Vector2& other), Null length vector!"<<std::endl;// COUTOK
        return Vector2(0,0);
    }
#endif
}
 
CLASS_DECLSPEC
inline Vector2 operator-(const Vector2& in)
{
#if GEOM_PSEUDO3D==GEOM_TRUE
    return Vector2(-in.x(),-in.y(),-in.z());
#else
    return Vector2(-in.x(),-in.y());
#endif
}
 
CLASS_DECLSPEC
inline Vector2 operator*(double d,const Vector2& vec)
{
#if GEOM_PSEUDO3D==GEOM_TRUE
    return Vector2(d*vec.x(),d*vec.y(),d*vec.z());
#else
    return Vector2(d*vec.x(),d*vec.y());
#endif
}
 
CLASS_DECLSPEC
inline Vector2 operator+(const Vector2& vec0,const Vector2& vec1)
{
#if GEOM_PSEUDO3D==GEOM_TRUE
    return Vector2(vec0.x()+vec1.x(), vec0.y()+vec1.y() , vec0.z()+vec1.z());
#else
    return Vector2(vec0.x()+vec1.x(), vec0.y()+vec1.y() );
#endif
}
 
CLASS_DECLSPEC
inline Vector2 operator-(const Vector2& vec0,const Vector2& vec1)
{
#if GEOM_PSEUDO3D==GEOM_TRUE
    return Vector2(vec0.x()-vec1.x(), vec0.y()-vec1.y() , vec0.z()-vec1.z());
#else
    return Vector2(vec0.x()-vec1.x(), vec0.y()-vec1.y() );
#endif
}
 
 
} // (namespace)