DGtal  1.4.beta
Functions
testInHalfPlane.cpp File Reference
DGtal Test Files
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/geometry/tools/determinant/Simple2x2DetComputer.h"
#include "DGtal/geometry/tools/determinant/AvnaimEtAl2x2DetSignComputer.h"
#include "DGtal/geometry/tools/determinant/COrientationFunctor2.h"
#include "DGtal/geometry/tools/determinant/InHalfPlaneBy2x2DetComputer.h"
#include "DGtal/geometry/tools/determinant/InHalfPlaneBySimple3x3Matrix.h"
#include "DGtal/geometry/tools/determinant/InGeneralizedDiskOfGivenRadius.h"
Include dependency graph for testInHalfPlane.cpp:

Go to the source code of this file.

Functions

template<typename OrientationFunctor >
bool testInHalfPlane (OrientationFunctor f)
 
bool testInGeneralizedDiskOfGivenRadius ()
 
int main (int argc, char **argv)
 

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program 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 this program. If not, see http://www.gnu.org/licenses/.

Author
Tristan Roussillon (trist.nosp@m.an.r.nosp@m.oussi.nosp@m.llon.nosp@m.@liri.nosp@m.s.cn.nosp@m.rs.fr ) Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), CNRS, France
Date
2013/11/22

Functions for testing models of COrientationFunctor2.

This file is part of the DGtal library.

Definition in file testInHalfPlane.cpp.

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 269 of file testInHalfPlane.cpp.

270 {
271  trace.beginBlock ( "Testing class InHalfPlane" );
272  trace.info() << "Args:";
273  for ( int i = 0; i < argc; ++i )
274  trace.info() << " " << argv[ i ];
275  trace.info() << endl;
276 
277  typedef PointVector<2, DGtal::int16_t> Point;
278 
279  bool res = true;
280 
281  typedef Simple2x2DetComputer<DGtal::int32_t, DGtal::int64_t> DetComputer;
282  typedef InHalfPlaneBy2x2DetComputer<Point, DetComputer> Functor1;
283  res = res && testInHalfPlane( Functor1() );
284 
285  typedef InHalfPlaneBySimple3x3Matrix<Point, DGtal::int32_t> Functor2;
286  res = res && testInHalfPlane( Functor2() );
287 
288  res = res && testInGeneralizedDiskOfGivenRadius();
289 
290  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
291  trace.endBlock();
292  return res ? 0 : 1;
293 }
DGtal::InHalfPlaneBy2x2DetComputer
Aim: Class that implements an orientation functor, ie. it provides a way to compute the orientation o...
Definition: InHalfPlaneBy2x2DetComputer.h:93
DGtal::InHalfPlaneBySimple3x3Matrix
Aim: Class that implements an orientation functor, ie. it provides a way to compute the orientation o...
Definition: InHalfPlaneBySimple3x3Matrix.h:92
DGtal::PointVector
Aim: Implements basic operations that will be used in Point and Vector classes.
Definition: PointVector.h:593
DGtal::Simple2x2DetComputer
Aim: Small class useful to compute the determinant of a 2x2 matrix from its four coefficients,...
Definition: Simple2x2DetComputer.h:93
DGtal::Trace::beginBlock
void beginBlock(const std::string &keyword="")
DGtal::Trace::emphase
std::ostream & emphase()
DGtal::Trace::info
std::ostream & info()
DGtal::Trace::endBlock
double endBlock()
DGtal::trace
Trace trace
Definition: Common.h:153
Point
MyPointD Point
Definition: testClone2.cpp:383
testInGeneralizedDiskOfGivenRadius
bool testInGeneralizedDiskOfGivenRadius()
Definition: testInHalfPlane.cpp:97
testInHalfPlane
bool testInHalfPlane(OrientationFunctor f)
Definition: testInHalfPlane.cpp:57

References DGtal::Trace::beginBlock(), DGtal::Trace::emphase(), DGtal::Trace::endBlock(), DGtal::Trace::info(), testInGeneralizedDiskOfGivenRadius(), testInHalfPlane(), and DGtal::trace.

◆ testInGeneralizedDiskOfGivenRadius()

bool testInGeneralizedDiskOfGivenRadius ( )

Example of a test. To be completed.

Definition at line 97 of file testInHalfPlane.cpp.

