Design of specially adapted reactive coordinates to economically compute potential and kinetic energy operators including geometry relaxation.