HermiteUnitTests.cpp

Go to the documentation of this file.
00001 // $Id$ 
00002 // $Source$ 
00003 // @HEADER
00004 // ***********************************************************************
00005 // 
00006 //                           Sacado Package
00007 //                 Copyright (2006) Sandia Corporation
00008 // 
00009 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
00010 // the U.S. Government retains certain rights in this software.
00011 // 
00012 // This library is free software; you can redistribute it and/or modify
00013 // it under the terms of the GNU Lesser General Public License as
00014 // published by the Free Software Foundation; either version 2.1 of the
00015 // License, or (at your option) any later version.
00016 //  
00017 // This library is distributed in the hope that it will be useful, but
00018 // WITHOUT ANY WARRANTY; without even the implied warranty of
00019 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00020 // Lesser General Public License for more details.
00021 //  
00022 // You should have received a copy of the GNU Lesser General Public
00023 // License along with this library; if not, write to the Free Software
00024 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00025 // USA
00026 // Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps
00027 // (etphipp@sandia.gov).
00028 // 
00029 // ***********************************************************************
00030 // @HEADER
00031 
00032 #include "HermiteUnitTests.hpp"
00033 
00034 CPPUNIT_TEST_SUITE_REGISTRATION( HermiteUnitTest );
00035 
00036 HermiteUnitTest::HermiteUnitTest() :
00037   urand(0.0, 1.0), tol_a(1.0e-15), tol_r(1.0e-14) {
00038   std::vector< Teuchos::RCP<const Stokhos::OneDOrthogPolyBasis<int,double> > > bases(1); 
00039   bases[0] = Teuchos::rcp(new basis_type(0));
00040   Teuchos::RCP<const Stokhos::CompletePolynomialBasis<int,double> > basis = 
00041     Teuchos::rcp(new Stokhos::CompletePolynomialBasis<int,double>(bases));
00042   expansion = Teuchos::rcp(new exp_type(basis));
00043 }
00044 
00045 HermiteUnitTest::HermiteUnitTest(double absolute_tolerance, 
00046          double relative_tolerance) :
00047   urand(0.0, 1.0), 
00048   tol_a(absolute_tolerance), 
00049   tol_r(relative_tolerance) {}
00050 
00051 void HermiteUnitTest::setUp() {
00052   double val;
00053 
00054   val = urand.number();
00055   ac = pce_type(expansion);
00056   ac.fastAccessCoeff(0) = val;
00057   a = val;
00058   
00059   val = urand.number();
00060   bc = pce_type(expansion);
00061   bc.fastAccessCoeff(0) = val;
00062   b = val;
00063 
00064   cc = pce_type(1.123);
00065   c = 1.123;
00066 }
00067 
00068 void HermiteUnitTest::tearDown() {}
00069 
00070 void HermiteUnitTest::comparePCEs(const pce_type& xc, double x) {
00071 
00072   // Compare sizes
00073   CPPUNIT_ASSERT(xc.size() == 1);
00074   
00075   // Compare hasFastAccess
00076   CPPUNIT_ASSERT(xc.hasFastAccess(0) == true);
00077   
00078   // Compare values
00079   compareDoubles(xc.coeff(0), x);
00080 }
00081 
00082 void HermiteUnitTest::compareDoubles(double x, double y) {
00083   CPPUNIT_ASSERT( fabs(x-y) < tol_a + tol_r*fabs(x) );
00084 }

Generated on Wed May 12 21:39:33 2010 for Sacado Package Browser (Single Doxygen Collection) by  doxygen 1.4.7