Oleg Nesterov improvement for la.diag.

grame-cncm · Nov 26, 2024 · 537bb1d · 537bb1d
1 parent 6f0e6e5
commit 537bb1d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/linearalgebra.lib b/linearalgebra.lib
@@ -245,6 +245,6 @@ declare identity copyright "MIT License";
 //
 // * `N`: The size of each axis of the matrix.
 //---------------------------------------------------------
-diag(N) = si.bus(N) <: par(i, N, par(j, N, _*(select2(i == j,0,1))));
+diag(N) = route(N, N*N, par(i,N, (i+1, (N+1)*i+1)));
 declare diag author "David Braun";
 declare diag copyright "MIT License";