Переписать программу с С на fortran

@Olga1994 · Регистрация: 12.03.2013

Студворк — интернет-сервис помощи студентам

C++
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
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
#include <stdio. h>
#include <stdlib. h>
#include <math. h>
#include <mpi. h>
#define N1 50
#define N2 50
#define eps 0.00001
double Y [N1 + 1] [N2 + 1], Ysh [N1 + 1] [N2 + 1];
double hx = 1. f / N1, hy = 2. f / N2;
/*Функия точного решения*/
double Fresh (double x, double y) {
return pow (x,2) * pow (y,2);
}
double RoFresh (double x, double y) {
return 2 * (pow (x,2) + pow (y,2));
}
/*Подпрограмма инициализации матрицы*/
void Inic () {
int i, j;
for (i = 0; i < N1 + 1; i++)
for (j = 0; j < N2 + 1; j++) {
if ( (i! = 0) && (j! = 0) && (i! = N1) && (j! = N2))
Y [i] [j] = 0;
else
Y [i] [j] = Fresh ( (i * hx), (j * hy));
}
}
int main (int argc, char **argv) {
int size, rank, flag = 1;
int i, j, f, it = 0;
double A = pow (hx,2),
B = A, 
D = pow (hy,2),
 C = - 2. f / A - 2. f / D, F1, Fi, pogr = 0;
double t_c = 0.0, time = 0.0, s_t = 0.0;
double max, m;
double alfa [N-1], beta [N-1];
// double Y [N1 + 1] [N2 + 1], Ysh [N1 + 1] [N2 + 1];
// double hx = 1. f / N1, hy = 2. f / N2;
MPI_Status stat;
MPI_Init (&argc, &argv);
MPI_Comm_size (MPI_COMM_WORLD, &size);
MPI_Comm_rank (MPI_COMM_WORLD, &rank);
time - = MPI_Wtime ();
t_c - = MPI_Wtime ();
if (rank == 0) {
printf ("%d \n", size);
Inic ();
}
MPI_Barrier (MPI_COMM_WORLD);
MPI_Bcast (Y, (N1 + 1) * (N2 + 1), MPI_DOUBLE, 0, MPI_COMM_WORLD);
if (rank == 0) {
do {
for (i = 0; i <= N1; i++)
for (j = 0; j <= N2; j++)
Ysh [i] [j] = Y [i] [j];
for (i = 1; i <= (N1 - 1) / size; i++) {
for (j = 1; j < N2; j++) {
Fi = ( - (Y [i+1] [j] + Y [i] [j+1])) / D + RoFresh;
}
}
for (i = 0; i <= (N - 1); i++) {
alfa += ( - D [i]) / (C [i] + Ai*alfa [i-1]);
beta += (F [i] - A [i] *beta [i-1]) / (C [i] + Ai*alfa [i-1]);
}
for (j = (N+1); j <= 1; j--)
Ysh += alfa [i] *Y [i+1] [j] + beta [i];
if (size == 1) {
do {
for (i = 1; i < N1; i++) {
for (j = 1; j < N2; j++) {
Fi = ( - (Y [i+1] [j] + Y [i] [j+1])) / D + RoFresh;
}
}
max = m = - 999;
for (i = 0; i <= N1; i++) {
for (j = 0; j <= N2; j++) {
pogr = fabs (Ysh [i] [j] - Y [i] [j]);
if (pogr > max) max = pogr;
pogr = fabs (Ysh [i] [j]);
if (pogr > m) m = pogr;
}
}
if (max / m < eps) {
t_c += MPI_Wtime ();
printf ("Bce 4etKo!!!! it = %d, time = %f\n", it,t_c);
flag = 0;
}
// Перезапись данных
for (i = 1; i < N1; i++)
for (j = 1; j < N2; j++)
Y [i] [j] = Ysh [i] [j];
it++;
} while (flag);
exit (1);
} else {
s_t - = MPI_Wtime ();
// Посылка следующему процессору верхнюю строку из основной матрицы в нижнюю границу
MPI_Send (&Ysh, (N1 + 1) * (N2 + 1), MPI_DOUBLE, 1, 1, MPI_COMM_WORLD);
MPI_Recv (&Y, (N1 + 1) * (N2 + 1), MPI_DOUBLE, 1, 1, MPI_COMM_WORLD, &stat);
s_t += MPI_Wtime ();
}
} while (1);
}
if ( (rank > 0) && (rank < size - 1)) {
do {
MPI_Recv (&Ysh, (N1 + 1) * (N2 + 1), MPI_DOUBLE, rank - 1, 1, MPI_COMM_WORLD, &stat);
for (i = (N1 - 1) / size * rank + 1; i < (N1 - 1) / size * (rank + 1) + 1; i++) {
for (j = 1; j < N2; j++) {
Fi = ( - (Y [i+1] [j] + Y [i] [j+1])) / D + RoFresh;
}
}
s_t - = MPI_Wtime ();
MPI_Send (&Ysh, (N1 + 1) * (N2 + 1), MPI_DOUBLE, rank - 1, 1, MPI_COMM_WORLD);
MPI_Send (&Ysh, (N1 + 1) * (N2 + 1), MPI_DOUBLE, rank + 1, 1, MPI_COMM_WORLD);
MPI_Recv (&Y, (N1 + 1) * (N2 + 1), MPI_DOUBLE, rank + 1, 1, MPI_COMM_WORLD, &stat);
s_t += MPI_Wtime ();
} while (1);
}
if ( (rank == size - 1) && (rank! = 0)) {
do {
s_t - = MPI_Wtime ();
MPI_Recv (&Ysh, (N1 + 1) * (N2 + 1), MPI_DOUBLE, rank - 1, 1, MPI_COMM_WORLD, &stat);
s_t += MPI_Wtime ();
// for (i = (N1-1) / size * (size - 1) + 1; i < N1; i++) {
for (i = (N1-1) / size * rank + 1; i < N1; i++) {
for (j = 1; j < N2; j++) {
Fi = ( - (Y [i+1] [j] + Y [i] [j+1])) / D + RoFresh;
}
}
max = m = - 999;
for (i = 0; i <= N1; i++) {
for (j = 0; j <= N2; j++) {
pogr = fabs (Ysh [i] [j] - Y [i] [j]);
if (pogr > max) max = pogr;
pogr = fabs (Ysh [i] [j]);
if (pogr > m) m = pogr;
}
}
if (max / m < eps) {
time += MPI_Wtime ();
printf ("Bce 4etKo!!!! it = %d, time = %f, s_t = %f, time_s4eta = %f \n", it, time, s_t, time-s_t);
exit (1);
}
s_t - = MPI_Wtime ();
MPI_Send (&Ysh, (N1 + 1) * (N2 + 1), MPI_DOUBLE, rank - 1, 1, MPI_COMM_WORLD);
s_t += MPI_Wtime ();
// Перезапись данных
for (i = 0; i <= N1; i++)
for (j = 0; j <= N2; j++)
Y [i] [j] = Ysh [i] [j];
it++;
} while (1);
}
MPI_Finalize ();
return 0;
}

