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
|
/*
Copyright 2012 University of Washington
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at
http://www.apache.org/licenses/LICENSE-2.0
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
*/
#ifndef _COMETPREPROCESS_H_
#define _COMETPREPROCESS_H_
#include "Common.h"
#include "ThreadPool.h"
struct PreprocessThreadData
{
Spectrum mstSpectrum;
int iAnalysisType;
int iFileLastScan;
bool *pbMemoryPool; //MH: Manages active memory pool
PreprocessThreadData()
{
}
PreprocessThreadData(Spectrum &spec_in,
int iAnalysisType_in,
int iFileLastScan_in)
{
mstSpectrum = spec_in;
iAnalysisType = iAnalysisType_in;
iFileLastScan = iFileLastScan_in;
}
~PreprocessThreadData()
{
//MH: Mark that the memory is no longer in use.
//DO NOT FREE MEMORY HERE. Just release pointer.
Threading::LockMutex(g_preprocessMemoryPoolMutex);
if(pbMemoryPool!=NULL)
*pbMemoryPool=false;
pbMemoryPool=NULL;
Threading::UnlockMutex(g_preprocessMemoryPoolMutex);
}
void SetMemory(bool *pbMemoryPool_in)
{
pbMemoryPool = pbMemoryPool_in;
}
};
class CometPreprocess
{
public:
CometPreprocess();
~CometPreprocess();
static void Reset();
static bool LoadAndPreprocessSpectra(MSReader &mstReader,
int iFirstScan,
int iLastScan,
int iAnalysisType,
int minNumThreads,
int maxNumThreads);
static void PreprocessThreadProc(PreprocessThreadData *pPreprocessThreadData);
static bool DoneProcessingAllSpectra();
static bool AllocateMemory(int maxNumThreads);
static bool DeallocateMemory(int maxNumThreads);
private:
// Private static methods
static bool PreprocessSpectrum(Spectrum &spec,
double *pdTmpRawData,
double *pdTmpFastXcorrData,
double *pdTmpCorrelationData,
double *pdTmpSmoothedSpectrum,
double *pdTmpPeakExtracted);
static bool CheckExistOutFile(int iCharge,
int iScanNum);
static bool AdjustMassTol(struct Query *pScoring);
static void PreloadIons(MSReader &mstReader,
Spectrum &spec,
bool bNext=false,
int scNum=0);
static bool CheckActivationMethodFilter(MSActivation act);
static bool CheckExit(int iAnalysisType,
int iScanNum,
int iTotalScans,
int iLastScan,
int iReaderLastScan,
int iNumSpectraLoaded);
static bool Preprocess(struct Query *pScoring,
Spectrum mstSpectrum,
double *pdTmpRawData,
double *pdTmpFastXcorrData,
double *pdTmpCorrelationData,
double *pdSmoothedSpectrum,
double *pdTmpPeakExtracted);
static bool LoadIons(struct Query *pScoring,
double *pdTmpRawData,
Spectrum mstSpectrum,
struct PreprocessStruct *pPre);
static void MakeCorrData(double *pdTmpRawData,
double *pdTmpCorrelationData,
struct Query *pScoring,
struct PreprocessStruct *pPre);
static bool Smooth(double *data,
int iArraySize,
double *pdSmoothedSpectrum);
static bool PeakExtract(double *data,
int iArraySize,
double *pdTmpPeakExtracted);
static void GetTopIons(double *pdTmpRawData,
struct msdata *pTmpSpData,
int iArraySize);
static int QsortByIon(const void *p0,
const void *p1);
static void StairStep(struct msdata *pTmpSpData);
static bool IsValidInputType(int inputType);
// Private member variables
static Mutex _maxChargeMutex;
static bool _bFirstScan;
static bool _bDoneProcessingAllSpectra;
//MH: Common memory to be shared by all threads during spectral processing
static bool *pbMemoryPool; //MH: Regulator of memory use
static double **ppdTmpRawDataArr; //MH: Number of arrays equals threads
static double **ppdTmpFastXcorrDataArr; //MH: Ditto
static double **ppdTmpCorrelationDataArr; //MH: Ditto
static double **ppdTmpSmoothedSpectrumArr; //MH: Ditto
static double **ppdTmpPeakExtractedArr; //MH: Ditto
};
#endif // _COMETPREPROCESS_H_
|
