Разговорчики в строю

[spamsink]

Из найденных на диске исходных текстов стандартной математической библиотеки АЛГОЛА-60 (библиотечные процедуры хранились в текстовом виде, по-видимому, в силу каких-то особенностей автора системы, и в кодировке, в точности совпадающей с кодировкой перфокарт, т. е. с битом четности, что было довольно неочевидно при поиске интересностей во всех возможных кодировках), одна содержит довольно забавную диагностику об ошибке. Сейчас в такой манере вычислительные библиотеки с пользователями не общаются.

_REAL _PRОCEDURE ЛАГРАНЖ2(A,Х,У,Х1,Х2,P,П,П1);
_REAL Х1,Х2;_INTEGER P,П,П1;_ARRAУ A,Х,У;
_BEGIN _REAL S,Z,Z1;_INTEGER I,J,J1,К;
_ARRAУ S1[0:20];
_INTEGER Я,Я1;
_IF Х2<A[1]_THEN _BEGIN Z:=Х2;_GО _TО М3;_END ;
_IF Х1<Х[1]_THEN _BEGIN Z:=Х1;_GО _TО М3;_END ;
S:=0;_FOR J1:=1_STEP 1_UNTIL П1_DО
_IF Х2≥A[J1]∧Х2≤A[J1+1]_THEN _GО _TО М1;Z:=Х2;_GО _TО М2;
М1:_FOR I:=1_STEP 1_UNTIL П_DО
_IF Х1≥Х[I]∧Х1≤Х[I+1]_THEN _GО _TО М;Z:=Х1;
_IF Х1≤Х[П]_THEN _GО _TО М;
М2:_IF Х2≤A[П1]_THEN _GО _TО М;
М3:
ОUTPUT(‘T’,‘ПРИ  ОБРАЩЕНИИ  К  ПРОЦЕДУРЕ  ЛАГРАНЖ2  АРГУМЕНТ=’,
‘E’,Z,‘T’,
‘ВЫШЕЛ  ЗА  ПРЕДЕЛЫ  ТАБЛИЦЫ. РАЗБЕРИСЬ  В  ЧЕМ  ДЕЛО, А  Я  ПРЕРЫВАЮ  СЧЕТ. ДО  ВСТРЕЧИ!’);
ОUTPUT(‘×’);
S:=S/0;
М:_IF I≤П-P_THEN J:=I_ELSE J:=П-P;
_IF J1>П1-P_THEN J1:=П1-P;
Я1:=J1;
_FOR Я:=
_STEP 1_UNTIL P_DО _BEGIN
S:=0;
_FOR I:=0_STEP 1_UNTIL P_DО _BEGIN
Z:=Z1:=1;
_FOR К:=0_STEP 1_UNTIL P_DО _BEGIN
_IF К≠I_THEN Z:=Z×(Х1-Х[К+J])/(Х[I+J]-Х[К+J]);
_END ;
S:=S+У[I+J+(J1-1)×П]×Z;_END ;
S1[Я]:=S;J1:=J1+1;_END ;
S:=0;
_FOR I:=0_STEP 1_UNTIL P_DО _BEGIN
Z:=Z1:=1;
_FOR К:=0_STEP 1_UNTIL P_DО _BEGIN
_IF К≠I_THEN Z:=Z×(Х2-A[К+Я1])/(A[I+Я1]-A[К+Я1]);
_END ;
S:=S+S1[I]×Z;_END ;
ЛАГРАНЖ2:=S;
_END 

Необходимое примечание: Название процедуры ЛАГРАНЖ2 вряд ли имеет отношение к точке Лагранжа L2, по направлению к которой сейчас летит телескоп имени Уэбба.

Comments

[vak]

Прикольный способ завершения программы.

S:=S/0;

[spamsink]

Эквивалентно kill(getpid(), SIGFPE);, а пишется короче. :)

[vak]

Что бы означало М3:0ОUTPUT?
Вроде тут пробел должен быть вместо нуля.

[spamsink]

Это перевод строки. В коде УПП байт 200 - это ноль, а 000 - пустышка. Я вчера руками на строки делил, а это проглядел, спасибо.

[juan_gandhi]

Вот эти две строки ценные тоже:

М:_IF I≤П-P_THEN J:=I_ELSE J:=П-P;
_IF J1>П1-P_THEN J1:=П1-P;

[spamsink]

Что не так?

[juan_gandhi]

О. Прочитал снова. Нет, все так. Я перепутал J и J1.

[dijifi]

Сейчас говорит «что-то пошло не так. Попробуйте ещё раз позже»

[spamsink]

По крайней мере не говорит "разберись, в чем дело".

[dijifi]

Это только персонально от эппл суппорт

[spamsink]

Я не удивлюсь, если формулировка "Вы это неправильно <глагол 2 л. мн. ч.>" впервые применили в СССР.

[vak]

Попытался было переписать на современный алгол и попробовать выполнить, но упёрся в непонимание, что же тут вычисляется. Какая-то интерполяция по Лагранжу.

[spamsink]

Я тоже не понимаю. Интерполяция по Лагранжу должна принимать два массива Х и Y одной длины и возвращать один массив коэффициентов полинома той же длины. А почему там три входных массива, где А - не той длины, что X и Y, и что там за значение функции вычисляется - неясно.

[vak]

Возможно, это вычисление определённого интеграла функции, представленной в виде полинома Лагранжа.

[spamsink]

Всё равно непонятно, зачем третий массив.

Мне тут объяснили - это интерполяция функции от двух переменных. В массиве S образуются интерполированные значения по одной координате по линии X1, а потом они используются для интерполяции по линии X2. Одно осталось не очень понятно, что там за трюки с индексами по массиву Y.