Форум на Исходниках.RU

FORTRAN: необычный GO TO

Pers_tmb — Wed, 05 Jan 2011 17:28:33 +0000

Pers_tmb: Спасибо! :) Всё достаточно логично, но необходимо было переубедиться.

AQ(n,m,14) используется как вспомогательный массив.

FORTRAN: необычный GO TO

amk — Wed, 05 Jan 2011 14:12:54 +0000

amk: COMMON в фортране кроме всего прочего гарантирует последовательное расположение своих элементов в памяти.
Поэтому (надо учесть, что в фортране быстрее всего меняются первые индексы, а в C - последние)
AQ(n,m,k) соответствует A(n,m,k) для 1 <= k <= 9
AQ(n,m,10) соответствует G1(n,m)
AQ(n,m,11) соответствует G2(n,m)
AQ(n,m,12) соответствует R0(n,m)
AQ(n,m,13) соответствует ZMU(n,m)
AQ(n,m,14) ничему не соответствует, возможно ошибка при подсчете индексов или в COMMON когда-то был еще один массив
везде 1 <= n, m <= 21

Поэтому ответ на последний ответ изменится, потому что это та же самая память

FORTRAN: необычный GO TO

Pers_tmb — Tue, 04 Jan 2011 23:50:04 +0000

Pers_tmb: Так, это понятно. И всё же массив G1 содержится в памяти следом за массивом A. При изменении элемента AQ(1,1,10) изменится ли элемент G1(1,1)?

FORTRAN: необычный GO TO

NetVir — Tue, 04 Jan 2011 22:11:53 +0000

NetVir: Pers_tmb все описано с стандарте

Цитата

An EQUIVALENCE statement is used to specify the sharing of storage units by two or more entities in a program unit. This causes association of the entities that share the storage units.

Т.е. изменение A(1,1,1) влечет извенение AQ(1,1,1) и наоборот.
На счет COMMON: здесь другая ситуация - к A(21,21,9),G1(21,21),G2(21,21),RO(21,21),ZMU(21,21) можно одновременно обратиться через CVP, но изменение какого либо элемента из заданого перечисления не влечет за собой изменение всех остальных элементов из перечисления. Собсно

Цитата

The COMMON statement provides a means of associating entities in different program units. This allows different program units to define and reference the same data without using arguments, and to share storage units.

FORTRAN: необычный GO TO

Pers_tmb — Tue, 04 Jan 2011 16:50:15 +0000

Pers_tmb: Решив некоторые проблемы продолжаю заниматься переводом с фортрана на Си...

В подпрограмме подключаются глобальные массивы, объявляется ещё один большой массив и используется функция EQUIVALENCE.

Если выкинуть лишнее, код представляет из себя вот что:

[ code ]
COMMON/CVP/A(21,21,9),G1(21,21),G2(21,21),RO(21,21),ZMU(21,21)
DIMENSION AQ(21,21,14)
EQUVALENCE (A(1,1,1),AQ(1,1,1))
[/ code ]

Мне интересно, как связаны теперь массивы...
Если мы изменим A(1,1,1) изменится также AQ(1,1,1) и наоборот при изменении AQ(1,1,1) изменится ли A(1,1,1)?
И, если предыдущее предположение верно, что соответствует элементу AQ(1,1,10) - G1(1,1) или это ни с чем не связанный элемент?

Добавлено 04.01.11, 16:52
Книга редкая. Наверняка издана небольшим тиражом в 1972 году...

FORTRAN: необычный GO TO

NetVir — Fri, 10 Dec 2010 11:50:48 +0000

NetVir: Pers_tmb Книгу не нашел, поэтому по-сути трудно давать комментарии.

FORTRAN: необычный GO TO

NetVir — Thu, 09 Dec 2010 20:42:02 +0000

NetVir:

Цитата andriano @ 09.12.10, 19:12

синтаксически - это присваивание X12-му элементу массива WP результата произведения "0.0" yf "X12"

Так и есть (по идее).

