МКЕ (перевод Fortran -> C++)

@мухобой · Регистрация: 02.04.2015

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

Имеется код на фортране, нужна помощь в переводе на с++

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
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
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
    implicit real*8 (a-h,o-z)
    parameter (NKU=13, NKE=20)
 
    dimension AK(4,4),ISSG(4),A(2*NKU,2*NKU),A1(2*NKU,2*NKU),B(2*NKU),
     * X(2*NKU),RE(2*NKU), XY(2,NKU),NC(2,NKE),NUZST(NKU),IPZST(2,NKU),
     * NZUZ(NKU), F(2,NKU), A2(2*NKU,2*NKU+1), B1(2*NKU), ST(NKE)
 
    Open (8,file='DATA.TXT')
    READ (8,*) KU,KUZST,KNU
    DO i=1,KU
        Read (8,*) XY(1,i),XY(2,i)
    End Do
    DO i=1,KUZST
        Read (8,*) NUZST(i), IPZST(1,i), IPZST(2,i)
    End Do
    DO i=1,KNU
        Read (8,*) NZUZ(I),F(1,I),F(2,I)
    End Do
    Read (8,*) KE
    DO  I=1,KE
        Read(8,*) NC(1,i),NC(2,i)
    End Do
    Read (8,*) E, AREA
    Read (8,*) UVEL
    CLOSE(8)
    
    Open (6, file='Picture')
    DO i=1,KU
        write(6,10) i,XY(1,i),XY(2,i)
    End Do
   10 format(i5, F15.4,F15.4,'   0.0')
    write(6,11)
   11 format('FACES')
    DO i=1,KE
        write(6,12) NC(1,i), NC(2,i)
    End Do
   12 Format(i8,i8,'.')
      CLOSE(6)
    
    KSS=2*KU
 
    DO i=1,KSS
        B(i)=0.0
        RE(i)=0.0
        DO j=1,KSS
            A(i,j)=0.0
        End Do
    End Do
    
    DO IJ=1,KE
 
        N1=NC(1,IJ)
        N2=NC(2,IJ)
        ALENGHT=SQRT( (XY(1,N2)-XY(1,N1))**2+(XY(2,N2)-XY(2,N1))**2)
        C=(XY(1,N2)-XY(1,N1))/ALENGHT
        S=(XY(2,N2)-XY(2,N1))/ALENGHT
 
        AK(1,1)=C*C
        AK(1,2)=S*C
        AK(1,3)=-C*C
        AK(1,4)=-S*C
        AK(2,1)=S*C
        AK(2,2)=S*S
        AK(2,3)=-S*C
        AK(2,4)=-S*S
        AK(3,1)=-C*C
        AK(3,2)=-S*C
        AK(3,3)=C*C
        AK(3,4)=S*C
        AK(4,1)=-S*C
        AK(4,2)=-S*S
        AK(4,3)=S*C
        AK(4,4)=S*S
 
        DO I=1,4
            DO j=1,4
                AK(I,J)=AK(I,J)*E*AREA/ALENGHT
            END DO
        END DO
 
        ISSG(1)=NC(1,IJ)*2-1
        ISSG(2)=NC(1,IJ)*2
        ISSG(3)=NC(2,IJ)*2-1
        ISSG(4)=NC(2,IJ)*2
 
        Do i=1,4
            Do j=1,4
                IG=ISSG(i)
                JG=ISSG(j)
                A(IG,JG)=A(IG,JG)+AK(I,J)
            End Do
        End Do 
    End Do
 
    Do IJ=1,KNU
        I1=NZUZ(IJ)*2-1
        I2=NZUZ(IJ)*2
        B(I1)=F(1,IJ)
        B(I2)=F(2,IJ)
    End Do
    
    DO i=1,KSS
        DO j=1,KSS
            A1(i,j)=A(i,j)
        End Do
        B1(i)=B(i)
    End Do
 
    Do IJ=1,KUZST
        I1=NUZST(IJ)*2-1
        I2=NUZST(IJ)*2
        IF (IPZST(1,IJ).EQ.1) THEN
            DO I=1,KSS
                A(I1,I)=0.0
            End Do
            DO I=1,KSS
                A(I,I1)=0.0
            End Do
            B(I1)=0.0
            A(I1,I1)=1
        END IF
        IF (IPZST(2,IJ).EQ.1) THEN
            DO I=1,KSS
                A(I2,I)=0.0
            END DO
            DO I=1,KSS
                A(I,I2)=0.0
            END DO
            B(I2)=0.0
            A(I2,I2)=1
        END IF
    END DO
 
    
      N=KSS
    DO i=1,N
        A2(i,N+1)=B(i)
        Do j=1,N
            A2(i,j)=A(i,j)
        End Do
    End Do
    N1=N+1
    DO 25 K=1,N
    K1=K+1
    S=A2(K,K)
    J=K
    DO 21 I=K1,N
    R=A2(I,K)
    IF (ABS(R).LE.ABS(S)) GOTO 21
    S=R
    J=I
   21   Continue
    IF (S.EQ.0) GOTO 99
    IF (J.EQ.K) GOTO 23
    DO 22 I=K,N1
    R=A2(K,I)
    A2(K,I)=A2(J,I)
   22   A2(J,I)=R
   23   DO 24 J=K1,N1
   24   A2(K,J)=A2(K,J)/S
    DO 25 I=K1,N
    R=A2(I,K)
    DO 25 J=K1,N1
   25   A2(I,J)=A2(I,J)-A2(K,J)*R
 
    X(N)=A2(N,N1)
    DO 27 I=N-1,1,-1
    S=A2(I,N1)
    DO 26 J=I+1,N
   26   S=S-A2(I,J)*X(J)
   27   X(I)=S
 
    Open (11,file='solve.txt')
    Do i=1,KSS
    Write(11,13) X(i) 
      End Do
    Close(11)   
   13 format(F15.6,'  ')    
   99 Continue
    Open (6, file='DPicture')
    DO i=1,KU
        write(6,100) i,XY(1,i)+Uvel*X(2*i-1),XY(2,i)+Uvel*X(2*i)
    End Do
  100 format(i5, F15.4,F15.4,'   0.0')
    write(6,111)
  111 format('FACES')
    DO i=1,KE
        write(6,112) NC(1,i), NC(2,i)
    End Do
  112 Format(i8,i8,'.')
      CLOSE(6)
 
    Do i=1,KSS
        s = 0
        Do j=1,KSS
            s = s + A1(i,j)*X(j)
        End Do
        RE(i) = s - B1(i)
    End Do
    Open (12, file='reactions.txt')
    Do i=1,KSS
        write(12,13) RE(i)
    End Do
    Close(12)
    
    Do IJ=1,KE
        
        N1=NC(1,IJ)
        N2=NC(2,IJ)
        ALENGHT=SQRT( (XY(1,N2)-XY(1,N1))**2+(XY(2,N2)-XY(2,N1))**2)
        C=(XY(1,N2)-XY(1,N1))/ALENGHT
        S=(XY(2,N2)-XY(2,N1))/ALENGHT
        
        IX=NC(1,IJ)*2-1
        IY=NC(1,IJ)*2
        JX=NC(2,IJ)*2-1
        JY=NC(2,IJ)*2
        
        ELENGHT=(X(JX)-X(IX))*C + (X(JY)-X(IY))*S
 
        ST(IJ)=ELENGHT/ALENGHT*E
    End do
 
    Open (14, file='stress.txt')
    Do i=1,KE
        write(14,130) ST(i)
    End Do
    Close(14)
  130 format(F20.6,'  ')
    STOP
    END

