File tree 26 files changed +133
-133
lines changed
26 files changed +133
-133
lines changed Original file line number Diff line number Diff line change 1
- alpha = 2
2
- beta = 3
1
+ alpha = 1.5
2
+ beta = 1.2
3
3
4
- ∀(i) | 0<=i<=100 {
5
- ∀(j) | 0<=j<=150 {
6
- temp[i][j] = Σ(alpha * a[i][k] * b[k][j]) | 0<=k<=200
4
+ ∀(i) | 0<=i<1600 {
5
+ ∀(j) | 0<=j<=1800 {
6
+ temp[i][j] = Σ(alpha * a[i][k] * b[k][j]) | 0<=k<=2200
7
7
}
8
8
}
9
9
10
- ∀(i) | 0<=i<=100 {
11
- ∀(j) | 0<=j<=250 {
10
+ ∀(i) | 0<=i<=1600 {
11
+ ∀(j) | 0<=j<=2400 {
12
12
d[i][j] = d[i][j] * beta
13
- d[i][j] = Σ(temp[i][k] * c[k][j]) | 0<=k<=150
13
+ d[i][j] = Σ(temp[i][k] * c[k][j]) | 0<=k<=1800
14
14
}
15
15
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<=100 {
2
- ∀(j) | 0<=j<=150 {
3
- ab[i][j] = Σ(a[i][k] * b[k][j]) | 0<=k<=200
1
+ ∀(i) | 0<=i<1600 {
2
+ ∀(j) | 0<=j<1800 {
3
+ ab[i][j] = Σ(a[i][k] * b[k][j]) | 0<=k<2000
4
4
}
5
5
}
6
6
7
- ∀(i) | 0<=i<=150 {
8
- ∀(j) | 0<=j<=250 {
9
- cd[i][j] = Σ(c[i][k] * d[k][j]) | 0<=k<=300
7
+ ∀(i) | 0<=i<1800 {
8
+ ∀(j) | 0<=j<2200 {
9
+ cd[i][j] = Σ(c[i][k] * d[k][j]) | 0<=k<2400
10
10
}
11
11
}
12
12
13
- ∀(i) | 0<=i<=100 {
14
- ∀(j) | 0<=j<=250 {
15
- abcd[i][j] = Σ(ab[i][k] * cd[k][j]) | 0<=k<=150
13
+ ∀(i) | 0<=i<1600 {
14
+ ∀(j) | 0<=j<2200 {
15
+ abcd[i][j] = Σ(ab[i][k] * cd[k][j]) | 0<=k<1800
16
16
}
17
17
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<20000 {
2
- tmp[i] = Σ(A[i][j]*x[j]) | 0<=j<24000
3
- y[k] = Σ(A[i][k]*tmp[i]) | 0<=k<24000
1
+ ∀(i) | 0<=i<1800 {
2
+ tmp[i] = Σ(A[i][j]*x[j]) | 0<=j<2200
3
+ y[k] = Σ(A[i][k]*tmp[i]) | 0<=k<2200
4
4
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<=200 {
2
- s[j] = Σ(A[i][j] * r[i]) | 0<=j<=150
3
- q[i] = Σ(A[i][j] * p[j]) | 0<=j<=150
1
+ ∀(i) | 0<=i<2200 {
2
+ s[j] = Σ(A[i][j] * r[i]) | 0<=j<1800
3
+ q[i] = Σ(A[i][j] * p[j]) | 0<=j<1800
4
4
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<200 {
1
+ ∀(i) | 0<=i<4000 {
2
2
∀(j) | 0<=j<i {
3
3
∀(k) | 0<=k<j {
4
4
A[i][j] = A[i][j] - A[i][k] * A[j][k]
Original file line number Diff line number Diff line change 1
- ∀(x) | 0<=x<=100 {
2
- mean[x] = Σ(data[k][x]/100 ) | 0<=k<=100
1
+ ∀(x) | 0<=x<2600 {
2
+ mean[x] = Σ(data[k][x]/3000 ) | 0<=k<3000
3
3
}
4
4
5
- ∀(i) | 0<=i<=100 {
6
- ∀(j) | 0<=j<=100 {
7
- cov[i][j] = Σ(((data[k][i] - mean[i])*(data[k][j] - mean[j]))/99) | 0<=k<=100
5
+ ∀(i) | 0<=i<2600 {
6
+ ∀(j) | 0<=j<2600 {
7
+ cov[i][j] = Σ(((data[k][i] - mean[i])*(data[k][j] - mean[j]))/99) | 0<=k<3000
8
8
}
9
9
}
10
10
11
- ∀(t) | 0<=t<=100 {
12
- std_dev[t] = Σ(((data[k][t] - mean[t])*(data[k][t] - mean[t]))/99) | 0<=k<=100
11
+ ∀(t) | 0<=t<2600 {
12
+ std_dev[t] = Σ(((data[k][t] - mean[t])*(data[k][t] - mean[t]))/99) | 0<=k<3000
13
13
}
14
14
15
- ∀(p) | 0<=p<=100 {
15
+ ∀(p) | 0<=p<2600 {
16
16
std_dev[p] = √(std_dev[p])
17
17
}
18
18
19
- ∀(i) | 0<=i<=100 {
20
- ∀(j) | 0<=j<=100 {
19
+ ∀(i) | 0<=i<2600 {
20
+ ∀(j) | 0<=j<2600 {
21
21
corr[i][j] = cov[i][j]/(std_dev[i]*std_dev[j])
22
22
}
23
23
}
Original file line number Diff line number Diff line change 1
- ∀(x) | 0<=x<=100 {
2
- mean[x] = Σ(data[k][x]/100 ) | 0<=k<=100
1
+ ∀(x) | 0<=x<2600 {
2
+ mean[x] = Σ(data[k][x]/3000 ) | 0<=k<3000
3
3
}
4
4
5
- ∀(i) | 0<=i<=100 {
6
- ∀(j) | i<=j<=100 {
7
- cov[i][j] = Σ(((data[k][i] - mean[i])*(data[k][j] - mean[j]))/99 ) | 0<=k<=100
5
+ ∀(i) | 0<=i<2600 {
6
+ ∀(j) | i<=j<2600 {
7
+ cov[i][j] = Σ(((data[k][i] - mean[i])*(data[k][j] - mean[j]))/2999 ) | 0<=k<3000
8
8
cov[j][i] = cov[i][j]
9
9
}
10
10
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<=150 {
2
- ∀(j) | 0<=j<=250 {
3
- ∀(k) | 0<=k<=350 {
4
- sum[k] = Σ(A[i][j][h] * C4[h][k]) | 0<=h<=350
1
+ ∀(i) | 0<=i<250 {
2
+ ∀(j) | 0<=j<220 {
3
+ ∀(k) | 0<=k<270 {
4
+ sum[k] = Σ(A[i][j][h] * C4[h][k]) | 0<=h<270
5
5
}
6
- ∀(k) | 0<=k<=350 {
6
+ ∀(k) | 0<=k<270 {
7
7
A[i][j][k] = sum[k]
8
8
}
9
9
}
Original file line number Diff line number Diff line change 1
- ∀(t) | 0<=t<=250 {
2
- ∀(j) | 1<=j<100 {
1
+ ∀(t) | 0<=t<40 {
2
+ ∀(j) | 1<=j<80 {
3
3
ey[0][j] = _fict_[t]
4
4
}
5
- ∀(i) | 1<=i<150 {
6
- ∀(j) | 0<=j<100 {
5
+ ∀(i) | 1<=i<60 {
6
+ ∀(j) | 0<=j<80 {
7
7
ey[i][j] = ey[i][j] - 0.5 * (hz[i][j] - hz[i-1][j])
8
8
}
9
9
}
10
- ∀(i) | 0<=i<150 {
11
- ∀(j) | 1<=j<100 {
10
+ ∀(i) | 0<=i<60 {
11
+ ∀(j) | 1<=j<80 {
12
12
ex[i][j] = ex[i][j] - 0.5 * (hz[i][j] - hz[i][j-1])
13
13
}
14
14
}
15
- ∀(i) | 0<=i<149 {
16
- ∀(j) | 0<=j<99 {
15
+ ∀(i) | 0<=i<59 {
16
+ ∀(j) | 0<=j<79 {
17
17
hz[i][j] = hz[i][j] - 0.7 * (ex[i][j+1] - ex[i][j] + ey[i+1][j] - ey[i][j])
18
18
}
19
19
}
Original file line number Diff line number Diff line change 1
- alpha = 10
2
- beta = 2
1
+ alpha = 1.5
2
+ beta = 1. 2
3
3
4
- ∀(i) | 0<=i<=50 {
4
+ ∀(i) | 0<=i<2000 {
5
5
6
- ∀(j) | 0<=j<=100 {
6
+ ∀(j) | 0<=j<2300 {
7
7
C[i][j] = C[i][j] * beta
8
8
}
9
9
10
- ∀(k) | 0<=k<=250 {
11
- ∀(j) | 0<=j<=100 {
12
- C[i][j] = C[i][j] + alpha * A[i][k] * B[k][j]
13
- }
10
+ ∀(k) | 0<=k<2600 {
11
+ C[i][j] = Σ(alpha * A[i][k] * B[k][j]) | 0<=j<2300
14
12
}
15
13
}
Original file line number Diff line number Diff line change 1
- alpha = 2
2
- beta = 3
1
+ alpha = 1.5
2
+ beta = 1.2
3
3
4
- ∀(i) | 0<=i<=150 {
5
- A[i][j] = Σ(u1[i] * v1[j] + u2[i] * v2[j]) | 0<=j<=150
4
+ ∀(i) | 0<=i<4000 {
5
+ A[i][j] = Σ(u1[i] * v1[j] + u2[i] * v2[j]) | 0<=j<4000
6
6
}
7
7
8
- ∀(i) | 0<=i<=150 {
9
- x[i] = Σ(beta * A[j][i] * y[j]) | 0<=j<=150
8
+ ∀(i) | 0<=i<4000 {
9
+ x[i] = Σ(beta * A[j][i] * y[j]) | 0<=j<4000
10
10
}
11
11
12
- ∀(i) | 0<=i<=150 {
12
+ ∀(i) | 0<=i<4000 {
13
13
x[i] = x[i] + z[i]
14
14
}
15
15
16
- ∀(i) | 0<=i<=150 {
17
- w[i][j] = Σ(alpha * A[i][j] * x[i]) | 0<=j<=150
16
+ ∀(i) | 0<=i<4000 {
17
+ w[i][j] = Σ(alpha * A[i][j] * x[i]) | 0<=j<4000
18
18
}
Original file line number Diff line number Diff line change 1
- alpha = 2
2
- beta = 3
1
+ alpha = 1.5
2
+ beta = 1.2
3
3
4
- ∀(i) | 0<=i<=150 {
5
- ∀(j) | 0<=j<=150 {
4
+ ∀(i) | 0<=i<2800 {
5
+ ∀(j) | 0<=j<2800 {
6
6
tmp[i] = x[j] * A[i][j] + tmp[i]
7
7
y[i] = x[j] * B[i][j] + y[i]
8
8
}
Original file line number Diff line number Diff line change 1
- ∀(k) | 0<=k<=150 {
2
- norm = Σ(A[i][k] * A[i][k]) | 0<=i<=250
1
+ ∀(k) | 0<=k<2600 {
2
+ norm = Σ(A[i][k] * A[i][k]) | 0<=i<2000
3
3
R[k][k] = √(norm)
4
- ∀(i) | 0<=k<=150 {
5
- Q[i][k] = A[j ][k] / R[k][k]
4
+ ∀(i) | 0<=i<2000 {
5
+ Q[i][k] = A[i ][k] / R[k][k]
6
6
}
7
- ∀(j) | k+1<=j<=150 {
8
- R[k][j] = Σ(Q[i][k] * A[i][j]) | 0<=i<=250
9
- A[i][j] = A[i][j] - Q[i][k] * R[k][j]
7
+ ∀(j) | k+1<=j<2600 {
8
+ R[k][j] = Σ(Q[i][k] * A[i][j]) | 0<=i<2000
9
+ ∀(j) | 0<=i<2000 {
10
+ A[i][j] = A[i][j] - Q[i][k] * R[k][j]
11
+ }
10
12
}
11
13
}
Original file line number Diff line number Diff line change 1
- ∀(t) | 0<=t<=250 {
2
- ∀(i) | 1<=i<100 {
3
- ∀(j) | 1<=j<100 {
4
- ∀(k) | 1<=k<100 {
1
+ ∀(t) | 0<=t<=1000 {
2
+ ∀(i) | 1<=i<199 {
3
+ ∀(j) | 1<=j<199 {
4
+ ∀(k) | 1<=k<199 {
5
5
B[i][j][k] = A[i][j][k] + 0.125 * (A[i+1][j][k] - 2.0 * A[i][j][k] + A[i-1][j][k]) + 0.125 * (A[i][j+1][k] - 2.0 * A[i][j][k] + A[i][j-1][k]) + 0.125 * (A[i][j][k+1] - 2.0 * A[i][j][k] + A[i][j][k-1])
6
6
}
7
7
}
8
8
}
9
- ∀(i) | 1<=i<100 {
10
- ∀(j) | 1<=j<100 {
11
- ∀(k) | 1<=k<100 {
9
+ ∀(i) | 1<=i<199 {
10
+ ∀(j) | 1<=j<199 {
11
+ ∀(k) | 1<=k<199 {
12
12
A[i][j][k] = B[i][j][k] + 0.125 * (B[i+1][j][k] - 2.0 * B[i][j][k] + B[i-1][j][k]) + 0.125 * (B[i][j+1][k] - 2.0 * B[i][j][k] + B[i][j-1][k]) + 0.125 * (B[i][j][k+1] - 2.0 * B[i][j][k] + B[i][j][k-1])
13
13
}
14
14
}
Original file line number Diff line number Diff line change 1
- ∀(t) | 0<=t<=250 {
2
- ∀(i) | 1<=i<=99 {
1
+ ∀(t) | 0<=t<1000 {
2
+ ∀(i) | 1<=i<3999 {
3
3
B[i] = 0.33333 * (A[i] + A[i-1] + A[i+1])
4
4
}
5
- ∀(i) | 1<=i<=99 {
5
+ ∀(i) | 1<=i<3999 {
6
6
A[i] = 0.33333 * (B[i] + B[i-1] + B[i+1])
7
7
}
8
8
}
Original file line number Diff line number Diff line change 1
- ∀(t) | 0<=t<=250 {
2
- ∀(i) | 1<=i<100 {
3
- ∀(j) | 1<=j<100 {
1
+ ∀(t) | 0<=t<1000 {
2
+ ∀(i) | 1<=i<2799 {
3
+ ∀(j) | 1<=j<2799 {
4
4
B[i][j] = 0.2 * (A[i][j] + A[i][j-1] + A[i][1+j] + A[1+i][j] + A[i-1][j])
5
5
}
6
6
}
7
- ∀(p) | 1<=p<100 {
8
- ∀(q) | 1<=q<100 {
7
+ ∀(p) | 1<=p<2799 {
8
+ ∀(q) | 1<=q<2799 {
9
9
A[p][q] = 0.2 * (B[p][q] + B[p][q-1] + B[p][1+q] + B[1+p][q] + B[p-1][q])
10
10
}
11
11
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<200 {
1
+ ∀(i) | 0<=i<4000 {
2
2
∀(j) | 0<=j<i {
3
3
∀(k) | 0<=k<j {
4
4
A[i][j] = A[i][j] - A[i][k] * A[k][j]
5
5
}
6
6
A[i][j] = A[i][j]/A[j][j]
7
7
}
8
- ∀(j) | i<=j<200 {
8
+ ∀(j) | i<=j<4000 {
9
9
∀(k) | 0<=k<i {
10
10
A[i][j] = A[i][j] - A[i][k] * A[k][j]
11
11
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<=250 {
1
+ ∀(i) | 0<=i<4000 {
2
2
∀(j) | 0<=j<i {
3
3
w = A[i][j]
4
4
∀(k) | 0<=k<j {
5
5
w = w - (A[i][k] * A[k][j])
6
6
}
7
7
A[i][j] = w / A[j][j]
8
8
}
9
- ∀(j) | i<=j<=250 {
9
+ ∀(j) | i<=j<4000 {
10
10
w = A[i][j]
11
11
∀(k) | 0<=k<i {
12
12
w = w - (A[i][k] * A[k][j])
15
15
}
16
16
}
17
17
18
- ∀(i) | 0<=i<=250 {
18
+ ∀(i) | 0<=i<4000 {
19
19
w = b[i]
20
20
∀(j) | 0<=j<i {
21
21
w = w - A[i][j] * y[j]
22
22
}
23
23
y[i] = w
24
24
}
25
25
26
- ∀(i) | 0<=i<250 {
27
- w = y[250 -1-i]
28
- ∀(j) | 250 -i<=j<250 {
29
- w = w - A[249 -i][j] * x[j]
26
+ ∀(i) | 0<=i<4000 {
27
+ w = y[4000 -1-i]
28
+ ∀(j) | 4000 -i<=j<4000 {
29
+ w = w - A[3999 -i][j] * x[j]
30
30
}
31
- x[250 -1-i] = w / A[250 -1-i][250 -1-i]
31
+ x[4000 -1-i] = w / A[4000 -1-i][4000 -1-i]
32
32
}
Original file line number Diff line number Diff line change 1
- ∀(i) | 0<=i<=150 {
2
- x1[i] = Σ(A[i][j] * y_1[j]) | 0<=j<=150
1
+ ∀(i) | 0<=i<4000 {
2
+ x1[i] = Σ(A[i][j] * y_1[j]) | 0<=j<4000
3
3
}
4
4
5
- ∀(i) | 0<=i<=150 {
6
- x2[i] = Σ(A[j][i] * y_2[j]) | 0<=j<=150
5
+ ∀(i) | 0<=i<4000 {
6
+ x2[i] = Σ(A[j][i] * y_2[j]) | 0<=j<4000
7
7
}
Original file line number Diff line number Diff line change 1
- ∀(t) | 0<=t<=250 {
2
- ∀(i) | 1<=i<100 {
3
- ∀(j) | 1<=j<100 {
1
+ ∀(t) | 0<=t<1000 {
2
+ ∀(i) | 1<=i<4000-2 {
3
+ ∀(j) | 1<=j<4000-2 {
4
4
A_t[i][j] = (A[i][j] + A[i-1][j-1] + A[i-1][j] + A[i-1][j+1] + A[i][j-1] + A[i][j+1] + A[i+1][j-1] + A[i+1][j] + A[i+1][j+1])/9.0
5
5
}
6
6
}
7
- ∀(i) | 1<=i<=100 {
8
- ∀(j) | 1<=j<=100 {
7
+ ∀(i) | 1<=i<4000-2 {
8
+ ∀(j) | 1<=j<4000-2 {
9
9
A[i][j] = A_t[i][j]
10
10
}
11
11
}
Original file line number Diff line number Diff line change 1
- alpha = 10
2
- beta = 2
1
+ alpha = 1.5
2
+ beta = 1. 2
3
3
4
- ∀(i) | 0<=i<=150 {
5
- ∀(j) | 0<=j<=250 {
4
+ ∀(i) | 0<=i<2000 {
5
+ ∀(j) | 0<=j<2600 {
6
6
temp2 = 0
7
7
C[k][j] = Σ(alpha *B[i][j] * A[i][k]) | 0<=k<i
8
8
temp2 = Σ(B[k][j] * A[i][k]) | 0<=k<i
You can’t perform that action at this time.
0 commit comments