98 {
99  unsigned int nbok = 0;
100  unsigned int nb = 0;
101 
102  typedef PointVector<2, DGtal::int16_t> Point;
103  typedef AvnaimEtAl2x2DetSignComputer<DGtal::int64_t> DetComputer;
104  typedef InGeneralizedDiskOfGivenRadius<Point, DetComputer> Functor;
105  BOOST_CONCEPT_ASSERT(( concepts::COrientationFunctor2<Functor> ));
106  typedef Functor::Value Value;
107 
108  Value res;
109 
110  trace.beginBlock ( "Infinite radius..." );
111 
112  Functor f_inf; //infinite radius by default
113  trace.info() << f_inf << " " << f_inf.isValid() << endl;
114 
115  f_inf.init( Point(0,0), Point(5,2) );
116  res = f_inf( Point(7,3) );
117  trace.info() << res << " > 0 " << std::endl;
118  if ( res > NumberTraits<Value>::ZERO )
119  nbok++;
120  nb++;
121  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
122 
123  f_inf.init( Point(0,0), Point(5,2) );
124  res = f_inf( Point(8,3) );
125  trace.info() << res << " < 0 " << std::endl;
126  if ( res < NumberTraits<Value>::ZERO )
127  nbok++;
128  nb++;
129  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
130 
131  f_inf.init( Point(0,0), Point(5,2) );
132  res = f_inf( Point(10,4) );
133  trace.info() << res << " == 0 " << std::endl;
134  if ( res == NumberTraits<Value>::ZERO )
135  nbok++;
136  nb++;
137  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
138 
139  Functor f_inf2(false); //infinite radius too, but with another orientation
140  trace.info() << f_inf2 << " " << f_inf2.isValid() << endl;
141 
142  f_inf2.init( Point(0,0), Point(5,2) );
143  res = f_inf2( Point(7,3) );
144  trace.info() << res << " < 0 " << std::endl;
145  if ( res < NumberTraits<Value>::ZERO )
146  nbok++;
147  nb++;
148  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
149 
150  f_inf2.init( Point(0,0), Point(5,2) );
151  res = f_inf2( Point(8,3) );
152  trace.info() << res << " > 0 " << std::endl;
153  if ( res > NumberTraits<Value>::ZERO )
154  nbok++;
155  nb++;
156  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
157 
158  f_inf2.init( Point(0,0), Point(5,2) );
159  res = f_inf2( Point(10,4) );
160  trace.info() << res << " == 0 " << std::endl;
161  if ( res == NumberTraits<Value>::ZERO )
162  nbok++;
163  nb++;
164  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
165 
166  trace.endBlock();
167 
168  trace.beginBlock ( "Finite radius..." );
169 
170  Functor f5(true, 25, 1); //radius 5 with positive orientation
171  trace.info() << f5 << " " << f5.isValid() << endl;
172 
173  f5.init( Point(5,0), Point(0,5) );
174  res = f5( Point(-4,1) );
175  trace.info() << res << " > 0 " << std::endl;
176  if ( res > NumberTraits<Value>::ZERO )
177  nbok++;
178  nb++;
179  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
180 
181  res = f5( Point(-5,1) );
182  trace.info() << res << " < 0 " << std::endl;
183  if ( res < NumberTraits<Value>::ZERO )
184  nbok++;
185  nb++;
186  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
187 
188  res = f5( Point(-3,4) );
189  trace.info() << res << " == 0 " << std::endl;
190  if ( res == NumberTraits<Value>::ZERO )
191  nbok++;
192  nb++;
193  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
194 
195  res = f5( Point(-1,6) );
196  trace.info() << res << " < 0 " << std::endl;
197  if ( res < NumberTraits<Value>::ZERO )
198  nbok++;
199  nb++;
200  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
201 
202  res = f5( Point(-1,7) );
203  trace.info() << res << " < 0 " << std::endl;
204  if ( res < NumberTraits<Value>::ZERO )
205  nbok++;
206  nb++;
207  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
208 
209  Functor f52(false, 25, 1); //radius 5 with negative orientation
210  trace.info() << f52 << " " << f52.isValid() << endl;
211 
212  f52.init( Point(-5,0), Point(0,5) );
213  res = f52( Point(4,1) );
214  trace.info() << res << " < 0 " << std::endl;
215  if ( res < NumberTraits<Value>::ZERO )
216  nbok++;
217  nb++;
218  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
219 
220  res = f52( Point(3,1) );
221  trace.info() << res << " < 0 " << std::endl;
222  if ( res < NumberTraits<Value>::ZERO )
223  nbok++;
224  nb++;
225  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
226 
227  res = f52( Point(5,1) );
228  trace.info() << res << " > 0 " << std::endl;
229  if ( res > NumberTraits<Value>::ZERO )
230  nbok++;
231  nb++;
232  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
233 
234  res = f52( Point(3,4) );
235  trace.info() << res << " == 0 " << std::endl;
236  if ( res == NumberTraits<Value>::ZERO )
237  nbok++;
238  nb++;
239  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
240 
241  res = f52( Point(4,3) );
242  trace.info() << res << " == 0 " << std::endl;
243  if ( res == NumberTraits<Value>::ZERO )
244  nbok++;
245  nb++;
246  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
247 
248  res = f52( Point(1,6) );
249  trace.info() << res << " > 0 " << std::endl;
250  if ( res > NumberTraits<Value>::ZERO )
251  nbok++;
252  nb++;
253  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
254 
255  res = f52( Point(1,7) );
256  trace.info() << res << " > 0 " << std::endl;
257  if ( res > NumberTraits<Value>::ZERO )
258  nbok++;
259  nb++;
260  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
261 
262  trace.endBlock();
263 
264 
265  return nbok == nb;
266 }
DGtal::AvnaimEtAl2x2DetSignComputer
Aim: Class that provides a way of computing the sign of the determinant of a 2x2 matrix from its four...
Definition: AvnaimEtAl2x2DetSignComputer.h:145
DGtal::InGeneralizedDiskOfGivenRadius
Aim: This class implements an orientation functor that provides a way to determine the position of ...
Definition: InGeneralizedDiskOfGivenRadius.h:127
DGtal::NumberTraits
Aim: The traits class for all models of Cinteger.
Definition: NumberTraits.h:564
DGtal::concepts::COrientationFunctor2
Aim: This concept is a refinement of COrientationFunctor, useful for simple algebraic curves that can...
Definition: COrientationFunctor2.h:85
Functor
InHalfPlaneBySimple3x3Matrix< Point, double > Functor
Definition: testConvexHull2DReverse.cpp:49

