// done with all structures along the way to deepest
int update_deepest(encoded &enc, int &deepest, short *min_pt, degen °, bool verbose)
{
// debug
/*if (deg.opt->f_point) {
fprintf(stderr, "UD: %s %.2f (%d, %d) (%d, %d)\n", pt_to_str(enc.pt).c_str(), deepest/100.0, enc.bp_left, enc.bp_right, enc.bp_left2, enc.bp_right2);
}*/
int tmp_en;
if (deg.opt->EOM) {
tmp_en = deg.current + energy_of_move_pt(enc.pt, enc.s0, enc.s1, enc.bp_left, enc.bp_right);
do_moves(enc, true, false);
if (enc.bp_left2 != 0) {
tmp_en += energy_of_move_pt(enc.pt, enc.s0, enc.s1, enc.bp_left2, enc.bp_right2);
do_moves(enc, false, true);
}
} else {
do_moves(enc);
tmp_en = energy_of_structure_pt(enc.seq, enc.pt, enc.s0, enc.s1, 0);
}
// our function
bool res = 0;
if (deg.opt->f_point) {
res = deg.opt->f_point(enc.pt, tmp_en);
undo_moves(enc);
return res;
}
// default behaviour
if (tmp_en <= deepest) {
deepest = tmp_en;
copy_arr(min_pt, enc.pt);
}
//printf("%s %7.2f\n", pt_to_str(enc.pt).c_str(), tmp_en/100.0);
undo_moves(enc);
return res;
}
PRIVATE int try_moves(intermediate_t c, int maxE, intermediate_t *next, int dist) {
int *loopidx, len, num_next=0, en, oldE;
move_t *mv;
short *pt;
len = c.pt[0];
loopidx = pair_table_to_loop_index(c.pt);
oldE = c.Sen;
for (mv=c.moves; mv->i!=0; mv++) {
int i,j;
if (mv->when>0) continue;
i = mv->i; j = mv->j;
pt = (short *) space(sizeof(short)*(len+1));
memcpy(pt, c.pt,(len+1)*sizeof(short));
if (j<0) { /*it's a delete move */
pt[-i]=0;
pt[-j]=0;
} else { /* insert move */
if ((loopidx[i] == loopidx[j]) && /* i and j belong to same loop */
(pt[i] == 0) && (pt[j]==0) /* ... and are unpaired */
) {
pt[i]=j;
pt[j]=i;
} else {
free(pt);
continue; /* llegal move, try next; */
}
}
#ifdef LOOP_EN
en = c.curr_en + energy_of_move_pt(c.pt, S, S1, i, j);
#else
en = energy_of_structure_pt(seq, pt, S, S1, 0);
#endif
if (en<maxE) {
next[num_next].Sen = (en>oldE)?en:oldE;
next[num_next].curr_en = en;
next[num_next].pt = pt;
mv->when=dist;
mv->E = en;
next[num_next++].moves = copy_moves(c.moves);
mv->when=0;
}
else free(pt);
}
free(loopidx);
return num_next;
}
/* done with all structures along the way to deepest*/
PRIVATE int
update_deepest(Encoded *Enc, struct_en *str, struct_en *min){
/* apply move + get its energy*/
int tmp_en;
tmp_en = str->energy + energy_of_move_pt(str->structure, Enc->s0, Enc->s1, Enc->bp_left, Enc->bp_right);
do_move(str->structure, Enc->bp_left, Enc->bp_right);
if (Enc->bp_left2 != 0) {
tmp_en += energy_of_move_pt(str->structure, Enc->s0, Enc->s1, Enc->bp_left2, Enc->bp_right2);
do_move(str->structure, Enc->bp_left2, Enc->bp_right2);
}
int last_en = str->energy;
str->energy = tmp_en;
/* use f_point if we have it */
if (Enc->funct) {
int end = Enc->funct(str, min);
/* undo moves */
if (Enc->bp_left2!=0) do_move(str->structure, -Enc->bp_left2, -Enc->bp_right2);
do_move(str->structure, -Enc->bp_left, -Enc->bp_right);
str->energy = last_en;
Enc->bp_left=0;
Enc->bp_right=0;
Enc->bp_left2=0;
Enc->bp_right2=0;
return (end?1:0);
}
if (Enc->verbose_lvl>1) { fprintf(stderr, " "); print_str(stderr, str->structure); fprintf(stderr, " %d\n", tmp_en); }
/* better deepest*/
if (tmp_en < min->energy) {
min->energy = tmp_en;
copy_arr(min->structure, str->structure);
/* delete degeneracy*/
free_degen(Enc);
/* undo moves*/
if (Enc->bp_left2!=0) do_move(str->structure, -Enc->bp_left2, -Enc->bp_right2);
do_move(str->structure, -Enc->bp_left, -Enc->bp_right);
str->energy = last_en;
Enc->bp_left=0;
Enc->bp_right=0;
Enc->bp_left2=0;
Enc->bp_right2=0;
return 1;
}
/* degeneracy*/
if ((str->energy == min->energy) && (Enc->current_en == min->energy)) {
int found = 0;
int i;
for (i=Enc->begin_pr; i<Enc->end_pr; i++) {
if (equals(Enc->processed[i], str->structure)) {
found = 1;
break;
}
}
for (i=Enc->begin_unpr; !found && i<Enc->end_unpr; i++) {
if (equals(Enc->unprocessed[i], str->structure)) {
found = 1;
break;
}
}
if (!found) {
/* print_stren(stderr, str); // fprintf(stderr, " %6.2f\n", str->energy); */
Enc->unprocessed[Enc->end_unpr]=allocopy(str->structure);
Enc->end_unpr++;
}
}
/* undo moves*/
if (Enc->bp_left2!=0) do_move(str->structure, -Enc->bp_left2, -Enc->bp_right2);
do_move(str->structure, -Enc->bp_left, -Enc->bp_right);
str->energy = last_en;
Enc->bp_left=0;
Enc->bp_right=0;
Enc->bp_left2=0;
Enc->bp_right2=0;
return 0;
}