@мухобой · 01.06.2017, 20:50 **[ТС]**

Помогите пожалуйста

Добавлено через 15 минут
Сложно я так понимаю))

@мухобой · 05.06.2017, 19:21 **[ТС]**

(

MrGluck · 05.06.2017, 19:24

мухобой, здесь мало мастеров фортрана, вы бы лучше задание выложили

@мухобой · 05.06.2017, 19:26 **[ТС]**

Нужно написать программу которая считала мке ферму. Вот пример кода на Фортране, ее по сути нужно в с++ перевести

@afront · 05.06.2017, 20:06

мухобой, давно бы сам все переписал, элементарная программа. Нужно заменить чтение файла и фортранные циклы DO от 1 до N заменить на c++ типа for (int i = 0; i < N; i++).
А что будет не понятно, конкретно спрашивайте

@мухобой · 06.06.2017, 13:09 **[ТС]**

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
210
211
212
 #include <iostream>
 #include <fstream>
 #include <math.h>
 using namespace std;
 const int NKU=13, NKE=20
 float AK(4,4),ISSG(4),A(2*NKU,2*NKU),A1(2*NKU,2*NKU),B(2*NKU), X(2*NKU),RE(2*NKU), XY(2,NKU),NC(2,NKE),NUZST(NKU),IPZST(2,NKU),NZUZ(NKU), F(2,NKU), A2(2*NKU,2*NKU+1), B1(2*NKU), ST(NKE);
 ifstream read("DATA.TXT");
  READ >> KU,KUZST,KNU;
 for (int i=0, i<KU; i++){
    Read >> XY[1][i],XY[2][i];
    }
    for (int i=0, i<KUZST; i++){
    Read >> NUZST[i], IPZST[1][i], IPZST[2][i];
    }
    for (int i=0, i<KNU; i++){
    Read >> NZUZ[I],F[1][I],F[2][I];
    }   
    Read >>KE;
    for (int i=0, i<KE; i++){
    Read >> NC[1][i],NC[2][i];
    }
    read.close();
    ifstream read("Picture.TXT");
    for (int i=0, i<KU; i++){
    write<< setw(8)<<i,XY[1][i],XY[2][i];
    }
    for (int i=0, i<KE; i++){
    write<< setw(8)<<NC[1][i], NC[2][i];
    }
    read.close();
        KSS=2*KU;
    for (int i=0, i<KSS; i++){
        B[i]=0.0;
        RE[i]=0.0;
        for (int j=0, j<KSS; j++){
        
            A[i][j]=0.0;
    }
}
        for (int IJ=0, IJ<KE; IJ++){
1 = NC[1, IJ];
N2 = NC[2, IJ];
ALENGHT = sqrt((XY[1, N2] - XY[1, N1])**2 + (XY[2, N2] - XY[2, N1])**2);
C = (XY[1, N2] - XY[1, N1])/ALENGHT;
S = (XY[2, N2] - XY[2, N1])/ALENGHT;
 
AK[1, 1] = C*C;
AK[1, 2] = S*C;
AK[1, 3] = -C*C;
AK[1, 4] = -S*C;
AK[2, 1] = S*C;
AK[2, 2] = S*S;
AK[2, 3] = -S*C;
AK[2, 4] = -S*S;
AK[3, 1] = -C*C;
AK[3, 2] = -S*C;
AK[3, 3] = C*C;
AK[3, 4] = S*C;
AK[4, 1] = -S*C;
AK[4, 2] = -S*S;
AK[4, 3] = S*C;
AK[4, 4] = S*S;
for (int i=0, i<4; i++){
for (int j=0, j<4; j++){
 
AK[I, J] = AK[I, J]*E*AREA/ALENGHT;     
}
}
ISSG[1] = NC[1, IJ]*2 - 1;
ISSG[2] = NC[1, IJ]*2;
ISSG[3] = NC[2, IJ]*2 - 1;
ISSG[4] = NC[2, IJ]*2;
 
for (int i=0, i<4; i++){
for (int j=0, j<4; j++){
    
    IG = ISSG[i];
JG = ISSG[j];
A[IG, JG] = A[IG, JG] + AK[I, J];
}
}
for (int IJ=0, IJ<KNU; IJ++){
    I1 = NZUZ[IJ]*2 - 1;
I2 = NZUZ[IJ]*2;
B[I1] = F[1, IJ];
B[I2] = F[2, IJ];
}
for (int i=0, i<KSS; i++){
    for (int j=0, j<KSS; j++){
        A1[i, j] = A[i, j];
    }
    B1[i] = B[i];
}
for (int IJ=0, IJ<KUZST; IJ++){
    I1 = NUZST[IJ]*2 - 1;
I2 = NUZST[IJ]*2;
if (IPZST[1, IJ] == 1){
    for (int i=0, i<KSS; i++){
        A[I1, I] = 0.0;
}
B[I1] = 0.0;
A[I1, I1] = 1;
}
if (IPZST[2, IJ] == 1) 
{
        for (int i=0, i<KSS; i++){
            A[I2, I] = 0.0;
        }
        for (int i=0, i<KSS; i++){
            A[I, I2] = 0.0;
        }
        B[I2] = 0.0;
        A[I2, I2] = 1;
}
}
N=KSS;
for (int i=0, i<n; i++){
    A2[i, N + 1] = B[i];
    for (int j=0, j<n; j++){
        A2[i, j] = A[i, j];
    }
}
N1=N+1;
for (int k=0, k<n; k++){
    K1=K+1;
    S=A2[K,K];
    J=K;
    for (int i=K1, i<n; i++){
        R = A2[I, K];
     if (ABS(R) <= ABS(S)) goto g21;
        S=R;
        J=I;
     g21: else continue; 
    }
    if (S == 0) goto g99;
if (J == K) goto g23;
for(I=K; I<=N1; I++)
{
R = A2[K, I];
A2[K, I] = A2[J, I];
A2[J, I] = R;
}
g23:
for(J=K1; J<=N1; J++)
{
A2[K, J] = A2[K, J]/S;
}
for(I=K1; I<=N; I++)
{
R = A2[I, K];
for(J=K1; J<=N1; J++)
{
A2[I, J] = A2[I, J] - A2[K, J]*R;
}
}
}
 
X[N] = A2[N, N1];
for( I = N - 1, 1, -1)\n{
S = A2[I, N1];
for( J = I + 1, N)\n{
S = S - A2[I, J]*X[J];
}
X[I] = S;
}
ofstream read("solve.txt");
for (int i=1, i<KSS; i++){
    write << X[i];
}
read.close();
g99:
    else continue;
        ifstream read("DPicture.TXT");
        for (int i=1, i<KU; i++){
        
        write << i, XY[1, i] + Uvel*X[2*i - 1], XY[2, i] + Uvel*X[2*i];
}
for (int i=1, i<KE; i++){
write << NC[1, i], NC[2, i];
}
read.close();
for (int i=1, i<KSS; i++){
    s=0;
for (int j=1, j<KSS; j++){  
s = s + A1[i, j]*X[j];
}
RE[i] = s - B1[i];
}
ofstream read("reactions.txt");
for (int i=1, i<KSS; i++){
write << Re[i];
}
for (int IJ=0, IJ<KE; IJ++){
    N1 = NC[1, IJ];
N2 = NC[2, IJ];
ALENGHT = sqrt((XY[1, N2] - XY[1, N1])**2 + (XY[2, N2] - XY[2, N1])**2);
C = (XY[1, N2] - XY[1, N1])/ALENGHT;
S = (XY[2, N2] - XY[2, N1])/ALENGHT;
IX = NC[1, IJ]*2 - 1;
IY = NC[1, IJ]*2;
JX = NC[2, IJ]*2 - 1;
JY = NC[2, IJ]*2;
ELENGHT = (X[JX] - X[IX])*C + (X[JY] - X[IY])*S;
 
ST[IJ] = ELENGHT/ALENGHT*E;
}
ofstream read("stress.txt");
for (int i=1, i<KE; i++){
    write << ST[i];
}
read.close();
}

Добавлено через 54 секунды
Попробовал переписать, но не работает... не знаток программирования

@afront · 06.06.2017, 13:25

Сообщение от мухобой

не знаток программирования

спасение утопающих дело рук самих утопающих и не тешьте себя иллюзиями что кто то будет разбираться с вашим кодом

@мухобой · 06.06.2017, 13:31 **[ТС]**

Смысл тогда всего этого если подсказать даже не могут?

@Antikl · 06.06.2017, 13:36

Сообщение от мухобой

Смысл тогда всего этого

ну не каждый же во всем разбираться плюс люди тоже учатся

@afront · 06.06.2017, 13:43

Сообщение от мухобой

подсказать даже не могут?

вы не правы, очень даже подсказывают, но как правило не конкретные вопросы, а лазить по чужому говнокоду никому удовольствия не доставляет. Добейтесь сначала что бы ваш код компилировался, затем в дебагере сравните ход выполнения фортранного и с++ кода.

@karaulov6 · 06.06.2017, 15:33

Это ищите раздел где понимают этот набор букв и цифр ))

