close

Вход

Забыли?

вход по аккаунту

?

Вычисление интегралов по многомерным областям на многопроцессорных вычислительных системах.

код для вставкиСкачать
Вычислительные технологии
Том 11, № 3, 2006
ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ
ПО МНОГОМЕРНЫМ ОБЛАСТЯМ
НА МНОГОПРОЦЕССОРНЫХ
ВЫЧИСЛИТЕЛЬНЫХ СИСТЕМАХ∗
Д. Я. Рахматуллин
Институт математики с вычислительным центром
Уфимского НЦ РАН, Россия
e-mail: nicknamej@rambler.ru
An algorithm for approximate integration of functions is considered. The domain is
a hyper prism with one smooth curved boundary. Results of the computations using the
proposed parallel algorithm are analyzed.
Введение
Вычисление интегралов по областям больших размерностей сдерживается отсутствием
достаточных вычислительных мощностей. Это NP-полная задача, и даже самые быстрые
персональные компьютеры вкупе с лучшими методами не могут справиться с этой задачей
за разумное время при удовлетворительной точности вычислений. Применение суперкомпьютеров — мощных многопроцессорных вычислительных систем — хотя и не устранило
полностью, но существенно ослабило ограничения на входные параметры задачи — размерность пространства и количество верных знаков в ответе. Если вычислительный алгоритм
достаточно хорошо распараллеливается, то предельная размерность пространства, в котором требуется найти интеграл, ограничивается лишь количеством процессоров конкретного суперкомпьютера, а оно на данный момент может измеряться тысячами. Мы используем
алгоритм, основанный на применении решетчатых кубатурных формул с ограниченным
пограничным слоем. Он почти идеально распараллеливается, т. е. исходная задача почти
поровну разбивается на ряд подзадач, независимо решаемых на каждом процессоре. При
этом обмен данными между процессорами сведен к минимуму, так что итоговое время вычислений уменьшается с увеличением числа процессоров практически пропорционально.
Суперкомпьютер МВС-15 000, на котором тестировалась программа, установлен в Межведомственном суперкомпьютерном центре РАН. Для вычислений предоставлялось 922
процессора с распределенной, т. е. отдельной для каждого двухпроцессорного вычислительного модуля, памятью.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований
(грант № 03-07-90077-в, № 06-01-00597-a) и Программы президиума РАН (№ 17).
c Институт вычислительных технологий Сибирского отделения Российской академии наук, 2006.
°
∗
117
118
Д. Я. Рахматуллин
1. Алгоритм
Рассмотрим пространство Соболева Wpm (Ω), Ω ⊂ Rn , функций, интегрируемых c p-й степенью (p > 1) вместе с производными до m-го порядка включительно. Здесь m ∈ N —
гладкость функций, а Ω — ограниченная область с гладкой границей Γ, лежащая вместе
со своим замыканием в некоторой гиперпризме (многомерном прямоугольном параллелепипеде) Q [1].
Мы берем это пространство в одной из эквивалентных норм:

¯p 1/p
¯
Z
¯
¯X
¯
p
2 m/2 2πixk ¯ 

