# convert to sgn(x); tosgn(abs(x_)) := sgn(x)*x; tosgn(theta(x_)) := sgn(x)/2+1/2;