getf SCI/util/testexamples.sci
 
reinit_for_test()
 
%U=mopen('SCI/tests/automatic_tests/chsolve_data.ref','rb');
 
A = sprand(20, 20, 0.1);
 
A = A * A' + eye();
 
spcho = chfact(A);
 
sol = (1:20)';rhs = A * sol;
 
spcho = chfact(A);
 
%ans = chsolve(spcho, rhs);
 
if load_ref('%ans') then   bugmes();quit;end,
 
 
xdel_run(winsid());
 
 
mclose(%U);