Sophie

Sophie

distrib > Mandriva > 9.1 > ppc > media > contrib > by-pkgid > 263386785cefb9ae5d63b926d214d809 > files > 1387

mpqc-2.1.2-4mdk.ppc.rpm

<!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> &nbsp; <a class="qindex" href="hierarchy.html">Class Hierarchy</a> &nbsp; <a class="qindex" href="annotated.html">Compound List</a> &nbsp; <a class="qindex" href="files.html">File List</a> &nbsp; <a class="qindex" href="functions.html">Compound Members</a> &nbsp; <a class="qindex" href="pages.html">Related Pages</a> &nbsp; </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 &lt;cljanss@limitpt.com&gt;</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 &lt;math/isosurf/tricoef.h&gt;</font>
00036 <font class="preprocessor">#include &lt;math/isosurf/edge.h&gt;</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&lt;Edge&gt; _edges[3];
00056     Ref&lt;Vertex&gt; *_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&lt;Edge&gt;&amp; v1, <font class="keyword">const</font> Ref&lt;Edge&gt;&amp; v2, <font class="keyword">const</font> Ref&lt;Edge&gt;&amp; v3,
00061              <font class="keywordtype">unsigned</font> <font class="keywordtype">int</font> orient0 = 0);
00062     Ref&lt;Edge&gt; 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&lt;Edge&gt;&amp;) <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&lt;Edge&gt;&amp;) <font class="keyword">const</font>;
00069     ~Triangle();
00070     <font class="keywordtype">void</font> add_edges(AVLSet&lt;Ref&lt;Edge&gt; &gt;&amp;);
00071     <font class="keywordtype">void</font> add_vertices(AVLSet&lt;Ref&lt;Vertex&gt; &gt;&amp;);
00072 
00073     <font class="comment">// returns the surface area element</font>
00074     <font class="comment">// 0&lt;=r&lt;=1, 0&lt;=s&lt;=1, 0&lt;=r+s&lt;=1</font>
00075     <font class="comment">// Ref&lt;Vertex&gt; is the intepolated vertex (both point and normal)</font>
00076     <font class="keywordtype">void</font> interpolate(<font class="keyword">const</font> Ref&lt;TriInterpCoef&gt;&amp;,
00077                      <font class="keywordtype">double</font> r,<font class="keywordtype">double</font> s,<font class="keyword">const</font> Ref&lt;Vertex&gt;&amp;v, SCVector3&amp; 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&lt;Vertex&gt;&amp;v, SCVector3&amp; 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&lt;Vertex&gt;&amp;v, SCVector3&amp; dA,
00080                      <font class="keyword">const</font> Ref&lt;Volume&gt; &amp;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&lt;Vertex&gt; 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&lt;Volume&gt;&amp;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&lt;TriInterpCoef&gt; **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&lt;KeyVal&gt;&amp;);
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&lt;TriInterpCoef&gt; 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&lt;KeyVal&gt;&amp;);
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>