kf kW
,
kgkW
fpm(Ω) = g∈Winf
fpm(Q) , kgkWfpm(Q) = |g0| + dx ¯¯ gk (1 + |k| ) e ¯¯
m
p (Q),
Q
g|Ω =f |Ω
где
gk =
Z
k6=0
dx g(x) e−2πixk .
Q
Будем использовать последовательность кубических решеток {hk}k∈Zn при h → 0. Наша
задача состоит в возможно более точном приближении интеграла
Z
Ω
I (f ) := dxf (x)
Ω
решетчатыми кубатурными формулами вида
X
KhΩ (f ) := hn
ck (h)f (hk),
h → 0.
hk∈Ω,
k∈Z n
Введем несколько определений.
Определение 1. Последовательность кубатурных формул (ПКФ) называется оптимальной, если каждый ее элемент минимизирует норму его разности с точным значением
интеграла по множеству коэффициентов этого элемента:
°
°
∀h KhΩ,opt := arg min °I Ω − KhΩ °(W M (Ω))∗ .
{ck (h)}
p
Определение 2. Последовательность кубатурных формул KhΩ,as называется асимптотически оптимальной, если она имеет оптимальные свойства в пределе, т. е. выполняется
следующее равенство:
°
°
° Ω
Ω,as °
°I − Kh ° m ∗
(W (Ω))
°
° p
lim °
= 1.
°
Ω,opt
h→0
°I Ω − Kh ° m ∗
(Wp (Ω))
Определение 3. Последовательность кубатурных формул называется оптимальной
по порядку, если существует такая, не зависящая от h константа C, что верна оценка
°
°
° Ω
Ω,ord °
°I − Kh ° m ∗
(W (Ω))
° p
6 C.
lim °
Ω,opt °
h→0 ° Ω
°I − Kh ° m ∗
(Wp (Ω))
ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ ПО МНОГОМЕРНЫМ ОБЛАСТЯМ ...
119
Определение 4. Последовательностью кубатурных формул с ограниченным пограничным слоем (ОПС) называется такая ПКФ, для которой выполняются два условия
1. Все ее коэффициенты равномерно ограничены по h и k:
sup |ck (h)| 6 L1 .
h, k
2. Всем узлам решетки, содержащимся в области Ω, за исключением пограничного слоя
толщины L2 h, соответствуют коэффициенты, равные единице:
ρ(kh, Rn \Ω) > L2 h ⇒ ck (h) = 1.
∀k, h :
Мы будем рассматривать именно ОПС-формулы. Верна следующая теорема:
n
Теорема 1. Пусть 1 6 p1 6 p2 < ∞ и
< m1 < m2 . Последовательность кубатурp1
ных формул KhΩ (f ) с ОПС асимптотически оптимальна на пространствах
n
o
f m (Ω)
W
p
m∈(m1 , m2 ), p∈(p1 , p2 )
тогда и только тогда, когда она оптимальна на них по порядку [2].
Определение 5. Последовательностью функционалов погрешности (ПФП) lhω (f ), соответствующей ПКФ Khω (f ), называется последовательность, каждый элемент которой
представляет собой разность точного интеграла I ω (f ) от функции f по области ω и соответствующего элемента ПКФ:
lhω (f ) := I ω (f ) − Khω (f ),
h → 0,
или в виде обобщенной функции:
lhω (x) := χω (x) − hn
X
ck (h)δ(x − hk),
h → 0,
hk∈ω,
k∈Z n
где χω (x) — характеристическая функция области ω.
Приведем теперь алгоритм построения ПКФ с нужными нам свойствами.
Для каждой точки
T x ∈ Ω можно указать окрестность U (x) такую, что попавшая в
нее часть границы Γ U (x) может быть гладко спроектирована на одну из координатных
плоскостей. Так как Ω ограничена, существует ее конечное покрытие {Ut }Tt=1 , Ut = U (x(t) ).
Пусть разбиение единицы {αt }Tt=1 подчинено этому покрытию, т. е.
∀αt ∈
C0∞
supp αt ⊂ Ut ,
T
X
αt (x) = 1.
t=1
T
Введем обозначение: Vt = Ω Ut . Последовательность функционалов погрешности для
всей области Ω получим суммированием по областям Vt локальных ПФП, помноженных
на функции разбиения единицы:
lhΩ (x) =
T
X
t=1
αt (x) lhVt (x),
120
Д. Я. Рахматуллин
где
lhVt (x) := χVt (x) − hn
X
ctk (h)δ(x − hk),
h → 0.
(1)
hk∈Vt ,
k∈Z n
Тогда lhΩ (x) будет иметь коэффициенты cak (h), вычисляемые по формуле
cak (h) =
T
X
αt (hk) ctk (h).
(2)
t=1
Таким образом, если уметь считать коэффициенты локальных ПФП lhVt (x) для любого t,
то задача будет решена. Далее мы приведем формулы, дающие явный вид коэффициентов
локальных ПФП.
Рассмотрим любую из локальных ПФП lhVt (x). Будем обозначать штрихом (n−1)-мерный
вариант того или иного обозначения, например
(x©1 , ..., xn−1 , xn ) = (x′ , xnª); (δ ′ (x′ ), ϕ(x)) =
T
′
ϕ(0 , xn ). Пусть для определенности Γt ≡ Γ Ut = x : xn = γ(x′ ), γ ∈ C M , т. е. последняя
координата любого вектора из множества Γt явно выражается через остальные координаты. Будем предполагать, что xn > γ(x′ ). Тогда для коэффициентов ПФП lhVt (x) можно
записать явное выражение [1]:

min{kn −ξ−1,M
M
+1 1 M
+1
P
P
P +1}

q−1

η
wsq

p=1 p q=1
s=1
ctk =



0,
min{kn −ξ−s,M
P +1}
wrp , kn > 1 + ξ,
r=1
kn 6 1 + ξ.
¸
½
¾
γ(hk ′ )
γ(hk ′ )
γ(hk ′ )
+1
, а {wi,j }M
, η=
— целая и дробная части числа
Здесь ξ =
i,j=1 —
h
h
h
элементы матрицы, обратной к матрице Вандермонда (M + 1) × (M + 1):
·

1

 1

 ...

1
...
1
2
...
M +1
...
...
...
1 2M ... (M + 1)M



.


Построенная таким образом ПФП lhΩ (x) является оптимальной по порядку на любом
пространстве Wpm (Ω) с m < M [1], а значит, по теореме 1, и асимптотически оптимальной.
В частности, для любого m < M выполняется оценка
° Ω°
°lh ° m ∗ 6 C hm , h → 0.
(Wp (Ω))
° °
¡
¢
Так как lhΩ (x), f (x) 6 °lhΩ °(W m (Ω))∗ kf kWpm (Ω) , рассчитанный программой ответ может
p
разниться с точным на величину произведения константы, нормы функции f и значение hm .
Z
Составим программу приближенного вычисления интеграла dxf (x)α(x) с использоVt
ванием формул (1) и (2).
ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ ПО МНОГОМЕРНЫМ ОБЛАСТЯМ ...
121
2. Программа
Для вычислений по изложенному алгоритму написана программа на языке C++ с использованием библиотеки параллельных функций MPI. На данный момент она реализует
основную часть алгоритма: при заданной области Vt с явно определенной границей Γt
она вычисляет интеграл от произведения функции f (x) и срезывающей функции αt (x),
вычисляя коэффициенты ctk локальной ПФП lhVt (x).
Мы предполагаем, что область Vt вложена в два соединенных в направлении n-й координаты единичных гиперкуба, в одном из которых находится гладкая граница Γt . Таким
образом, Vt — вытянутая гиперпризма, в которой одна из плоских границ заменена произвольной гладкой поверхностью.
В качестве параметров программы следует указывать:
— размерность пространства n;
— подынтегральную функцию f (x);
— срезывающую функцию αt (x);
— γ(x′ ) — функцию, задающую поверхность Γt ;
— параметр гладкости M ;
— параметр разбиения h. Задается как N −1 , где N — количество точек решетки на ребре
единичного куба. Общее количество узлов в гиперпризме Q вычисляется как
b := |Q|N n ;
N
— коэффициенты растяжения области Vt по всем координатам;
— количество процессоров P.
Следует отметить, что используемые нами кубатурные формулы являются условноненасыщаемыми. Хотя параметр гладкости M задается заранее, наш алгоритм применим
и для функций, реальная гладкость m которых меньше или больше M . В первом случае
теоретическая погрешность вычислений имеет порядок hm , во втором — hM . Так как наши эксперименты проводились с бесконечно гладкими функциями f , мы имели второй
случай, т. е. гладкость таких функций искусственно понижалась до порядка гладкости M ,
фактически учитываемого в программе. Ясно, что теоретически в этом случае лучше
брать возможно большее M . Однако, так как от параметра M зависит размер матрицы
Вандермонда, с обратной к которой мы имеем дело, при увеличении M усиливаются ее
плохие свойства. Кроме того, поскольку параметр M лежит в основе нескольких вложенных циклов нашей программы, при его увеличении, как будет показано позже, значительно возрастает время вычислений. В связи с этим величина M варьировалась в не очень
широком диапазоне — от 2 до 6.
3. Вычислительные эксперименты
Программа тестировалась при следующих параметрах:
1) n от 2 до 10;
¶
µn
P
bi
ai xi , a = (1, 2, 3, 1, 1..., 1), b = (1, 2, 3, 1, 1..., 1);
2) f (x) = sin
i=1
3) α(x1 , ..., xn ) =
n−1
Q
i=1
ϕ(xi ) ψ(xn ), где ϕ(t) = ξ(2εt)ξ(2ε(1 − t)), ψ(t) = ξ(1 + ε − 2t),
122
Д. Я. Рахматуллин
аε=1и

