Rotate a point by another point in 2D

If you rotate point (px, py) around point (ox, oy) by angle theta you’ll get:

p'x = cos(theta) * (px-ox) - sin(theta) * (py-oy) + ox
p'y = sin(theta) * (px-ox) + cos(theta) * (py-oy) + oy