References DGtal::Trace::beginBlock(), DGtal::Trace::endBlock(), DGtal::Trace::info(), and DGtal::trace.

Referenced by main().

◆ testInHalfPlane()

template<typename OrientationFunctor >
bool testInHalfPlane ( OrientationFunctor  f)

Example of a test. To be completed.

Parameters
fany orientation functor
Template Parameters
OrientationFunctora model of COrientationFunctor2

Definition at line 57 of file testInHalfPlane.cpp.

58 {
59  BOOST_CONCEPT_ASSERT(( concepts::COrientationFunctor2<OrientationFunctor> ));
60 
61  unsigned int nbok = 0;
62  unsigned int nb = 0;
63 
64  trace.beginBlock ( "Testing block ..." );
65  trace.info() << f << endl;
66 
67  typedef typename OrientationFunctor::Value Value;
68  typedef typename OrientationFunctor::Point Point;
69  Point a(0,0);
70  Point b(5,2);
71 
72  //first quadrant
73  f.init(a, b);
74  if (f( Point(2,1) ) == NumberTraits<Value>::ONE)
75  nbok++; //a, b, (2,1) are CCW oriented
76  nb++;
77  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
78 
79  if (f( Point(3,1) ) == -NumberTraits<Value>::ONE)
80  nbok++; //a, b, (3,1) are CW oriented
81  nb++;
82  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
83 
84  if (f( Point(10,4) ) == NumberTraits<Value>::ZERO)
85  nbok++; //a, b, (10,4) belong to the same line
86  nb++;
87  trace.info() << "(" << nbok << "/" << nb << ") " << endl;
88 
89  trace.endBlock();
90 
91  return nbok == nb;
92 }

References DGtal::Trace::beginBlock(), DGtal::Trace::endBlock(), DGtal::Trace::info(), and DGtal::trace.

Referenced by main().