@@ -138,10 +138,10 @@ void print_mat(const Matrix& A) {
138138}
139139
140140template <class Matrix>
141- void init_mat(Matrix& A, typename Matrix::value_type geta = 1 ) {
141+ void init_mat(Matrix& A) {
142142 for(int i = 0; i < A.extent(0); ++i) {
143143 for(int j = 0; j < A.extent(1); ++j) {
144- A[i, j] = i * A.extent(1) + j + geta ;
144+ A[i, j] = i * A.extent(1) + j;
145145 }
146146 }
147147}
@@ -173,7 +173,7 @@ int main()
173173 std::mdspan X(X_vec.data(), N, N);
174174 std::mdspan B(B_vec.data(), N, N);
175175
176- init_mat (A)
176+ init_tria_mat (A);
177177 init_mat(B);
178178
179179 // (1)
@@ -197,7 +197,7 @@ int main()
197197 B,
198198 X,
199199 std::divides<void>{});
200- print (X);
200+ print_mat (X);
201201
202202 // (3)
203203 std::cout << "(3)\n";
@@ -207,7 +207,7 @@ int main()
207207 std::linalg::implicit_unit_diagonal,
208208 B,
209209 X);
210- print (X);
210+ print_mat (X);
211211
212212 // (4)
213213 std::cout << "(4)\n";
@@ -218,13 +218,12 @@ int main()
218218 std::linalg::implicit_unit_diagonal,
219219 B,
220220 X);
221- print (X);
221+ print_mat (X);
222222
223223 return 0;
224224}
225225```
226226* A.extent[ link /reference/mdspan/extents/extent.md]
227- * v.extent[ link /reference/mdspan/extents/extent.md]
228227* std::mdspan[ link /reference/mdspan/mdspan.md]
229228* std::extents[ link /reference/mdspan/extents.md]
230229* std::linalg::layout_blas_packed[ link /reference/linalg/layout_blas_packed.md]
0 commit comments