может там https://www.cyberforum.ru/fortran/

@nmcf · 06.06.2017, 19:54

Откуда там запятые?

C++
1
Read >> NUZST[i] >> IPZST[1][i] >> IPZST[2][i];

@мухобой · 06.06.2017, 20:20 **[ТС]**

Отсутствуют вроде как

Добавлено через 7 минут
Исправил везде где есть запятые

@afront · 06.06.2017, 20:43

исправил несколько строк в начале, дальше лень, в фортране индексация массивов с 1 а в с++ с 0 и это нужно везде смотреть

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
#include <iostream>
#include <fstream>
#include <math.h>
using namespace std;
int main()
{
const int NKU=13, NKE=20;
float AK[4][4],ISSG[4],A[2*NKU][2*NKU],A1[2*NKU][2*NKU],B[2*NKU], X[2*NKU],RE[2*NKU], XY[2][NKU],NC[2][NKE],NUZST[NKU],IPZST[2][NKU],NZUZ[NKU], F[2][NKU], A2[2*NKU][2*NKU+1], B1[2*NKU], ST[NKE];
int KU=0,KUZST,KNU;
ifstream read ("DATA.TXT");
read  >> KU >> KUZST >> KNU;
for (int i=0; i<KU; i++){
    read >> XY[0][i],XY[1][i];
}
for (int i=0; i<KUZST; i++){
    read >> NUZST[i], IPZST[0][i], IPZST[1][i];
}
for (int i=0; i<KNU; i++){
    read >> NZUZ[i],F[0][i],F[1][i];
}
int KE;
read >>KE;
for (int i=0; i<KE; i++){
    read >> NC[0][i],NC[1][i];
}
read.close();

Добавлено через 1 минуту
да еще в фортране массивы real*8 а вы поставили float

@мухобой · 06.06.2017, 20:55 **[ТС]**

И ещё у меня сомнения где есть операторы перехода go to и метки, если не затруднит посмотрите, пожалуйста)

