<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> <html><head><meta name="robots" content="noindex"> <meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1"> <title>triangle.h Source File</title> <link href="doxygen.css" rel="stylesheet" type="text/css"> </head><body bgcolor="#ffffff"> <!-- Generated by Doxygen 1.2.5 on Mon Oct 14 14:16:38 2002 --> <center> <a class="qindex" href="index.html">Main Page</a> <a class="qindex" href="hierarchy.html">Class Hierarchy</a> <a class="qindex" href="annotated.html">Compound List</a> <a class="qindex" href="files.html">File List</a> <a class="qindex" href="functions.html">Compound Members</a> <a class="qindex" href="pages.html">Related Pages</a> </center> <hr><h1>triangle.h</h1><div class="fragment"><pre>00001 <font class="comment">//</font> 00002 <font class="comment">// triangle.h</font> 00003 <font class="comment">//</font> 00004 <font class="comment">// Copyright (C) 1996 Limit Point Systems, Inc.</font> 00005 <font class="comment">//</font> 00006 <font class="comment">// Author: Curtis Janssen <cljanss@limitpt.com></font> 00007 <font class="comment">// Maintainer: LPS</font> 00008 <font class="comment">//</font> 00009 <font class="comment">// This file is part of the SC Toolkit.</font> 00010 <font class="comment">//</font> 00011 <font class="comment">// The SC Toolkit is free software; you can redistribute it and/or modify</font> 00012 <font class="comment">// it under the terms of the GNU Library General Public License as published by</font> 00013 <font class="comment">// the Free Software Foundation; either version 2, or (at your option)</font> 00014 <font class="comment">// any later version.</font> 00015 <font class="comment">//</font> 00016 <font class="comment">// The SC Toolkit is distributed in the hope that it will be useful,</font> 00017 <font class="comment">// but WITHOUT ANY WARRANTY; without even the implied warranty of</font> 00018 <font class="comment">// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</font> 00019 <font class="comment">// GNU Library General Public License for more details.</font> 00020 <font class="comment">//</font> 00021 <font class="comment">// You should have received a copy of the GNU Library General Public License</font> 00022 <font class="comment">// along with the SC Toolkit; see the file COPYING.LIB. If not, write to</font> 00023 <font class="comment">// the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.</font> 00024 <font class="comment">//</font> 00025 <font class="comment">// The U.S. Government is granted a limited license as per AL 91-7.</font> 00026 <font class="comment">//</font> 00027 00028 <font class="preprocessor">#ifndef _math_isosurf_triangle_h</font> 00029 <font class="preprocessor"></font><font class="preprocessor">#define _math_isosurf_triangle_h</font> 00030 <font class="preprocessor"></font> 00031 <font class="preprocessor">#ifdef __GNUC__</font> 00032 <font class="preprocessor"></font><font class="preprocessor">#pragma interface</font> 00033 <font class="preprocessor"></font><font class="preprocessor">#endif</font> 00034 <font class="preprocessor"></font> 00035 <font class="preprocessor">#include <math/isosurf/tricoef.h></font> 00036 <font class="preprocessor">#include <math/isosurf/edge.h></font> 00037 00038 <font class="keyword">namespace </font>sc { 00039 00040 <font class="keyword">class </font>Triangle: <font class="keyword">public</font> RefCount { 00041 <font class="keyword">protected</font>: 00042 <font class="comment">// these break gcc 2.5.8</font> 00043 <font class="comment">//unsigned int _order:5;</font> 00044 <font class="comment">//unsigned int _orientation0:1;</font> 00045 <font class="comment">//unsigned int _orientation1:1;</font> 00046 <font class="comment">//unsigned int _orientation2:1;</font> 00047 <font class="comment">//unsigned char _order;</font> 00048 <font class="comment">//unsigned char _orientation0;</font> 00049 <font class="comment">//unsigned char _orientation1;</font> 00050 <font class="comment">//unsigned char _orientation2;</font> 00051 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> _order; 00052 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> _orientation0; 00053 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> _orientation1; 00054 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> _orientation2; 00055 Ref<Edge> _edges[3]; 00056 Ref<Vertex> *_vertices; 00057 <font class="keyword">public</font>: 00058 <font class="keyword">enum</font> {max_order = 10}; 00059 00060 Triangle(<font class="keyword">const</font> Ref<Edge>& v1, <font class="keyword">const</font> Ref<Edge>& v2, <font class="keyword">const</font> Ref<Edge>& v3, 00061 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> orient0 = 0); 00062 Ref<Edge> edge(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> _edges[i]; }; 00063 <font class="keywordtype">int</font> contains(<font class="keyword">const</font> Ref<Edge>&) <font class="keyword">const</font>; 00064 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> orientation(<font class="keywordtype">int</font> i)<font class="keyword"> const</font> 00065 <font class="keyword"> </font>{ 00066 <font class="keywordflow">return</font> i==0?_orientation0:i==1?_orientation1:_orientation2; 00067 } 00068 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> orientation(<font class="keyword">const</font> Ref<Edge>&) <font class="keyword">const</font>; 00069 ~Triangle(); 00070 <font class="keywordtype">void</font> add_edges(AVLSet<Ref<Edge> >&); 00071 <font class="keywordtype">void</font> add_vertices(AVLSet<Ref<Vertex> >&); 00072 00073 <font class="comment">// returns the surface area element</font> 00074 <font class="comment">// 0<=r<=1, 0<=s<=1, 0<=r+s<=1</font> 00075 <font class="comment">// Ref<Vertex> is the intepolated vertex (both point and normal)</font> 00076 <font class="keywordtype">void</font> interpolate(<font class="keyword">const</font> Ref<TriInterpCoef>&, 00077 <font class="keywordtype">double</font> r,<font class="keywordtype">double</font> s,<font class="keyword">const</font> Ref<Vertex>&v, SCVector3& dA); 00078 <font class="keywordtype">void</font> interpolate(<font class="keywordtype">double</font> r,<font class="keywordtype">double</font> s,<font class="keyword">const</font> Ref<Vertex>&v, SCVector3& dA); 00079 <font class="keywordtype">void</font> interpolate(<font class="keywordtype">double</font> r,<font class="keywordtype">double</font> s,<font class="keyword">const</font> Ref<Vertex>&v, SCVector3& dA, 00080 <font class="keyword">const</font> Ref<Volume> &vol, <font class="keywordtype">double</font> isovalue); 00081 00082 <font class="comment">// returns a corner vertex from the triangle</font> 00083 <font class="comment">// i = 0 is the (0,0) vertex (or L1 = 1, L2 = 0, L3 = 0)</font> 00084 <font class="comment">// i = 1 is the (r=1,s=0) vertex (or L1 = 0, L2 = 1, L3 = 0)</font> 00085 <font class="comment">// i = 2 is the (r=0,s=1) vertex (or L1 = 0, L2 = 0, L3 = 1)</font> 00086 Ref<Vertex> vertex(<font class="keywordtype">int</font> i); 00087 00088 <font class="keywordtype">double</font> flat_area(); 00089 00090 <font class="comment">// flip the orientation</font> 00091 <font class="keywordtype">void</font> flip(); 00092 00093 <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> order()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> _order; } 00094 00095 <font class="keywordtype">void</font> set_order(<font class="keywordtype">int</font> order, <font class="keyword">const</font> Ref<Volume>&vol,<font class="keywordtype">double</font> isovalue); 00096 }; 00097 00098 00099 00100 <font class="keyword">class </font>TriangleIntegrator: <font class="keyword">public</font> DescribedClass { 00101 <font class="keyword">private</font>: 00102 <font class="keywordtype">int</font> _n; 00103 <font class="keywordtype">double</font>* _r; 00104 <font class="keywordtype">double</font>* _s; 00105 <font class="keywordtype">double</font>* _w; 00106 <font class="comment">// precomputed interpolation coefficients for triangles of various orders</font> 00107 Ref<TriInterpCoef> **coef_; <font class="comment">// max_order by _n</font> 00108 <font class="keyword">protected</font>: 00109 <font class="keywordtype">void</font> set_r(<font class="keywordtype">int</font> i,<font class="keywordtype">double</font> r); 00110 <font class="keywordtype">void</font> set_s(<font class="keywordtype">int</font> i,<font class="keywordtype">double</font> s); 00111 <font class="keywordtype">void</font> set_w(<font class="keywordtype">int</font> i,<font class="keywordtype">double</font> w); 00112 <font class="keywordtype">void</font> init_coef(); 00113 <font class="keywordtype">void</font> clear_coef(); 00114 <font class="keyword">public</font>: 00115 TriangleIntegrator(<font class="keyword">const</font> Ref<KeyVal>&); 00116 TriangleIntegrator(<font class="keywordtype">int</font> n); 00117 <font class="keyword">virtual</font> ~TriangleIntegrator(); 00118 <font class="keyword">inline</font> <font class="keywordtype">double</font> w(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> _w[i]; } 00119 <font class="keyword">inline</font> <font class="keywordtype">double</font> r(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> _r[i]; } 00120 <font class="keyword">inline</font> <font class="keywordtype">double</font> s(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> _s[i]; } 00121 <font class="keyword">inline</font> <font class="keywordtype">int</font> n()<font class="keyword"> </font>{ <font class="keywordflow">return</font> _n; } 00122 <font class="keyword">virtual</font> <font class="keywordtype">void</font> set_n(<font class="keywordtype">int</font> n); 00123 Ref<TriInterpCoef> coef(<font class="keywordtype">int</font> order, <font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> coef_[order-1][i]; } 00124 }; 00125 00126 00127 <font class="keyword">class </font>GaussTriangleIntegrator: <font class="keyword">public</font> TriangleIntegrator { 00128 <font class="keyword">private</font>: 00129 <font class="keywordtype">void</font> init_rw(<font class="keywordtype">int</font> order); 00130 <font class="keyword">public</font>: 00131 GaussTriangleIntegrator(<font class="keyword">const</font> Ref<KeyVal>&); 00132 GaussTriangleIntegrator(<font class="keywordtype">int</font> order); 00133 ~GaussTriangleIntegrator(); 00134 <font class="keywordtype">void</font> set_n(<font class="keywordtype">int</font> n); 00135 }; 00136 00137 } 00138 00139 <font class="preprocessor">#endif</font> 00140 <font class="preprocessor"></font> 00141 <font class="comment">// Local Variables:</font> 00142 <font class="comment">// mode: c++</font> 00143 <font class="comment">// c-file-style: "CLJ"</font> 00144 <font class="comment">// End:</font> </div></pre><hr> <address> <small> Generated at Mon Oct 14 14:16:38 2002 for <a href="http://aros.ca.sandia.gov/~cljanss/mpqc">MPQC</a> 2.1.2 using the documentation package <a href="http://www.stack.nl/~dimitri/doxygen/index.html">Doxygen</a> 1.2.5. </small> </address> </body> </html>