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

UptoLike

ВУЗ: 

ТГУ | Томск

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

Рубрика: 

186
c Определение числа элементов, обрабатываемых одним
c процессорным элементом
m=Ny/size
If (rank.lt.(Nx-size*m)) then
m=m+1
End if
c Задание начальных значений
Do i=0,Nx+1
Do j=0,m+1
T(i,j)=100.0d0
Tnew(i,j)=100.0d0
End do
End do
c Задание граничных условий
If (rank.eq.0) then
Do i=0,Nx+1
T(i,0)=10.0d0
End do
End if
If (rank.eq.size-1) then
Do i=0,Nx+1
T(i,m+1)=10.0d0
End do
End if
Do j=0,m+1
T(0,j)=1.0d0
T(Nx+1,j)=10.0d0
End do
c Подготовка к расчетам
tau = 0.2*(hx*hx*hy*hy)/(alfa*(hx**2+hy**2))
ae=alfa*tau/(hx*hx)
aw=alfa*tau/(hx*hx)
an=alfa*tau/(hy*hy)
as=alfa*tau/(hy*hy)
ap=-(ae+aw+an+as)
time = 0
If (rank.eq.0) then
open (2,file='rezult1.dat')
timeStart = MPI_Wtime()
timeStop1=0.0

Страницы