Changeset 3802

Timestamp:

May 12, 2019, 6:31:06 AM (4 years ago)

Author:

Peter

Message:

closes #918. implement function Pileup::size

Location:

trunk

Files:

• test/pileup.cc (modified) (1 diff)
• yat/omic/Pileup.h (modified) (2 diffs)

trunk/test/pileup.cc

@@ -158 +158 @@
   Pileup<BamReadIterator>::const_reverse_iterator rend = pileup.rend();
   std::distance(pileup.rbegin(), rend);
+  pileup.size();
 }
 

trunk/yat/omic/Pileup.h

@@ -210 +210 @@
      */
     void increment(void);
+
+    /**
+       \return number of reads in pileup. Reads with a deletion of the
+       current position is also counted.
+
+       \since new in yat 0.17
+     */
+    size_t size(void) const;
 
     /**
@@ -380 +388 @@
 
   template<typename Iterator>
+  size_t Pileup<Iterator>::size(void) const
+  {
+    return data_.size();
+  }
+
+
+  template<typename Iterator>
   Pileup<Iterator>::Entry::Entry(const BamRead& b)
     : bam_(new BamRead(b)), qpos_(0), cigar_(bam_->cigar()), skip_(0)