a_tetrahedron.h

Go to the documentation of this file.

/*
Copyright 2016 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 _A_TETRAHEDRON_H_
#define _A_TETRAHEDRON_H_
 
#include <iostream>
#include <stdlib.h>
#include "a_point.h"
#include "a_segment.h"
#include "a_triangle.h"
 
class a_tetrahedron: public a_geom_base
{
public:
        a_tetrahedron(const a_point& p1, const a_point& p2, const a_point& p3, const a_point& p4) {p_[0]=p1;p_[1]=p2;p_[2]=p3;p_[3]=p4;}
        a_tetrahedron(const a_tetrahedron& p) {for (unsigned short i=0; i<4; i++) p_[i] = p.p_[i];}
 
        virtual const std::string classname() {return "a_tetrahedron";}
        static const std::string help();
        //: get points
        inline a_point  p1() const { return p_[0]; }
        inline a_point  p2() const { return p_[1]; }
        inline a_point  p3() const { return p_[2]; }
        inline a_point  p4() const { return p_[3]; }
        inline a_point  operator[](const int i) const {return p_[i%4];}
        inline a_point  p(const int ref) const {return p_[ref%4];}
 
        //: set points
        inline void    p1(const a_point& v) { p_[0] = v; }
        inline void    p2(const a_point& v) { p_[1] = v; }
        inline void    p3(const a_point& v) { p_[2] = v; }
        inline void    p4(const a_point& v) { p_[3] = v; }
        inline void     set(const a_point& p1, const a_point& p2, const a_point& p3, const a_point& p4) {p_[0] = p1; p_[1] = p2; p_[2] = p3; p_[3] = p4;}
 
        a_segment s(const int ref) const;
        a_triangle f(const int ref) const;
        inline a_point  c() const {return (p_[0]+p_[1]+p_[2]+p_[3])/4.; }
        a_tetrahedron& translate(double x, double y, double z);
        a_tetrahedron& rotate(const a_point& x_axis, const a_point& y_axis, const a_point& z_axis);
        double V() const;
        double S() const;
        bool contains(const a_point& p) const;
        double quality() const;
        bool operator==(const a_tetrahedron& p);
        a_tetrahedron& operator=(const a_tetrahedron& p);
 
        //: input/output
        virtual void read(std::istream &i);
        virtual void write(std::ostream &o) const;
 
protected:
        a_point   p_[4];
};
 
#endif