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



#if GEOM_PSEUDO3D==GEOM_TRUE

#else

#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)