hmmer-3.3/easel/esl_msaweight.c
/* As above, but for BLOSUM weights with fractional identity threshold <maxid> */
As above, but for BLOSUM weights with fractional identity threshold 
As above but for BLOSUM weights with fractional identity threshold
static void
do_BLOSUM(ESL_ALPHABET *abc, ESL_MSA *msa, double maxid, double *expect, char *msg)
{
  if (esl_msaweight_BLOSUM(msa, maxid)                     != eslOK) esl_fatal(msg);
  if (esl_vec_DCompare(msa->wgt, expect, msa->nseq, 0.001) != eslOK) esl_fatal(msg);
  if (abc)
    {
      if (esl_msa_Digitize(abc, msa, NULL)                     != eslOK) esl_fatal(msg);
      if (esl_msaweight_BLOSUM(msa, maxid)                     != eslOK) esl_fatal(msg);
      if (esl_vec_DCompare(msa->wgt, expect, msa->nseq, 0.001) != eslOK) esl_fatal(msg);
      if (esl_msa_Textize(msa)                                 != eslOK) esl_fatal(msg);
    }
}
static void do_BLOSUM ( ESL_ALPHABET * abc , ESL_MSA * msa , double maxid , double * expect , char * msg ) { if ( esl_msaweight_BLOSUM ( msa , maxid ) != eslOK ) esl_fatal ( msg ) ; if ( esl_vec_DCompare ( msa -> wgt , expect , msa -> nseq , 0.001 ) != eslOK ) esl_fatal ( msg ) ; if ( abc ) { if ( esl_msa_Digitize ( abc , msa , NULL ) != eslOK ) esl_fatal ( msg ) ; if ( esl_msaweight_BLOSUM ( msa , maxid ) != eslOK ) esl_fatal ( msg ) ; if ( esl_vec_DCompare ( msa -> wgt , expect , msa -> nseq , 0.001 ) != eslOK ) esl_fatal ( msg ) ; if ( esl_msa_Textize ( msa ) != eslOK ) esl_fatal ( msg ) ; } }
As above, but for BLOSUM weights with fractional identity threshold <maxid>
As above but for BLOSUM weights with fractional identity threshold <maxid >