Цитата andriano @ 09.12.10, 19:12

Смущает, правда, бессмысленность вычисления этого произведения.

Не должно так уж сильно смущать. Я и не такое видел. Как показывает практика, "великие мужи" хоть и великие в разрезе своего дела (области науки по своему профилю), но когда доходит дело до прикладного эксперимента (хоть и численного), то реализации оставляют желать лучшего, тем более, когда работы оформляются самими разработчиками численных методов. Обычно, такую работу перекладывают на плечи Ph.D. или чаще Ph.D. Candidate, которые уже более грамотно "прикладывают" работы своего начальства со знанием ЯП как минимум.

Добавлено 09.12.10, 20:44

Цитата Pers_tmb @ 09.12.10, 13:48

Там рядом есть строчка WVW ( DELF,WQ,ROP,DELRO,ETA2 ) = -( DELF/((ETA*DX1)**2)+WQ*(ROP/6.+DELRO/8.))/(ROP/3.+5.*DELRO/24.). Сомневаюсь, что это массив.

И? В чем повод сомнения "WQ" - это не "WP".
Если найду книгу, то смогу дать более подробно консультацию.

FORTRAN: необычный GO TO

andriano — Thu, 09 Dec 2010 19:12:01 +0000

andriano:

Цитата NetVir @ 09.12.10, 12:29

Цитата Pers_tmb @ 09.12.10, 12:23

WP(X12)=0.0*X12

Забивается массив нулями, если X12 - это переменная цикла.

Меня мучают смутные сомнения:
синтаксически - это присваивание X12-му элементу массива WP результата произведения "0.0" yf "X12". Смущает, правда, бессмысленность вычисления этого произведения.
Но, насколько я помню, присваивание значений всему массиву производится либо в заголовке оператором
DATA WP /0.0*X12/
либо в теле без указания индекса в скобках. Т.е. массив записывается как скалярная переменная. Ну еще есть всякие FORALL, но они синтаксически должны отличаться от присваивания элементу массива с заданным индексом.
В общем, хотелось бы видеть фрагмент кода побольше.

FORTRAN: необычный GO TO

Pers_tmb — Thu, 09 Dec 2010 13:48:38 +0000

Pers_tmb: Госмен А.Д., Пан В.М., Ранчел А.К. "Численные методы исследования течений вязкой жидкости". В оригинале издана в 1969 ACADEMIC PRESS - LONDON AND NEW YORK. У нас - 1972 издательство "Мир".

Там рядом есть строчка WVW ( DELF,WQ,ROP,DELRO,ETA2 ) = -( DELF/((ETA*DX1)**2)+WQ*(ROP/6.+DELRO/8.))/(ROP/3.+5.*DELRO/24.). Сомневаюсь, что это массив.

FORTRAN: необычный GO TO

NetVir — Thu, 09 Dec 2010 13:20:13 +0000

NetVir:

Цитата Pers_tmb @ 09.12.10, 12:48

WP(X12)=0.0*X12

Это не функция. Это глупое забивание массива. ИМХО.
Что за книга?

FORTRAN: необычный GO TO

Pers_tmb — Thu, 09 Dec 2010 13:00:05 +0000

Pers_tmb: Нет... сначала SUBROUTINE BOUNCT(K)
потом переменные COMMON
сразу после них вот эти "функции"... причём некоторые гораздо правдоподобнее выглядят: и входных значений побольше, и выражение внутри посложнее.

Перепечатывать все долго. С книги пишу.

FORTRAN: необычный GO TO

NetVir — Thu, 09 Dec 2010 12:52:24 +0000

NetVir: в начале подпрограммы есть implicit none?

FORTRAN: необычный GO TO

Pers_tmb — Thu, 09 Dec 2010 12:48:21 +0000

Pers_tmb: В начале подпрограммы после описания глобальных переменных есть несколько строчек. Одна из них:

WP(X12)=0.0*X12

В коде же непосредственно есть строчка:

A(1,J,NW)=WP(X2(J))/R(J)

