По сути изначально задача ставилась как счет грвитационных сил Земли, с точночтью О(h^2). Координаты брались простые декартовы, а не сферические. Алгоритм прост. Но требует слишком много операций в R^3. Надо думать
думать, думать, думать 
Вроде там возникает Теплицева структура, расширяется до циркулянта... туда сюда, а толку пока мало. Все дело сводится к дискр. преобр. Фурье и там должен быть реальный выигрышь по времени 
Нужен совет, желательно бы и про ДПФ в R^3 и может вообще мысли по всему этому поводу
может кто с этим сталкивался?