Оценка прогресса расчета частот

[Гесс]

Считаю я себе частоты в гауссиане.
И тянется этот процесс мучительно медленно.
И задался я вопросом есть ли там какой либо спопоб посмотреть когда ж это мучение закончится (по логике вроде должен быть).
Мне посоветовали grepать "EnFreq"
Я думал чтоу меня при 8 атомной системе будет 3*8-6=18 циклов, однако сейчас идет уже двадцатый.
Выглядят они "расширяющимся треугольником", ну либо "половиной матрицы":
[cut]

...
Leave EnFreq: IXYZ=17 JXYZ=15 IStep= 1.
Leave EnFreq: IXYZ=17 JXYZ=15 IStep= 2.
Leave EnFreq: IXYZ=17 JXYZ=16 IStep= 1.
Leave EnFreq: IXYZ=17 JXYZ=16 IStep= 2.
Leave EnFreq: IXYZ=18 JXYZ= 1 IStep= 1.
Leave EnFreq: IXYZ=18 JXYZ= 1 IStep= 2.
Leave EnFreq: IXYZ=18 JXYZ= 2 IStep= 1.
Leave EnFreq: IXYZ=18 JXYZ= 2 IStep= 2.
...
Leave EnFreq: IXYZ=18 JXYZ=16 IStep= 1.
Leave EnFreq: IXYZ=18 JXYZ=16 IStep= 2.
Leave EnFreq: IXYZ=18 JXYZ=17 IStep= 1.
Leave EnFreq: IXYZ=18 JXYZ=17 IStep= 2.
Leave EnFreq: IXYZ=19 JXYZ= 1 IStep= 1.
Leave EnFreq: IXYZ=19 JXYZ= 1 IStep= 2.
Leave EnFreq: IXYZ=19 JXYZ= 2 IStep= 1.
Leave EnFreq: IXYZ=19 JXYZ= 2 IStep= 2.
...
Leave EnFreq: IXYZ=19 JXYZ=17 IStep= 1.
Leave EnFreq: IXYZ=19 JXYZ=17 IStep= 2.
Leave EnFreq: IXYZ=19 JXYZ=18 IStep= 1.
Leave EnFreq: IXYZ=19 JXYZ=18 IStep= 2.
Leave EnFreq: IXYZ=20 JXYZ= 1 IStep= 1.
Leave EnFreq: IXYZ=20 JXYZ= 1 IStep= 2.
Leave EnFreq: IXYZ=20 JXYZ= 2 IStep= 1.
Leave EnFreq: IXYZ=20 JXYZ= 2 IStep= 2.

[/cut]
Общее число таких точек уже выше 400. Есть ли какая то оценка хотя бы сверху?
Как это делать правильно?

[amge]

Я бы предположил не 3*8-6=18, а 3*8=24.

[Гесс]

Я бы предположил не 3*8-6=18, а 3*8=24.

То есть еще 19+20+21+22+23 чуть больше 120, каждый из них в два степа, итого еще 240 итераций. меньше чем по итерации в час... Как говорила русская мышь: "Ладно, мы подождем"...

[amge]

Кстати, для удобства слежения (за Природой) я использую скрипт prir (аналог комбинации tail -f и grep).
Для вашего случая будет что-то вроде

Код: Выделить всё

#!/usr/bin/perl -w

die "Usage: $0 output.file\n" unless @ARGV;

open FILE, '<', $ARGV[0] or die "Can't read $ARGV[0]: $!\n";
for (;;) {
  for ($curpos = tell(FILE); $_ = <FILE>; $curpos = tell(FILE)) {
    # for numeric hessian
    /Leave EnFreq: / && print;
    /Normal termination of Gaussian/ && exit;
  }
  sleep 1;
  seek(FILE, $curpos, 0);
}

[surius]

Я бы предположил не 3*8-6=18, а 3*8=24.

каждый из них в два степа, итого еще 240 итераций. меньше чем по итерации в час... Как говорила русская мышь: "Ладно, мы подождем"...

вопрос из чисто спортивного интереса: у вас симметрия системы не Cs случайно? не до конца понятно почему в 2 шага? почему не в 1 или в 6 или....

[VTur]

Гесс, если считаете частоты в попловском базисе, то считайте в шаровых функциях, если в даннинговских - может не протянуть, переходите в сферические гармоники, но каждый раз добавляйте int=ultrafine.

Верить мне не надо. сделайте расчеты и так, и так. Когда перейдете от сферических гармоник к шаровым функциям, сравнивайте верхние частоты, когда от prund сетки к ultrafine нижние.

[Гесс]

вопрос из чисто спортивного интереса: у вас симметрия системы не Cs случайно? не до конца понятно почему в 2 шага? почему не в 1 или в 6 или....

Незнаю почему. Симметрия задавалась как Td, через Z-matrix, однако гауссиан считает что она С1 и вроде как ведет расчет без симметрии.

[Гесс]

Гесс, если считаете частоты в попловском базисе, то считайте в шаровых функциях, если в даннинговских - может не протянуть, переходите в сферические гармоники, но каждый раз добавляйте int=ultrafine.

Верить мне не надо. сделайте расчеты и так, и так. Когда перейдете от сферических гармоник к шаровым функциям, сравнивайте верхние частоты, когда от prund сетки к ultrafine нижние.

Дык у меня вопрос не про качество полученных частот, а про то получу я их в обозримом будущем или нет.
Частоты считаются численно, проверить их аналитически в СFOURE несложилось, возможно когда нибудь мои ручки долезут до дальтона, посмотрю там.
При том что мне счас светит изучение Амбера и Дэша которые меня ни разу не греют - шанс что я полезу в Дальтон достаточно неплохи. Инпут у него выглядит так себе, но если будет падать реже чем СFOUR (a аутпут будет содержательнее чем у последнего) - может и подружусь.