Форум программистов, компьютерный форум, киберфорум
Fortran
Войти
Регистрация
Восстановить пароль
 
0 / 0 / 0
Регистрация: 16.01.2021
Сообщений: 1
1

Прошу Вас перевести 2-а немаленьких кода с C++ на fortran90

16.01.2021, 15:27. Просмотров 253. Ответов 1

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
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <imsl.h>
#include <malloc.h>
#define nu2 4
double F_r(double);
double K_r(double);
double Q_r(double);
int main()
{
 int N = 20, i, k;
 double* C, * B, * r, * rh, * A, * phi, * v, * dv, *vt, dvmax[5]; // массивы С и В -
верхняя и главная диагонали матрицы А, r - узлы, rh - вспомогательные узлы, vt - u* точное
решение, v - решение, которое ищем с помощью функции, dv - их разность
 double a = 0, b = 2;
 double h;
 FILE* outf;
 outf = fopen("output1.txt", "wt");
 for (i=0; i < 5; i++)
 {
 fprintf(outf, "N=%i\n", N);
 h = (b - a) / N;
 C = (double*)malloc(N * sizeof(double));
 phi = (double*)malloc(N * sizeof(double));
 B = (double*)malloc(N * sizeof(double));
 r = (double*)malloc((N + 1) * sizeof(double));
 rh = (double*)malloc(N * sizeof(double));
 A = (double*)malloc((2 * N) * sizeof(double));
 v = (double*)malloc(N * sizeof(double));
 dv = (double*)malloc(N * sizeof(double));
 vt = (double*)malloc(N * sizeof(double));
 for (k = 0; k < (N+1); k++)
 {
 r[k] = k * h;
 }
 for (k = 0; k < N; k++)
 {
 vt[k] = -r[k]*r[k]+8; // vt=u*=-r^2+8
 /*vt[k] = -r[k]*r[k]*r[k]+12;*/ // для теста с ненулевой погрешностью
аппроксимации
 }
 for (k = 0; k < N; k++)
 {
 rh[k] = (r[k + 1] + r[k]) / 2;
 }
 B[0] = 0;
 C[0] = (rh[0] * K_r(rh[0])) / h + (rh[0] * h * Q_r(r[0])) / 4;
 phi[0] = rh[0] * h * F_r(r[0])/ 4 ;
 for (k = 1; k < N; k++)
 {
 B[k] = -rh[k-1] * K_r(rh[k-1]) / h;
 C[k] = (rh[k] * K_r(rh[k]) + rh[k - 1] * K_r(rh[k - 1])) / h + r[k] * h *
Q_r(r[k]);
16
 phi[k] = F_r(r[k]) * r[k] * h;
 }
 phi[N - 1] = F_r(r[N - 1]) * r[N - 1] * h + (rh[N - 1] * K_r(rh[N - 1]) * nu2) / h;
 for (k = 0; k < (2 * N); k++)
 {
 if (k < N)
 A[k] = B[k];
 else if (k < 2 * N)
 A[k] = C[k - N];
 }
 imsl_d_lin_sol_posdef_band(N, A, 1, phi, IMSL_RETURN_USER, v, 0);
 for (k = 0; k < N; k++)
 {
 dv[k] = fabs(v[k]-vt[k]);
 printf("%e\n", dv[k]);
 }
 dvmax[i] = dv[0];
 for (k = 1; k < N; k++)
 {
 if (dv[k] > dvmax[i]) dvmax[i] = dv[k];
 }
 fprintf(outf, "raznost=%e\n", dvmax[i]);
 if (i!=0) fprintf(outf, "otnoshenie=%e\n", dvmax[i-1]/dvmax[i]);
 
 free(C);
 free(B);
 free(r);
 free(rh);
 free(A);
 free(v);
 free(dv);
 N = N * 2;
 }
 return 0;
}
double F_r(double r)
{
 double F;
 F = -r*r+12;
 /*F = -pow(r,4)-pow(r,3)+21*r + 12;*/ // это для задачи с ненулевой погрешностью
аппроксимации
 return F;
}
double K_r(double r)
{
 double K;
 return 1;
}
double Q_r(double r)
{
 double Q;
 Q = r + 1; // это для задачи с ненулевой погрешностью аппроксимации
 return 1;
}
 
 
=============================================================
 
 
 
