[med-svn] [psortb] 07/08: Adapt to hmmer2 interface
Andreas Tille
tille at debian.org
Tue Apr 18 20:55:53 UTC 2017
This is an automated email from the git hooks/post-receive script.
tille pushed a commit to branch master
in repository psortb.
commit b91dba4b157ec52eba705c32ee58e4395e2d72ce
Author: Andreas Tille <tille at debian.org>
Date: Wed Apr 12 22:24:03 2017 +0200
Adapt to hmmer2 interface
---
debian/patches/fix_hmmer_interface.patch | 24 +++++++++++++++++++++++-
1 file changed, 23 insertions(+), 1 deletion(-)
diff --git a/debian/patches/fix_hmmer_interface.patch b/debian/patches/fix_hmmer_interface.patch
index 6b65eb7..4dba247 100644
--- a/debian/patches/fix_hmmer_interface.patch
+++ b/debian/patches/fix_hmmer_interface.patch
@@ -9,7 +9,21 @@
if(! HMMFileRead(hmmfp, &tmp_model)) {
HMMFileClose(hmmfp);
-@@ -107,7 +107,7 @@ HMMReport *HMM::search(char *seq) {
+@@ -95,25 +95,26 @@ int HMM::load(char *filename) {
+
+ HMMReport *HMM::search(char *seq) {
+ struct tophit_s *ghit, *dhit;
++ struct dpmatrix_s *mx;
+ struct p7trace_s *trace;
+ double pvalue, evalue;
+ int numhits;
+ HMMReport *hmmrep;
+ float score;
+ int seqlen;
+- char *dseq;
++ unsigned char *dseq;
+
+ // Ensure we have a valid model and sequence.
if((model == NULL) ||(seq == NULL) || (!(seqlen = strlen(seq)))) return NULL;
// Convert the sequence to the format the HMMER library expects.
@@ -18,3 +32,11 @@
// Calculate the raw scores for the sequence.
if (P7ViterbiSize(seqlen, model->M) <= RAMLIMIT)
+- score = P7Viterbi(dseq, seqlen, model, &trace);
++ score = P7Viterbi(dseq, seqlen, model, mx, &trace);
+ else
+- score = P7SmallViterbi(dseq, seqlen, model, &trace);
++ score = P7SmallViterbi(dseq, seqlen, model, mx, &trace);
+
+
+ // Calculate the score using the forward algorithm if we need.
--
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-med/psortb.git
More information about the debian-med-commit
mailing list