function out=veclength(v)
out=sqrt(realdot(v,v));