Just assume P1 is at the origin, and the body is not rotated. Then the split operation is easy, P2 will be at -l/4 and P3 at l/4 along the split axis.
Next you apply the rotation matrix to the new points, and then you add the real position of P1. This should give you the correct points.
For example, let’s assume you want to split along the x axis, M is the rotation matrix, P1 is the original center and l is the length along the x axis:
P2 = M * (-l/4,0,0) + P1
P3 = M * (+l/4,0,0) + P1
Note that the matrix multiplication is the same as just taking the first column of the matrix and multiplying it by ±l/4, because the other vector components are zero.