int main(int argc, const char *argv[]) {
OptArgs opts;
opts.ParseCmdLine(argc, argv);
string queryFile, goldFile;
double epsilon;
bool help = false;
bool version = false;
int allowedWrong = 0;
double maxAbsVal = 0;
double minCorrelation = 1;
opts.GetOption(queryFile, "", 'q', "query-wells");
opts.GetOption(goldFile, "", 'g', "gold-wells");
opts.GetOption(epsilon, "0.0", 'e', "epsilon");
opts.GetOption(allowedWrong, "0", 'm', "max-mismatch");
opts.GetOption(minCorrelation, "1", 'c', "min-cor");
opts.GetOption(maxAbsVal, "1e3", '-', "max-val");
opts.GetOption(help, "false", 'h', "help");
opts.GetOption(version, "false", 'v', "version");
opts.CheckNoLeftovers();
if (version) {
fprintf (stdout, "%s", IonVersion::GetFullVersion("RawWellsEquivalent").c_str());
exit(0);
}
if (queryFile.empty() || goldFile.empty() || help) {
cout << "RawWellsEquivalent - Check to see how similar two wells files are to each other" << endl
<< "options: " << endl
<< " -g,--gold-wells trusted wells to compare against." << endl
<< " -q,--query-wells new wells to check." << endl
<< " -e,--epsilon maximum allowed difference to be considered equivalent." << endl
<< " -m,--max-mixmatch maximum number of non-equivalent entries to allow." << endl
<< " -c,--min-cor minimum correlation allowed to be considered equivalent." << endl
<< " --max-val maximum absolute value considered (avoid extreme values)." << endl
<< " -h,--help this message." << endl
<< "" << endl
<< "usage: " << endl
<< " RawWellsEquivalent -e 10 --query-wells query.wells --gold-wells gold.wells " << endl;
exit(1);
}
NumericalComparison<double> compare = CompareWells(queryFile, goldFile, epsilon, maxAbsVal);
cout << compare.GetCount() << " total values. " << endl
<< compare.GetNumSame() << " (" << (100.0 * compare.GetNumSame())/compare.GetCount() << "%) are equivalent. " << endl
<< compare.GetNumDiff() << " (" << (100.0 * compare.GetNumDiff())/compare.GetCount() << "%) are not equivalent. " << endl
<< "Correlation of: " << compare.GetCorrelation() << endl;
if((compare.GetCount() - allowedWrong) >= compare.GetNumSame() ||
compare.GetCorrelation() < minCorrelation) {
cout << "Wells files not equivalent for allowed mismatch: " << allowedWrong
<< " minimum correlation: " << minCorrelation << endl;
return 1;
}
cout << "Wells files equivalent for allowed mismatch: " << allowedWrong
<< " minimum correlation: " << minCorrelation << endl;
return 0;
}
int main(int argc, const char *argv[]) {
OptArgs opts;
opts.ParseCmdLine(argc, argv);
string queryFile, goldFile;
double epsilon;
bool help = false;
bool version = false;
int allowedWrong = 0;
double maxAbsVal = 0;
double minCorrelation = 1;
bool dumpMisMatch = false;
opts.GetOption(queryFile, "", 'q', "query-wells");
opts.GetOption(goldFile, "", 'g', "gold-wells");
opts.GetOption(epsilon, "0.0", 'e', "epsilon");
opts.GetOption(allowedWrong, "0", 'm', "max-mismatch");
opts.GetOption(minCorrelation, "1", 'c', "min-cor");
opts.GetOption(maxAbsVal, "1e3", '-', "max-val");
opts.GetOption(help, "false", 'h', "help");
opts.GetOption(version, "false", 'v', "version");
opts.GetOption(dumpMisMatch, "false", 'o', "dump-mismatch");
opts.CheckNoLeftovers();
if (version) {
fprintf (stdout, "%s", IonVersion::GetFullVersion("RawWellsEquivalent").c_str());
exit(0);
}
if (queryFile.empty() || goldFile.empty() || help) {
printUsage();
exit(1);
}
DumpMismatches dump(dumpMisMatch);
NumericalComparison<double> compare = CompareWells(queryFile, goldFile, epsilon, maxAbsVal, dump);
cout << compare.GetCount() << " total values. " << endl
<< compare.GetNumSame() << " (" << (100.0 * compare.GetNumSame())/compare.GetCount() << "%) are equivalent. " << endl
<< compare.GetNumDiff() << " (" << (100.0 * compare.GetNumDiff())/compare.GetCount() << "%) are not equivalent. " << endl
<< "Correlation of: " << compare.GetCorrelation() << endl;
if((compare.GetCount() - allowedWrong) > compare.GetNumSame() ||
(compare.GetCorrelation() < minCorrelation && compare.GetCount() != compare.GetNumSame())) {
cout << "Wells files not equivalent for allowed mismatch: " << allowedWrong
<< " minimum correlation: " << minCorrelation << endl;
return 1;
}
cout << "Wells files equivalent for allowed mismatch: " << allowedWrong
<< " minimum correlation: " << minCorrelation << endl;
return 0;
}
