source: branches/0.10-stable/test/utility.cc @ 3093

Last change on this file since 3093 was 3093, checked in by Peter, 8 years ago

implement replace(3), which fixes #768

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 17.6 KB
Line 
1// $Id: utility.cc 3093 2013-10-21 04:07:37Z peter $
2
3/*
4  Copyright (C) 2006 Jari Häkkinen, Markus Ringnér
5  Copyright (C) 2007, 2008 Jari Häkkinen, Peter Johansson
6  Copyright (C) 2009, 2010, 2011 Peter Johansson
7  Copyright (C) 2012 Jari Häkkinen, Peter Johansson
8
9  This file is part of the yat library, http://dev.thep.lu.se/yat
10
11  The yat library is free software; you can redistribute it and/or
12  modify it under the terms of the GNU General Public License as
13  published by the Free Software Foundation; either version 3 of the
14  License, or (at your option) any later version.
15
16  The yat library is distributed in the hope that it will be useful,
17  but WITHOUT ANY WARRANTY; without even the implied warranty of
18  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
19  General Public License for more details.
20
21  You should have received a copy of the GNU General Public License
22  along with yat. If not, see <http://www.gnu.org/licenses/>.
23*/
24
25#include <config.h>
26
27#include "Suite.h"
28
29#include "yat/utility/Exception.h"
30#include "yat/utility/Matrix.h"
31#include "yat/utility/utility.h"
32#include "yat/utility/Vector.h"
33#include "yat/utility/VectorConstView.h"
34#include "yat/utility/sort_index.h"
35#include "yat/utility/stl_utility.h"
36
37#include <boost/concept_archetype.hpp>
38
39#include <cassert>
40#include <cerrno>
41#include <fstream>
42#include <iostream>
43#include <limits>
44#include <list>
45#include <map>
46#include <string>
47#include <vector>
48
49using namespace theplu::yat;
50void test_errno_error(test::Suite& suite);
51void test_errno_error_func(void);
52void test_get_map(test::Suite&);
53void test_load(test::Suite&);
54void test_load2(test::Suite&);
55void test_inverse(test::Suite&);
56
57void test_dirname(test::Suite& suite);
58void test_basename(test::Suite& suite);
59void test_fnmatch(test::Suite& suite);
60bool test_fnmatch(bool, std::string, std::string);
61void test_replace(test::Suite& suite);
62
63template<typename InputIterator, typename Key>
64void test_inverse_validate(InputIterator first, InputIterator last,
65                           const std::map<Key, std::vector<size_t> >& m,
66                           test::Suite& suite);
67template<typename InputIterator, typename Key>
68void test_inverse_validate(InputIterator first,
69                           const std::multimap<Key, size_t>& m,
70                           test::Suite& suite);
71
72void test_sort_index(test::Suite& suite);
73void test_less_nan(test::Suite& suite);
74void test_ptr_compare(test::Suite& suite);
75void test_compose_functors(test::Suite& suite);
76void test_dereferencer(test::Suite& suite);
77
78int main(int argc, char* argv[])
79{
80  using namespace theplu::yat;
81  test::Suite suite(argc, argv);
82  suite.err() << "testing utility ... " << std::endl;
83
84  // test float/double
85  std::string s("1.2");
86  if (!utility::is<double>(s)){
87    suite.add(false);
88  }
89  else if (!utility::is<float>(s)) {
90    suite.add(false);
91  }
92  else if (utility::is<int>(s)) {
93    suite.add(false);
94  }
95  else if (utility::is_nan(s)) {
96    suite.add(false);
97  }
98
99  // test int
100  s="23";
101  if (!utility::is<double>(s)){
102    suite.add(false);
103  }
104  else if (!utility::is<float>(s)) {
105    suite.add(false);
106  }
107  else if (!utility::is<int>(s)) {
108    suite.add(false);
109  }
110  else if (utility::is_nan(s)) {
111    suite.add(false);
112  }
113
114  // test nan
115  s=" nAn  ";
116  if (!utility::is<double>(s)){
117    suite.add(false);
118    suite.err() << "error: " << s << " is a double\n";
119  }
120  /* we don't require NaN for float and int
121  else if (!utility::is_float(s)) {
122    suite.add(false);
123    suite.err() << "error: " << s << " is a float\n";
124  }
125  else if (!utility::is_int(s)) {
126    suite.add(false);
127    suite.err() << "error: " << s << " is a int\n";
128  }
129  */
130  else if (!utility::is_nan(s)) {
131    suite.add(false);
132    suite.err() << "error: " << s << " is nan\n";
133
134  }
135
136  // testing trailing values
137  s=" 23 23   ";
138  if (utility::is<double>(s)){
139    suite.add(false);
140  }
141  else if (utility::is<float>(s)) {
142    suite.add(false);
143  }
144  else if (utility::is<int>(s)) {
145    suite.add(false);
146  }
147  else if (utility::is_nan(s)) {
148    suite.add(false);
149  }
150
151  if (utility::convert<double>("1.23")!=1.23)
152    suite.add(false);
153  utility::convert<double>("-inf");
154  utility::convert<double>("inf");
155  utility::convert<double>("NaN");
156
157  suite.add(utility::convert(1.23)=="1.23");
158  suite.add(utility::convert(1)=="1");
159
160  suite.add(utility::is<unsigned int>("1"));
161  suite.add(!utility::is<unsigned int>("   "));
162  suite.add(!utility::is<unsigned int>("   -1"));
163  // test -1 cannot be assigned to unsigned int
164  try {
165    unsigned int x = utility::convert<unsigned int>("-1");
166    suite.err() << "convert<unsigned>(\"-1\") did not throw\n";
167    suite.err() << "x: " << x << "\n";
168    suite.add(false);
169    exit(1);
170  }
171  catch (std::runtime_error& e) {
172    suite.out() << "expected error: " << e.what() << "\n";
173  }
174
175  if (!suite.add(utility::is<double>("-inf")))
176    suite.err() << "is<double>(\"-inf\") should return true\n";
177  if (!suite.add(utility::is<double>("inf")))
178    suite.err() << "is<double>(\"inf\") should return true\n";
179  if (!suite.add(utility::is<double>("NaN")))
180    suite.err() << "is<double>(\"NaN\") should return true\n";
181  if (!suite.add(utility::is<double>("1.23")))
182    suite.err() << "is<double>(\"1.23\") should return true\n";
183  if (!suite.add(!utility::is<double>("1.23.2")))
184    suite.err() << "is<double>(\"1.23.2\") should return false\n";
185  if (!suite.add(!utility::is<int>("1.23")))
186    suite.err() << "is<int>(\"1.23\") should return false\n";
187  if (!suite.add(!utility::is<int>("")))
188    suite.err() << "is<int>(\"\") should return false\n";
189
190  {
191    utility::Log<long double> f;
192    utility::Log<double> f2(1.3);
193    f(2.0);
194    utility::Exp<double> e;
195    e(3.2);
196  }
197
198  test_get_map(suite);
199  test_inverse(suite);
200  test_sort_index(suite);
201  test_less_nan(suite);
202  test_errno_error(suite);
203  test_ptr_compare(suite);
204  test_load(suite);
205  test_load2(suite);
206
207  test_compose_functors(suite);
208
209  double x = 8.0;
210  suite.add(suite.equal(utility::log2(x), 3.0));
211  float xf = x;
212  suite.add(suite.equal(utility::log2(xf), 3.0));
213  long double xld = 8.0;
214  suite.add(suite.equal(utility::log2(xld), 3.0));
215
216  test_basename(suite);
217  test_dirname(suite);
218  test_fnmatch(suite);
219  test_replace(suite);
220
221  return suite.return_value();
222}
223
224
225void test_errno_error_func(void)
226{
227  errno=1;
228  throw utility::errno_error("");
229}
230
231
232void test_errno_error(test::Suite& suite)
233{
234  suite.out() << "testing errno_error\n";
235  try {
236    test_errno_error_func();
237  }
238  catch ( utility::errno_error& e) {
239    suite.out() << "catching expected exception with what():\n";
240    suite.out() << e.what() << "\n";
241    return;
242  }
243  suite.err() << "error: expected thrown exception\n";
244  suite.add(false);
245}
246
247
248void test_compose_functors(test::Suite& suite)
249{
250  typedef utility::abs<double> U1;
251  // we use float here to check conversion from float to double in functors
252  typedef utility::Exp<float> U2;
253  typedef std::plus<double> B;
254  U1 u1;
255  U2 u2;
256  B b;
257
258  utility::compose_f_gx_hy<B, U1, U2> binary =
259    utility::make_compose_f_gx_hy(b, u1, u2);
260
261  utility::compose_f_gxy<U2, B> binary2 =
262    utility::make_compose_f_gxy(u2, b);
263
264  utility::compose_f_gx<U1, U2> unary =
265    utility::make_compose_f_gx(u1, u2);
266}
267
268
269void test_dereferencer(test::Suite& suite)
270{
271  using utility::Dereferencer;
272  Dereferencer<double*> deref;
273  double x = 0;
274  double* px = &x;
275  deref(px) = 1.662;
276  if (!suite.add(x==1.662)) {
277    suite.out() << "test_dereferencer failed: x: " << x << "\n";
278  }
279}
280
281
282void test_get_map(test::Suite& suite)
283{
284  std::map<std::string, int> m;
285  m["one"] = 1;
286  m["two"] = 2;
287  const std::map<std::string, int> m2(m);
288  int i = utility::get(m2, static_cast<std::string>("one"));
289  i = 0; // avoid compiler warning
290  try {
291    utility::get(m2, static_cast<std::string>("three"));
292    suite.add(false);
293  }
294  catch (std::runtime_error& e) {
295    suite.out() << "expected exception thrown with what() = "
296                << e.what() << "\n";
297  }
298}
299
300void test_less_nan(test::Suite& suite)
301{
302  suite.out() << "test less_nan\n";
303  bool prior_ok=suite.ok();
304  utility::less_nan<double> f;
305  suite.add(f(2.7,3.14));
306  suite.add(!f(2.7,-3.14));
307  suite.add(!f(2.7,2.7));
308  suite.add(f(3.14, std::numeric_limits<double>::quiet_NaN()));
309  suite.add(!f(std::numeric_limits<double>::quiet_NaN(), 2.7));
310  suite.add(!f(std::numeric_limits<double>::quiet_NaN(),
311               std::numeric_limits<double>::quiet_NaN()));
312  utility::less_nan<utility::DataWeight> fw;
313  suite.add(fw(utility::DataWeight(2), utility::DataWeight(3)));
314
315  bool post_ok=suite.ok();
316  if (!post_ok && prior_ok)
317    suite.err() << "test_less_nan failed\n";
318}
319
320
321void test_load(test::Suite& suite)
322{
323  std::istringstream iss(std::string("1.69 3.14"));
324  std::vector<double> double_vec;
325  utility::load(iss, double_vec);
326  suite.add(double_vec.size()==2);
327
328  std::istringstream iss2(std::string("1 2"));
329  std::vector<double> uint_vec;
330  utility::load(iss2, uint_vec);
331  suite.add(uint_vec.size()==2);
332
333  std::istringstream iss3(std::string("1.69 3.14"));
334  std::vector<std::string> str_vec;
335  utility::load(iss3, str_vec);
336  if (str_vec.size()==2) {
337    if (str_vec[1]!="3.14") {
338      suite.out() << "error: load<string>:" << "\n"
339                  << "  str_vec[1] = " << str_vec[1] << "\n"
340                  << "  expected:  3.14n";
341      suite.add(false);
342    }
343  }
344  else {
345    suite.out() << "str_vec.size() is not 2\n";
346    suite.add(false);
347  }
348}
349
350
351void test_load2(test::Suite& suite)
352{
353  suite.out() << "test load<double>(matrix)\n";
354  std::string str("1.69 3.14\n1.23 1.22\n");
355  std::istringstream iss(str);
356  std::vector<std::vector<double> > double_m;
357  utility::load(iss, double_m, '\0', '\n');
358  suite.add(double_m.size()==2);
359
360  suite.out() << "test load<std::string>(matrix)\n";
361  suite.out() << str;
362  std::istringstream iss2(str);
363  std::vector<std::vector<std::string> > str_m;
364  utility::load(iss2, str_m, '\0', '\n');
365  std::cout << str_m.size() << std::endl;
366  for (size_t i=0; i<str_m.size(); ++i)
367    for (size_t j=0; j<str_m[i].size(); ++j)
368      suite.out() << str_m[i][j] << "\n";
369  suite.add(str_m.size()==2);
370}
371
372
373void test_inverse(test::Suite& suite)
374{
375  suite.err() << "Testing inverse\n";
376  std::vector<std::string> vec;
377  vec.push_back("No");
378  vec.push_back("one");
379  vec.push_back("shall");
380  vec.push_back("be");
381  vec.push_back("subjected");
382  std::map<std::string, std::vector<size_t> > map;
383  utility::inverse(vec.begin(), vec.end(), map);
384  test_inverse_validate(vec.begin(), vec.end(), map, suite);
385  utility::inverse(vec.begin()+1, vec.end(), map);
386  test_inverse_validate(vec.begin()+1, vec.end(), map, suite);
387  const std::vector<std::string> vec2(vec);
388  utility::inverse(vec2.begin(), vec2.end(), map);
389  test_inverse_validate(vec2.begin(), vec2.end(), map, suite);
390  std::multimap<std::string, size_t> mmap;
391  utility::inverse(vec.begin(), vec.end(), mmap);
392  test_inverse_validate(vec.begin(), mmap, suite);
393  utility::inverse(vec.begin()+1, vec.end(), mmap);
394  test_inverse_validate(vec.begin()+1, mmap, suite);
395  // do not run compile tests
396  if (false) {
397    std::map<std::string, std::vector<size_t> > m;
398    boost::input_iterator_archetype<std::string> start;
399    boost::input_iterator_archetype<std::string> end;
400    utility::inverse(start, end, m);
401    std::map<std::string, std::vector<size_t>, std::greater<std::string> > m2;
402    utility::inverse(start, end, m2);
403    std::multimap<std::string, size_t> m3;
404    utility::inverse(start, end, m3);
405    std::multimap<std::string, size_t, std::greater<std::string> > m4;
406    utility::inverse(start, end, m4);
407
408    std::map<std::string, size_t> unique_map;
409    utility::inverse(start, end, unique_map);
410  }
411}
412
413
414template<typename InputIterator, typename Key>
415void test_inverse_validate(InputIterator first, InputIterator last,
416                           const std::map<Key, std::vector<size_t> >& m,
417                           test::Suite& suite)
418{
419  typedef typename std::map<Key, std::vector<size_t> >::const_iterator Iter;
420  for ( InputIterator iter=first; iter != last; ++iter) {
421    Iter map_iter = m.find(*iter);
422    if (!suite.add(map_iter!=m.end())) {
423      suite.err() << "test_inverse_validate() failed\n";
424      suite.err() << "  could not find: " << *first << " in map m\n";
425    }
426    for (size_t i=0; i<map_iter->second.size(); ++i)
427      if (!suite.add(map_iter->second[i] ==
428                     static_cast<size_t>(distance(first, iter))) ) {
429        suite.err() << "test_inverse_validate() failed\n";
430        suite.err() << "  comparing: " << map_iter->second[i]
431                    << " and " << distance(first, iter)
432                    << " expected them to be equal\n";
433      }
434  }
435}
436
437
438template<typename InputIterator, typename Key>
439void test_inverse_validate(InputIterator first,
440                           const std::multimap<Key, size_t>& m,
441                           test::Suite& suite)
442{
443  for (typename std::multimap<Key, size_t>::const_iterator iter=m.begin();
444       iter!=m.end(); ++iter) {
445    suite.add(*(first + iter->second) == iter->first);
446  }
447
448}
449
450void test_sort_index(test::Suite& suite)
451{
452  suite.err() << "testing sort_index" << std::endl;
453  utility::Vector a(10);
454  for (size_t i=0; i<a.size(); ++i)
455    a(i) = std::pow(i-4.2,2);
456  std::vector<size_t> vec;
457  utility::sort_index(vec, a);
458
459  std::vector<size_t> vec2;
460  utility::sort_index(a.begin(), a.end(), vec2);
461  if (vec.size()==vec2.size()) {
462    if (!suite.equal_range(vec.begin(), vec.end(), vec2.begin())) {
463      suite.add(false);
464    }
465  }
466  else {
467    suite.add(false);
468    suite.err() << "size mismatch: vec.size()=" << vec.size()
469                << "vec2.size()=" << vec2.size() << "\n";
470  }
471  const utility::VectorConstView b(a, 0, 5, 2);
472
473  std::vector<size_t> vec3;
474  utility::sort_index(vec3, b);
475  std::vector<size_t> vec4;
476  utility::sort_index(b.begin(), b.end(), vec4);
477  if (vec3.size()!=vec4.size()) {
478    suite.add(false);
479    suite.err() << "size mismatch: vec3.size()=" << vec3.size()
480                << " vec4.size()=" << vec4.size() << "\n";
481  }
482  else {
483    if (!suite.equal_range(vec3.begin(), vec3.end(), vec4.begin())){
484      suite.add(false);
485    }
486  }
487
488  std::list<double> list;
489  for (size_t i=0; i<a.size(); ++i)
490    list.push_back(a(i));
491  std::vector<size_t> vec5;
492  utility::sort_index(list.begin(), list.end(), vec5);
493  if (vec.size()!=vec5.size()) {
494    suite.add(false);
495    suite.err() << "size mismatch: vec.size()=" << vec.size()
496                << " vec5.size()=" << vec5.size() << "\n";
497  }
498  else {
499    if (!suite.equal_range(vec.begin(), vec.end(), vec5.begin())){
500      suite.add(false);
501    }
502  }
503
504  // do not run compiler tests
505  if (false) {
506    std::vector<size_t> vec;
507    utility::sort_index(boost::forward_iterator_archetype<double>(),
508                        boost::forward_iterator_archetype<double>(),
509                        vec);
510  }
511}
512
513void test_ptr_compare(test::Suite& suite)
514{
515  using utility::make_ptr_compare;
516  std::vector<double*> vec(10);
517  std::greater<double> d_greater;
518  for (size_t i=0; i<vec.size(); ++i)
519    vec[i] = new double(i);
520  std::sort(vec.begin(), vec.end(), make_ptr_compare(*vec.begin(), d_greater));
521  for (size_t i=1; i<vec.size(); ++i) {
522    suite.add( *vec[i-1] > *vec[i]);
523  }
524  std::sort(vec.begin(), vec.end(), make_ptr_compare(*vec.begin()));
525  for (size_t i=1; i<vec.size(); ++i) {
526    suite.add( *vec[i-1] < *vec[i]);
527  }
528  for (size_t i=0; i<vec.size(); ++i)
529    delete vec[i];
530
531}
532
533
534void test_fnmatch(test::Suite& suite)
535{
536  bool ok=true;
537  ok &= test_fnmatch(true,"peter", "peter");
538  ok &= test_fnmatch(false,"peter", "peterj");
539
540  ok &= test_fnmatch(true,"*", "peterj");
541  ok &= test_fnmatch(true,"p*", "peterj");
542  ok &= test_fnmatch(true, "p*", "peter");
543  ok &= test_fnmatch(false, "p*j", "peter");
544  ok &= test_fnmatch(true, "p*j", "peterj");
545  ok &= test_fnmatch(true, "*peter", "peter");
546
547  ok &= test_fnmatch(true, "p?ter", "peter");
548  ok &= test_fnmatch(false, "p?er", "peter");
549  ok &= test_fnmatch(false, "p?eter", "peter");
550
551  ok &= test_fnmatch(true, "filename", "filename");
552  ok &= test_fnmatch(true, "*name", "filename");
553  ok &= test_fnmatch(true, "[fa]il?name", "filename");
554  ok &= test_fnmatch(true, "[fa]*il?name", "ffilename");
555
556  ok &= test_fnmatch(true, "[fa]*il?name", "fafafailename");
557  ok &= test_fnmatch(false, "[fa]?il?name", "ilename");
558  ok &= test_fnmatch(false, "?[fa]il?name", "ilename");
559  ok &= test_fnmatch(true, "[fa]il?name", "filename");
560  ok &= test_fnmatch(false, "[fa]?il?name", "fafafailename");
561
562  ok &= test_fnmatch(true, "*name", "/path/to/filename");
563  ok &= test_fnmatch(true, "*name", "file.name");
564  ok &= test_fnmatch(true, "*.txt", ".file.txt");
565  suite.add(ok);
566}
567
568
569void test_dirname(test::Suite& suite)
570{
571  std::vector<std::string> path;
572  path.push_back("/usr/lib");
573  path.push_back("/usr/");
574  path.push_back("usr");
575  path.push_back("/");
576  path.push_back(".");
577  path.push_back("..");
578  std::vector<std::string> dir;
579  dir.push_back("/usr");
580  dir.push_back("/");
581  dir.push_back(".");
582  dir.push_back("/");
583  dir.push_back(".");
584  dir.push_back(".");
585  assert(dir.size()==path.size());
586  using utility::dirname;
587  for (size_t i=0; i<dir.size(); ++i) {
588    if (dir[i] != dirname(path[i])) {
589      suite.add(false);
590      suite.out() << "error:\n";
591      suite.out() << "path:           " << path[i] << "\n";
592      suite.out() << "directory_name: " << dirname(path[i]) << "\n";
593      suite.out() << "expected:       " << dir[i] << "\n";
594    }
595  }
596}
597
598
599void test_basename(test::Suite& suite)
600{
601  std::vector<std::string> path;
602  path.push_back("/usr/lib");
603  path.push_back("/usr/");
604  path.push_back("usr");
605  path.push_back("/");
606  path.push_back(".");
607  path.push_back("..");
608  std::vector<std::string> file;
609  file.push_back("lib");
610  file.push_back("usr");
611  file.push_back("usr");
612  file.push_back("/");
613  file.push_back(".");
614  file.push_back("..");
615  assert(file.size()==path.size());
616  using utility::basename;
617  for (size_t i=0; i<file.size(); ++i) {
618    if (file[i] != basename(path[i])) {
619      suite.add(false);
620      suite.out() << "error:\n";
621      suite.out() << "path:           " << path[i] << "\n";
622      suite.out() << "basename:       " << basename(path[i]) << "\n";
623      suite.out() << "expected:       " << file[i] << "\n";
624    }
625  }
626}
627
628
629bool test_fnmatch(bool answ, std::string a, std::string b)
630{
631  using namespace utility;
632  bool res = fnmatch(a, b);
633  // check that fnmatch and regexp agree
634  if (res == answ)
635    return true;
636  if (res!=answ)
637    std::cerr << "fnmatch(" << a << ", " << b << ") results "
638              << res
639              << ". Expects " << answ << std::endl;
640  return false;
641}
642
643
644void test_replace(test::Suite& suite)
645{
646  std::string s = "Some magic string!!!";
647  utility::replace(s, "magic", "arbitrary");
648  if (!suite.add(s=="Some arbitrary string!!!"))
649    suite.err() << "error: " << s << "\n";
650}
Note: See TracBrowser for help on using the repository browser.