1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
|
/*++
Module Name:
AlignerContext.h
Abstract:
Common parameters for running single & paired alignment.
Authors:
Ravi Pandya, May, 2012
Environment:
`
User mode service.
Revision History:
Integrated from SingleAligner.cpp & PairedAligner.cpp
--*/
#pragma once
#include "stdafx.h"
#include "Genome.h"
#include "RangeSplitter.h"
#include "AlignerOptions.h"
#include "AlignerStats.h"
#include "ParallelTask.h"
#include "GenomeIndex.h"
class AlignerExtension;
/*++
Common context state shared across threads during alignment process
--*/
class AlignerContext : public TaskContextBase
{
public:
AlignerContext(int i_argc, const char **i_argv, const char *i_version, AlignerExtension* i_extension = NULL);
~AlignerContext();
// running alignment
void runAlignment(int argc, const char **argv, const char *version, unsigned *nArgsConsumed);
// ParallelTask template
void initializeThread();
void runThread();
void finishThread(AlignerContext* common);
void printStatsHeader();
void printStats();
void beginIteration();
void finishIteration();
// advance to next iteration in range, return false when past end
bool nextIteration();
// overrideable by concrete single/paired alignment subclasses
// parse options from the command line
AlignerOptions* parseOptions(int argc, const char **argv, const char *version, unsigned *argsConsumed, bool paired);
// initialize from options
virtual bool initialize();
// new stats object
virtual AlignerStats* newStats() = 0;
// instantiate and run a parallel task
virtual void runTask() = 0;
// run single thread within single iteration
virtual void runIterationThread() = 0;
virtual void typeSpecificBeginIteration() = 0;
virtual void typeSpecificNextIteration() = 0;
virtual bool isPaired() = 0;
friend class AlignerContext2;
// common state across all threads
GenomeIndex *index;
ReadWriterSupplier *writerSupplier;
ReaderContext readerContext;
_int64 alignStart;
_int64 alignTime;
AlignerOptions *options;
AlignerStats *stats;
AlignerExtension *extension;
unsigned maxDist;
unsigned numSeedsFromCommandLine;
double seedCoverage;
unsigned minWeightToCheck;
int maxHits;
bool detailedStats;
ReadClippingType clipping;
unsigned extraSearchDepth;
int argc;
const char **argv;
const char *version;
FILE *perfFile;
bool noUkkonen;
bool noOrderedEvaluation;
bool noTruncation;
int maxSecondaryAlignmentAdditionalEditDistance;
int maxSecondaryAlignments;
int maxSecondaryAlignmentsPerContig;
unsigned minReadLength;
// iteration variables
int maxHits_;
int maxDist_;
// Per-thread context state used during alignment process
ReadWriter *readWriter;
};
// abstract class for extending base context
class AlignerExtension
{
public:
virtual ~AlignerExtension() {}
virtual AbstractOptions* extraOptions() { return NULL; }
virtual AbstractStats* extraStats() { return NULL; }
virtual bool skipAlignment() { return false; }
virtual void initialize() {}
virtual void beginIteration() {}
virtual AlignerExtension* copy() { return new AlignerExtension(); }
virtual void beginThread() {}
virtual bool runIterationThread(PairedReadSupplier* supplier, AlignerContext* threadContext) { return false; }
virtual bool runIterationThread(ReadSupplier* supplier, AlignerContext* threadContext) { return false; }
virtual void finishThread() {}
virtual void finishIteration() {}
virtual void printStats() {}
virtual void finishAlignment() {}
};
|
