R/norm.lp.R

norm.lp <-
function(x,lp=2){
	if (lp<0) stop('p>0');
	if (lp==0) return(sum(abs(x)>10*.Machine$double.eps));
	sum(x^lp)^(1/lp)
}
yfyang86/paftoga documentation built on May 4, 2019, 2:33 p.m.