d1 = r*(cos(a), sin(a),0)-p0;
in_out = step(0,zcomp(d0^d1));