Главная страница  Алгоритмы 

[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [ 42 ] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57]

,,дующем «Замечании и подтверждении к алгоритму 50» ;П. Наура.

I

Подтверждение к алгоритму 50

Б. Рандел (Randel В. «САСМ», 1962, № 1)

Процедура INVHILBERT была проверена с использованием транслятора DEUCE ALGOL. Необходимо было внести следующие исправления... *.

Транслированная программа с использованием ~20-разрядной мантиссы и плавающей запятой дала следующий сегмент:

16.0 120.0000 240.0002 -140.0000 -120.0 1200.0000 -2700.0000 1680.0019 240.0 -2700.0000 6480.0000 -4200.0000 - 140.0 1680.0019 -4200.0000 2800.0039

Замечание и подтверждение к алгоритму 50 П. Наур (Naur Р. «САСМ», 1963, № 1)

Кроме внесения исправлений, указанных Б. Ранделом («САСМ», 1962, № 1), мы модифицировали и упростили . алгоритм следующим образом ...**.

В сумме эти изменения, кроме сок)ращ€ния длины программы, повысили скорость в 1.6 раза. В результате получен следующий алгоритм***.

Как исходный, так и улучшенный алгоритм успещно были проверены в системе GIER-ALGOL (30-разряд-ная мантисса). Тест-программ а включала следующее.

1. Вывод матрицы 4X4 для сравнения с результатами Рандела. (приведенными выше. - Прим. ред. перев.). Результаты:

16.000000 -120.000000 240.0000000 -140.000000 -120.000000 1200.000000 -2700.0000000 1680.000000 240.000000 -2700.000000 6480.0000000 -4200.000000 -140.000000 1680.000000 -4200.0000000 2799.999977

2. Вычисление обращенного сегмента для п-\, .., п с помощью процедуры INVHILBERT, перемножение на сегмент гильбертовой матрицы. Результаты

t * Указываются две поправки к алгоритму 50. (Прим. ред.) Jf ** Указываются четыре изменения, внесенных в алгоритм 50а. (Прим. ред.).

*** "Приводится новый вид процедуры INVHILBERT, из которой ординарной переработкой получен алгоритм 50а. (Прим ред.)



сравнивались с единичной матрицей. Для получения от, носительной ошибки максимальное отклонение делилось на наибольший элемент обращенной 1матрицы. Некоторые из результатов, которые в целом были удовлетворительными, приводятся ниже (табл. 23).

Т а б л и ца 23

Порядок п

Элемент с максимальной ошиокой

Абсолютное значение максимальной ошибки

Наибольший элемент результата

Относительная ошибка

* 3

9 12 15

S [3,3 S 2,4 S 2,8 S 5,9 s[1.12]

2.38X10- 4.39X10- 1.24X102 1.54X10» 1.06X10"

1.92X10"

4.41X10»

1.22X10"

3.66X10»

1.15Х1№»

1.24X10-» 9.96X10-10 1.01X10-» 4.21X10-» 9.22X10-»

Время выполнения обращения к пересмотренной процедуре было следующим:

п Время, .с.

5 0.2 10 0.6 15 -1.3



Подтверждения и замечания к алгоритмам опубликованным в предыдущей серии

Подтверждение к алгоритмам 4а, 7а, 18а, 19а, 23а, 24а, 26а, 31а, 41а, 46а, 48а, 51а, 55а, 56а, 60а, 66а, 70а, 75а, 85а, 105а, 120а, 123а, 126а, 140а, 167а, 168а, 169а, 180а, 185а, 197а, 203а, 204а, 205а, 209а, 218а, 220а, 230а, 231а, 233а, 238а Л. С. Кривонос, 3. А. Шиншинова, Москва май 1970

Вышеперечисленные алгоритмы были транслированы на машине БЭСМ-6 в системе БЭСМ-АЛГОЛ {Щ и для примеров, приведенных в «Свидетельствах» к этим алгоритмам, дали правильные результаты. Предварительно в алгоритмы 48а, 105а и 126а были внесены исправления, опубликованные в «Приложениях 1» к выпускам алгоритмов данной серии [27, 28]. В табл. 24 приведены значения времени трансляции и времени счета одного примера для некоторых из Вышеуказанных алгоритмов.

В «Свидетельстве к алгоритму 105а» {26, с. 13, 9-я стр. сверху была замечена одна опечатка. Вместо х2=2 должно быть х2=-2.

Таблица 24

Номер , алгоритма

Время, с

Номер алгоритма

Время, с

Номер алгоритма

Время, с

трансляции

счета

трансляции

счета

транс-лягции

счета

0.94

0.12

0.90

0.16

1.82

0.08

1.18

0.02

0.76

0.06

1.16

0.10

0.56

0.06

0.52

0.12

2.14

9.54

1.00

0.04

0.58

0.04

1.88

2.42

0.78

0.08

0.86

0.08

2.40

8.01

0.82

0.04

1.38

0.20

1.80

0.14

Подтверждение и замечания к алгоритмам 29а, 33а, 36а, 37а, 38а, 67а, 77а, 80а, 84а, 122а, 124а, 128а, 131а, 134а, 188а, 189а, 191а, 199а, 227а 1 М. И. Агеев, Москва, апрель 1971

Вышеперечисленные алгоритмы были транслированы на машине БЭСМ-6 в системе БЭСМ-АЛГОЛ и для примеров, приведенных в «Свидетельствах» к этим алгоритмам, дали правильные результаты.




[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [ 42 ] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57]

0.0087