Fullwave design of cm-scale cylindrical metasurfaces via fast direct solvers