Skip to content

Commit

Permalink
Expose expectedCoverage parameter with option -xC; no canu support ye…
Browse files Browse the repository at this point in the history
…t. Replace atoi() and atof().
  • Loading branch information
brianwalenz committed Jun 28, 2021
1 parent 98f4cd2 commit e192966
Showing 1 changed file with 10 additions and 7 deletions.
17 changes: 10 additions & 7 deletions src/correction/generateCorrectionLayouts.C
Original file line number Diff line number Diff line change
Expand Up @@ -181,7 +181,6 @@ main(int argc, char **argv) {
bool doLogging = false;

uint32 expectedCoverage = 40; // How many overlaps per read to save, global filter
uint32 minEvidenceOverlap = 40;

uint32 iidMin = 1;
uint32 iidMax = UINT32_MAX;
Expand Down Expand Up @@ -212,27 +211,29 @@ main(int argc, char **argv) {


} else if (strcmp(argv[arg], "-b") == 0) { // READ SELECTION
iidMin = atoi(argv[++arg]);
iidMin = strtouint32(argv[++arg]);

} else if (strcmp(argv[arg], "-e") == 0) {
iidMax = atoi(argv[++arg]);
iidMax = strtouint32(argv[++arg]);

} else if (strcmp(argv[arg], "-eL") == 0) { // EVIDENCE SELECTION
minEvidenceLength = atoi(argv[++arg]);
minEvidenceLength = strtouint32(argv[++arg]);

} else if (strcmp(argv[arg], "-eE") == 0) {
maxEvidenceErate = atof(argv[++arg]);
maxEvidenceErate = strtodouble(argv[++arg]);

} else if (strcmp(argv[arg], "-eC") == 0) {
maxEvidenceCoverage = atof(argv[++arg]);
maxEvidenceCoverage = strtodouble(argv[++arg]);

} else if (strcmp(argv[arg], "-xC") == 0) {
expectedCoverage = strtouint32(argv[++arg]);

} else if (strcmp(argv[arg], "-V") == 0) {
doLogging = true;

} else if (strcmp(argv[arg], "-D") == 0) {
dumpScores = true;


} else {
fprintf(stderr, "ERROR: unknown option '%s'\n", argv[arg]);
err++;
Expand Down Expand Up @@ -268,6 +269,8 @@ main(int argc, char **argv) {
fprintf(stderr, " -eE erate maximum error rate of evidence overlaps\n");
fprintf(stderr, " -eC coverage maximum coverage of evidence reads to emit\n");
fprintf(stderr, "\n");
fprintf(stderr, " -xC coverage estimated coverage in input reads\n");
fprintf(stderr, "\n");

if (seqName == NULL)
fprintf(stderr, "ERROR: no input seqStore (-S) supplied.\n");
Expand Down

0 comments on commit e192966

Please # to comment.