axisL2norm(x,normAxis)

% function [xnorm, xnormalized] = axisL2norm(x,normAxis)
% Normalize an N-d tensor along a given axis using the L2 (euclidean) norm

function [xnormalized, xnorm] = axisL2norm(x,normAxis)

dims = size(x);
Ndims = length(dims);

rpt = ones(1,Ndims);
rpt(normAxis) = dims(normAxis);
xnorm = sqrt(sum(x.^2,normAxis));
xnormalized = x ./ repmat(xnorm,rpt);