<!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>gaussbas.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:36 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>gaussbas.h</h1><div class="fragment"><pre>00001 <font class="comment">//</font> 00002 <font class="comment">// gaussbas.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 _chemistry_qc_basis_gaussbas_h</font> 00029 <font class="preprocessor"></font><font class="preprocessor">#define _chemistry_qc_basis_gaussbas_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 <iostream></font> 00036 00037 <font class="preprocessor">#include <util/state/state.h></font> 00038 <font class="preprocessor">#include <util/state/array.h></font> 00039 <font class="preprocessor">#include <util/keyval/keyval.h></font> 00040 <font class="preprocessor">#include <math/scmat/matrix.h></font> 00041 <font class="preprocessor">#include <math/scmat/vector3.h></font> 00042 <font class="preprocessor">#include <chemistry/molecule/molecule.h></font> 00043 00044 <font class="keyword">namespace </font>sc { 00045 00046 <font class="keyword">class </font>GaussianShell; 00047 <font class="keyword">class </font>BasisFileSet; 00048 <font class="keyword">class </font>Integral; 00049 00050 <font class="keyword">class </font>CartesianIter; 00051 <font class="keyword">class </font>SphericalTransformIter; 00052 <a name="l00114"></a><a class="code" href="class_sc__GaussianBasisSet.html">00114</a> <font class="keyword">class </font>GaussianBasisSet: <font class="keyword">public</font> SavableState 00115 { 00116 <font class="keyword">private</font>: 00117 <font class="keywordtype">char</font>* name_; 00118 GaussianShell** shell_; 00119 Arrayint shell_to_function_; 00120 Arrayint function_to_shell_; 00121 00122 Ref<Molecule> molecule_; 00123 00124 Ref<SCMatrixKit> matrixkit_; 00125 Ref<SCMatrixKit> so_matrixkit_; 00126 RefSCDimension basisdim_; 00127 00128 <font class="keywordtype">int</font> ncenter_; 00129 SSBArrayint shell_to_center_; 00130 SSBArrayint center_to_shell_; 00131 SSBArrayint center_to_nshell_; 00132 SSBArrayint center_to_nbasis_; 00133 00134 <font class="keywordtype">int</font> nshell_; 00135 <font class="keywordtype">int</font> nbasis_; 00136 <font class="keywordtype">int</font> nprim_; 00137 00138 <font class="keywordtype">void</font> recursively_get_shell(<font class="keywordtype">int</font>&,Ref<KeyVal>&, 00139 <font class="keyword">const</font> <font class="keywordtype">char</font>*,<font class="keyword">const</font> <font class="keywordtype">char</font>*,BasisFileSet&, 00140 <font class="keywordtype">int</font>,<font class="keywordtype">int</font>,<font class="keywordtype">int</font>); 00141 00142 <font class="keywordtype">void</font> init(Ref<Molecule>&,Ref<KeyVal>&, 00143 BasisFileSet&, 00144 <font class="keywordtype">int</font> have_userkeyval, 00145 <font class="keywordtype">int</font> pure); 00146 <font class="keywordtype">void</font> init2(<font class="keywordtype">int</font> skip_ghosts=0); 00147 00148 <font class="keyword">protected</font>: 00149 <a class="code" href="class_sc__GaussianBasisSet.html#a0">GaussianBasisSet</a>(<font class="keyword">const</font> GaussianBasisSet&); 00150 <font class="keyword">virtual</font> <font class="keywordtype">void</font> set_matrixkit(<font class="keyword">const</font> Ref<SCMatrixKit>&); 00151 00152 <font class="keyword">public</font>: <a name="l00154"></a><a class="code" href="class_sc__GaussianBasisSet__ValueData.html">00154</a> <font class="keyword">class </font>ValueData { 00155 <font class="keyword">protected</font>: 00156 CartesianIter **civec_; 00157 SphericalTransformIter **sivec_; 00158 <font class="keywordtype">int</font> maxam_; 00159 <font class="keyword">public</font>: 00160 ValueData(<font class="keyword">const</font> Ref<GaussianBasisSet> &, <font class="keyword">const</font> Ref<Integral> &); 00161 ~ValueData(); 00162 CartesianIter **civec()<font class="keyword"> </font>{ <font class="keywordflow">return</font> civec_; } 00163 SphericalTransformIter **sivec()<font class="keyword"> </font>{ <font class="keywordflow">return</font> sivec_; } 00164 }; 00165 00276 <a class="code" href="class_sc__GaussianBasisSet.html#a0">GaussianBasisSet</a>(<font class="keyword">const</font> Ref<KeyVal>&); 00277 <a class="code" href="class_sc__GaussianBasisSet.html#a0">GaussianBasisSet</a>(StateIn&); 00278 <font class="keyword">virtual</font> ~GaussianBasisSet(); 00279 00280 <font class="keywordtype">void</font> <a class="code" href="class_sc__GaussianBasisSet.html#a3">save_data_state</a>(StateOut&); 00281 <a name="l00283"></a><a class="code" href="class_sc__GaussianBasisSet.html#a4">00283</a> <font class="keyword">const</font> <font class="keywordtype">char</font>* <a class="code" href="class_sc__GaussianBasisSet.html#a4">name</a>()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> name_; } 00284 <a name="l00286"></a><a class="code" href="class_sc__GaussianBasisSet.html#a5">00286</a> Ref<Molecule> <a class="code" href="class_sc__GaussianBasisSet.html#a5">molecule</a>()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> molecule_; } <a name="l00288"></a><a class="code" href="class_sc__GaussianBasisSet.html#a6">00288</a> Ref<SCMatrixKit> <a class="code" href="class_sc__GaussianBasisSet.html#a6">matrixkit</a>()<font class="keyword"> </font>{ <font class="keywordflow">return</font> matrixkit_; } <a name="l00290"></a><a class="code" href="class_sc__GaussianBasisSet.html#a7">00290</a> Ref<SCMatrixKit> <a class="code" href="class_sc__GaussianBasisSet.html#a7">so_matrixkit</a>()<font class="keyword"> </font>{ <font class="keywordflow">return</font> so_matrixkit_; } <a name="l00292"></a><a class="code" href="class_sc__GaussianBasisSet.html#a8">00292</a> RefSCDimension <a class="code" href="class_sc__GaussianBasisSet.html#a8">basisdim</a>()<font class="keyword"> </font>{ <font class="keywordflow">return</font> basisdim_; } 00293 00295 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a9">ncenter</a>() <font class="keyword">const</font>; <a name="l00297"></a><a class="code" href="class_sc__GaussianBasisSet.html#a10">00297</a> <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a10">nshell</a>()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> nshell_; } 00299 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a11">nshell_on_center</a>(<font class="keywordtype">int</font> icenter) <font class="keyword">const</font>; 00302 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a12">shell_on_center</a>(<font class="keywordtype">int</font> icenter, <font class="keywordtype">int</font> shell) <font class="keyword">const</font>; <a name="l00304"></a><a class="code" href="class_sc__GaussianBasisSet.html#a13">00304</a> <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a13">shell_to_center</a>(<font class="keywordtype">int</font> ishell)<font class="keyword"> const </font>{ <font class="keywordflow">return</font> shell_to_center_(ishell); } <a name="l00306"></a><a class="code" href="class_sc__GaussianBasisSet.html#a14">00306</a> <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a14">nbasis</a>()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> nbasis_; } 00308 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a15">nbasis_on_center</a>(<font class="keywordtype">int</font> icenter) <font class="keyword">const</font>; <a name="l00310"></a><a class="code" href="class_sc__GaussianBasisSet.html#a16">00310</a> <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a16">nprimitive</a>()<font class="keyword"> const </font>{ <font class="keywordflow">return</font> nprim_; } 00311 00313 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a17">max_nfunction_in_shell</a>() <font class="keyword">const</font>; 00316 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a18">max_ncartesian_in_shell</a>(<font class="keywordtype">int</font> aminc=0) <font class="keyword">const</font>; 00318 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a19">max_angular_momentum</a>() <font class="keyword">const</font>; 00320 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a20">max_ncontraction</a>() <font class="keyword">const</font>; 00323 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a21">max_am_for_contraction</a>(<font class="keywordtype">int</font> con) <font class="keyword">const</font>; 00325 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a22">max_cartesian</a>() <font class="keyword">const</font>; 00326 <a name="l00328"></a><a class="code" href="class_sc__GaussianBasisSet.html#a23">00328</a> <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a23">shell_to_function</a>(<font class="keywordtype">int</font> i)<font class="keyword"> const </font>{ <font class="keywordflow">return</font> shell_to_function_(i); } 00330 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a24">function_to_shell</a>(<font class="keywordtype">int</font> i) <font class="keyword">const</font>; 00331 <a name="l00333"></a><a class="code" href="class_sc__GaussianBasisSet.html#a25">00333</a> <font class="keyword">const</font> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a25">operator()</a>(<font class="keywordtype">int</font> i)<font class="keyword"> const </font>{ <font class="keywordflow">return</font> *shell_[i]; } <a name="l00335"></a><a class="code" href="class_sc__GaussianBasisSet.html#a26">00335</a> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a25">operator()</a>(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> *shell_[i]; } <a name="l00337"></a><a class="code" href="class_sc__GaussianBasisSet.html#a27">00337</a> <font class="keyword">const</font> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a27">operator[]</a>(<font class="keywordtype">int</font> i)<font class="keyword"> const </font>{ <font class="keywordflow">return</font> *shell_[i]; } <a name="l00339"></a><a class="code" href="class_sc__GaussianBasisSet.html#a28">00339</a> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a27">operator[]</a>(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> *shell_[i]; } <a name="l00341"></a><a class="code" href="class_sc__GaussianBasisSet.html#a29">00341</a> <font class="keyword">const</font> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a29">shell</a>(<font class="keywordtype">int</font> i)<font class="keyword"> const </font>{ <font class="keywordflow">return</font> *shell_[i]; } <a name="l00343"></a><a class="code" href="class_sc__GaussianBasisSet.html#a30">00343</a> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a29">shell</a>(<font class="keywordtype">int</font> i)<font class="keyword"> </font>{ <font class="keywordflow">return</font> *shell_[i]; } 00344 00346 <font class="keyword">const</font> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a25">operator()</a>(<font class="keywordtype">int</font> icenter,<font class="keywordtype">int</font> ishell) <font class="keyword">const</font>; 00348 GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a25">operator()</a>(<font class="keywordtype">int</font> icenter,<font class="keywordtype">int</font> ishell); <a name="l00350"></a><a class="code" href="class_sc__GaussianBasisSet.html#a33">00350</a> <font class="keyword">const</font> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a29">shell</a>(<font class="keywordtype">int</font> i,<font class="keywordtype">int</font> j)<font class="keyword"> const </font>{ <font class="keywordflow">return</font> <a class="code" href="class_sc__GaussianBasisSet.html#a25">operator()</a>(i,j); } <a name="l00352"></a><a class="code" href="class_sc__GaussianBasisSet.html#a34">00352</a> GaussianShell& <a class="code" href="class_sc__GaussianBasisSet.html#a29">shell</a>(<font class="keywordtype">int</font> i,<font class="keywordtype">int</font> j)<font class="keyword"> </font>{ <font class="keywordflow">return</font> <a class="code" href="class_sc__GaussianBasisSet.html#a25">operator()</a>(i,j); } 00353 00356 <font class="keywordtype">double</font> <a class="code" href="class_sc__GaussianBasisSet.html#a35">r</a>(<font class="keywordtype">int</font> icenter,<font class="keywordtype">int</font> xyz) <font class="keyword">const</font>; 00357 00360 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a36">values</a>(<font class="keyword">const</font> SCVector3& r, ValueData *, <font class="keywordtype">double</font>* basis_values) <font class="keyword">const</font>; 00365 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a37">grad_values</a>(<font class="keyword">const</font> SCVector3& r, ValueData *, 00366 <font class="keywordtype">double</font>*g_values,<font class="keywordtype">double</font>* basis_values=0) <font class="keyword">const</font>; 00371 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a38">hessian_values</a>(<font class="keyword">const</font> SCVector3& r, ValueData *, <font class="keywordtype">double</font> *h_values, 00372 <font class="keywordtype">double</font>*g_values=0,<font class="keywordtype">double</font>* basis_values=0) <font class="keyword">const</font>; 00375 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a39">shell_values</a>(<font class="keyword">const</font> SCVector3& r, <font class="keywordtype">int</font> sh, 00376 ValueData *, <font class="keywordtype">double</font>* basis_values) <font class="keyword">const</font>; 00380 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a40">grad_shell_values</a>(<font class="keyword">const</font> SCVector3& r, <font class="keywordtype">int</font> sh, 00381 ValueData *, 00382 <font class="keywordtype">double</font>*g_values, <font class="keywordtype">double</font>* basis_values=0) <font class="keyword">const</font>; 00386 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a41">hessian_shell_values</a>(<font class="keyword">const</font> SCVector3& r, <font class="keywordtype">int</font> sh, 00387 ValueData *, <font class="keywordtype">double</font> *h_values, 00388 <font class="keywordtype">double</font>*g_values=0,<font class="keywordtype">double</font>* basis_values=0) <font class="keyword">const</font>; 00389 00391 <font class="keywordtype">int</font> <a class="code" href="class_sc__GaussianBasisSet.html#a42">equiv</a>(<font class="keyword">const</font> Ref<GaussianBasisSet> &b); 00392 00394 <font class="keywordtype">void</font> <a class="code" href="class_sc__GaussianBasisSet.html#a43">print_brief</a>(std::ostream& =ExEnv::out0()) <font class="keyword">const</font>; 00396 <font class="keywordtype">void</font> <a class="code" href="class_sc__GaussianBasisSet.html#a44">print</a>(std::ostream& =ExEnv::out0()) <font class="keyword">const</font>; 00397 }; 00398 00399 } 00400 00401 <font class="preprocessor">#endif</font> 00402 <font class="preprocessor"></font> 00403 <font class="comment">// Local Variables:</font> 00404 <font class="comment">// mode: c++</font> 00405 <font class="comment">// c-file-style: "CLJ"</font> 00406 <font class="comment">// End:</font> </div></pre><hr> <address> <small> Generated at Mon Oct 14 14:16:36 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>