Compact numerical solutions to the two-dimensional repulsive Hubbard model obtained via nonunitary similarity transformations