ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE-4/src/math/Triangle.cpp
Revision: 3482
Committed: Fri Dec 5 16:20:39 2008 UTC (15 years, 7 months ago) by chuckv
File size: 3176 byte(s)
Log Message:
Different method again for Langevin Dynamics.

File Contents

# Content
1 /* Copyright (c) 2008 The University of Notre Dame. All Rights Reserved.
2 *
3 * The University of Notre Dame grants you ("Licensee") a
4 * non-exclusive, royalty free, license to use, modify and
5 * redistribute this software in source and binary code form, provided
6 * that the following conditions are met:
7 *
8 * 1. Acknowledgement of the program authors must be made in any
9 * publication of scientific results based in part on use of the
10 * program. An acceptable form of acknowledgement is citation of
11 * the article in which the program was described (Matthew
12 * A. Meineke, Charles F. Vardeman II, Teng Lin, Christopher
13 * J. Fennell and J. Daniel Gezelter, "OOPSE: An Object-Oriented
14 * Parallel Simulation Engine for Molecular Dynamics,"
15 * J. Comput. Chem. 26, pp. 252-271 (2005))
16 *
17 * 2. Redistributions of source code must retain the above copyright
18 * notice, this list of conditions and the following disclaimer.
19 *
20 * 3. Redistributions in binary form must reproduce the above copyright
21 * notice, this list of conditions and the following disclaimer in the
22 * documentation and/or other materials provided with the
23 * distribution.
24 *
25 * This software is provided "AS IS," without a warranty of any
26 * kind. All express or implied conditions, representations and
27 * warranties, including any implied warranty of merchantability,
28 * fitness for a particular purpose or non-infringement, are hereby
29 * excluded. The University of Notre Dame and its licensors shall not
30 * be liable for any damages suffered by licensee as a result of
31 * using, modifying or distributing the software or its
32 * derivatives. In no event will the University of Notre Dame or its
33 * licensors be liable for any lost revenue, profit or data, or for
34 * direct, indirect, special, consequential, incidental or punitive
35 * damages, however caused and regardless of the theory of liability,
36 * arising out of the use of or inability to use software, even if the
37 * University of Notre Dame has been advised of the possibility of
38 * such damages.
39 *
40 *
41 * Triangle.cpp
42 *
43 * Purpose: Provide basic triangle object for OOPSE
44 *
45 * Created by Charles F. Vardeman II on 29 July 2008.
46 * @author Charles F. Vardeman II
47 * @version $Id: Triangle.cpp,v 1.2 2008-12-05 16:20:38 chuckv Exp $
48 *
49 */
50
51 #include "math/Triangle.hpp"
52
53 using namespace oopse;
54
55
56 Triangle::Triangle() : HaveNormal_(false), HaveCentroid_(false),HaveArea_(false), area_(0.0), normal_(V3Zero),
57 centroid_(V3Zero),facetVelocity_(V3Zero), mass_(0.0),
58 a_(V3Zero),b_(V3Zero),c_(V3Zero){
59 }
60
61 void Triangle::addVertices(Vector3d P1, Vector3d P2, Vector3d P3){
62 vertices_[0] = P1;
63 vertices_[1] = P2;
64 vertices_[2] = P3;
65
66 // Compute some quantites like a,b,c
67 a_ = P1-P2;
68 b_ = P1-P3;
69 c_ = P2-P3;
70 }
71
72
73 RealType Triangle::computeArea(){
74 HaveArea_ = true;
75 area_ = getNormal().length() * 0.5;
76 return area_;
77 }
78
79 Vector3d Triangle::computeNormal(){
80 HaveNormal_ = true;
81 normal_ = cross(a_,b_);
82 return normal_;
83 }
84
85 Vector3d Triangle::computeCentroid(){
86 HaveCentroid_ = true;
87 centroid_ = (vertices_[0] + vertices_[1] + vertices_[2])/3.0;
88 return centroid_;
89 }