NB. dtype is here an argument to asmatrix, not to ones as above