Плохо знаю фортран, не могу разобраться как правильно написать массив задала и переменные переписала, а дальше ступор...Помогите , ПОЖАЛУЙСТА, бедной студентке

Добавлено через 4 минуты
Решение двумерного уравнения Пуассона итерационным методом Зейделя
Найти численное решение задачи Дирихле для уравнения Пуассона

@vua72 · 10.04.2015, 07:47

Одномерный массив: real(8), dimension (0:n+1) :: y, как варианты real*8 y(0:n+1), double precision y(0:n+1) - это если компилятор старый (до фортрана 90).
Двумерный массив double Y [N1 + 1] [N2 + 1] - real(8), dimension (0:N1 + 1, 0:N2 + 1) :: y.
Если не писать левую границу (0

, то индексация с 1 по умолчанию.
Если у вас компилятор поддерживает фортран 90 и выше, вашу задачу просто решить или через внутренние подпрограммы, или с использованием модулей. Одно и другое дает возможность рассматривать массивы как глобальные.
Плюс посмотрите http://www.opennet.ru/docs/RUS/MPI_intro/ - это поможет переписать участки кода с MPI.
Начало вашей программы могло бы быть

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
41
42
43
44
45
46
47
48
49
program puason
!include 'mpi. h'
integer, parameter :: N1=50            !#define N1 50
integer, parameter :: N2=50             !#define N2 50
real(8), parameter :: eps=1d-5        !#define eps 0.00001
real(8), dimension (0:N2, 0:N2):: Y 
real(8), dimension (0:N1, 0:N1):: Ysh 
real(8) ::hx = 1_8 / N1, hy = 1_8 / N2
! здесь идет основная программа
    y=0; ysh=0
    call Inic()
    call print_matrix()
    
contains
  
       !/*Функия точного решения*/
    real(8) function Fresh (x, y) 
        real(8), intent (in) :: x, y
        Fresh=x**2*y**2
    end function Fresh
    real(8) function RoFresh (x, y) 
        real(8), intent (in) :: x, y
        RoFresh=2*(x**2+y**2)
    end function RoFresh
    !/*Подпрограмма инициализации матрицы*/
    subroutine Inic()
        implicit none
        integer :: i,j
        do i=0, N1
            do j=0, N2
                if((i/= 0).and.(j/= 0).and.(i/= N1).and.(j/= N2)) then
                    y(i,j)=0
                else
                    y(i,j)=Fresh(i*hx, j*hy)
                end if
            enddo
        enddo
    end subroutine Inic
    subroutine print_matrix()
        integer :: i,j
        do i=0, N1
            do j=0, N2
                print('(e12.4,2x,$)'), y(i,j)
            enddo
            print*
        enddo
        
    end subroutine print_matrix
end program

Новые блоги и статьи Все статьи Все блоги /
Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .	My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .	Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .
Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .	Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .