Советское - значит, отличное

[spamsink]

...ракета-носитель отклонилась от намеченного курса и была взорвана по команде с Земли

Знакомо? А что ж вы хотели:

F О R Е Х ИПМ АН СССР 4.12 ОТ 25.06.85 ВРЕМЯ:  00.00

   1                РRОGRАМ МАIN
   2                DОUВLЕ РRЕСISIОN Х,У
   3                Х=1.D0
   4                У=1.0D0
   5                Х=Х-У
   6                РRINТ 10,У,Х
   7           10   FОRМАТ (1Х,2D50.40)
   8                SТОР
   9                ЕND


               ТАБЛИЦА ОТНОСИТЕЛЬНЫХ АДРЕСОВ

        1    2    3    4    5    6    7    8    9   10
    0 :    :    :   1:   3:   6:  12:    :  21:    :  22
КОМАНД 00030    КОНСТАНТ 00006    ПАМЯТЬ ДЛЯ ПЕРЕМЕННЫХ 00004
ДЛИНА ПОДПРОГРАММЫ МАIN   00042

           *ЕХЕСUТЕ
     9.9999999999999999999999735302203982320601-001    1.6543612251060553497428173841399257071316-024

Comments

[raydac.livejournal.com]

ну зато денег было столько что как говаривал один выделяющий денег на пуски и всю эту хрень про каждый старт - "мы стреляем городами"

[mochalkina.livejournal.com]

а почему так?

[spamsink]

Я сам удивляюсь. Но, как говорится, "в программах бывают ошибки, и очень серьезные".

[raydac.livejournal.com]

а может там написано было в лицензии - "не для использования в ракетной и атомной промышленности"? :)

[lupus-lupusum.livejournal.com]

так флоатинг поинт же
с тех давних пор ничего не изменилось кстати, см. Кнута

[dvv.livejournal.com]

Можно ещё см. IEEE. Там тоже много интересного пишут.

[spamsink]

Представление степеней двойки даже в плавающей точке должно быть точное (да и не только степеней двойки, а всех целых от 1 до 2^{разрядность мантиссы}-1). А видим - неточное.

[potan.livejournal.com]

На каком компьютере?
Слышал что на БЭСМ-6 арифметика с двойной точностью работала очень плохо и ей ни кто не пользовался.

[dvv.livejournal.com]

Очевидно, на БЭСМ-6. С компилятором FOREX.

[potan.livejournal.com]

Тогда это известный аппаратный баг. Типа fdiv на некоторых интелловских процессорах.

[spamsink]

Очень медленно она работала - это да. А так вроде все более или менее правильно, например, мой любимый тест: x=(cos(x)+x)/2 точно - без дребезга в младшем бите - сходится к
0.739085133215160641655308,
в то время как истинное значение -
0.739085133215160641655312

[reedcat.livejournal.com]

Ей не пользовались не поэтому...
Ей не пользвались потому, что арифмтетика с одинарной точностью на БЭСМ-6 была совсем не вдвое менее точной, чем арифметика с двойной. Смысла не было, удваивать расход памяти и замедлять счет на порядок, ради 5 десятичных разрядов.

[dvv.livejournal.com]

А что говорят *FORTRAN и *FTN?

[spamsink]

*FTN - то же самое. *FORTRAN - правильно!

[dvv.livejournal.com]

(В.П. Иванников: „Ракета ушла за бугор”.)

[dvv.livejournal.com]

Кстати, 80*log10(2.) чему у нас равно? То–то же. И дайте, дайте мне ОТЛАДЧИК ОС. И напомните, как им пользоваться.

ОТЛАДЧИК ОС

[potan.livejournal.com]

gdb

[spamsink]

Правильно, ошибка на один младший разряд мантиссы.

А чего отладчик ОС? Пишешь TЕЛЕ 77770 (или ему 011 слов было нужно? - тогда 77767) - получаешь результат.

[dvv.livejournal.com]

  00000            УИИ   ’00010’(’07’)
                   ПВ    ВЕGIN (’11’)
  00001   /0001/   СЧ    =В’0000000000000000’
                   ЗП    Х     +’00001’
  00002            ЗП    АСС*
                   СЧ    =В’4050000000000000’
  00003            ЗП    Х
                   СЧ    =В’0017777777777776’
  00004            ЗП    У     +’00001’
                   ЗП    АСС*
  00005            СЧ    =В’4017777777777777’

Шо за поэбень?!

[spamsink]

My sentiment exactly.

[reedcat.livejournal.com]

А что не нравится? :)

[dvv.livejournal.com]

И нех на советское гнать:

Ф О Р Т Р А Н
   /16.07.73/
                РRОGRАМ МАIN
                DОUВLЕ РRЕСISIОN Х,У
                Х=1.D0
       2        У=1.0D0
       3        Х=Х-У
       4        РRINТ 10,У,Х
          10    FОRМАТ (1Х,2D50.40)
       5        SТОР
                ЕND



                                      ТАБЛИЦА  ОТНОСИТЕЛЬНЫХ АДРЕСОВ

         1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16
17    18    19    20
 000  :00001:00003:00005:00010 00016
           *ЕХЕСUТЕ
     1.0000000000000000000000000000000000000000+000    0.0000000000000000000000000000000000000000+000

ИПМ, ИПМ…

[spamsink]

Как раз этот Фортран - калька с CDC-шного:

Весь транслятор был разбит на блоки, и каждый докладывал свой блок... Списывали, по правде говоря, с транслятора на CDC-1604. Но для нас, конечно, не нюхавших программирования, тем более, системного, это была работа творческая...

[reedcat.livejournal.com]

Ну и что? :)

Вот не удивил ни разу... :)

Что же касается комментариев:

- к ракетам "это" не имеет никакого отношения :)

- то, что на БЭСМ-6 тип REAL имел диапазон 10 в 19, а не 10 в 12, как было в IBM / CDC c которых драли Фортран, так это плюс для БЭСМ-6, а не минус... А арифметика DOUBLE PRECISION как имела точность 10 в 24, так и имеет... И на тех и на тех машинах :)
- на дураков дипломников и аспирантов, которые думали, что если поставить DOUBLE PRECISION то "оно будет считать точнее", я, в свое время, насмотрелся
- физическая задача, требующая для своего решения DOUBLE PRECISION в принципе не может выдавать осмысленный результат, оно будет считать бред.
- и вообще, "а у вас негров линчуют!" :)

[spamsink]

Насчет ракет это была шутка. Но, как говорится, "если вы услышали тринадцатый удар часов, возникают сомнения и в правильности первых двенадцати". С такими ошибками в тривиальных местах я удивляюсь, как вообще что-то работало.