PlotSparse does not works anymore
Reported by Clement DAVID (@davidcl)
ERROR LOG:
----------
--> PlotSparse(U,"x");
at line 2 of function is_digit
at line 17 of function ana_style ( /home/davidcl/work/branches/master-review/scilab/modules/umfpack/macros/ana_style.sci line 31 )
at line 59 of function PlotSparse ( /home/davidcl/work/branches/master-review/scilab/modules/umfpack/macros/PlotSparse.sci line 73 )
Undefined variable: _str2code
HOW TO REPRODUCE THE BUG:
-------------------------
// from ordmmd.xml help page
A = [
4. 1. 2. 0.5 2.
1. 0.5 0. 0. 0.
2. 0. 3. 0. 0.
0.5 0. 0. 5./8. 0.
2. 0. 0. 0. 16.
];
A = sparse(A);
// See the sparsity pattern of the Cholesky factors of A
U = sparse(chol(full(A)));
scf();
subplot(2,1,1);
PlotSparse(U,"x");
xtitle("Sparsity pattern of U, such that A=U''*U");
// Get the adjacency structure
[xadj,iadj,val]=sp2adj(A);
// Compute multiple minimum degree ordering
n = size(A,"r");
[perm,invp,nofsub]=ordmmd(xadj,iadj,n);
// Convert the permutation vector into matrix.
P=spzeros(n,n);
for i=1:n
P(perm(i),i)=1;
end
// See the sparsity pattern of the Cholesky factors
// of P'*A*P
U = sparse(chol(full(P'*A*P)));
subplot(2,1,2);
PlotSparse(U,"x");
xtitle("Sparsity pattern of U, such that P''*A*P=U''*U");
OTHER INFORMATION:
------------------