matthiasm@8: function [outmat, scale] = qnormalise(inmat, varargin)
matthiasm@8: % qnormalise(inmat, q, dim) returns the field inmat normalised by the q
matthiasm@8: % norm over dimension dim
matthiasm@8: dim = 1;
matthiasm@8: if nargin > 1
matthiasm@8:     q = varargin{1};
matthiasm@8:     if nargin > 2
matthiasm@8:         dim = varargin{2};
matthiasm@8:     end
matthiasm@8: else
matthiasm@8:     q = 2;
matthiasm@8: end
matthiasm@8: sizeInmat = size(inmat);
matthiasm@8: nDim = length(sizeInmat);
matthiasm@8: repeatpattern = [ones(1,dim-1) sizeInmat(dim) ones(1,nDim-dim)];
matthiasm@8: scale = qnorm(inmat, q, dim);
matthiasm@8: outmat = inmat ./ repmat(scale, repeatpattern);