I have been looking into the generalisation of rotation about an axis in 3D. Which I have found to be is rotation of a vector in N dimensional space about an N-2 dimesnional subspace. The paper General n-Dimensional Rotations introduces a very interesting algorithm to compute such a matrix and do such transformation. But during the change of basis(the initial rotations in order to embed the N-2 subspace into $X_1 X_2 X_3 ... X_{N-2}$) they use the coordinates of the points that are produced after the applications of each rotation matrix.
A matrix is formed which contains the coordinates of all the (N-1) elements.
The change of basis here makes the columns of the last two elements equal to zero due to which the said matrix can be embedded in N-2 dimensional space. The transformation can be written as
For making the $v^{k}_{r,c}$ equal to 0 we use the matrix given by
Now as you can observe the rotation matrix takes in as input $v^{k-1}_{r,c}$ which is the coordinates after the application of the last rotation matrix. This procedure is repeated until most of the coordinates of the starting plane(the plane we wanted to do rotation about) is zero and then we apply the rotation about this plane. And finally we apply the inverse of the basis change.
My Question - Is there any method in which we do not have to check the coordinates obtained at the intermediate steps, rather we can at the beginning look at the coordinates of the plane and decide the input to all the rotation matrices, to do the basis change and finally the rotation about the subspace.
Perhaps this is somehow possible by using spherical coordinates, I have been trying but not able to find any such method.


