zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Cache pileup bases and mapping quals after first call to getBases() and getMappingQuals(), respectively. Subsequent calls to these method will return cached arrays. 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@2525 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
asivache 2010-01-06 18:05:00 +00:00
parent d6fb19bb67
commit a17d725c35
1 changed files with 13 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
java/src/org/broadinstitute/sting/utils/pileup/ReadBackedPileup.java
View File

@ -20,7 +20,9 @@ import net.sf.samtools.SAMRecord;
public class ReadBackedPileup implements Iterable<PileupElement> {
private GenomeLoc loc = null;
private ArrayList<PileupElement> pileup = null;
byte [] bases = null; // will be used to cache pileup bases after first call to getBases
byte [] mapquals = null; // will be used to cache pileup mapping quals after first call to getMappingQuals
private int size = 0; // cached value of the size of the pileup
private int nDeletions = 0; // cached value of the number of deletions
private int nMQ0Reads = 0; // cached value of the number of MQ0 reads
@ -395,9 +397,11 @@ public class ReadBackedPileup implements Iterable<PileupElement> {
* @return
*/
public byte[] getBases() {
byte[] v = new byte[size()];
for ( ExtendedPileupElement pile : this.extendedForeachIterator() ) { v[pile.getPileupOffset()] = pile.getBase(); }
return v;
if ( bases == null ) {
bases = new byte[size()];
for ( ExtendedPileupElement pile : this.extendedForeachIterator() ) { bases[pile.getPileupOffset()] = pile.getBase(); }
}
return bases;
}
/**
@ -425,9 +429,11 @@ public class ReadBackedPileup implements Iterable<PileupElement> {
* @return
*/
public byte[] getMappingQuals() {
byte[] v = new byte[size()];
for ( ExtendedPileupElement pile : this.extendedForeachIterator() ) { v[pile.getPileupOffset()] = (byte)pile.getRead().getMappingQuality(); }
return v;
if ( mapquals == null ) {
mapquals = new byte[size()];
for ( ExtendedPileupElement pile : this.extendedForeachIterator() ) { mapquals[pile.getPileupOffset()] = (byte)pile.getRead().getMappingQuality(); }
}
return mapquals;
}