При этом, переменная X12 нигде не встречается.

FORTRAN: необычный GO TO

NetVir — Thu, 09 Dec 2010 12:39:52 +0000

NetVir:

Цитата Pers_tmb @ 09.12.10, 12:36

Иными словами, если на фортране написано I=5., то это float I=5, а если I=.5, то - float I=0.5?

да

Добавлено 09.12.10, 12:40

Цитата Pers_tmb @ 09.12.10, 12:36

Но просто смысла в такой функции нет.

такое часто бывает. Я лично уже привык.

Добавлено 09.12.10, 12:41
покажи описание переменной

Цитата Pers_tmb @ 09.12.10, 12:36

и блок кода, где она учавствует

FORTRAN: необычный GO TO

Pers_tmb — Thu, 09 Dec 2010 12:36:22 +0000

Pers_tmb:

Цитата

WP(X12)=0.0*X12

Я так полагаю, что это что-то вроде функции... потому что эта строчка соседствует с оными. Но просто смысла в такой функции нет.

Добавлено 09.12.10, 12:39

Цитата NetVir @ 09.12.10, 12:29

Цитата Pers_tmb @ 09.12.10, 12:23

I=0. и I=.0 ?

Присвоение вещественного нуля переменной.

Иными словами, если на фортране написано I=5., то это float I=5, а если I=.5, то - float I=0.5?

FORTRAN: необычный GO TO

NetVir — Thu, 09 Dec 2010 12:29:12 +0000

NetVir:

Цитата Pers_tmb @ 09.12.10, 12:23

I=0. и I=.0 ?

Присвоение вещественного нуля переменной.

Цитата Pers_tmb @ 09.12.10, 12:23

WP(X12)=0.0*X12

Забивается массив нулями, если X12 - это переменная цикла.

FORTRAN: необычный GO TO

Pers_tmb — Thu, 09 Dec 2010 12:23:59 +0000

Pers_tmb: Я пробовал использовать f2c, но что-то не получилось. Думаю, надёжнее самому перевести код. Тут всего тысячи полторы-две строчек.. :)

Дабы не засорять раздел новой темой, задам вопрос тут.

Что значат записи вроде: I=0. и I=.0 ?

И самое интересное: WP(X12)=0.0*X12 ?

FORTRAN: необычный GO TO

tulebaev — Tue, 30 Nov 2010 08:25:17 +0000

tulebaev: Существует утилита Fortran2C (f2c) -- преобразует Fortran-код в C-листинг.

А самое быстрое -- скомпилировать Fortran-код в obj и подключить к проекту на C

FORTRAN: необычный GO TO

Pers_tmb — Mon, 29 Nov 2010 08:24:00 +0000

Pers_tmb: Примерно так и предполагал, но надо было удостовериться.
Спасибо! :)

FORTRAN: необычный GO TO

volvo877 — Sun, 28 Nov 2010 23:30:59 +0000

volvo877: Нормальное использование, "вычислимый goto", если мне память не изменяет.

GO TO (список_меток), n
Сначала вычисляется выражение n, потом происходит переход к n-ой метке из списка. В твоем случае, только при K = 2 произойдет переход ко второй метке, при K = 1, или 3 .. 8 - переход к метке 1.

FORTRAN: необычный GO TO

Pers_tmb — Sun, 28 Nov 2010 19:19:07 +0000

Pers_tmb: Здравствуйте, уважаемые! :)

Вынужден переводить на понятный мне язык Си некоторую программку на фортране из книги, изданной в Лондоне аж в 1969 году. Фортран этот настолько бородат, что ни в одном русскоязычном учебнике я такого применения GO TO не увидел. Признаюсь, просматривал очень бегло.

Вот, собственно, код, вызвавший затруднение.

[ code ]

GO TO (1,2,1,1,1,1,1,1),K
1 SOURCE=0.
RETURN

2 SOURCE=A(I,J,NW)
RETURN
END

[ /code ]

Кто-нибудь знает, как это работает?