samer@4: function M=data2maps(RFn,X,N)
samer@4: % data2maps - build linear maps for covering data ranges
samer@4: %
samer@4: % data2maps :: 
samer@4: %    ([[N,M]] -> [[2,M]]) ~ 'function to compute min and max of range',
samer@4: %    [[N,M]]              ~ 'the data'
samer@4: %    L:natural            ~ 'number of bins in output maps'
samer@4: % -> { [M] -> dmap(L) }   ~ 'one dmap per column in a cell array'.
samer@4: %
samer@4: % Two argument version computes L=sqrt(N)
samer@4: %
samer@4: % data2maps :: 
samer@4: %    ([[N,M]] -> [[2,M]]) ~ 'function to compute min and max of range',
samer@4: %    [[N,M]]              ~ 'the data'
samer@4: % -> { [M] -> dmap(L) }   ~ 'one dmap per column in a cell array'.
samer@4: 
samer@4: if nargin<3, N=sqrt(size(X,1)); end
samer@4: M=maprows(@(r){linmap(r(1),r(2),N)}, closed2hopen(feval(RFn,X)'));
samer@4: