GlobiPack Package Browser (Single Doxygen Collection) Version of the Day
TestLagrPolyMeritFunc1D_UnitTests.cpp
Go to the documentation of this file.
00001 /*
00002 // @HEADER
00003 // ***********************************************************************
00004 // 
00005 //    GlobiPack: Collection of Scalar 1D globalizaton utilities
00006 //                 Copyright (2009) Sandia Corporation
00007 // 
00008 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00009 // license for use of this work by or on behalf of the U.S. Government.
00010 // 
00011 // Redistribution and use in source and binary forms, with or without
00012 // modification, are permitted provided that the following conditions are
00013 // met:
00014 //
00015 // 1. Redistributions of source code must retain the above copyright
00016 // notice, this list of conditions and the following disclaimer.
00017 //
00018 // 2. Redistributions in binary form must reproduce the above copyright
00019 // notice, this list of conditions and the following disclaimer in the
00020 // documentation and/or other materials provided with the distribution.
00021 //
00022 // 3. Neither the name of the Corporation nor the names of the
00023 // contributors may be used to endorse or promote products derived from
00024 // this software without specific prior written permission.
00025 //
00026 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00027 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00028 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00029 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00030 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00031 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00032 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00033 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00034 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00035 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00036 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00037 //
00038 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov) 
00039 // 
00040 // ***********************************************************************
00041 // @HEADER
00042 */
00043 
00044 
00045 #include "GlobiPack_TestLagrPolyMeritFunc1D.hpp"
00046 #include "Teuchos_UnitTestHarness.hpp"
00047 #include "Teuchos_Tuple.hpp"
00048 
00049 
00050 namespace {
00051 
00052 
00053 //
00054 // Helper code and declarations
00055 //
00056 
00057 
00058 using Teuchos::as;
00059 using Teuchos::outArg;
00060 using Teuchos::null;
00061 using Teuchos::RCP;
00062 using Teuchos::Array;
00063 using Teuchos::tuple;
00064 using GlobiPack::TestLagrPolyMeritFunc1D;
00065 using GlobiPack::computeValue;
00066 
00067 
00068 template<class Scalar>
00069 inline Scalar sqr(const Scalar &x) { return x*x; }
00070 
00071 
00072 double g_tol = Teuchos::ScalarTraits<double>::eps()*100.0;
00073 
00074 
00075 TEUCHOS_STATIC_SETUP()
00076 {
00077   Teuchos::UnitTestRepository::getCLP().setOption(
00078     "tol", &g_tol, "Floating point tolerance" );
00079 }
00080 
00081 
00082 //
00083 // Unit tests
00084 //
00085 
00086 
00087 //
00088 // Here we represent the simple quadratic merit function:
00089 //
00090 //   phi(alpha) = 0.5 * (alpha - 2.0)^2 + 2.0
00091 //
00092 // with the derivative:
00093 //
00094 //   Dphi(alpha) = alpha - 2.0
00095 //
00096 
00097 template<class Scalar>
00098 Scalar phi_quad1(const Scalar &alpha)
00099 {
00100   return as<Scalar>(0.5)*sqr(alpha-as<Scalar>(2.0)) + as<Scalar>(2.0);
00101 }
00102 
00103 template<class Scalar>
00104 Scalar Dphi_quad1(const Scalar &alpha)
00105 {
00106   return alpha-as<Scalar>(2.0);
00107 }
00108 
00109 TEUCHOS_UNIT_TEST_TEMPLATE_1_DECL( TestLagrPolyMeritFunc1D, basic, Scalar )
00110 {
00111   
00112   typedef Teuchos::ScalarTraits<Scalar> ST;
00113   
00114   ECHO(Array<Scalar> alphaPoints = tuple<Scalar>(0.0, 2.0, 4.0));
00115   ECHO(Array<Scalar> phiPoints = tuple<Scalar>(4.0, 2.0, 4.0));
00116   ECHO(TestLagrPolyMeritFunc1D<Scalar> meritFunc(alphaPoints, phiPoints));
00117   TEST_ASSERT(meritFunc.supportsDerivEvals());
00118 
00119   Array<Scalar> alphaTestPoints = tuple<Scalar>(0.0, 1.0, 2.0, 3.0, 4.0);
00120   for (int test_i = 0; test_i < as<int>(alphaTestPoints.size()); ++test_i) {
00121     out << "\ntest_i="<<test_i<<"\n\n";
00122     Teuchos::OSTab tab(out);
00123     ECHO(const Scalar alpha = alphaTestPoints[test_i]);
00124     out << "alpha="<<alpha<<"\n";
00125     ECHO(Scalar phi = as<Scalar>(-1.0));
00126     ECHO(Scalar Dphi = as<Scalar>(-1.0));
00127     ECHO(meritFunc.eval(alpha, outArg(phi), outArg(Dphi)));
00128     TEST_FLOATING_EQUALITY(phi, phi_quad1(alpha), g_tol);
00129     TEST_FLOATING_EQUALITY(Dphi, Dphi_quad1(alpha), g_tol);
00130   }
00131 }
00132 
00133 TEUCHOS_UNIT_TEST_TEMPLATE_1_INSTANT_REAL_SCALAR_TYPES( TestLagrPolyMeritFunc1D, basic )
00134 
00135 
00136 } // namespace
00137 
00138 
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Friends