| 3077 | */ |
| 3078 | |
| 3079 | void parameter_set4search(const int xlen, const int ylen, |
| 3080 | double &D0_MIN, double &Lnorm, |
| 3081 | double &score_d8, double &d0, double &d0_search, double &dcu0) |
| 3082 | { |
| 3083 | //parameter initialization for searching: D0_MIN, Lnorm, d0, d0_search, score_d8 |
| 3084 | D0_MIN=0.5; |
| 3085 | dcu0=4.25; //update 3.85-->4.25 |
| 3086 | |
| 3087 | Lnorm=getmin(xlen, ylen); //normalize TMscore by this in searching |
| 3088 | if (Lnorm<=19) //update 15-->19 |
| 3089 | d0=0.168; //update 0.5-->0.168 |
| 3090 | else d0=(1.24*pow((Lnorm*1.0-15), 1.0/3)-1.8); |
| 3091 | D0_MIN=d0+0.8; //this should be moved to above |
| 3092 | d0=D0_MIN; //update: best for search |
| 3093 | |
| 3094 | d0_search=d0; |
| 3095 | if (d0_search>8) d0_search=8; |
| 3096 | if (d0_search<4.5) d0_search=4.5; |
| 3097 | |
| 3098 | score_d8=1.5*pow(Lnorm*1.0, 0.3)+3.5; //remove pairs with dis>d8 during search & final |
| 3099 | } |
| 3100 | |
| 3101 | void parameter_set4final_C3prime(const double len, double &D0_MIN, |
| 3102 | double &Lnorm, double &d0, double &d0_search) |
no test coverage detected