Форум программистов, компьютерный форум, киберфорум
palva
Войти
Регистрация
Восстановить пароль
Оценить эту запись

Пример рекурсии в фортране

Запись от palva размещена 01.01.2014 в 19:46
Обновил(-а) palva 03.01.2014 в 23:58

Современный фортран допускает рекурсивные подпрограммы и функции. Вот пример вычисления определителя при помощи разложения по строке. Определитель всегда разлагается по верхней строке, при этом вычеркивается верхняя строка и один из столбцов. Номера невычернутых столбцов передаются на нижний уровень рекурсии в массиве ia. Первоначально этот массив содержит номера всех столбцов матрицы. Параметр k содержит размер массива ia. Номер первой невычеркнутой строки равен n-k+1. Именно по этой строке будет сделано разложение на очередном уровне рекурсии.
Fortran
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
! Д. К. Фаддеев, И. С. Соминский.
! Сборник задач по высшей алгебре. 1977.
! № 279. Вычислить определитель.
! Ответ: 900.
 
    integer, parameter :: n=4
    real*8 :: a(n,n) = &
! Обратите внимание, что матрица в фортране записывается по столбцам.
    (/ 1, -2,  3,  4, &
       2,  1, -4,  3, &
       3, -4, -1, -2, &
       4,  3,  2, -1  /)
    integer :: ia(n) = (/ (i, i=1,n) /)
    real*8 det
    print *, det(a, n, ia, n)
    end
    
    recursive real*8 function det(a, n, ia, k)
!    Функция вычисляет минор матрицы a размера n x n, стоящий на пересечении
!    k последних строк и k столбцов, номера которых перечислены в массиве ia.
    integer k,n
    real*8, dimension (n,n) :: a
    integer, dimension(k) :: ia
    integer i, j, sign
    integer, dimension(k-1) :: iar
    real*8 :: s
    if(k==1) then
        det=a(n,ia(1))
        return
    end if
    if(k==n) ia = (/ (i, i=1,n) /)
    s=0.0
    sign = 1
    do i = 1,k
        iar=(/ (ia(j), j=1,i-1), (ia(j), j=i+1,k) /)
        s = s + sign * a(n-k+1,ia(i)) * det(a, n, iar, k-1)
        sign = -sign
    end do
    det = s
    end function det
Размещено в Без категории
Просмотров 3613 Комментарии 0
Всего комментариев 0
Комментарии
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.