Высокопроизводительные вычисления на кластерах. Беликов Д.А - 75 стр.

UptoLike

ВУЗ: 

ТГУ | Томск

Составители: 

Рубрика: 

75
менить любую квадратурную формулу, например формулу трапе-
ций.
Положим
, , 0,1,...,
i
b a
h x a i h i n
n
и
0 1 2 1
( ) ( 2 2 2 )
2
b
n n
a
h
I F x dx F F F F F
,
где
( ) ( , ) , 0,1,...,
d
i i i
c
.
В конечном итоге поставленная задача сводится к вычислению
( 1)
n
-интеграла
( , )
d
i i
c
F f x y dy
.
Описанный метод можно применить для криволинейной области.
Для этого запишем интеграл (4.1) в виде
2 2
( , ) , 0, , 0, sin , ( , ) 3 sin ( ),
b d
a c
I dx f x y dy a b c d x f x y y x
здесь верхний предел интегрирования
d
зависит от
x
.
4.4 Параллельная программа расчета двойного интеграла
методом повторного применения квадратурной формулы
трапеции
Program povtor
C Вычисление кратного интеграла с помощью повторного
C применения квадратурной формулы трапеции
Implicit None
Include 'mpif.h'
Integer Size, Rank, Ierr, I,k, N,m
Double Precision GSum,s,s1,eps,df,pi,
$ time, a1, b1,hh, h1,a,b,c,d,x,y,f
C Задание числа интервалов разбиения по x, пределов
С интегрирования и точности вычисления однократного
С интеграла
Parameter (n=100000, pi=3.14159265358979d0, a=0d0, b=pi,
$ c=0d0, eps=1d-9)

Страницы