function A=pnormalise(p,A,i)
% pnormalise - normalise A relative to p-norm in given array dimension
%
% pnormalise :: nonneg, [D:[[1 E]]], [E] -> [D].
if nargin<3, i=1; end
A=divnorm(@(x)pnorm(p,x,i),A);
	

