---

test_distance_constraint.cpp

//$Header: /home/ben/Mapper/c++/test/RCS/test_distance_constraint.cpp,v 6.4 2002/06/23 12:38:16 ben Exp $
//Unit tests of Mapper distance_constraint classes.
// 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 "distance_constraint.h"
#include "test/constraints.h"
#include "test/test.h"
#include "through_constraint.h"



void test::test_distance_constraint(void)
{
   test_start_section("distance_constraint");
   {
      test_start_section("distance_constraint::distance_constraint");
      const double u = 1.0;
      const double d = 1.0;
      distance_constraint dc( u, d );
      require( dc.uncertainty() == u );
      require( dc.wanted_distance() == d );
      require( dc.points().size() == 2 );
   };
   {
      test_start_section("distance_constraint::wanted_distance");
      const double u = 1.0;
      const double d0 = 0.0;
      const double d1 = 2.0;
      distance_constraint dc( u, d0 );
      assert( dc.wanted_distance() == d0 );
      dc.wanted_distance( d1 );
      require( dc.wanted_distance() == d1 );
   };
   {
      test_start_section("distance_constraint::energy");
      const double u = 1.0;
      const double d = 1.0;
      distance_constraint dc( u, d );
      vector< vector2 > p(2), dedp(2);
      double e;
      p[0] = vector2( 0, 0 );
      {
         p[1] = vector2( 1.0, 0.0 );
         dc.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] = vector2( 0.0, 1.0 );
         dc.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] = vector2( 0.5, 0.5 );
         dc.energy( tol, p, e, dedp );
         require( 0.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] = vector2( 2.0, 0.0 );
         dc.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] = vector2( 1.0, 1.0 );
         dc.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 );
      };
   };
}

---