Solved – Pseudo-random orthogonal matrix generation

A little while back, J.M. suggested using the Stewart algorithm for generating $n$ by $n$ pseudo random orthogonal matrices in $O(n^2)$ time. He further noted that this methodology is implemented in Nick Higham's matrix computation toolbox (matlab). Now this package contains a bunch of .m files. Could anyone direct me to which one implements the pseudorandom orthogonal matrix generator? Alternatively, would anyone know of a R implementation of this algorithm?

It's in the Test Matrix Toolbox, not the Matrix Computation Toolbox. The M-file qmult.m (premultiplication by a Haar-distributed pseudorandom orthogonal matrix) can be found here or here.

Similar Posts:

Rate this post

Leave a Comment