Расчетное задание / КМ3 / ТФ-11-19 Рыжков М.В. Вариант 16

IPVT2, RJ, IPATH2, AINV(1,J), RES) RJ(J) = 0.0E0 20 CONTINUE

!Вывод результатов

!CALL WRRRN ('AINV', N, N, AINV, N, 0)

! 99999 FORMAT (' RCOND = ',F5.3,/,' L1 Condition number = ',F6.3)

! Решение исходного уравнения

CALL LSLRB (N, A, LDA, NLCA3, NUCA3, B, IPATH3, X) ! Вывод результатов

!CALL WRRRN ('X', 1, N, X, 1, 0) END DO

T=Timef()

PRINT *,' TIME: ',T end program one

21

10. Листинг программы двойной точности

program double implicit none

INTEGER*4 IPATH, LDA, LDFAC, N, NLCA, NUCA, IPATH2, LDFAC2, NLCA2, NUCA2, NOUT, IPATH3, NLCA3, NUCA3, j, k, i

PARAMETER (IPATH=1, LDA=7, LDFAC=10, N=8, NLCA=3, NUCA=3, IPATH2=1, LDFAC2=10, NLCA2=3, NUCA2=3, IPATH3=1, NLCA3=3, NUCA3=3)

REAL T,Timef

EXTERNAL DLFTRB, DLFSRB, DWRRRN, DLFCRB, DLFIRB, DSET, UMACH, DLSLRB T=Timef()

do i=1,11000000 do j=1,8

p(i,j) = i+j END DO

END DO

do i = 1,10500000,7 do j = 1,8,1

A(1,j)=p(i,j)

A(2,j)=p(i+1,j)

A(3,j)=p(i+2,j)

A(4,j)=p(i+3,j)

A(5,j)=p(i+4,j)

A(6,j)=p(i+5,j)

A(7,j)=p(i+6,j) END DO

!Зададим правую часть линейной системы

DO k=1,N

DO j=1,4

B(k,j)= j+k END DO

END DO

!решение исходного уравнения

CALL DLFTRB (N, A, LDA, NLCA, NUCA, FAC, LDFAC, IPVT)

!решение для правых сторон

DO 10 J=1, 4

CALL DLFSRB (N, FAC, LDFAC, NLCA, NUCA, IPVT, B(1,J), IPATH,

X(1,J))

10 CONTINUE

!вывод результатов

!CALL DWRRRN ('X', N, 2, X, N, 0) !решение исходного уравнения

CALL DLFCRB (N, A, LDA, NLCA2, NUCA2, FAC2, LDFAC2, IPVT2, RCOND)

!Вывод оценки числа обусловленности

!CALL UMACH (2, NOUT)

!WRITE (NOUT,99999) RCOND, 1.0E0/RCOND

CALL DSET (N, 0.0E0, RJ, 1)

DO 20 J=1, N

22

RJ(J) = 1.0E0

CALL DLFIRB (N, A, LDA, NLCA2, NUCA2, FAC2, LDFAC2, IPVT2, RJ, IPATH2, AINV(1,J), RES)

RJ(J) = 0.0E0 20 CONTINUE

!вывод результатов

!CALL DWRRRN ('AINV', N, N, AINV, N, 0)

!99999 FORMAT (' RCOND = ',F5.3,/,' L1 Condition number = ',F6.3) !решение исходного уравнения

CALL DLSLRB (N, A, LDA, NLCA3, NUCA3, B, IPATH3, X)

!вывод результатов

!CALL DWRRRN ('X', 1, N, X, 1, 0)

END DO

T=Timef()

PRINT *,' TIME: ',T

end program double

11.Результаты измерений времени выполнения расчета

Параметры системы:

Операционная система: Windows 10 Корпоративная LTSC (10.0, Build 17763) Тип системы: 64-разрядная ОС, процессор х64

Модель ПК: Lenovo G70-80 (LENOVO 80FF)

Процессор: Название: Intel(R) Core(TM) i3-5020U CPU @ 2.20GHz

Установленная память (ОЗУ): 16,00 Гб Количество ядер: 2 Количество логических процессоров: 4

12.Список литературы

1.Ю. Б. Воробьев Методические указания по проведению лабораторных работ по дисциплине "Компьютерные и сетевые технологии ЯЭ".

2.IMSL. Math Library. Volumes 1 and 2.

3.О. В. Бартеньев «Фортран для профессионалов. Математическая библиотека

IMSL»

4.Л. С. Шихобалов «Матрицы и определители».

28