@afront · 08.06.2017, 08:39

завтра

Добавлено через 44 минуты
кое что сегодня

описание массивов вначале нужно разбить на double и int

C++
1
2
double AK[4][4],A[2*NKU][2*NKU],A1[2*NKU][2*NKU],B[2*NKU], X[2*NKU],RE[2*NKU], XY[2][NKU],F[2][NKU], A2[2*NKU][2*NKU+1], B1[2*NKU], ST[NKE];
int ISSG[4],NC[2][NKE],NUZST[NKU],IPZST[2][NKU],NZUZ[NKU];

переходы по меткам в с++ такие же как в фортране поэтому пишите просто

C++
1
 g99:  continue;

и нужно было все оставить в верхнем регистре как в фортране, меньше путаницы
да и в вашем фортране есть несколько символов n в нижнем регистре, фортран не делает разницы между верхним и нижнем регистрами а с++ делает

Добавлено через 12 часов 52 минуты

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
210
211
212
213
214
215
216
217
218
219
220
221
222
#include <iostream>
#include <fstream>
#include <math.h>
using namespace std;
int main()
{
    const int nku=13, nke=20;
    double ak[4][4],a[2*nku][2*nku],a1[2*nku][2*nku],b[2*nku], x[2*nku],re[2*nku], xy[2][nku],f[2][nku], a2[2*nku][2*nku+1], b1[2*nku], st[nke];
    int issg[4],nc[2][nke],nuzst[nku],ipzst[2][nku],nzuz[nku];
    int ku,kuzst,knu;
 
    ifstream read ("data.txt");
    read  >> ku >> kuzst >> knu;
    for (int i=0; i<ku; i++) 
        read >> xy[0][i] >> xy[1][i];
     
    for (int i=0; i<kuzst; i++) 
        read >> nuzst[i] >> ipzst[0][i] >> ipzst[1][i];
     
    for (int i=0; i<knu; i++) 
        read >> nzuz[i] >> f[0][i] >> f[1][i];
     
    int ke;
    read >> ke;
    for (int i=0; i<ke; i++) 
        read >> nc[0][i] >> nc[1][i];
     
    double e,area,uvel;
    read >> e >> area >> uvel; 
    read.close();
 
    ofstream write("picture.txt");
    for (int i=0; i<ku; i++) 
        write <<  i << xy[0][i] << xy[1][i];
     
    for (int i=0; i<ke; i++) 
        write << nc[0][i] << nc[1][i];
     
    write.close();
 
    int kss=2*ku;
    for (int i=0; i<kss; i++){
        b[i]=0.0;
        re[i]=0.0;
        for (int j=0; j<kss; j++)
            a[i][j]=0.0;
    }
 
    for (int ij=0; ij<ke; ij++){
        int n1 = nc[0][ij];
        int n2 = nc[1][ij];
        double alenght = sqrt((xy[0][n2] - xy[0][n1])*(xy[0][n2] - xy[0][n1]) + (xy[1][n2] - xy[1][n1])*(xy[1][n2] - xy[2][n1]));
        double c = (xy[0][n2] - xy[0][n1])/alenght;
        double s = (xy[1][n2] - xy[1][n1])/alenght;
 
        ak[0][0] = c*c;
        ak[0][1] = s*c;
        ak[0][2] = -c*c;
        ak[0][3] = -s*c;
        ak[1][0] = s*c;
        ak[1][1] = s*s;
        ak[1][2] = -s*c;
        ak[1][3] = -s*s;
        ak[2][0] = -c*c;
        ak[2][1] = -s*c;
        ak[2][2] = c*c;
        ak[2][3] = s*c;
        ak[3][0] = -s*c;
        ak[3][1] = -s*s;
        ak[3][2] = s*c;
        ak[3][3] = s*s;
 
        for (int i=0; i<4; i++)
            for (int j=0; j<4; j++)
                ak[i][j] = ak[i][j]*e*area/alenght;     
 
        issg[0] = nc[1][ij]*2 - 1;
        issg[1] = nc[1][ij]*2;
        issg[2] = nc[2][ij]*2 - 1;
        issg[3] = nc[2][ij]*2;
 
        for (int i=0; i<4; i++)
            for (int j=0; j<4; j++){
                int ig = issg[i];
                int jg = issg[j];
                a[ig][jg] = a[ig][jg] + ak[i][j];
            }
        
        for (int ij=0; ij<knu; ij++){
            int i1 = nzuz[ij]*2 - 1;
            int i2 = nzuz[ij]*2;
            b[i1] = f[0][ij];
            b[i2] = f[1][ij];
        }
 
        for (int i=0; i<kss; i++){
            for (int j=0; j<kss; j++)
                a1[i][j] = a[i][j];
            b1[i] = b[i];
        }
        for (int ij=0; ij<kuzst; ij++){
            int i1 = nuzst[ij]*2 - 1;
            int i2 = nuzst[ij]*2;
            if (ipzst[1][ij] == 1){
                for (int i=0; i<kss; i++) 
                    a[i1][i] = 0.0;
                 
                b[i1] = 0.0;
                a[i1][i1] = 1;
            }
            if (ipzst[2][ij] == 1) 
            {
                for (int i=0; i<kss; i++) 
                    a[i2][i] = 0.0;
                 
                for (int i=0; i<kss; i++) 
                    a[i][i2] = 0.0;
                 
                b[i2] = 0.0;
                a[i2][i2] = 1;
            }
        }
        int n=kss;
        for (int i=0; i<n; i++){
            a2[i][n + 1] = b[i];
            for (int j=0; j<n; j++) 
                a2[i][j] = a[i][j];
             
        }
        n1=n+1;
        for (int k=0; k<n; k++){
            int k1=k+1;
            s=a2[k][k];
            int j=k;
            for (int i=k1; i<n; i++){
                double r = a2[i][k];
                if (abs(r) <= abs(s)) goto g21;
                s=r;
                j=i;
g21: continue; 
            }
            if (s == 0) goto g99;
            if (j == k) goto g23;
            for(int i=k; i<=n1; i++)
            {
                double r = a2[k][i];
                a2[k][i] = a2[j][i];
                a2[j][i] = r;
            }
g23:
            for(j=k1; j<=n1; j++)
                a2[k][j] = a2[k][j]/s;
 
            for(int i=k1; i<=n; i++)
            {
                double r = a2[i][k];
                for(int j=k1; j<=n1; j++)
                    a2[i][j] = a2[i][j] - a2[k][j]*r;
            }
        }
 
        x[n] = a2[n][n1];
        for( int i = n - 1; 1; -1){
            s = a2[i][n1];
            for( int j = i + 1; n;j++)
                s = s - a2[i][j]*x[j];
 
            x[i] = s;
        }
        {
            ofstream write1("solve.txt");
            for (int i=1; i<kss; i++) 
                write1 << x[i];
             
            write1.close();
        }
g99:    continue;
 
        ofstream write2("dpicture.txt");
        for (int i=0; i<ku; i++){
            write2 << i << xy[0][i] + uvel * x[2 * i - 1] << xy[1][i] + uvel * x[2 * i];
        }
        for (int i=0; i<ke; i++){
            write2 << nc[0][i], nc[2][i];
        }
        write2.close();
 
        for (int i=0; i<kss; i++){
            double s=0;
            for (int j=1; j<kss; j++)   
                s = s + a1[i][j]*x[j];
             
            re[i] = s - b1[i];
        }
 
        ofstream write3("reactions.txt");
        for (int i=0; i < kss; i++) 
            write3 << re[i];
        write3.close();
 
        for (int ij=0; ij<ke; ij++){
            n1 = nc[0][ij];
            n2 = nc[1][ij];
            alenght = sqrt(xy[0][n2] - xy[0][n1])*(xy[0][n2] - xy[0][n1]) + (xy[1][n2] - xy[1][n1])*(xy[1][n2] - xy[1][n1]);
            c = (xy[0][n2] - xy[0][n1])/alenght;
            s = (xy[1][n2] - xy[1][n1])/alenght;
            int ix = nc[0][ij]*2 - 1;
            int iy = nc[0][ij]*2;
            int jx = nc[1][ij]*2 - 1;
            int jy = nc[1][ij]*2;
            double elenght = (x[jx] - x[ix])*c + (x[jy] - x[iy])*s;
 
            st[ij] = elenght/alenght*e;
        }
 
        ofstream write4("stress.txt");
        for (int i=0; i<ke; i++) 
            write4 << st[i];
         
        write4.close();
    }
}

