質問編集履歴
2
タイプミス
test
CHANGED
File without changes
|
test
CHANGED
@@ -96,9 +96,9 @@
|
|
96
96
|
|
97
97
|
|
98
98
|
|
99
|
-
def cg(
|
99
|
+
def cg(A, b, x0):
|
100
100
|
|
101
|
-
x = x
|
101
|
+
x = x0
|
102
102
|
|
103
103
|
|
104
104
|
|
@@ -150,7 +150,7 @@
|
|
150
150
|
|
151
151
|
|
152
152
|
|
153
|
-
|
153
|
+
A = mmread('bcsstm01.mtx').todense()
|
154
154
|
|
155
155
|
|
156
156
|
|
@@ -164,7 +164,7 @@
|
|
164
164
|
|
165
165
|
|
166
166
|
|
167
|
-
ans = cg(
|
167
|
+
ans = cg(A, b, x0)
|
168
168
|
|
169
169
|
```
|
170
170
|
|
1
質問をより具体的に改良する
test
CHANGED
File without changes
|
test
CHANGED
@@ -104,6 +104,20 @@
|
|
104
104
|
|
105
105
|
r0 = (b - (mydot(A,x)))
|
106
106
|
|
107
|
+
|
108
|
+
|
109
|
+
#r0 = [e1 - e2 for e1, e2 in zip(b, (mydot(A,x)))]
|
110
|
+
|
111
|
+
|
112
|
+
|
113
|
+
#r2 = set(b) - set(mydot(A,x))
|
114
|
+
|
115
|
+
#r0 = list(r2)
|
116
|
+
|
117
|
+
|
118
|
+
|
119
|
+
r0 =
|
120
|
+
|
107
121
|
p = r0
|
108
122
|
|
109
123
|
k = 0
|