hmmer-3.3/src/eweight.c
/* Evaluate fx = rel entropy - etarget, which we want to be = 0,
 * for effective sequence number <x>.
 */
Evaluate fx = rel entropy - etarget, which we want to be = 0,
for effective sequence number .
Evaluate fx = rel entropy - etarget which we want to be = 0 for effective sequence number .
static int
eweight_target_f(double Neff, void *params, double *ret_fx)
{
  struct ew_param_s *p = (struct ew_param_s *) params;
  p7_hmm_CopyParameters(p->hmm, p->h2);
  p7_hmm_Scale(p->h2, Neff / (double) p->h2->nseq);
  p7_ParameterEstimation(p->h2, p->pri);
  *ret_fx = p7_MeanMatchRelativeEntropy(p->h2, p->bg) - p->etarget;
  return eslOK;
}
static int eweight_target_f ( double Neff , void * params , double * ret_fx ) { struct ew_param_s * p = ( struct ew_param_s * ) params ; p7_hmm_CopyParameters ( p -> hmm , p -> h2 ) ; p7_hmm_Scale ( p -> h2 , Neff / ( double ) p -> h2 -> nseq ) ; p7_ParameterEstimation ( p -> h2 , p -> pri ) ; * ret_fx = p7_MeanMatchRelativeEntropy ( p -> h2 , p -> bg ) - p -> etarget ; return eslOK ; }
Evaluate fx = rel entropy - etarget, which we want to be = 0,
 for effective sequence number <x>.
Evaluate fx = rel entropy - etarget which we want to be = 0 for effective sequence number <x > .