0, t < 0,


 Rt
Rε
ξ(t) =
z(t)dt/ z(t)dt, z(t) = (t(ε − t))M , 0 6 t < ε,

0

 0
1, ε 6 t;
¶
µn−1
P
1
1
d
4) γ(x′ ) = sin
ci xi i + , где c = (1, 2, 3, 1, 1..., 1), d = (1, 2, 3, 1, 1..., 1);
4
2
i=1
5) M от 2 до 6;
6) N от 5 до 100 000;
7) коэффициенты растяжения области Vt : (1, ...1, 2) — n-мерный вектор;
8) P от 1 до 900.
Заметим, что независимыми параметрами являются лишь n, M, N и P, поэтому в
процессе тестирования следует варьировать лишь их.
Точность вычислений рассчитывалась по устойчивости десятичных знаков в ответе
при использовании последовательности уменьшающегося параметра h. Верными при этом
считались совпадающие десятичные знаки в ответах при текущем и следующем — меньшем значении h (см., например, [3]). В результате экспериментов установлено, что при
фиксированных значениях параметров n, M, N и изменении P точность вычислений не
изменяется, т. е. она не зависит от числа процессоров. Далее, если фиксировать n, M и
увеличивать N, то точность повышается. То же самое происходит при постоянных параметрах n, N и возрастающем M. Например, при n = 2 получаем табл. 1.
Сравнивая при каждом M вышестоящие ответы с нижестоящими, мы нашли верные
(устойчивые) знаки и убрали несовпадающие (округление учитывалось). Количество верных знаков для каждого результата показано в табл. 2.
Теоретически ожидаемые порядки точности представлены в табл. 3. Как видно, при
теоретической точности до 16 знаков практические результаты не менее точны, тогда как
Таблица 1. Результаты экспериментов, n = 2
2
0.03539
0.03539291
0.03539291153
0.035392911532
N \M
100
1000
10 000
100 000
N \M
100
1000
10 000
100 000
3
0.032350
0.0323504440
0.0323504440474
0.03235044404743
5
0.028926864
0.02892686406822
0.028926864068224
0.0289268640682242
Таблица 2. Количество верных знаков
N \M
100
1000
10 000
2
5
8
11
3
6
10
13
4
8
13
15
5
9
14
15
6
10
15
16
4
0.03034590
0.0303459075339
0.030345907533967
0.0303459075339675
6
0.0278703404
0.027870340463228
0.0278703404632276
0.0278703404632276
Таблица 3. Теоретическая точность
N \M
100
1000
10 000
2
4
6
8
3
6
9
12
4
8
12
16
5
10
15
20
6
12
18
24
123
ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ ПО МНОГОМЕРНЫМ ОБЛАСТЯМ ...
далее сказываются ошибки округления. В программе используются числа типа long double
с 15 (в ОС UNIX) значащими цифрами после точки.
Согласно теории, при увеличении размерности n и cохранении величин N и M постоянными погрешность ПФП имеет порядок hm = N −m , т. е. не зависит от размерности
пространства. Однако не нужно забывать, что на практике мы должны учитывать важный
дополнительный параметр — предоставленное для вычислений время. Следует подбирать
число N так, чтобы счет не оказался слишком долгим. Для этого нужно поддерживать
b := |Q|N n , а значит, при увелипостоянным общее количество узлов в гиперпризме Q : N
чении n уменьшать N.
Практика показала, что при максимальном на МВС-15 000 количестве процессоров оптимальное по затратам времени общее количество узлов — 1010 точек. Поэтому, если для
размерности n = 2 мы можем взять до 105 точек в одном измерении, то для n = 10 вынуждены ограничиться значением N = 10. В этом случае пограничный слой кубатурной
формулы, т. е. приграничная часть области, где пересчитываются коэффициенты, сравним по объему со всей областью, что плохо сказывается на точности вычислений. В самом
деле, счет с параметрами n = 10, N = 10, M = 2 дает лишь четыре верных знака.
Теперь проанализируем практическую значимость применения большого количества
процессоров в нашей задаче. Время вычислений при варьировании параметров n, M, N
и числа процессоров P изменялось следующим образом.
При фиксированных n, M, N и увеличивающемся числе процессоров P время вычислений уменьшалось почти пропорционально возрастанию P. Пусть, например, n = 7,
M = 5, N = 15.
Для лучшего анализа введем понятия ускорения и эффективности распараллеливания
программы:
T1
SP =
(ускорение),
TP
SP
(эффективность),
EP =
P
где TP — время, за которое задача выполняется на P процессорах. Будем варьировать P
от 100 до 900 с шагом 100. При этом положим T1 = 100T100 . Дело в том, что реальное
время T1 слишком велико, поэтому фактически идеальным для нас будет случай, когда
P = 100. Имеем табл. 4.
В идеальном (пропорциональном) случае было бы так, как в табл. 5.
Таблица 4. Реальный случай (n = 7, M = 5, N = 15)
P
TP
SP
EP
100
262.0
100.0
1.0
200
135.0
194.07
0.97
300
89.0
294.38
0.98
400
68.0
385.29
0.96
500
53.0
494.34
0.99
600
47.0
557.45
0.93
700
39.0
671.79
0.96
800
34.0
770.59
0.96
900
30.0
873.33
0.97
Таблица 5. Идеальный случай
P
TP
SP
EP
100
262.0
100.0
1.0
200
131.0
200.0
1.0
300
87.33
300.0
1.0
400
65.50
400.0
1.0
500
52.40
500.0
1.0
600
43.67
600.0
1.0
700
37.43
700.0
1.0
800
32.75
800.0
1.0
900
29.11
900.0
1.0
124
Д. Я. Рахматуллин
70
800
60
600
T (c) 50
Sp
400
40
200
30
100 200 300 400 500 600 700 800 900
P
Рис. 1. Зависимость ускорения от числа процессоров.
2
3
4
5
6
M
Рис. 2. Зависимость времени вычислений от
гладкости.
На графике (рис. 1) наглядно показано отклонение экспериментального ускорения SP
(ломаная кривая) от идеального ускорения (прямая).
Зависимость времени вычислений от параметра гладкости M при прочих равных условиях продемонстрируем на примере (рис. 2). Пусть n = 10, N = 5, P = 20. Как видно,
при увеличении M время счета возрастает. Это связано с тем, что в алгоритме активно
используется матрица, обратная к матрице Вандермонда размера (M + 1) × (M + 1).
Список литературы
[1] Соболев С.Л., Васкевич В.Л. Кубатурные формулы. Новосибирск: Изд-во ИМ СО РАН,
1996. 484 с.
[2] Ramazanov M.D. To the Lp -theory of Sobolev formulas // Siberian Advances in Mathematics.
1999. Vol. 9, N 1. P. 99–125.
[3] Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы: Учеб. пособие. М.:
Наука. Гл. ред. физ.-мат. лит., 1987. 600 с.
Поступила в редакцию 23 декабря 2005 г.,
в переработанном виде — 7 февраля 2006 г.
Правила для авторов
<http://www.ict.nsc.ru/mathpub/comp-tech/>
1. Статья должна быть представлена в редакцию в двух экземплярах в виде рукописи, отпечатанной
на одной стороне листа стандартного формата A4 и подписанной авторами, файла рукописи в формате
LATEX 2ε и файлов рисунков на дискете.
2. Все файлы предоставляются на дискете 3.5” формата 1440 Кбайт. Предпочтительнее пересылка
файлов по электронной почте jct@ict.nsc.ru в виде *.zip архива.
3. На отдельной странице на русском и английском языках прилагаются: название статьи, имена
авторов, аннотация (не более 300 знаков) и ключевые слова (в электронном виде — в файле рукописи, в
конце)
4. Статья должна сопровождаться разрешением на опубликование от учреждения, в котором выполнена данная работа. В сопроводительном письме необходимо указать почтовый адрес, телефоны, e-mail
автора, с которым будет проводиться переписка.
5. Для каждого автора должна быть представлена (на русском и английском языках) в виде отдельного файла следующая информация:
◦ Фамилия, имя, отчество
◦ место работы и должность
◦ ученая степень и звание
◦ почтовый адрес
◦ телефоны с кодом города (дом. и служебный), факс, e-mail, URL домашней страницы
◦ область научных интересов (краткое резюме)
6. Материалы следует направлять по адресу: редакция журнала “Вычислительные технологии”, Институт вычислительных технологий СО РАН, просп. Акад. Лаврентьева, 6, 630090, Новосибирск, 90, Россия,
Игорю Алексеевичу Пестунову (отв. секретарь) — тел.: +7(383)3308785, e-mail: jct@ict.nsc.ru; Галине Григорьевне Митиной (зав. РИО).
Рекомендации по оформлению статьи в LATEX
В редакцию следует направлять исходный файл, подготовленный в формате LATEX 2ε в классе jctart
(допускается использование стандартного класса article).
Файл класса jctart.cls можно скачать с сайта ЖВТ: http://www.ict.nsc.ru/win/mathpub/comp-tech/.
1. Структура файла в формате LATEX 2ε :
\documentclass{jctart}
\usepackage{amsmath}
...
\begin{document}
\setcounter{page}{1}\pagestyle{myheadings}
\markboth{<И.О. Фамилия автора(ов)>}{<КРАТКОЕ НАЗВАНИЕ СТАТЬИ (ДО 40 СИМВОЛОВ)>}
\title{<НАЗВАНИЕ СТАТЬИ>\footnote{<Ссылка на поддержку (факультативно)>.}}
\author{\sc{<И.О. Фамилия первого автора>}\\
\it{<Место работы первого автора>}\\
e-mail: \tt{<Адрес первого автора>}\\[2mm]
\sc{<И.О. Фамилия второго автора>}\\
\it{<Место работы второго автора (отличное от первого)>}\\ ...}
\date{}
\maketitle
\begin{abstract}
<Текст аннотации>
\end{abstract}
<Текст статьи>
\begin{thebibliography}{9}
<Библиография (\bibitem-список)>
{\small
\bibitem{} {\sc Иванов~И.И., Иванова~И.И.} К вопросу о вычислительных технологиях //
Вычисл. технологии. 1999. Т.~11, №~11. С.~1123–1135.
...}
\end{thebibliography}
\end{document}
(В конце файла даются:
<Перевод названия статьи на английский язык (или на русский, если статья на английском)>
<аннотации на английский язык (или на русский, если статья на английском)>)
2. Список литературы составляется по ходу упоминания работы в тексте и оформляется по образцу:
Книга
Шокин Ю.И. Метод дифференциального приближения. Новосибирск: Наука, 1979. 222 с.
Бренстед А. Введение в теорию выпуклых многогранников: Пер. с англ. М.: Мир, 1988.
Рояк М.Э., Соловейчик Ю.Г., Шурина Э.П. Сеточные методы решения краевых задач математической физики: Учеб. пособие. Новосибирск: Изд-во НГТУ, 1998.
Finlayson B.A. The Method of Weighted Resuduals and Variational Principles. N.Y.: Acad. Press, 1972.
Книга четырех авторов
Проблемы вычислительной математики / А.Ф. Воеводин, В.В. Остапенко, В.В. Пивоваров, С.М. Шургин. Новосибирск: Изд-во СО РАН, 1995.
Статья из продолжающегося тематического сборника
Федорова А.А., Черных Г.Г. О численном моделировании струйных течений вязкой несжимаемой
жидкости // Моделирование в механике: Сб. науч. тр. / РАН. Сиб. отд-ние. Вычисл. центр. Ин-т теор. и
прикл. механики. 1992. Т. 6 (23). С. 129–140.
Статья из журнала
Игнатьев Н.А. Выбор минимальной конфигурации нейронных сетей // Вычисл. технологии. 2001.
Т. 6, № 1. С. 23–28.
Venkatakrishnan V. Newton solution of inviscid and viscous problems // AIAA J. 1989. Vol. 27, N 7.
P. 285–291.
Труды конференции
Ivanov I.I. Problems in computational techologies // Intern. Conf. Comput. Techs. Novosibirsk, 1988.
P. 225–229.
Препринт
Гуськов А.Е., Федотов А.М., Молородов Ю.И. Информационная система“Конференции”. Новосибирск, 2003 (Препр. РАН. Сиб. отд-ние. ИВТ. № 1–03).
Диссертация
Деменков А.Г. Численное моделирование турбулентных следов в однородной жидкости: Дис. ...
канд. физ.-мат. наук. Новосибирск, 1997. 123 с.
3. Иллюстрации вставляются в текст статьи с помощью команды \includegraphics, например:
\begin{figure}[htbp]
\centering
\includegraphics{fig1.eps}
\caption{<Подрисуночная подпись.>}
\end{figure}
Наиболее предпочтительной формой представления иллюстраций являются файлы рисунков в векторном формате PostScript (.eps) или черно-белых растровых в форматах .bmp, .pcx, .tif с разрешением
300 dpi.
Все надписи на рисунках (обозначение осей и т.п.) должны быть выполнены в том же начертании
(гарнитура “Roman”), что и в тексте статьи. Латинские символы — курсивом, из математической моды
(x[k], z × 10−3 , ψ, P, ...), цифровые обозначения на графиках — наклонно (№ кривой — 1, 2,...), цифры
по осям — прямо (10, 15,...), единицы измерения — по-русски (кг, м,...).
Instructions for Authors
<http://www.ict.nsc.ru/mathpub/comp-tech/>
1. Papers may be submitted to the editorial board as two copies of the manuscript typed on one side of the
standard A4 sheet (297x210 mm) and files of the manuscript in LATEX 2ε format and files of the figures on a
diskette.
2. All files should be submitted on a 3.5"floppy disc (1440 Kbytes) or sent to jct@ict.nsc.ru as a
*.zip - archive.
3. A separate page should contain a title, names of the authors, an abstract (not more than 300 characters)
and keywords.
4. The paper should be accompanied by the publication permission from the organization, where the work
was done. The enclosed letter should contain the postal address, phone numbers and e-mail of the corresponding
author.
5. A separate file should contain the following information on each author:
◦ First name, second name, last name
◦ Affiliation, position
◦ Academic degree and title
◦ Postal address
◦ Office and home phone numbers (including area code), fax number, e-mail address, homepage URL
◦ Scientific interests (brief curriculum vitae)
6. All materials should be sent to the following address: Dr. Igor A. Pestunov (executive secretary), Journal
of Computational Technologies, Institute of Computational Technologies SB RAS, Academician Lavrentyev
Ave. 6, Novosibirsk, 630090, Russia. Phone +7(383)3308785, E-mail: jct@ict.nsc.ru; Galina G. Mitina (publishing
department manager).
Recommendations on submitting paper in LATEX
The source file should be submitted in LATEX 2ε format using jctart class file (standard article class can also
be used).
The files of appropriate jctart.cls class file can be downloaded from JCT web site:
http://www.ict.nsc.ru/win/mathpub/comp-tech/.
1. The file structure in LATEX 2ε format:
\documentclass[english]{jctart}
\usepackage{amsmath}
...
\begin{document}
\setcounter{page}{1}\pagestyle{myheadings}
\markboth{<Name(s) of author(s)>}{<SHORT TITLE (LESS THAN 40 CHARACTERS)>}
\title{<TITLE OF PAPER>\footnote{<Reference to supporting organization (optional)>.}}
\author{\sc{<Name of the first author>}\\
\it{<Affiliation of the first author>}\\
e-mail: \tt{<Address of the first author>}\\[2mm]
\sc{<Name of the second author>}\\
\it{<Affiliation of the second author>}\\ ...}
\date{}
\maketitle
\begin{abstract}
<Abstract>
\end{abstract}
<Text of paper>
\begin{thebibliography}{9}
<References (\bibitem-list)>
{\small
\bibitem{} {\sc Ivanov~I.I., Ivanova~I.I.} On computational technologies //
Computational technologies. 1999. Vol.~11, No.~11. P.~1123–1135.
...}
\end{thebibliography}
\end{document}
2. A list of the references should be sorted according to the order of citations in the text and it should
be written as in the following example:
Book
Finlayson B.A. The method of weighted residuals and variational principles. N.Y.: Acad. Press, 1972.
Book by four authors
Problems of computational mathematics / A.F. Voevodin, V.V. Ostapenko, V.V. Pivovarov, S.M. Shurgin.
Novosibirsk: SB RAS Publishing House, 1995.
Paper from continued subject transactions
Fedorova A.A., Chernykh G.G. On numerical modelling of viscous incompressible jet fluid flows //
Modelling in mechamics: Scientific transactions / RAS. Siberian branch. Computing Center. Institute of Theoretical and Applied Mechanics. 1992. Vol. 6 (23). P. 129–140.
Paper from journal
Venkatakrishnan V. Newton solution of inviscid and viscous problems // AIAA J. 1989. Vol. 27, N 7.
P. 285–291.
Conference proceedings
Ivanov I.I. Problems in computational techologies // Intern. Conf. Comput. Techs. Novosibirsk, 1988.
P. 225–229.
Dissertation
Demenkov A.G. Numerical modelling of turbulent wakes in homogeneous fluid: Dissertation for degree of
candidate of physical and mathematical sciences. Novosibirsk, 1997. 123 p.
3. Figures should be included into the text using command \includegraphics{<figure file name>}, for
example:
\begin{figure}[htbp]
\centering
\includegraphics{fig1.eps}
\caption{<Figure caption.>}
\end{figure}
The preferred presentation form for illustrations is a figure file in vector format PostScript (.eps) or black
and white bitmap formats .pcx, .bmp, .tif with 300 dpi resolution.
All figure inscriptions (axes definitions, etc.) should be done by the same font as in the text of paper (“Roman” type family). Latin characters should be done in italics in mathematical mode (x[k], z × 10 −3 , ψ, P, ...),
figures on axes — by straight font.
In papers, which are written in Russian, the units of measurement should be written in Russian.
В ближайших номерах/Forthcoming papers
Liseikin V.D., Likhanova Yu.V., Patrakhin D.V., Vaseva I.A. Generation of block
structured smooth grids
Лисейкин В.Д., Лиханова Ю.В., Патрахин Д.В., Васева И.А. Построение блочных, гладких сеток
Ащепков Л.Т., Давыдов Д.В. Показатель интервального неравенства: свойства
и применение
Aschepkov L.T., Davydov D.V. Numerical characterization for interval inequalities: properties
and applications
Башуров В.В. Применение методов геометрической оптики к решению задач
безопасности объекта
Bashurov V.V. Application of optics methods to security problems
Горелов Д.Н., Редреев Д.Г. Построение квадратурной формулы для сингулярного интеграла с ядром Коши по контуру крылового профиля
Gorelov D.N., Redreev D.G. Construction of a quadrature formula for a singular integral
with the Cauchy kernel over wing profile contour
Данаев Н.Т., Ергалиев Е.К. Об одном итерационном методе решения стационарных уравнений Навье — Стокса
Danaev N.T., Ergaliev E.K. On an iteration method for stationary Navie — Stocks equations
Картошкина А.Е. Влияние динамики на термодиффузию в плоском слое со
свободными границами
Kartoshkina A.E. Influence of dynamics on thermodiffusion in a 2D layer with free boundaries
Клочков Ю.В., Николаев А.П., Джабраилов А.Ш. К вопросу о неадекватности
изопараметрической параметризации в МКЭ
Klochkov J.V., Nikolaev А.P., Dzhabrailov A.Sh. On inadequacy of isoparametrical
parametrization in finite element method
Новиков Е.А., Шорников Ю.В. Контроль устойчивости метода Фельберга седьмого порядка точности
Novikov E.A., Shornikov Yu.V. Stability control of Felberg method of the 7th order of
accuracy
Червов В.В. Моделирование трехмерной конвекции в мантии Земли с применением неявного метода расщепления по физическим процессам
Chervov V.V. Modelling of 3-D convection in the Earth mantle using an implicit splitting
method over the physical processes
Чикина Л.Г. Двухпараметрические итерационные методы
Chikina L.G. Two-parameter iterative methods
Шапеев В.П., Черепанов А.Н. Конечно-разностный алгоритм для численного
моделирования процесса лазерной сварки металлических пластин
Shapeev V.P., Cherepanov A.N. Finite-difference algorithm for numerical simulation of
process of metal plates laser welding
Карбановский В.В., Каиров Т.В. Замечания к статье В.В. Альчикова “Решение
уравнений магнитостатики для ферромагнетиков различной формы”
Karbanovskiy V.V., Kairov T.V. Remarks on the article “A solutions of magnetostatic
equations for ferromagnetics of various shape”
Документ
Категория
Без категории
Просмотров
4
Размер файла
324 Кб
Теги
интеграл, областям, многопроцессорных, вычисления, система, вычислительной, многомерная
1/--страниц
Пожаловаться на содержимое документа