Changes to ilund:
An option should allow variable limits on the fill (the maximum number
of elements).  One possible value for this array would be idx=0;
idx(nxn)++; idx = maxn * max(1,idx); (that is, the number of incident
nodes in the flow ordering).  

For David's test problems, we need to add a source and/or sink problem.
