Sophie: mpqc-2.1.2-4mdk ppc

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>tricoef.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>tricoef.h</h1><div class="fragment"><pre>00001 <font class="comment">//</font>
00002 <font class="comment">// tricoef.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_tricoef_h</font>
00029 <font class="preprocessor"></font><font class="preprocessor">#define _math_isosurf_tricoef_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;util/ref/ref.h&gt;</font>
00036 
00037 <font class="keyword">namespace </font>sc {
00038 
00039 <font class="keyword">class </font>TriInterpCoefKey {
00040   <font class="keyword">private</font>:
00041     <font class="keywordtype">int</font> order_;
00042     <font class="keywordtype">double</font> L2_;
00043     <font class="keywordtype">double</font> L3_;
00044   <font class="keyword">public</font>:
00045     TriInterpCoefKey(<font class="keywordtype">int</font> order, <font class="keywordtype">double</font> L2, <font class="keywordtype">double</font> L3):
00046       order_(order), L2_(L2), L3_(L3) {}
00047     <font class="keywordtype">int</font> order()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> order_; }
00048     <font class="keywordtype">double</font> L1()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> 1.0 - L2_ - L3_; }
00049     <font class="keywordtype">double</font> L2()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> L2_; }
00050     <font class="keywordtype">double</font> L3()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> L3_; }
00051     <font class="keywordtype">int</font> cmp(<font class="keyword">const</font> TriInterpCoefKey&amp;t)<font class="keyword"> const </font>{
00052         <font class="keywordflow">if</font> (order_ &lt; t.order_) <font class="keywordflow">return</font> -1;
00053         <font class="keywordflow">if</font> (order_ &gt; t.order_) <font class="keywordflow">return</font> 1;
00054         <font class="keywordflow">if</font> (L2_ &lt; t.L2_) <font class="keywordflow">return</font> -1;
00055         <font class="keywordflow">if</font> (L2_ &gt; t.L2_) <font class="keywordflow">return</font> 1;
00056         <font class="keywordflow">if</font> (L3_ &lt; t.L3_) <font class="keywordflow">return</font> -1;
00057         <font class="keywordflow">if</font> (L3_ &gt; t.L3_) <font class="keywordflow">return</font> 1;
00058         <font class="keywordflow">return</font> 0;
00059       }
00060 };
00061 
00062 <font class="preprocessor">#define TriInterpCoefKeyEQ(k1,k2) ((k1).cmp(k2)==0)</font>
00063 <font class="preprocessor"></font><font class="preprocessor">#define TriInterpCoefKeyCMP(k1,k2) ((k1).cmp(k2))</font>
00064 <font class="preprocessor"></font>
00065 <font class="keyword">class </font>TriInterpCoef: <font class="keyword">public</font> RefCount {
00066     <font class="keywordtype">double</font> *coef_;
00067     <font class="keywordtype">double</font> *rderiv_;
00068     <font class="keywordtype">double</font> *sderiv_;
00069   <font class="keyword">public</font>:
00070     TriInterpCoef(<font class="keyword">const</font> TriInterpCoefKey&amp; key);
00071     ~TriInterpCoef();
00072     <font class="keywordtype">double</font>&amp; coef(<font class="keywordtype">int</font> i, <font class="keywordtype">int</font> j, <font class="keywordtype">int</font> k)<font class="keyword"> </font>{<font class="keywordflow">return</font> coef_[ijk_to_index(i,j,k)];}
00073     <font class="keywordtype">double</font>&amp; rderiv(<font class="keywordtype">int</font> i, <font class="keywordtype">int</font> j, <font class="keywordtype">int</font> k)<font class="keyword"> </font>{<font class="keywordflow">return</font> rderiv_[ijk_to_index(i,j,k)];}
00074     <font class="keywordtype">double</font>&amp; sderiv(<font class="keywordtype">int</font> i, <font class="keywordtype">int</font> j, <font class="keywordtype">int</font> k)<font class="keyword"> </font>{<font class="keywordflow">return</font> sderiv_[ijk_to_index(i,j,k)];}
00075 
00076     <font class="keyword">static</font> <font class="keywordtype">int</font>
00077     ijk_to_index(<font class="keywordtype">int</font> i, <font class="keywordtype">int</font> j, <font class="keywordtype">int</font> k)<font class="keyword"></font>
00078 <font class="keyword">    </font>{
00079       <font class="keywordtype">int</font> n = i + j + k;
00080       <font class="keywordtype">int</font> ir = n - i;
00081       <font class="keywordflow">return</font> (ir*(ir+1)&gt;&gt;1) + j;
00082     }
00083 
00084     <font class="keyword">static</font> <font class="keywordtype">int</font>
00085     order_to_nvertex(<font class="keywordtype">int</font> order)<font class="keyword"></font>
00086 <font class="keyword">    </font>{
00087       <font class="keywordflow">return</font> ((order+1)*(order+2)&gt;&gt;1);
00088     }
00089 };
00090 
00091 }
00092 
00093 <font class="preprocessor">#endif</font>
00094 <font class="preprocessor"></font>
00095 <font class="comment">// Local Variables:</font>
00096 <font class="comment">// mode: c++</font>
00097 <font class="comment">// c-file-style: "CLJ"</font>
00098 <font class="comment">// End:</font>
</div></pre><hr>
<address>
<small>

Generated at Mon Oct 14 14:16:39 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>