/home/dorival/mechsys/lib/dem/cylinder.h

Go to the documentation of this file.

/************************************************************************
 * MechSys - Open Library for Mechanical Systems                        *
 * Copyright (C) 2005 Dorival M. Pedroso, Raul Durand                   *
 * Copyright (C) 2009 Sergio Galindo                                    *
 *                                                                      *
 * 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 3 of the License, or    *
 * 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 MECHSYS_DEM_CYLINDER_H
#define MECHSYS_DEM_CYLINDER_H

// Std lib
#include <cmath>

// MechSys
#include <mechsys/dem/torus.h>

class Cylinder
{
public:
    // Constructor
    Cylinder (Torus * T0, Torus * T1, Vec3_t const * Y0, Vec3_t const * Y1);  
    Cylinder (Torus & T0, Torus & T1, Vec3_t const & Y0, Vec3_t const & Y1);  

    // Methods
    void UpdatedL ();                                                                                     
    void Draw      (std::ostream & os, double Radius=1.0, char const * Color="Blue", bool BPY=false);

    // Data
    Torus       * T0; 
    Torus       * T1; 
    Vec3_t const *Y0; 
    Vec3_t const *Y1; 
    Vec3_t        X0; 
    Vec3_t        X1; 
};




inline Cylinder::Cylinder (Torus * TheT0, Torus * TheT1, Vec3_t const * TheY0, Vec3_t const * TheY1)
    : T0(TheT0), T1(TheT1)
{
    X0 = *T0->X0;
    X1 = *T1->X0;
    Y0 = TheY0;
    Y1 = TheY1;
    T0->X0 = &X0;
    T1->X0 = &X1;

}

inline Cylinder::Cylinder (Torus & TheT0, Torus & TheT1, Vec3_t const & TheY0, Vec3_t const & TheY1)
    : T0(&TheT0), T1(&TheT1)
{
    X0 = *T0->X0;
    X1 = *T1->X0;
    Y0 = &TheY0;
    Y1 = &TheY1;
    T0->X0 = &X0;
    T1->X0 = &X1;
}

inline void Cylinder::UpdatedL ()
{
    X0 = 0.5*(*T0->X1 + *Y0);
    X1 = 0.5*(*T1->X1 + *Y1);
}

inline void Cylinder::Draw (std::ostream & os, double Radius, char const * Color, bool BPY)
{
    if (BPY)
    {
        throw new Fatal("Cylinder shape not implemented for blender");
    }
    else
    {
        double R0 = norm((*T0->X0)-(*T0->X1));
        double R1 = norm((*T1->X0)-(*T1->X1));
        os << "cone { \n <" << (*T0->X0)(0) << "," << (*T0->X0)(1) << "," << (*T0->X0)(2) << "> , " << R0 + Radius << "\n"
                     << "<" << (*T1->X0)(0) << "," << (*T1->X0)(1) << "," << (*T1->X0)(2) << "> , " << R1 + Radius << "\n"
                            << " pigment { color " << Color        <<" } }\n"
           << "cone { \n <" << (*T0->X0)(0) << "," << (*T0->X0)(1) << "," << (*T0->X0)(2) << "> , " << R0 - Radius << "\n"
                     << "<" << (*T1->X0)(0) << "," << (*T1->X0)(1) << "," << (*T1->X0)(2) << "> , " << R1 - Radius << "\n"
                            << " pigment { color " << Color        <<" } }\n";
    }
}

#endif // MECHSYS_DEM_CYLINDER_H