htslib 0.4.1 → 0.4.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 86838244ad390e04fb293d3db9fd639c837dcd7b2309cc0c83da0b5925e7185c
-  data.tar.gz: 547063279b934706efcf919d287217f4231fe76d2e8fb3ae6df8e4dd7b6040ad
+  metadata.gz: 69e910e95cfe1bef250416f8bbe480a39a1bc54b53103435cf15811860f3b92d
+  data.tar.gz: 060c7e5af8127c496af3b5d5bec5eca5e7476cb34ccc86a75546fd1a80d6dd8a
 SHA512:
-  metadata.gz: 787f9167abef88ddbd758d3155f1f3568bb99d8a5a340086a1cc677be46e0e77ec5a0ef56d201610f06143d0f54a20a1f2563ae3e080996564b5b5545e3e8a2a
-  data.tar.gz: d07ae5e437c9729a39e32adb12ff694560eaa33e44094679199bc8b742ccaab4c5a111c7e6391a314a8b6a19fae82acea549123b57497690678e1b2b1746f2c7
+  metadata.gz: 62adc672a3606265d231efa7f568232d55279634c8fb650190aee01c750b0dbac2421a1cf8fb2507bf5b21e16afbe4e4425354cd3dd2b1398d4df7b3f24e6655
+  data.tar.gz: 2a6bd47a73ff2695a515e934fa1919d8fcdf1108e10178404af31da769b4e7c40a8580f2382ab8f7e6ac8697261344fb32ecebc138807d5cd03724745457f4d4

data/lib/hts/bcf/header.rb

@@ -104,7 +104,6 @@ module HTS
         subset_samples = normalize_subset_samples(samples)
         validate_subset_samples!(subset_samples)
 
-        nil
         imap_pointer = nil
         if subset_samples.empty?
           subset_hdr = LibHTS.bcf_hdr_subset(@bcf_hdr, 0, ::FFI::Pointer::NULL, ::FFI::Pointer::NULL)

data/lib/hts/bcf.rb

@@ -322,7 +322,7 @@ module HTS
       record = Record.new(header, bcf1)
       begin
         loop do
-          slen = LibHTS.hts_itr_next(@hts_file[:fp][:bgzf], qiter, bcf1, ::FFI::Pointer::NULL)
+          slen = LibHTS.bcf_itr_next(@hts_file, qiter, bcf1)
           break if slen == -1
           raise if slen < -1
 
@@ -371,7 +371,7 @@ module HTS
       bcf1 = LibHTS.bcf_init
       record = Record.new(header, bcf1)
       loop do
-        slen = LibHTS.hts_itr_next(@hts_file[:fp][:bgzf], qiter, bcf1, ::FFI::Pointer::NULL)
+        slen = LibHTS.bcf_itr_next(@hts_file, qiter, bcf1)
         break if slen == -1
         raise if slen < -1
 

data/lib/hts/libhts/hfile.rb

@@ -64,6 +64,18 @@ module HTS
       [:string, :int, HFile],
       :string
 
+    # Read a line from a stream and append it to a kstring
+    attach_function \
+      :khgetline,
+      [KString, HFile],
+      :int
+
+    # Convenience function to call kgetline with a FILE *
+    attach_function \
+      :kfgetline,
+      [KString, :pointer],
+      :int
+
     # Peek at characters to be read without removing them from buffers
     attach_function \
       :hpeek,

data/lib/hts/libhts/tbx.rb

@@ -14,6 +14,21 @@ module HTS
       [Tbx, :string],
       :int
 
+    attach_function \
+      :tbx_itr_querys1,
+      [Tbx, :string],
+      HtsItr.by_ref
+
+    attach_function \
+      :tbx_itr_regarray,
+      [Tbx, :pointer, :uint],
+      HtsItr.by_ref
+
+    attach_function \
+      :tbx_itr_next1,
+      [HtsFile, Tbx, HtsItr, :pointer],
+      :int
+
     # Internal helper function used by tbx_itr_next()
     attach_function \
       :hts_get_bgzfp,

data/lib/hts/libhts/tbx_funcs.rb

@@ -16,11 +16,15 @@ module HTS
       end
 
       def tbx_itr_querys(tbx, s)
+        return tbx_itr_querys1(tbx, s) if respond_to?(:tbx_itr_querys1)
+
         hts_itr_querys(tbx[:idx], s, @@tbx_name2id, tbx, @ffi_functions[:hts_itr_query],
                        @ffi_functions[:tbx_readrec])
       end
 
       def tbx_itr_next(htsfp, tbx, itr, r)
+        return tbx_itr_next1(htsfp, tbx, itr, r) if respond_to?(:tbx_itr_next1)
+
         hts_itr_next(hts_get_bgzfp(htsfp), itr, r, tbx)
       end
 

data/lib/hts/libhts/vcf.rb

@@ -145,6 +145,16 @@ module HTS
       [BGZF, :pointer, :pointer, :pointer, :hts_pos_t, :hts_pos_t],
       :int
 
+    attach_function \
+      :bcf_itr_querys1,
+      [HtsIdx, BcfHdr, :string],
+      HtsItr.by_ref
+
+    attach_function \
+      :bcf_itr_regarray,
+      [HtsIdx, BcfHdr, :pointer, :uint],
+      HtsItr.by_ref
+
     # Write a line to a VCF file
     attach_function \
       :vcf_write_line,

data/lib/hts/libhts/vcf_funcs.rb

@@ -309,9 +309,22 @@ module HTS
       end
 
       def bcf_itr_querys(idx, hdr, s)
+        return bcf_itr_querys1(idx, hdr, s) if respond_to?(:bcf_itr_querys1)
+
         hts_itr_querys(idx, s, @@bcf_hdr_name2id, hdr, @ffi_functions[:hts_itr_query], @ffi_functions[:bcf_readrec])
       end
 
+      # Ruby implementation of htslib's static inline bcf_itr_next().
+      def bcf_itr_next(htsfp, itr, r)
+        return -2 unless htsfp[:is_bgzf] == 1
+
+        if itr[:multi] == 1
+          hts_itr_multi_next(htsfp, itr, r)
+        else
+          hts_itr_next(htsfp[:fp][:bgzf], itr, r, FFI::Pointer::NULL)
+        end
+      end
+
       # Load a BCF index
       def bcf_index_load(fn)
         hts_idx_load(fn, HTS_FMT_CSI)

data/lib/hts/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module HTS
-  VERSION = "0.4.1"
+  VERSION = "0.4.2"
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: htslib
 version: !ruby/object:Gem::Version
-  version: 0.4.1
+  version: 0.4.2
 platform: ruby
 authors:
 - kojix2