Cuspatial

import cuspatial

cgs = cuspatial.GeoSeries.from_points_xy(cdf[["longitude", "latitude"]].interleave_columns())

cgdf = cuspatial.from_geopandas(gdf)

cgdf["geometry"].within(cgeom)