void revers ( int ivec[], int kdim );
int simdo ( int qind, int qfor, int iprod[], int kdim, int *jsub, int ivec[] );