numtest.cc

Go to the documentation of this file.

//
// Copyright (c) 2005-2007 University of Southern California
// Rob Peters <rjpeters at usc dot edu>
//
// created: Tue Feb 22 11:24:07 2005
// commit: $Id: numtest.cc 10065 2007-04-12 05:54:56Z rjpeters $
// $HeadURL: file:///lab/rjpeters/svnrepo/code/trunk/groovx/src/pkgs/whitebox/numtest.cc $
//
// --------------------------------------------------------------------
//
// This file is part of GroovX.
//   [http://ilab.usc.edu/rjpeters/groovx/]
//
// GroovX is free software; you can redistribute it and/or modify it
// under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// GroovX is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with GroovX; if not, write to the Free Software Foundation,
// Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
//

#ifndef GROOVX_PKGS_WHITEBOX_NUMTEST_CC_UTC20050626084022_DEFINED
#define GROOVX_PKGS_WHITEBOX_NUMTEST_CC_UTC20050626084022_DEFINED

#include "pkgs/whitebox/numtest.h"

#include "mtx/mathspecial.h"

#include "tcl/pkg.h"

#include "rutz/unittest.h"

#include "rutz/trace.h"

namespace
{
  void testErfc()
  {
    // Compare our erfc() function with the one from cmath
    TEST_REQUIRE_APPROX(dash::erfc(0), ::erfc(0), 1e-5);
    TEST_REQUIRE_APPROX(dash::erfc(1), ::erfc(1), 1e-5);
    TEST_REQUIRE_APPROX(dash::erfc(-1), ::erfc(-1), 1e-5);

    // Compare our erfc() function with known values
    TEST_REQUIRE_APPROX(dash::erfc(0), 1, 1e-5);
    TEST_REQUIRE_APPROX(dash::erfc(1), 0.1573, 1e-5);
    TEST_REQUIRE_APPROX(dash::erfc(-1), 1.8427, 1e-5);

    // Check the symmetry properties of our erfc function
    //  erfc(x)+erfc(-x) == 2
    TEST_REQUIRE_APPROX(dash::erfc(1), 2-dash::erfc(-1), 1e-5);
    TEST_REQUIRE_APPROX(dash::erfc(2), 2-dash::erfc(-2), 1e-5);


    static rutz::prof p1("testprof/erfc/custom", __FILE__, __LINE__);
    static rutz::prof p2("testprof/erfc/builtin", __FILE__, __LINE__);

    {
      rutz::trace t(p1, false);
      double x = 0.0;
      for (int i = 0; i < 2500; ++i)
        for (int f = 1; f <= 100; ++f)
          x += dash::erfc(double(i)/1000.0);
    }

    {
      rutz::trace t(p2, false);
      double x = 0.0;
      for (int i = 0; i < 2500; ++i)
        for (int f = 1; f <= 100; ++f)
          x += ::erfc(double(i)/1000.0);
    }
  }

  void testGammaln()
  {
    // Compare custom gammaln() with builtin lgamma() function
    for (int f = 1; f <= 100; ++f)
      TEST_REQUIRE_APPROX(dash::gammaln(f), ::lgamma(f), 1e-5);

    static rutz::prof p1("testprof/gammaln/custom", __FILE__, __LINE__);
    static rutz::prof p2("testprof/gammaln/builtin", __FILE__, __LINE__);

    {
      rutz::trace t(p1, false);
      double x = 0.0;
      for (int i = 0; i < 2500; ++i)
        for (int f = 1; f <= 100; ++f)
          x += dash::gammaln(f);
    }

    {
      rutz::trace t(p2, false);
      double x = 0.0;
      for (int i = 0; i < 2500; ++i)
        for (int f = 1; f <= 100; ++f)
          x += ::lgamma(f);
    }
  }
}

extern "C"
int Numtest_Init(Tcl_Interp* interp)
{
GVX_TRACE("Numtest_Init");

  GVX_PKG_CREATE(pkg, interp, "Numtest", "4.$Revision: 10065 $");

  DEF_TEST(pkg, testErfc);
  DEF_TEST(pkg, testGammaln);

  GVX_PKG_RETURN(pkg);
}

// Need these to avoid dyld errors on Mac OS X
extern "C" int Numtest_SafeInit(Tcl_Interp*) { return 1; }

extern "C" int Numtest_Unload(Tcl_Interp* interp, int /*flags*/)
{
GVX_TRACE("Numtest_Unload");
  return tcl::pkg::destroy_on_unload(interp, "Numtest");
}

extern "C" int Numtest_SafeUnload(Tcl_Interp*, int /*flags*/) { return 1; }

static const char __attribute__((used)) vcid_groovx_pkgs_whitebox_numtest_cc_utc20050626084022[] = "$Id: numtest.cc 10065 2007-04-12 05:54:56Z rjpeters $ $HeadURL: file:
#endif // !GROOVX_PKGS_WHITEBOX_NUMTEST_CC_UTC20050626084022_DEFINED

---

The software described here is Copyright (c) 1998-2005, Rob Peters.
This page was generated Wed Dec 3 06:49:40 2008 by Doxygen version 1.5.5.