view general/numerical/pnormalise.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

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);