Форум программистов, компьютерный форум, киберфорум
Наши страницы
C для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.80/5: Рейтинг темы: голосов - 5, средняя оценка - 4.80
i_love_es
0 / 0 / 0
Регистрация: 22.05.2015
Сообщений: 1
1

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

22.05.2015, 07:14. Просмотров 870. Ответов 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
#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);
}
0
Лучшие ответы (1)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
22.05.2015, 07:14
Ответы с готовыми решениями:

Перевод кода с C++
помогите пожалуйста перевести #include &lt;iostream&gt; #include &lt;conio.h&gt; using...

Перевод кода из C++
помогите переписать с С++ на С int main(int argc, char* argv) { int...

Перевод кода из с++ в с
Господа, помогите перевести из с++ в си фрагменты кода, очень надеюсь) //...

Перевод кода с C++
Помогите,перевести,пожалуйста. #include &lt;iostream&gt; #define SIZE 10 int...

Перевод кода с C++ в C
Уже долго маюсь с этим кодом, вроде ничего сложного, но все равно не могу...

1
HighPredator
5692 / 2011 / 723
Регистрация: 10.12.2010
Сообщений: 5,778
Записей в блоге: 3
26.05.2015, 11:27 2
Лучший ответ Сообщение было отмечено Памирыч как решение

Решение

Что-то вроде такого:
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;
}
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
26.05.2015, 11:27

Перевод кода С++ на Си
Здравствуйте! Прошу помощи в переводе данного кода с языка Cpp на язык Си. Буду...

Перевод кода на С
#include &lt;stdio.h&gt; #include &lt;time.h&gt; #include &lt;stdlib.h&gt; int main(int...

Перевод кода из С++ в С
Помогите пожалуйста перевести код из С++ и С #include &lt;iostream&gt; using...


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

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

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