This function simplifies the benchmarking of ILU algorithms.
Examples
marray=spilu_iluhubavail()// A small 3-by-3 matrix// nnz(A)=7A=[120345067];A=sparse(A);// Test all methodsform=spilu_iluhubavail()mprintf("Method: %s\n",m);[L,U,perm]=spilu_iluhub(A,m);mprintf(" nnz(LU)= %d\n",nnz(L)+nnz(U));end