Mercurial > forge

changeset 11693:213fe8b0c3d7 octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
bug fix to csaps_sel (Binv had not been made sparse)
author nir-krakauer
date Thu, 09 May 2013 19:37:51 +0000
parents 4f0ae470d973
children c79a6e375d83
files main/splines/inst/csaps_sel.m
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/main/splines/inst/csaps_sel.m	Thu May 09 15:40:14 2013 +0000
+++ b/main/splines/inst/csaps_sel.m	Thu May 09 19:37:51 2013 +0000
@@ -195,7 +195,7 @@
 
 function Binv = banded_matrix_inverse(d, U, m) #given a (2m+1)-banded, symmetric n x n matrix B = U'*inv(diag(d))*U, where U is unit upper triangular with bandwidth (m+1), returns Binv, a sparse symmetric matrix containing the central 2m+1 bands of the inverse of B
 #Reference: Hutchinson and de Hoog 1985
-  Binv = diag(d);
+  Binv = sparse(diag(d));
   n = rows(U);
   for i = n:(-1):1
     p = min(m, n - i);