Is that possible to add a norm function to use the L1 and L infinite norms instead of L2 (euclidian norm) ?
These norms have often the advantage of being faster (no float product and square root)

Something like norm1, norminf or length1, lengthinf
Then for coherence, it would be nice to define an alias to the current length function named along length2 or norm2