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 Notessection 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