Procedure from library gaussman.lib (see gaussman_lib).
gaussman.lib
ncols(a)==size(m)==size(V); typeof(V[i])=="int"
order (a[i][,V[i]]) with multiplicity m[i] lexicographically