[Enhancement] The default value of maxsites of meme is confusing

[Yusuke Takahashi]

Dear All,

Under some conditions, "maxsites" is cut off at 1000.

This behavior, I think, is not written in the manual and is confusing.

I finally have found that the behavior seems to be intended and written at the line 898 of init.c of meme version 5.1.1.

    893   } else { // ANR/TCM model 

    894     if (min_nsites<2) min_nsites = 2; // default 

    895     if (!max_nsites) {

    896       if (objfun == Classic || objfun == NC) {

    897         //max_nsites = MIN(5*n_primary_samples, 50);

    898         max_nsites = MIN(5*n_primary_samples, classic_max_nsites);

    899       } else {

    900         max_nsites = 5*n_primary_samples;

    901       }

    902     }

 

I would like the behavior to be added to the manual.

In addition, I would like to know the reason of the behavior.