a_point2.cxx

Go to the documentation of this file.

/*
Copyright 2000-19 Pierre SMARS (smars@yuntech.edu.tw)
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 2 of the License, or
(at your option) any later version.
 
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "a_point2.h"
#include "a_point.h"
#include <math.h>
#include <sstream>
 
//---------------------------------------------------------------------------
const std::string a_point2::help()
{
        std::ostringstream o;
        return o.str();
}
//---------------------------------------------------------------------------
a_point2& a_point2::operator-()
{
        x_ = -x_;
        y_ = -y_;
        return *this;
}
//---------------------------------------------------------------------------
a_point2& a_point2::operator+=(const a_point2& p)
{
        x_ += p.x();
        y_ += p.y();
        return *this;
}
//---------------------------------------------------------------------------
a_point2& a_point2::operator-=(const a_point2& p)
{
        x_ -= p.x();
        y_ -= p.y();
        return *this;
}
//---------------------------------------------------------------------------
double a_point2::operator*=(const a_point2& p)
{
        return x_*p.x()+y_*p.y();
}
//---------------------------------------------------------------------------
bool a_point2::operator==(const a_point2& p)
{
        return ((x_ == p.x())&&(y_ == p.y()));
}
//---------------------------------------------------------------------------
a_point2 & a_point2::operator=(const a_point2& p)
{
        x_ = p.x();
        y_ = p.y();
        return *this;
}
//-operator*=----------------------------------------------------------------
a_point2 & a_point2::operator*=(double v)
{
        x_ *= v;
        y_ *= v;
        return *this;
}
//-operator/=----------------------------------------------------------------
a_point2 & a_point2::operator/=(double v)
{
        x_ /= v;
        y_ /= v;
        return *this;
}
//-sumsq---------------------------------------------------------------------
double a_point2::sumsq() const
{
        return x_*x_+y_*y_;
}
//-norm----------------------------------------------------------------------
double a_point2::norm() const
{
        return sqrt(sumsq());
}
//-normalise-----------------------------------------------------------------
a_point2& a_point2::normalise()
{
        double n = norm();
        if (n != 0) operator/=(n);
        return *this;
}
//-rotate--------------------------------------------------------------------
a_point2& a_point2::rotate(const double angle)
{
        double xn = x_*cos(angle)+y_*sin(angle);
        double yn = -x_*sin(angle)+y_*cos(angle);
        x_ = xn; y_ = yn;
        return *this;
}
//***************************************************************************
//-operator+-----------------------------------------------------------------
a_point2 operator+(const a_point2& a, const a_point2& b)
{
        a_point2 r = a;
        return r += b;
}
//-operator------------------------------------------------------------------
a_point2 operator-(const a_point2& a, const a_point2& b)
{
        a_point2 r = a;
        return r -= b;
}
//-operator*-----------------------------------------------------------------
double operator*(const a_point2& a, const a_point2& b)
{
        a_point2 r = a;
        return r *= b;
}
//-operator*-----------------------------------------------------------------
a_point2 operator*(const double v, const a_point2& a)
{
        a_point2 r = a;
        return r *= v;
}
//-operator*-----------------------------------------------------------------
a_point2 operator*(const a_point2& a, double v)
{
        a_point2 r = a;
        return r *= v;
}
//-operator/-----------------------------------------------------------------
a_point2 operator/(const a_point2& a, double v)
{
        a_point2 r = a;
        return r /= v;
}
//---------------------------------------------------------------------------
a_point2 circle_centre(const a_point2& p1, const a_point2& p2, const a_point2& p3)
{
exit(-1);
}
//***************************************************************************
//---------------------------------------------------------------------------
void a_point2::read(std::istream &in)    
{    
  in >> x_;    
  in >> y_;    
}    
//---------------------------------------------------------------------------    
void a_point2::write(std::ostream &o) const    
{    
  o << this->x() << " ";    
  o << this->y() << " ";    
}    
//---------------------------------------------------------------------------
double angle(const a_point2& a, const a_point2& b)
{
        a_point a3(a.x(),a.y(),0.);
        a_point b3(b.x(),b.y(),0.);
        return angle(a3,b3);
}
//---------------------------------------------------------------------------
a_point2 average_rot(const a_point2& a, const a_point2& b, double f)
{
        a_point a3(a.x(),a.y(),0.);
        a_point b3(b.x(),b.y(),0.);
        a_point r = average_rot(a3,b3);
        return a_point2(r.x(),r.y());
}