https://www.intel.com/content/www/us/en/docs/onemkl/developer-reference-fortran/2025-0/lags2.html
?lags2
Computes 2-by-2 orthogonal matrices U, V, and Q, and applies them to matrices A and B such that the rows of the transformed A and B are parallel.