Добавлено через 21 час 36 минут
мухобой, если вы скините вариант исходных данных попроще я попробую отладить

@мухобой · 09.06.2017, 00:46 **[ТС]**

data.txt

@afront · 10.06.2017, 09:18

в матрице а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
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
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
#include <iostream>
#include <fstream>
#include <math.h>
#include <iomanip>
using namespace std;
int main()
{
    const int nku=13, nke=20;
    double ak[4][4]={},a[2*nku][2*nku]={},a1[2*nku][2*nku]={},b[2*nku]={}, x[2*nku]={},re[2*nku]={};
    double xy[2][nku]={},f[2][nku]={}, a2[2*nku][2*nku+1]={}, b1[2*nku]={}, st[nke]={};
    int issg[4],nc[2][nke],nuzst[nku],ipzst[2][nku],nzuz[nku];
    int ku,kuzst,knu;
 
    ifstream read ("data.txt");
    read  >> ku >> kuzst >> knu;
    for (int i=0; i<ku; i++) 
        read >> xy[0][i] >> xy[1][i];
 
    for (int i=0; i<kuzst; i++) 
        read >> nuzst[i] >> ipzst[0][i] >> ipzst[1][i];
 
    for (int i=0; i<knu; i++) 
        read >> nzuz[i] >> f[0][i] >> f[1][i];
 
    int ke;
    read >> ke;
    for (int i=0; i<ke; i++) 
        read >> nc[0][i] >> nc[1][i];
 
    double e,area,uvel;
    read >> e >> area >> uvel; 
    read.close();
    cout << fixed << setprecision(4);
    ofstream write("Picture");
    //  cout << fixed << setprecision(4);
    for (int i=0; i<ku; i++){ 
        write << setw(5)<< i+1 <<  setw(15) << xy[0][i] << setw(15) << xy[1][i] <<"    0.0" << endl;
    }
    write <<  "FACES" << endl; 
    for (int i=0; i<ke; i++) 
        write << setw(8)<<nc[0][i] << setw(8)<<nc[1][i] << "." << endl;;
 
    write.close();
 
    int kss=2*ku;
    for (int i=0; i<kss; i++){
        b[i]=0.0;
        re[i]=0.0;
        for (int j=0; j<kss; j++)
            a[i][j]=0.0;
    }
 
    for (int ij=0; ij<ke; ij++){
        int n1 = nc[0][ij]-1;
        int n2 = nc[1][ij]-1;
        double al1 = xy[0][n2] - xy[0][n1];
        double al2 = xy[1][n2] - xy[1][n1];
        double alenght = sqrt(al1*al1 + al2*al2);
        double c = (xy[0][n2] - xy[0][n1])/alenght;
        double s = (xy[1][n2] - xy[1][n1])/alenght;
 
        ak[0][0] = c*c;
        ak[0][1] = s*c;
        ak[0][2] = -c*c;
        ak[0][3] = -s*c;
        ak[1][0] = s*c;
        ak[1][1] = s*s;
        ak[1][2] = -s*c;
        ak[1][3] = -s*s;
        ak[2][0] = -c*c;
        ak[2][1] = -s*c;
        ak[2][2] = c*c;
        ak[2][3] = s*c;
        ak[3][0] = -s*c;
        ak[3][1] = -s*s;
        ak[3][2] = s*c;
        ak[3][3] = s*s;
 
        for (int i=0; i<4; i++)
            for (int j=0; j<4; j++)
                ak[i][j] = ak[i][j]*e*area/alenght;     
 
        issg[0] = nc[0][ij]*2 - 1;
        issg[1] = nc[0][ij]*2;
        issg[2] = nc[1][ij]*2 - 1;
        issg[3] = nc[1][ij]*2;
 
        for (int i=0; i<4; i++)
            for (int j=0; j<4; j++){
                int ig = issg[i]-1;
                int jg = issg[j]-1;
                a[ig][jg] = a[ig][jg] + ak[i][j];
            }
    }
            for (int ij=0; ij<knu; ij++){
                int i1 = nzuz[ij]*2 - 1;
                int i2 = nzuz[ij]*2;
                b[i1] = f[0][ij];
                b[i2] = f[1][ij];
            }
 
            for (int i=0; i<kss; i++){
                for (int j=0; j<kss; j++)
                    a1[i][j] = a[i][j];
                b1[i] = b[i];
            }
            for (int ij=0; ij<kuzst; ij++)
            {
                int i1 = nuzst[ij]*2 - 1-1;
                int i2 = nuzst[ij]*2-1;
                if (ipzst[0][ij] == 1)
                {
                    for (int i=0; i<kss; i++) 
                        a[i1][i] = 0.0;
 
                    b[i1] = 0.0;
                    a[i1][i1] = 1;
                }
                if (ipzst[1][ij] == 1) 
                {
                    for (int i=0; i<kss; i++) 
                        a[i2][i] = 0.0;
 
                    for (int i=0; i<kss; i++) 
                        a[i][i2] = 0.0;
 
                    b[i2] = 0.0;
                    a[i2][i2] = 1;
                }
            }
            int n=kss;
            for (int i=0; i<n; i++)
            {
                a2[i][n] = b[i];
                for (int j=0; j<n; j++) 
                    a2[i][j] = a[i][j];
 
            }
 
            int n1=n+1;
            for (int k=0; k<n; k++)
            {
                int k1=k;
                double s=a2[k][k];
                int j=k;
                for (int i=k1; i<n; i++)
                {
                    double r = a2[i][k];
                    if (abs(r) > abs(s)) 
                    {
                        s=r;
                        j=i;
                    }
                }
                if (s > 0.00001) 
                {
                    if (j != k) 
                    {
                        for(int i=k; i<n1; i++)
                        {
                            double r = a2[k][i];
                            a2[k][i] = a2[j][i];
                            a2[j][i] = r;
                        }
                    }
                    for(j=k1; j<n1; j++)
                        a2[k][j] = a2[k][j]/s;
 
                    for(int i=k1; i<n; i++)
                    {
                        double r = a2[i][k];
                        for(int j=k1; j<n1; j++)
                            a2[i][j] = a2[i][j] - a2[k][j]*r;
                    }
                }
            }
             ofstream write1("a.txt");
             write1 << setprecision(4) << fixed;
             for (int i=0; i<n; i++){ 
                 for (int j=0; j<n1; j++) 
                        write1 << setw(15) << a[i][j];
                write1 << endl;
             }
            write1.close();
 
                x[n-1] = a2[n-1][n1-1];
                for( int i = n - 1;i>= 0; i--)
                {
                    double s = a2[i][n1-1];
                    for( int j = i; j<n;j++)
                        s = s - a2[i][j]*x[j];
 
                    x[i] = s;
                }
                {
                    ofstream write1("solve.txt");
                    write1 << setprecision(4) << fixed;
                    for (int i=0; i<kss; i++) 
                        write1 << setw(15) << x[i] << " "<< endl;
 
                    write1.close();
                }
            
 
            ofstream write2("Dpicture");
            write2 << setprecision(4) << fixed;
            for (int i=0; i<ku; i++){
                write2 << setw(5) << i <<setw(15) <<  xy[0][i] + uvel * x[2 * i - 1] << setw(15) << xy[1][i] + uvel * x[2 * i]<<"   0.0"<<endl;
            }
            for (int i=0; i<ke; i++){
                write2 << setw(8) <<  nc[0][i] <<setw(8) <<  nc[2][i]<<"."<<endl;
            }
            write2.close();
 
            for (int i=0; i<kss; i++){
                double s=0;
                for (int j=1; j<kss; j++)   
                    s = s + a1[i][j]*x[j];
 
                re[i] = s - b1[i];
            }
 
            ofstream write3("reactions.txt");
            write3 << setprecision(6) << fixed;
            for (int i=0; i < kss; i++) 
                write3 << setw(15)<< re[i]<<endl;
            write3.close();
 
            for (int ij=0; ij<ke; ij++){
                n1 = nc[0][ij];
                int n2 = nc[1][ij];
                double alenght = sqrt(xy[0][n2] - xy[0][n1])*(xy[0][n2] - xy[0][n1]) + (xy[1][n2] - xy[1][n1])*(xy[1][n2] - xy[1][n1]);
                double c = (xy[0][n2] - xy[0][n1])/alenght;
                double s = (xy[1][n2] - xy[1][n1])/alenght;
                int ix = nc[0][ij]*2 - 1;
                int iy = nc[0][ij]*2;
                int jx = nc[1][ij]*2 - 1;
                int jy = nc[1][ij]*2;
                double elenght = (x[jx] - x[ix])*c + (x[jy] - x[iy])*s;
 
                st[ij] = elenght/alenght*e;
            }
 
            ofstream write4("stress.txt");
            write4 << setprecision(6) << fixed;
            for (int i=0; i<ke; i++) 
                write4 << setw(20)<< st[i]<<"."<<endl;
 
            write4.close();
    }

