Mercurial > hg > medcouple

comparison jmedcouple.c++ @ 28:65f6fad65f49

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
jmedcouple: get rid of unused trace_lev parameter
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Sun, 18 Jan 2015 20:41:04 -0500
parents 40cbfb64e952
children cf8c8c825c6d
comparison
equal deleted inserted replaced
27:40cbfb64e952 28:65f6fad65f49
2 #include <limits> 2 #include <limits>
3 #include <iostream> 3 #include <iostream>
4 #include <string> 4 #include <string>
5 #include <fstream> 5 #include <fstream>
6 #include <iomanip> 6 #include <iomanip>
7 #include <sstream>
8 #include <algorithm> 7 #include <algorithm>
9 #include <utility> 8 #include <utility>
10 9
11 using namespace std; 10 using namespace std;
12 11
66 return 0; 65 return 0;
67 } 66 }
68 67
69 double medcouple(const std::vector<double>& X, 68 double medcouple(const std::vector<double>& X,
70 double eps1 = numeric_limits<double>::epsilon(), 69 double eps1 = numeric_limits<double>::epsilon(),
71 double eps2 = numeric_limits<double>::min(), 70 double eps2 = numeric_limits<double>::min())
72 int trace_lev = 0)
73 { 71 {
74 size_t n = X.size(), n2 = (n - 1)/2; 72 size_t n = X.size(), n2 = (n - 1)/2;
75 73
76 if (n < 3) 74 if (n < 3)
77 return 0; 75 return 0;
224 else { 222 else {
225 cerr << "Please provide input file." << endl; 223 cerr << "Please provide input file." << endl;
226 return 1; 224 return 1;
227 } 225 }
228 226
229 double trace_lev = 0;
230 if(argc > 2) {
231 stringstream ss;
232 ss << argv[2];
233 ss >> trace_lev;
234 }
235
236 ifstream ifs(fname); 227 ifstream ifs(fname);
237 if (not ifs) { 228 if (not ifs) {
238 cerr << "File " << fname << " not found." << endl; 229 cerr << "File " << fname << " not found." << endl;
239 return 1; 230 return 1;
240 } 231 }
244 while (ifs >> datum) { 235 while (ifs >> datum) {
245 data.push_back(datum); 236 data.push_back(datum);
246 } 237 }
247 238
248 cout << setprecision(16) 239 cout << setprecision(16)
249 << medcouple(data, eps1, eps2, trace_lev) << endl; 240 << medcouple(data, eps1, eps2) << endl;
250 241
251 return 0; 242 return 0;
252 } 243 }
253 244