a_3dh.h

Go to the documentation of this file.

/*
Copyright 2010-2019 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/>.
*/
#ifndef _VE_3D_H_
#define _VE_3D_H_
 
#include "a_mat_c.h"
#include "a_point.h"
 
// .NAME a_3dh
// .SECTION Author
// Pierre Smars
 
#include <iostream>
 
class a_3dh : public a_mat_c
{
public:
        inline a_3dh() : a_mat_c(4) {(*this)(3)=1;};
        inline a_3dh(double x, double y, double z, double t=1.) : a_mat_c(4) {(*this)(0)=x; (*this)(1)=y; (*this)(2)=z; (*this)(3)=t;};
        inline a_3dh(double * x) : a_mat_c(4) {(*this)(0)=x[0]; (*this)(1)=x[1]; (*this)(2)=x[2]; (*this)(3)=1.;};
        inline a_3dh(const a_point& p) : a_mat_c(4) {(*this)(0)=p.x(); (*this)(1)=p.y(); (*this)(2)=p.z(); (*this)(3)=1.;};
        a_3dh(const a_mat& m);
        static const std::string help();
        inline void set(double x, double y, double z, double t=1.) {(*this)(0)=x; (*this)(1)=y; (*this)(2)=z; (*this)(3)=t;};
        inline void set(double * x) {(*this)(0)=x[0]; (*this)(1)=x[1]; (*this)(2)=x[2]; (*this)(3)=1.;};
        inline void set(const a_point& p) {(*this)(0)=p.x(); (*this)(1)=p.y(); (*this)(2)=p.z(); (*this)(3)=1.;};
        inline a_point point() const {return a_point((*this)(0)/(*this)(3),(*this)(1)/(*this)(3),(*this)(2)/(*this)(3));};
        inline void X(double v) {(*this)(0) = v;}
        inline void Y(double v) {(*this)(1) = v;}
        inline void Z(double v) {(*this)(2) = v;}
        inline void T(double v) {(*this)(3) = v;}
        inline double X() const {return (*this)(0);}
        inline double Y() const {return (*this)(1);}
        inline double Z() const {return (*this)(2);}
        inline double T() const {return (*this)(3);}
        inline double x() const {return (*this)(0)/(*this)(3);}
        inline double y() const {return (*this)(1)/(*this)(3);}
        inline double z() const {return (*this)(2)/(*this)(3);}
 
        //: input/output
        friend   std::istream&    operator>> (std::istream& i, a_3dh& m);
        friend   std::ostream&    operator<< (std::ostream& o, const a_3dh& m);
};
 
#endif