result.push_back( std::polar(rho[i], theta[i] * cctbx::constants::pi_180));
result.push_back(std::polar(rho[i], theta[i]));
  py_polar_complex(    cctbx::af::shared<double> rho,    cctbx::af::shared<double> theta,    bool deg = false)  {    cctbx_assert(rho.size() == theta.size());    cctbx::af::shared<std::complex<double> > result;    result.reserve(rho.size());    if (deg) {      for(std::size_t i=0;i<rho.size();i++) {        result.push_back(std::polar(rho[i], theta[i]));      }    }    else {      for(std::size_t i=0;i<rho.size();i++) {        result.push_back(          std::polar(rho[i], theta[i] * cctbx::constants::pi_180));      }    }    return result;  }