The documentation for the Sparse BLAS routines is somewhat flawed:
For example for ?gthr (http://software.intel.com/sites/products/documentation/hpc/mkl/webhelp/bla/functn_gthr.html), the
Fortran 95 Interface Notes
section says "Holds the vector with the number of elements nz" for every argument.
This isn't really helpful.
And my question for gthr is, if size(x) has to be equal to size(indx) or if the the gathering stop if the smaller size limit is reached?
Fortran 95 Interface NotesDF