МКЕ (перевод 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 **[ТС]**

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

Новые блоги и статьи Все статьи Все блоги /
Модульный подход на примере 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 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .

@мухобой 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,631
	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