1 #include "Quantifications.h"
2 
Quantifications(uint32 nGeIn)3 Quantifications::Quantifications (uint32 nGeIn) {
4 
5     geneCounts.nType=3;
6     geneCounts.cAmbig = new uintQ[geneCounts.nType];
7     geneCounts.cNone = new uintQ[geneCounts.nType];
8 
9     geneCounts.nGe=nGeIn;
10     geneCounts.gCount = new uintQ* [geneCounts.nType];
11 
12     geneCounts.cMulti=0;
13     for (int itype=0; itype<geneCounts.nType; itype++)
14     {
15         geneCounts.cAmbig[itype]=0;
16         geneCounts.cNone[itype]=0;
17         geneCounts.gCount[itype] = new uintQ [geneCounts.nGe];
18         for (uint32 ii=0; ii<geneCounts.nGe; ii++)
19         {
20             geneCounts.gCount[itype][ii]=0;
21         };
22     };
23 };
24 
addQuants(const Quantifications & quantsIn)25 void Quantifications::addQuants(const Quantifications & quantsIn)
26 {
27     geneCounts.cMulti += quantsIn.geneCounts.cMulti;
28     for (int itype=0; itype<geneCounts.nType; itype++)
29     {
30         geneCounts.cAmbig[itype] += quantsIn.geneCounts.cAmbig[itype];
31         geneCounts.cNone[itype] += quantsIn.geneCounts.cNone[itype];
32         for (uint32 ii=0; ii<geneCounts.nGe; ii++)
33         {
34             geneCounts.gCount[itype][ii] += quantsIn.geneCounts.gCount[itype][ii];
35         };
36     };
37 };