#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <imsl.h>
#include <malloc.h>
#define nu2 4
double F_r(double);
double K_r(double);
double Q_r(double);
int main()
{
 int N = 20, i, k;
 double* C, * B, * r, * rh, * A, * phi, * v; // массивы С и В - верхняя и главна
диагонали матрицы А, r - узлы, rh - вспомогательные узлы, vt - u* точное решение, v -
решение, которое ищем с помощью функции, dv - их разность
 double a = 0, b = 2;
 double h;
 FILE* outf;
 outf = fopen("output.txt", "wt");
 for (i = 0; i < 3; i++)
 {
 fprintf(outf, "N=%i\n", N);
 h = (b - a) / N;
 C = (double*)malloc(N * sizeof(double));
 phi = (double*)malloc(N * sizeof(double));
 B = (double*)malloc(N * sizeof(double));
 r = (double*)malloc((N + 1) * sizeof(double));
 rh = (double*)malloc(N * sizeof(double));
 A = (double*)malloc((2 * N) * sizeof(double));
 v = (double*)malloc(N * sizeof(double));
 for (k = 0; k < (N + 1); k++)
 {
 r[k] = k * h;
 }
 for (k = 0; k < N; k++)
 {
 rh[k] = (r[k + 1] + r[k]) / 2;
 }
 B[0] = 0;
 C[0] = (rh[0] * K_r(rh[0])) / h + (rh[0] * h * Q_r(r[0])) / 4;
 phi[0] = rh[0] * h * F_r(r[0]) / 4;
 for (k = 1; k < N; k++)
 {
 B[k] = -rh[k - 1] * K_r(rh[k - 1]) / h;
 C[k] = (rh[k] * K_r(rh[k]) + rh[k - 1] * K_r(rh[k - 1])) / h + r[k] * h *
Q_r(r[k]);
 phi[k] = F_r(r[k]) * r[k] * h;
 }
 phi[N - 1] = F_r(r[N - 1]) * r[N - 1] * h + (rh[N - 1] * K_r(rh[N - 1]) * nu2) / h;
 for (k = 0; k < (2 * N); k++)
 {
 if (k < N)
 A[k] = B[k];
 else if (k < 2 * N)
 A[k] = C[k - N];
 }
18
 imsl_d_lin_sol_posdef_band(N, A, 1, phi, IMSL_RETURN_USER, v, 0);
 for (k = 0; k < N; )
 {
 fprintf(outf, "%e\n", v[k]);
 if (i == 2) k = k + i + 2;
 else k = k + i + 1;
 }
 free(C);
 free(B);
 free(r);
 free(rh);
 free(A);
 free(v);
 N = N * 2;
 }
 return 0;
}
double F_r(double r)
{
 double F;
 F = pow(r, 4) - 25 * r * r - 15 * r + 12;
 return F;
}
double K_r(double r)
{
 double K;
 K = 2 * r + 1;
 return K;
}
double Q_r(double r)
{
 double Q;
 Q = r + 1;
 return Q;
}
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
16.01.2021, 15:27
Ответы с готовыми решениями:

Прошу Вас подсказать по 2-м небольшим фрагментам кода
Друзья, помогите, пожалуйста мне вот с таким кодом: 1) array = ; function...

Матлаб.С 3 по 10 задание.Прошу выполнить до Вторника,а то отчислят с универа.Прошу вас!
я Скинул 2 фотки.Там задания Правила: 4.3. Создавайте темы с осмысленными и понятными...

Прошу вас разобраться
// KR_6.cpp: определяет точку входа для консольного приложения. // Информация об участниках...

Помогите по массивам, прошу вас.
Записать элементы массива Х=(х1, х2, …, х15) в массив Y=(y1, y2, ..., y15), сдвинув элементы...

1
5423 / 3959 / 1557
Регистрация: 02.02.2014
Сообщений: 10,628
16.01.2021, 15:48 2
Цитата Сообщение от tworog Посмотреть сообщение
немаленьких
ключевое слово...
вопрос: кто и почему этим должен заниматься?
да еще разбираться в коде без постановки самой задачи...
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
16.01.2021, 15:48

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь.

Вычислите значение выражения! Прошу Вас!
((|x-5|-sinx)/3)+(√x^2+2014)cos2x-3 при x=-2.34. Ответ: -1.76911 Как это расписать? Помогите...

Прошу у вас помощи! с доработкой программы!
Здравствуйте! Тема: Редактор блок-схем. У меня есть исходник работающей программы которая только...

Прошу Вас, посмотреть конфигурацию моего РС
Прошу Вас, посмотреть конфигурацию моего РС (согласно гарантийного талона) Покупал 19.07.2001г.,...

Курсовая работа!прошу ответить вас)
вобщем вот я взяла курсовую работу по теме &quot;Система тестирования&quot; это по дисциплине...

Уважаемые сайтостроители, прошу вас о помощи!
Я делаю сайт просмотра онлайн фильмов и делаю полностью с нуля...подскажите сначала мне плиз вот...

Прошу Вас посодействовать в решении небольшой задачи
Уважаемые коллеги! Прошу Вас помочь разобраться в задаче: Необходимо сделать так, чтобы текст...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.