g_point.cxx

Go to the documentation of this file.

/*
         Copyright 2002-2012 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 "g_point.h"
#include <math.h>
 
const double pi = 3.14159265358979;
const double r2d = 180./pi;
const double r2g = 200./pi;
 
//---------------------------------------------------------------------------
void g_point::print(std::ostream& o)
{
        g_point a(*this);
        a.normalise();
        o << "direction vector : " << a << std::endl;
        o << "distances" << std::endl;
        o << "d  " << this->norm() << std::endl;
        o << "dx " << this->x() << std::endl;
        o << "dy " << this->y() << std::endl;
        o << "dz " << this->z() << std::endl;
        o << "orientation" << std::endl;
        o << "cx " << r2d*acos(a.x()) << " deg = " << r2g*acos(a.x()) << " gon " << std::endl;
        o << "cy " << r2d*acos(a.y()) << " deg = " << r2g*acos(a.y()) << " gon " << std::endl;
        o << "cz " << r2d*acos(a.z()) << " deg = " << r2g*acos(a.z()) << " gon" << std::endl;
        g_point b(a.x(),a.y(),0.);
        b.normalise();
        if (b.y() < 0)
                o << "azimuth   " << 360.-r2d*acos(b.x()) << " deg = " << 400.-r2g*acos(b.x()) << " gon" << std::endl;
        else
                o << "azimuth   " << r2d*acos(b.x()) << " deg = " << r2g*acos(b.x()) << " gon" << std::endl;
        o << "elevation " << 90.-r2d*acos(a.z()) << " deg = " << 100.-r2g*acos(a.z()) << " gon" << std::endl;
}
//---------------------------------------------------------------------------
void g_point::print(a_text& o)
{
        g_point a(*this);
        a.normalise();
        o << "direction vector : " << a << std::endl;
        bool state0 = o.use_user_unit();
        o << "distances" << std::endl;
        o << "d  " << this->norm() << std::endl;
        o << "dx " << this->x() << std::endl;
        o << "dy " << this->y() << std::endl;
        o << "dz " << this->z() << std::endl;
        o.use_user_unit(0);
        o << "orientation" << std::endl;
        o << "cx " << r2d*acos(a.x()) << " deg = " << r2g*acos(a.x()) << " gon " << std::endl;
        o << "cy " << r2d*acos(a.y()) << " deg = " << r2g*acos(a.y()) << " gon " << std::endl;
        o << "cz " << r2d*acos(a.z()) << " deg = " << r2g*acos(a.z()) << " gon" << std::endl;
        g_point b(a.x(),a.y(),0);
        b.normalise();
        if (b.y() < 0)
                o << "azimuth   " << 360.-r2d*acos(b.x()) << " deg = " << 400.-r2g*acos(b.x()) << " gon" << std::endl;
        else
                o << "azimuth   " << r2d*acos(b.x()) << " deg = " << r2g*acos(b.x()) << " gon" << std::endl;
        o << "elevation " << 90.-r2d*acos(a.z()) << " deg = " << 100.-r2g*acos(a.z()) << " gon" << std::endl;
        if (state0) o.use_user_unit(1);
}