| 1 | /* | 
| 2 | * Copyright (C) 2000-2004  Object Oriented Parallel Simulation Engine (OOPSE) project | 
| 3 | * | 
| 4 | * Contact: oopse@oopse.org | 
| 5 | * | 
| 6 | * This program is free software; you can redistribute it and/or | 
| 7 | * modify it under the terms of the GNU Lesser General Public License | 
| 8 | * as published by the Free Software Foundation; either version 2.1 | 
| 9 | * of the License, or (at your option) any later version. | 
| 10 | * All we ask is that proper credit is given for our work, which includes | 
| 11 | * - but is not limited to - adding the above copyright notice to the beginning | 
| 12 | * of your source code files, and to any copyright notice that you may distribute | 
| 13 | * with programs based on this work. | 
| 14 | * | 
| 15 | * This program is distributed in the hope that it will be useful, | 
| 16 | * but WITHOUT ANY WARRANTY; without even the implied warranty of | 
| 17 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the | 
| 18 | * GNU Lesser General Public License for more details. | 
| 19 | * | 
| 20 | * You should have received a copy of the GNU Lesser General Public License | 
| 21 | * along with this program; if not, write to the Free Software | 
| 22 | * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA. | 
| 23 | * | 
| 24 | */ | 
| 25 |  | 
| 26 | /** | 
| 27 | * @file DirectionalAtom.hpp | 
| 28 | * @author    tlin | 
| 29 | * @date  10/23/2004 | 
| 30 | * @version 1.0 | 
| 31 | */ | 
| 32 |  | 
| 33 | #ifndef PRIMITIVES_DIRECTIONALATOM_HPP | 
| 34 | #define PRIMITIVES_DIRECTIONALATOM_HPP | 
| 35 |  | 
| 36 | #include "primitives/Atom.hpp" | 
| 37 | #include "types/DirectionalAtomType.hpp" | 
| 38 | namespace oopse{ | 
| 39 | class DirectionalAtom : public Atom { | 
| 40 | public: | 
| 41 | DirectionalAtom(DirectionalAtomType* dAtomType); | 
| 42 | /** | 
| 43 | * Returns the inertia tensor of this stuntdouble | 
| 44 | * @return the inertia tensor of this stuntdouble | 
| 45 | */ | 
| 46 | virtual Mat3x3d getI(); | 
| 47 |  | 
| 48 |  | 
| 49 | /** | 
| 50 | * Sets  the previous rotation matrix of this stuntdouble | 
| 51 | * @param a  new rotation matrix | 
| 52 | */ | 
| 53 | virtual void setPrevA(const RotMat3x3d& a); | 
| 54 |  | 
| 55 | /** | 
| 56 | * Sets  the current rotation matrix of this stuntdouble | 
| 57 | * @param a  new rotation matrix | 
| 58 | */ | 
| 59 | virtual void setA(const RotMat3x3d& a); | 
| 60 |  | 
| 61 | /** | 
| 62 | * Sets  the rotation matrix of this stuntdouble in specified snapshot | 
| 63 | * @param a rotation matrix to be set | 
| 64 | * @param snapshotNo | 
| 65 | * @see #getA | 
| 66 | */ | 
| 67 | virtual void setA(const RotMat3x3d& a, int snapshotNo); | 
| 68 |  | 
| 69 | /** | 
| 70 | * Left multiple rotation matrix by another rotation matrix | 
| 71 | * @param m a rotation matrix | 
| 72 | */ | 
| 73 | void rotateBy(const RotMat3x3d& m); | 
| 74 |  | 
| 75 | /** Sets the internal unit frame of this stuntdouble by three euler angles */ | 
| 76 | void setUnitFrameFromEuler(double phi, double theta, double psi); | 
| 77 |  | 
| 78 | /** | 
| 79 | * Returns the gradient of this stuntdouble | 
| 80 | * @return the gradient of this stuntdouble | 
| 81 | */ | 
| 82 | virtual std::vector<double> getGrad(); | 
| 83 |  | 
| 84 | virtual void accept(BaseVisitor* v); | 
| 85 |  | 
| 86 | protected: | 
| 87 | Mat3x3d inertiaTensor_;   /**< inertial tensor */ | 
| 88 | RotMat3x3d sU_;               /**< body fixed standard unit vector */ | 
| 89 | }; | 
| 90 |  | 
| 91 | }//namepace oopse | 
| 92 |  | 
| 93 | #endif //PRIMITIVES_DIRECTIONALATOM_HPP |