Повернуть вектор на 90 градусов через перекрестное произведение

Question

Я пытаюсь понять код, в котором авторы повернули векторы в 3D на 90 градусов. Для этого они использовали перекрестный продукт. Как кросс-продукт может включать вращение на 90 градусов? Любое объяснение будет полезно.

    % Inputs:
    %   V  #vertices by dim list of mesh vertex positions
    %   F  #faces by simplex-size list of mesh face indices
    % Gradient of a scalar function defined on piecewise linear elements (mesh)
    % is constant on each triangle i,j,k:
    % grad(Xijk) = (Xj-Xi) * (Vi - Vk)^R90 / 2A + (Xk-Xi) * (Vj - Vi)^R90 / 2A
    % grad(Xijk) = Xj * (Vi - Vk)^R90 / 2A + Xk * (Vj - Vi)^R90 / 2A + 
    %             -Xi * (Vi - Vk)^R90 / 2A - Xi * (Vj - Vi)^R90 / 2A
    % where Xi is the scalar value at vertex i, Vi is the 3D position of vertex
    % i, and A is the area of triangle (i,j,k). ^R90 represent a rotation of 
    % 90 degrees
    %
    % renaming indices of vertices of triangles for convenience
    i1 = F(:,1); i2 = F(:,2); i3 = F(:,3); 
    % #F x 3 matrices of triangle edge vectors, named after opposite vertices
    v32 = V(i3,:) - V(i2,:);  v13 = V(i1,:) - V(i3,:); v21 = V(i2,:) - V(i1,:);

    % area of parallelogram is twice area of triangle
    % area of parallelogram is || v1 x v2 || 
    n  = cross(v32,v13,2); 

    % This does correct l2 norm of rows so that it contains #F list of twice
    % triangle areas
    dblA = normrow(n);

    % now normalize normals to get unit normals
    u = normalizerow(n);

    % rotate each vector 90 degrees around normal ---- ????
    eperp21 = bsxfun(@times,cross(u,v21),1./dblA);
    eperp13 = bsxfun(@times,cross(u,v13),1./dblA);

Answer 1

В 3d перекрестное произведение a x b двух векторов a и b приводит к вектору p := a x b, который перпендикулярен к a и b.

Это означает, что если вы умножите вектор с единичным вектором u, представляющим ось вращения, вы получите вектор, повернутый на 90 градусов вокруг поворота. ось.