https://github.com/rapidsai/cuspatial
import cuspatial
cgs = cuspatial.GeoSeries.from_points_xy(cdf[["longitude", "latitude"]].interleave_columns())
cgdf = cuspatial.from_geopandas(gdf)
cgdf["geometry"].within(cgeom)