Перевод кода с C++

@i_love_es · Регистрация: 22.05.2015

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

Доброго времени суток! помоги пожалуйста переделать код задачи из С++ в С

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
#include <iostream>
#include <omp.h>
 
#define M 5000
 
#define E 0.00001 
#define T 0.01 
 
 
int main(){
    int  i, j, v;
    long int dt1;
    double alf, bet, mf;
    double Spr, Spr1, Spz;
//  double A[M][M], F[M], Xk[M], Zk[M], Rk[M], Sz[M];
    double ** A = new double*[M];
    double * F = new double[M];
    double * Xk = new double[M];
    double * Zk = new double[M];
    double * Rk = new double[M];
    double * Sz = new double[M];
 
    omp_set_num_threads(2);
 
    double ts = omp_get_wtime();
 
    for (int i=0;i<M;i++)
        A[i] = new double[M];
 
    for (mf=0,i = 0; i < M; i++) {
        #pragma omp parallel
        {
            #pragma omp for
            for (j = 0; j < M; j++) {
                if(i == j)
                    A[i][j] = 2.0;
                else
                    A[i][j] = 8.;
            }
        }
        F[i] = M + 1;
        mf += F[i] * F[i];
    }
 
    /* Задаем начальное приближение корней. В Хk хранятся значения корней
    * к-й итерации. */
    #pragma omp parallel
    {
        #pragma omp for
        for (i = 0; i < M; i++)
            Xk[i] = 0.2;
    }
/* Задаем начальное значение r0 и z0. */
    for (i = 0; i < M; i++) {
        Sz[i]=0;
        double sum = 0;
        #pragma omp parallel
        {
            #pragma omp for reduction(+:sum)
            for (int j = 0; j < M; j++) {
                sum += A[i][j] * Xk[j];
            }
        }
        Sz[i] = sum;
        Rk[i] = F[i] - Sz[i];
        Zk[i] = Rk[i];  
    }
 
 
    do {
        Spz = 0;
        Spr = 0;
        for(i = 0; i < M; i++) {
            Sz[i]=0;
            double sum = 0;
            #pragma omp parallel
            {
                #pragma omp for reduction(+:sum)
                for(int j = 0; j < M; j++)
                    sum += A[i][j] * Zk[j];
            }
            Sz[i] = sum;
            Spz += Sz[i] * Zk[i];
            Spr += Rk[i] * Rk[i];
        }
        alf = Spr/Spz;
    
        Spr1 = 0;
        #pragma omp parallel
        {
            #pragma omp for
            for(int i = 0; i < M; i++){
                Xk[i] += alf*Zk[i];
                Rk[i] -= alf*Sz[i];
                #pragma omp critical
                Spr1 += Rk[i]*Rk[i];
            }
        }
        bet = Spr1/Spr;
    
        #pragma omp parallel
        {
            #pragma omp for
            for(i = 0; i < M; i++)
                Zk[i] = Rk[i] + bet*Zk[i];
        }
    
    } while(Spr1/mf > E*E);
 
 
//  Xk[0]...Xk[M-1] - корни
 
    double tf = omp_get_wtime();
    std::cout<<"\nruntime: "<<tf - ts<<" s.\n";
    std::cout<<Xk[0]<<" "<<Xk[1];
    return(0);
}

HighPredator · 26.05.2015, 11:27

Что-то вроде такого:

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
#include <stdio.h>
#include <stdlib.h>
#include <omp.h>
 
#define M 5000
 
#define E 0.00001 
#define T 0.01 
 
 
int main(void){
    int  i, j, v;
    long int dt1;
    double sum, tf, ts;
    double alf, bet, mf;
    double Spr, Spr1, Spz;
//  double A[M][M], F[M], Xk[M], Zk[M], Rk[M], Sz[M];
    double ** A = malloc(sizeof(double*) * M);
    double * F = malloc(sizeof(double) * M);
    double * Xk = malloc(sizeof(double) * M);
    double * Zk = malloc(sizeof(double) * M);
    double * Rk = malloc(sizeof(double) * M);
    double * Sz = malloc(sizeof(double) * M);
 
    omp_set_num_threads(2);
 
    ts = omp_get_wtime();
 
    for (i=0;i<M;i++)
        A[i] = malloc(sizeof(double) * M);
 
    for (mf=0,i = 0; i < M; i++) {
        #pragma omp parallel
        {
            #pragma omp for
            for (j = 0; j < M; j++) {
                if(i == j)
                    A[i][j] = 2.0;
                else
                    A[i][j] = 8.;
            }
        }
        F[i] = M + 1;
        mf += F[i] * F[i];
    }
 
    /* Задаем начальное приближение корней. В Хk хранятся значения корней
    * к-й итерации. */
    #pragma omp parallel
    {
        #pragma omp for
        for (i = 0; i < M; i++)
            Xk[i] = 0.2;
    }
/* Задаем начальное значение r0 и z0. */
    for (i = 0; i < M; i++) {
        Sz[i]=0;
        sum = 0;
        #pragma omp parallel
        {
            #pragma omp for reduction(+:sum)
            for (j = 0; j < M; j++) {
                sum += A[i][j] * Xk[j];
            }
        }
        Sz[i] = sum;
        Rk[i] = F[i] - Sz[i];
        Zk[i] = Rk[i];  
    }
 
 
    do {
        Spz = 0;
        Spr = 0;
        for(i = 0; i < M; i++) {
            Sz[i]=0;
            sum = 0;
            #pragma omp parallel
            {
                #pragma omp for reduction(+:sum)
                for(j = 0; j < M; j++)
                    sum += A[i][j] * Zk[j];
            }
            Sz[i] = sum;
            Spz += Sz[i] * Zk[i];
            Spr += Rk[i] * Rk[i];
        }
        alf = Spr/Spz;
    
        Spr1 = 0;
        #pragma omp parallel
        {
            #pragma omp for
            for(i = 0; i < M; i++){
                Xk[i] += alf*Zk[i];
                Rk[i] -= alf*Sz[i];
                #pragma omp critical
                Spr1 += Rk[i]*Rk[i];
            }
        }
        bet = Spr1/Spr;
    
        #pragma omp parallel
        {
            #pragma omp for
            for(i = 0; i < M; i++)
                Zk[i] = Rk[i] + bet*Zk[i];
        }
    
    } while(Spr1/mf > E*E);
 
 
//  Xk[0]...Xk[M-1] - корни
 
    tf = omp_get_wtime();
    printf("\nruntime: %f s.\n", tf - ts);
    printf("%f %f", Xk[0], Xk[1]);
    return 0;
}

Новые блоги и статьи Все статьи Все блоги /
Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .

Перевод кода с C++

Решение