/* unsafe function that will be replaced by a threadsafe companion in the future */
PRIVATE void init_alifold(int length){
unsigned int n;
if (length<1) nrerror("initialize_fold: argument must be greater 0");
if (init_length>0) free_alifold_arrays();
get_arrays((unsigned) length);
make_pair_matrix();
init_length=length;
for (n = 1; n <= (unsigned) length; n++)
indx[n] = (n*(n-1)) >> 1; /* n(n-1)/2 */
update_fold_params();
}
/* unsafe function that will be replaced by a threadsafe companion in the future */
PRIVATE void init_alifold(int length){
#ifdef _OPENMP
/* Explicitly turn off dynamic threads */
omp_set_dynamic(0);
#endif
if (length < 1) nrerror("initialize_fold: argument must be greater 0");
free_alifold_arrays();
get_arrays((unsigned) length);
init_length = length;
indx = get_indx((unsigned)length);
update_alifold_params();
}
