#define X(i, j) i - 1 + *n * (j - 1)
-#define DINDEX(i, j) n * (i - 1) - i*(i - 1)/2 + j - i - 1
-
-int give_index(int i, int j, int n)
-{
- if (i > j) return(DINDEX(j, i));
- else return(DINDEX(i, j));
-}
-
void distDNA_indelblock(unsigned char *x, int *n, int *s, double *d)
{
int i1, i2, s1, s2, target, N, start_block, end_block;