Добавлено через 16 часов 22 минуты
еще с некоторыми исправлениями

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
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
#include <iostream>
#include <fstream>
#include <math.h>
#include <iomanip>
using namespace std;
int main()
{
    const int nku=13, nke=20;
    double ak[4][4]={},a[2*nku][2*nku]={},a1[2*nku][2*nku]={},b[2*nku]={}, x[2*nku]={},re[2*nku]={};
    double xy[2][nku]={},f[2][nku]={}, a2[2*nku][2*nku+1]={}, b1[2*nku]={}, st[nke]={};
    int issg[4],nc[2][nke],nuzst[nku],ipzst[2][nku],nzuz[nku];
    int ku,kuzst,knu;
 
    ifstream read ("data.txt");
    read  >> ku >> kuzst >> knu;
    for (int i=0; i<ku; i++) 
        read >> xy[0][i] >> xy[1][i];
 
    for (int i=0; i<kuzst; i++) 
        read >> nuzst[i] >> ipzst[0][i] >> ipzst[1][i];
 
    for (int i=0; i<knu; i++) 
        read >> nzuz[i] >> f[0][i] >> f[1][i];
 
    int ke;
    read >> ke;
    for (int i=0; i<ke; i++) 
        read >> nc[0][i] >> nc[1][i];
 
    double e,area,uvel;
    read >> e >> area >> uvel; 
    read.close();
     
    ofstream write("Picture");
    write << fixed << setprecision(4);
    for (int i=0; i<ku; i++){ 
        write << setw(5)<< i+1 <<  setw(15) << xy[0][i] << setw(15) << xy[1][i] <<"    0.0" << endl;
    }
    write <<  "FACES" << endl; 
    for (int i=0; i<ke; i++) 
        write << setw(8)<<nc[0][i] << setw(8)<<nc[1][i] << "." << endl;;
 
    write.close();
 
    int kss=2*ku;
    for (int i=0; i<kss; i++)
    {
        b[i]=0.0;
        re[i]=0.0;
        for (int j=0; j<kss; j++)
            a[i][j]=0.0;
    }
 
    for (int ij=0; ij<ke; ij++)
    {
        int n1 = nc[0][ij]-1;
        int n2 = nc[1][ij]-1;
        double al1 = xy[0][n2] - xy[0][n1];
        double al2 = xy[1][n2] - xy[1][n1];
        double alenght = sqrt(al1*al1 + al2*al2);
        double c = al1/alenght;
        double s = al2/alenght;
 
        ak[0][0] = c*c;
        ak[0][1] = s*c;
        ak[0][2] = -c*c;
        ak[0][3] = -s*c;
        ak[1][0] = s*c;
        ak[1][1] = s*s;
        ak[1][2] = -s*c;
        ak[1][3] = -s*s;
        ak[2][0] = -c*c;
        ak[2][1] = -s*c;
        ak[2][2] = c*c;
        ak[2][3] = s*c;
        ak[3][0] = -s*c;
        ak[3][1] = -s*s;
        ak[3][2] = s*c;
        ak[3][3] = s*s;
 
        for (int i=0; i<4; i++)
            for (int j=0; j<4; j++)
                ak[i][j] = ak[i][j]*e*area/alenght;     
 
        issg[0] = nc[0][ij]*2 - 1;
        issg[1] = nc[0][ij]*2;
        issg[2] = nc[1][ij]*2 - 1;
        issg[3] = nc[1][ij]*2;
 
        for (int i=0; i<4; i++)
            for (int j=0; j<4; j++)
            {
                int ig = issg[i] - 1;
                int jg = issg[j] - 1;
                a[ig][jg] = a[ig][jg] + ak[i][j];
            }
    }
 
    for (int ij=0; ij<knu; ij++)
    {
        int i1 = nzuz[ij]*2 - 1 - 1;
        int i2 = nzuz[ij]*2 - 1 ;
        b[i1] = f[0][ij];
        b[i2] = f[1][ij];
    }
 
    for (int i=0; i<kss; i++)
    {
        for (int j=0; j<kss; j++)
            a1[i][j] = a[i][j];
        b1[i] = b[i];
    }
    for (int ij=0; ij<kuzst; ij++)
    {
        int i1 = nuzst[ij]*2 - 1 - 1;
        int i2 = nuzst[ij]*2 - 1;
        if (ipzst[0][ij] == 1)
        {
            for (int i=0; i<kss; i++) 
                a[i1][i] = 0.0;
 
            b[i1] = 0.0;
            a[i1][i1] = 1;
        }
        if (ipzst[1][ij] == 1) 
        {
            for (int i=0; i<kss; i++) 
                a[i2][i] = 0.0;
 
            for (int i=0; i<kss; i++) 
                a[i][i2] = 0.0;
 
            b[i2] = 0.0;
            a[i2][i2] = 1;
        }
    }
    int n=kss;
    for (int i=0; i<n; i++)
    {
        a2[i][n] = b[i];
        for (int j=0; j<n; j++) 
            a2[i][j] = a[i][j];
 
    }
 
    int n1=n+1;
    for (int k=0; k<n; k++)
    {
        int k1=k;
        double s=a2[k][k];
        int j=k;
        for (int i=k1; i<n; i++)
        {
            double r = a2[i][k];
            if (abs(r) > abs(s)) 
            {
                s=r;
                j=i;
            }
        }
        if (s > 0.00001) 
        {
            if (j != k) 
            {
                for(int i=k; i<n1; i++)
                {
                    double r = a2[k][i];
                    a2[k][i] = a2[j][i];
                    a2[j][i] = r;
                }
            }
            for(j=k1; j<n1; j++)
                a2[k][j] = a2[k][j]/s;
 
            for(int i=k1; i<n; i++)
            {
                double r = a2[i][k];
                for(int j=k1; j<n1; j++)
                    a2[i][j] = a2[i][j] - a2[k][j]*r;
            }
        }
    }
 
    ofstream write10("a.txt");
    write10 << setprecision(4) << fixed;
    for (int i=0; i<n; i++)
    { 
        for (int j=0; j<n1; j++) 
            write10 << setw(15) << a2[i][j];
        write10 << endl;
    }
    write10.close();
 
    x[n-1] = a2[n-1][n1-1];
    for( int i = n - 1;i>= 0; i--)
    {
        double s = a2[i][n1-1];
        for( int j = i; j<n;j++)
            s = s - a2[i][j]*x[j];
        x[i] = s;
    }
 
    ofstream write1("solve.txt");
    write1 << setprecision(4) << fixed;
    for (int i=0; i<kss; i++) 
        write1 << setw(15) << x[i] << " "<< endl;
    write1.close();
 
 
 
    ofstream write2("Dpicture");
    write2 << setprecision(4) << fixed;
    for (int i=0; i<ku; i++)
        write2 << setw(5) << i <<setw(15) <<  xy[0][i] + uvel * x[2 * i] << setw(15) << xy[1][i] + uvel * x[2 * i+1]<<"   0.0"<<endl;
 
    for (int i=0; i<ke; i++)
        write2 << setw(8) <<  nc[0][i] <<setw(8) <<  nc[2][i]<<"."<<endl;
 
    write2.close();
 
    for (int i=0; i<kss; i++)
    {
        double s=0;
        for (int j=1; j<kss; j++)   
            s = s + a1[i][j]*x[j];
        re[i] = s - b1[i];
    }
 
    ofstream write3("reactions.txt");
    write3 << setprecision(6) << fixed;
    for (int i=0; i < kss; i++) 
        write3 << setw(15)<< re[i]<<endl;
    write3.close();
 
    for (int ij=0; ij<ke; ij++)
    {
        n1 = nc[0][ij];
        int n2 = nc[1][ij];
        double al1=xy[0][n2] - xy[0][n1];
        double al2=xy[1][n2] - xy[1][n1];
        double alenght = sqrt(al1*al1 + al2*al2);
        double c = al1/alenght;
        double s = al2/alenght;
        int ix = nc[0][ij]*2 - 1;
        int iy = nc[0][ij]*2;
        int jx = nc[1][ij]*2 - 1;
        int jy = nc[1][ij]*2;
        double elenght = (x[jx] - x[ix])*c + (x[jy] - x[iy])*s;
 
        st[ij] = elenght/alenght*e;
    }
 
    ofstream write4("stress.txt");
    write4 << setprecision(6) << fixed;
    for (int i=0; i<ke; i++) 
        write4 << setw(20)<< st[i]<<"."<<endl;
 
    write4.close();
}

