質問編集履歴
2
コードの修正
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -17,6 +17,17 @@
|
|
|
17
17
|
|
|
18
18
|
```
|
|
19
19
|
|
|
20
|
+
using namespace std;
|
|
21
|
+
#include <iostream>
|
|
22
|
+
#include <math.h>
|
|
23
|
+
#include <stdio.h>
|
|
24
|
+
#include <stdlib.h>
|
|
25
|
+
#include <string.h>
|
|
26
|
+
#include <cmath>
|
|
27
|
+
#include <limits>
|
|
28
|
+
|
|
29
|
+
|
|
30
|
+
|
|
20
31
|
void sf_F1_ijk(double vecRij[3],double vecRik[3], double Rij, double Rik, double cost, double lambda, double zeta,double *Fijk, double *dFijk_dRj[3], double *dFijk_dRk[3]) ;
|
|
21
32
|
|
|
22
33
|
void sf_cut(double Rij,double Rc, double *fc,double *dfc) ;
|
|
@@ -32,21 +43,18 @@
|
|
|
32
43
|
|
|
33
44
|
|
|
34
45
|
double dGi[3][3]={};
|
|
35
|
-
|
|
46
|
+
double Fijk,dFijk_dRj[3],dFijk_dRk[3], Fijk_dRk[3],vecRij[3]={ } ,vecRik[3]={ } ,vecRjk[3]={ } ,Rij, Rik;
|
|
36
|
-
double vecRij[3],vecRik[3], Rij, Rik,cost,lambda, zeta,
|
|
37
|
-
Fijk,dFijk_dRj[3], dFijk_dRk[3];
|
|
38
47
|
|
|
39
48
|
double s;
|
|
40
49
|
|
|
41
|
-
|
|
50
|
+
double cost,lambda, zeta ;
|
|
42
51
|
|
|
43
52
|
|
|
44
53
|
|
|
45
|
-
sf_F1_ijk(vecRij[3],vecRik[3], Rij, Rik, cost, lambda, zeta,&Fijk,&dFijk_dRj[3],&dFijk_dRk[3]) ;
|
|
46
54
|
|
|
47
55
|
|
|
56
|
+
sf_F1_ijk(vecRij[3],vecRik[3], Rij, Rik, cost, lambda, zeta,&Fijk,&dFijk_dRj[3],&dFijk_dRk[3]) ;
|
|
48
57
|
|
|
49
|
-
|
|
50
58
|
s=Fijk ;
|
|
51
59
|
|
|
52
60
|
cout<<s<<"\n" ;
|
|
@@ -55,18 +63,6 @@
|
|
|
55
63
|
}
|
|
56
64
|
|
|
57
65
|
|
|
58
|
-
/*
|
|
59
|
-
void sf_F1_ijk(double cost, double lambda, double zeta,double *Fijk
|
|
60
|
-
|
|
61
|
-
{
|
|
62
|
-
|
|
63
|
-
double arg ;
|
|
64
|
-
|
|
65
|
-
arg = 0.5000000*(1.0000000 + lambda*cost) ;
|
|
66
|
-
*Fijk =pow(arg,zeta) ;
|
|
67
|
-
}
|
|
68
|
-
*/
|
|
69
|
-
|
|
70
66
|
|
|
71
67
|
|
|
72
68
|
|
|
@@ -110,10 +106,10 @@
|
|
|
110
106
|
}
|
|
111
107
|
|
|
112
108
|
|
|
113
|
-
C++
|
|
114
109
|
|
|
115
110
|
|
|
116
111
|
|
|
112
|
+
|
|
117
113
|
void sf_G4_update(double nG4, double vecRij[3], double vecRik[3], double vecRjk[3], double Rij, double Rik, double Rjk,
|
|
118
114
|
double cost, double n[100], double G[100], int iG, double dGi[3][100], double dGj[3][100], double dGk[3][100])
|
|
119
115
|
|
|
@@ -129,9 +125,7 @@
|
|
|
129
125
|
//sf_F1_ijk
|
|
130
126
|
double F1, F1k,dF1k[3],dF1j[3];
|
|
131
127
|
//sf_F2
|
|
132
|
-
double F|列1|列2|列3|
|
|
133
|
-
|:--|:--:|--:|
|
|
134
|
-
|
|
128
|
+
double F2ij,dF2ij,F2ik ,dF2ik,F2jk,dF2jk ;
|
|
135
129
|
|
|
136
130
|
|
|
137
131
|
for( int iG4 = 0 ; iG4<nG4; iG4++){
|
|
@@ -215,4 +209,6 @@
|
|
|
215
209
|
|
|
216
210
|
}
|
|
217
211
|
|
|
212
|
+
|
|
213
|
+
|
|
218
214
|
```
|
1
コードとしてページに表示されるようにした
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -1,5 +1,5 @@
|
|
|
1
1
|
言語はC++，コンパイラはgcc,osはwslを使用しています．
|
|
2
|
-
下記のエラーが出てきて
|
|
2
|
+
下記のエラーが出てきてコンパイルが通らず困っています．どなたか原因をご教授お願い致します。
|
|
3
3
|
|
|
4
4
|
|
|
5
5
|
|
|
@@ -14,15 +14,9 @@
|
|
|
14
14
|
|
|
15
15
|
|
|
16
16
|
|
|
17
|
-
|
|
17
|
+
|
|
18
|
-
コード
|
|
19
18
|
```
|
|
20
19
|
|
|
21
|
-
|
|
22
|
-
|
|
23
|
-
|
|
24
|
-
|
|
25
|
-
//void sf_F1_ijk(double cost, double lambda, double zeta,double *Fijk) ;
|
|
26
20
|
void sf_F1_ijk(double vecRij[3],double vecRik[3], double Rij, double Rik, double cost, double lambda, double zeta,double *Fijk, double *dFijk_dRj[3], double *dFijk_dRk[3]) ;
|
|
27
21
|
|
|
28
22
|
void sf_cut(double Rij,double Rc, double *fc,double *dfc) ;
|
|
@@ -219,4 +213,6 @@
|
|
|
219
213
|
*zeta=coef_zeta[1];
|
|
220
214
|
*eta=coef_eta[1];
|
|
221
215
|
|
|
222
|
-
}
|
|
216
|
+
}
|
|
217
|
+
|
|
218
|
+
```
|