NoDiag on large arrays

[dhutchis]

The NoDiag function throws an error in Octave when passed a large Associative Array. The culprit is line 15: ij = sub2ind(size(A),i,j);. I ran into this error while processing a SCALE 18 power law graph.

Replacing this line with AA = A - diag(A); solves the problem. Anyone have an objection to using this implementation?

[dhutchis]

Somehow this doesn't always work. I observed differences at SCALE 10 with this change; some diagonal elements slipped through. I wonder why, and what a better change would be.