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

UptoLike

ВУЗ: 

ТГУ | Томск

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

Рубрика: 

118
30 continue
do 50 j = 1,n
a(j,j)=0.0
do 40 i = 1,n
! задание элементов на главной диагонали матрицы
if(i.ne.j) a(j,j) = a(j,j) + abs(a(i,j))*1.1d0
40 continue
50 continue
return
end
6.6 OpenMP-программа решения СЛАУ методом Якоби
Данная программа является записью представленного выше ал-
горитма с использованием стандарта OpenMP.
Program Jacobi
implicit none
!
integer
$ i, j, m_size, count1
!
double precision
$ summ, eps, max_err,
$ time1, dr, drmax
!
parameter
$ (m_size = 5000, !размер матрицы
$ eps = 0.000005) !точность вычисления
!
double precision
$ a(m_size, m_size+1),
$ x(m_size), xx(m_size)
!
!$ integer omp_get_num_threads
!$ double precision omp_get_wtime, tm0
!
! фиксируется время начала счета
!$ tm0=omp_get_wtime()

Страницы