| 3117 | } |
| 3118 | |
| 3119 | void parameter_set4final(const double len, double &D0_MIN, double &Lnorm, |
| 3120 | double &d0, double &d0_search, const int mol_type) |
| 3121 | { |
| 3122 | if (mol_type>0) // RNA |
| 3123 | { |
| 3124 | parameter_set4final_C3prime(len, D0_MIN, Lnorm, |
| 3125 | d0, d0_search); |
| 3126 | return; |
| 3127 | } |
| 3128 | D0_MIN=0.5; |
| 3129 | |
| 3130 | Lnorm=len; //normalize TMscore by this in searching |
| 3131 | if (Lnorm<=21) d0=0.5; |
| 3132 | else d0=(1.24*pow((Lnorm*1.0-15), 1.0/3)-1.8); |
| 3133 | if (d0<D0_MIN) d0=D0_MIN; |
| 3134 | d0_search=d0; |
| 3135 | if (d0_search>8) d0_search=8; |
| 3136 | if (d0_search<4.5) d0_search=4.5; |
| 3137 | } |
| 3138 | |
| 3139 | void parameter_set4scale(const int len, const double d_s, double &Lnorm, |
| 3140 | double &d0, double &d0_search) |
no test coverage detected