| 1 | < | /* Copyright (c) 2008 The University of Notre Dame. All Rights Reserved. | 
| 1 | > | /* Copyright (c) 2008, 2010 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 | 
| 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.4 2009-11-25 20:02:01 gezelter Exp $ | 
| 47 | > | *  @version $Id$ | 
| 48 |  | * | 
| 49 |  | */ | 
| 50 |  |  | 
| 53 |  | using namespace OpenMD; | 
| 54 |  |  | 
| 55 |  |  | 
| 56 | < | Triangle::Triangle() : HaveNormal_(false), HaveCentroid_(false), | 
| 56 | > | Triangle::Triangle() : HaveNormal_(false), HaveUnitNormal_(false), HaveCentroid_(false), | 
| 57 |  | HaveArea_(false), area_(0.0), normal_(V3Zero), | 
| 58 |  | centroid_(V3Zero), facetVelocity_(V3Zero), mass_(0.0), | 
| 59 |  | a_(V3Zero), b_(V3Zero), c_(V3Zero){ | 
| 76 |  | area_ = getNormal().length() * 0.5; | 
| 77 |  | return area_; | 
| 78 |  | } | 
| 79 | < |  | 
| 79 | > | // This should return the normal for our calculations. | 
| 80 |  | Vector3d Triangle::computeNormal(){ | 
| 81 |  | HaveNormal_ = true; | 
| 82 |  | normal_ = cross(a_,b_); | 
| 83 |  | return normal_; | 
| 84 |  | } | 
| 85 | + | // This should return the normal for our calculations. | 
| 86 | + | Vector3d Triangle::computeUnitNormal(){ | 
| 87 | + | HaveUnitNormal_ = true; | 
| 88 | + | unitnormal_ = cross(a_,b_); | 
| 89 | + | unitnormal_.normalize(); | 
| 90 | + | return unitnormal_; | 
| 91 | + | } | 
| 92 |  |  | 
| 93 |  | Vector3d Triangle::computeCentroid(){ | 
| 94 |  | HaveCentroid_ = true; |