|
Язык: Фортран
Программный комплекс SYSINT(SYSINTM) предназначен для решения задачи
на собственные значения для системы интегральных уравнений.
Модифицированная итерационная схема реализована на основе непрерывного
аналога метода Ньютона (НАМН) [2]. В программе SYSINTM реализована
модификация с одновременным обращением оператора производной Фреше,
предложенная в [3] и обеспечивающая существенный выигрыш во времени
при раcчетах на векторных вычислительных системах [3,1]. Более подробное
описание пограммного пакета дано в [1].
Обращение к комплексам SYSINT и SYSINTM осуществляется операторами:
CALL SYSINT (N, LM, LMN, X, EV, FI, TAU0, NTAU, HT5, GNORM, EPS,
NMAX, F, V, V1, R1, R2, R3, TKSI, B)
CALL SYSINTM (N, LM, LMN, X, EV, FI, TAU0, NTAU, GNORM, EPS, NMAX,
F, V, V1, R1, R2, R3, TKSI, B, A, DAT, NMOD)
где:
| N |
- |
число узлов сетки по аргументу x; |
| LM |
- |
число уравнений системы; |
| LMN |
- |
размерность рабочих массивов, LMN=N*LM; |
| X |
- |
массив узлов сетки по x размерности N; |
| EV |
- |
собственное значение lambda. При обращении к программе
этому параметру присваивается начальное приближение для
собственного значения lambda_0, после окончания работы
программы здесь находится k-e приближение lambda_k; |
| FI |
- |
двумерный массив размерности LM*N. При обращении к
программе в нем задается начальное приближение,
на выходе здесь находится полученное решение; |
| TAU0 |
- |
начальное значение шага; |
| HT5 |
- |
шаг сетки; |
| GNORM |
- |
значение G из условия нормировки; |
| EPS |
- |
заданное малое число; |
| NMAX |
- |
максимально допустимое число итераций k.
При его превышении происходит выход из программы; |
| F, V, V1, R1, R2, R3, B, TKSI |
- |
рабочие массивы;
массивы F, V, V1, R1, R2, R3 имеют размерность LMN;
массив TKSI имеет размерность N;
двумерный массив B имеет размерность LMN*LMN; |
| A, DAT |
- |
дополнительные двумерные рабочие массивы программы SYSINTM,
также имеющие размерность LMN*LMN; |
| NMOD |
- |
параметр модифицированного алгоритма.
|
Литература:
- Земляная Е.В. SYSINT(SYSINTM) - комплекс программ для численного решения
задачи на собственные значения для системы интегральных уравнений.
Сообщение ОИЯИ. P11-94-120, Дубна, 1994.
- Пузынин И.В., Амирханов И.В., Земляная Е.В., Первушин В.Н., Пузынина Т.П.,
Стриж Т.А., Лахно В.Д. Обобщенный непрерывный аналог метода Ньютона
для численного исследования некоторых квантово-полевых моделей.
ЭЧАЯ, 1999, т.30, вып.1, стр.210-265.
- Puzynin I.V., Amirkhanov I.V., Puzynina T.P., Zemlyanaya E.V.
The Newtonian Iterative Scheme with Simultaneous Calculating
the Inverse Operator for the Derivative of Nonlinear Function.
JINR Rapid Comm., No.5[62]-93, Dubna, 1993, p.63-73;
In: Proc. Intern. Conf. Programming and Mathematical Methods
for Solving Physical Problems. P&MM'93, Ed. by Yu.Yu.Lobanov and
E.P.Zhidkov, Dubna, 1993, p.30-34.
- Давиденко Д.Ф. Препринт ИАЭ им. И.В.Курчатова, ИАЭ-1963, Москва,1970.
Архив программы включает в себя описание программы (в формате .pdf),
фортранные тексты программ.
|
|
