Cool.  Some suggestions for this routine:

1) get rid of the maximum number of non-zeros per row arguments, they are stupid and just make this calling sequence too damn complicated looking.

2) the manual page doesn't make clear these are nonzero BLOCKS per row, not nonzero entries per row.

  Barry


  Input Arguments:
+  A - matrix being preallocated
.  bs - block size
.  dnz - maximum number of nonzero entries per row of diagonal block of parallel matrix
.  dnnz - number of nonzeros per row of diagonal block of parallel matrix
.  onz - maximum number of nonzero entries per row of off-diagonal block of parallel matrix
.  onnz - number of nonzeros per row of off-diagonal block of parallel matrix
.  dnzu - maximum number of nonzero entries per row of upper-triangular part of diagonal block of parallel matrix
.  dnnzu - number of nonzeros per row of upper-triangular part of diagonal block of parallel matrix
.  onzu - maximum number of nonzero entries per row of upper-triangular part of off-diagonal block of parallel matrix
-  onnzu - number of nonzeros per row of upper-triangular part of off-diagonal block of parallel matrix