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

UptoLike

164
Call VecView(u, PETSC_VIEWER_STDOUT_WORLD, ierr)
c Уничтожаем все использовавшиеся PETSc объекты
Call KSPDestroy(ksp, ierr)
Call VecDestroy(u, ierr)
Call VecDestroy(x, ierr)
Call VecDestroy(b, ierr)
Call MatDestroy(A, ierr)
c Завершаем работу PETSc
Call PetscFinalize(ierr)
end
Запустим программу на 4 процессорах, используя разные методы:
mpirun -np 4 ./laplas -ksp_type gmres
Norm of error 0.3097E-04,iterations 10
mpirun -np 4 ./laplas -ksp_type bcgs
Norm of error 0.5074E-04,iterations 6
mpirun -np 4 ./laplas -ksp_type cg
Norm of error 0.3140E-04,iterations 10
8.7 Задания
Решить задачу Дирихле для уравнения Пуассона. Численное ре-
шение получить в области G с границей Г, которая задается соот-
ветствующим образом в каждом варианте.
1.
2 2
2 2
6
u u
x y
x y
;
( , ) ( , ), 0 1, 0 1 ;
x y G x y x y
3
(0, )
u y y
,
y
;
3
( ,0)
u x x
,
0 1
x
;
3
(1, ) 1
u y y
,
y
;
3
,1 1
u x x
,
0 1
x
.
2.
2 2
2 2
2
x y
u u
e
x y
;
( , ) ( , ), 0 0.5, 0 1 ;
x y G x y x y
3
(0, )
u y y
,
y
; ( ,0)
x
u x e
,
0 0.5
x
;
(0.5, )
y
u y e e
,
y
;
1
,1
x
u x e
,
0 0.5
x
.