@мухобой · 11.06.2017, 21:51 **[ТС]**

Данные получение на с++ не сходятся с фортраном...

Новые блоги и статьи Все статьи Все блоги /
Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Нашел на реддите интересную статью под названием «Кто-нибудь знает, где получить бесплатный компьютер или. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro
Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	01.06.2017, 20:50 [ТС]
	Помогите пожалуйста Добавлено через 15 минут Сложно я так понимаю)) 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	05.06.2017, 19:21 [ТС]
	( 0

MrGluck Форумчанин 8216 / 5047 / 1437 Регистрация: 29.11.2010 Сообщений: 13,453
	05.06.2017, 19:24
	мухобой, здесь мало мастеров фортрана, вы бы лучше задание выложили 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	05.06.2017, 19:26 [ТС]
	Нужно написать программу которая считала мке ферму. Вот пример кода на Фортране, ее по сути нужно в с++ перевести 0

@afront 1498 / 1213 / 821 Регистрация: 29.02.2016 Сообщений: 3,630
	05.06.2017, 20:06
	мухобой, давно бы сам все переписал, элементарная программа. Нужно заменить чтение файла и фортранные циклы DO от 1 до N заменить на c++ типа for (int i = 0; i < N; i++). А что будет не понятно, конкретно спрашивайте 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	06.06.2017, 13:31 [ТС]
	Смысл тогда всего этого если подсказать даже не могут? 0

@karaulov6 23 / 23 / 6 Регистрация: 23.03.2013 Сообщений: 245
	06.06.2017, 15:33
	Это ищите раздел где понимают этот набор букв и цифр )) может там https://www.cyberforum.ru/fortran/ 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	06.06.2017, 20:20 [ТС]
	Отсутствуют вроде как Добавлено через 7 минут Исправил везде где есть запятые 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	06.06.2017, 20:55 [ТС]
	И ещё у меня сомнения где есть операторы перехода go to и метки, если не затруднит посмотрите, пожалуйста) 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	09.06.2017, 00:46 [ТС]
	data.txt 0

@мухобой 0 / 0 / 0 Регистрация: 02.04.2015 Сообщений: 61
	11.06.2017, 21:51 [ТС]
	Данные получение на с++ не сходятся с фортраном... 0