---

test_vector_constraint.cpp

//$Header: /home/ben/Mapper/c++/test/RCS/test_vector_constraint.cpp,v 6.4 2002/06/23 12:38:17 ben Exp $
//Unit tests of Mapper vector_constraint class.
// Copyright Benedict Adamson 2002.
// This file is part of Mapper.

// Mapper 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.

// Mapper 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 Mapper; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

#include <cmath>

#include "test/constraints.h"
#include "test/test.h"
#include "vector_constraint.h"



void test::test_vector_constraint(void)
{
   test_start_section("vector_constraint");
   {
      test_start_section("vector_constraint::vector_constraint");
      const double u = 1.0;
      const double x1 = 0.0;
      const double y1 = 2.0;
      vector_constraint vc( u, x1, y1 );
      require( vc.uncertainty() == u );
      require( vc.x == x1 );
      require( vc.y == y1 );
      require( vc.points().size() == 2 );
   };
   {
      test_start_section("vector_constraint::energy");
      const double u = 1.0;
      const double x1 = 0.0;
      const double y1 = 2.0;
      vector_constraint vc( u, x1, y1 );
      vector< vector2 > p(2), dedp(2);
      double e;
      p[0] = vector2( 0.0, 0.0 );
      {
         p[1] = p[0] + vector2( x1, y1 );
         vc.energy( tol, p, e, dedp );
         require( 0 == e );
         require( 0 == dedp[0].x && 0 == dedp[0].y );
         require( 0 == dedp[1].x && 0 == dedp[1].y );
      };
      {
         p[1] = p[0] + vector2( x1+1, y1+1 );
         vc.energy( tol, p, e, dedp );
         require( 0 < e );
         require( dedp[0].x == -dedp[1].x && dedp[0].y == -dedp[1].y );
         require( 0 < dedp[1].x && 0 < dedp[1].y );
      };
      {
         p[1] = p[0] + vector2( x1+1, y1-1 );
         vc.energy( tol, p, e, dedp );
         require( 0 < e );
         require( dedp[0].x == -dedp[1].x && dedp[0].y == -dedp[1].y );
         require( 0 < dedp[1].x && dedp[1].y < 0 );
      };
      {
         p[1] = p[0] + vector2( x1-1, y1-1 );
         vc.energy( tol, p, e, dedp );
         require( 0 < e );
         require( dedp[0].x == -dedp[1].x && dedp[0].y == -dedp[1].y );
         require( dedp[1].x < 0 && dedp[1].y < 0 );
      };
      {
         p[1] = p[0] + vector2( x1-1, y1+1 );
         vc.energy( tol, p, e, dedp );
         require( 0 < e );
         require( dedp[0].x == -dedp[1].x && dedp[0].y == -dedp[1].y );
         require( dedp[1].x < 0 && 0 < dedp[1].y );
      };
   };
}

---