def in_domain(self, x):
x = x.coordinates if isinstance(x, Mesh) else np.asarray(x)<EOL>if pimms.is_vector(x): return np.asarray(self.in_domain([x])[<NUM_LIT:0>], dtype=np.bool)<EOL>if x.shape[<NUM_LIT:0>] != <NUM_LIT:3>: x = x.T<EOL>if self.radius is None: return np.ones(x.shape[<NUM_LIT:1>], dtype=np.bool)<EOL>x = self.alignment_matrix.dot(np.concatenate((x, np.ones((<NUM_LIT:1>,x.shape[<NUM_LIT:1>])))))<EOL>x = np.clip(x[<NUM_LIT:2>] / self._sphere_radius, -<NUM_LIT:1>, <NUM_LIT:1>)<EOL>th = np.arccos(x)<EOL>rad = pimms.mag(self.radius, '<STR_LIT>')<EOL>return (th < rad)<EOL>
proj.in_domain(x) yields a boolean array whose elements indicate whether the coordinates in
  x are part of the domain of the given map projection. The argument x may be either a
  coordinate matrix, a coordinate vector, or a mesh (in which case this is equivalent to
  proj.in_domain(x.coordinates).