Reversible Karhunen-Loeve Transform

I have a problem with implementing RKLT. I've read through papers by P开发者_如何学Cengwei Hao and Qingyun Shi but i still have some questions. Has anybody seen an implementation of this transform? I don't know how the permutation matrix should look like and how to obtain SERM factorization from the TERM factorization.

---

Here's a C implementation od PLUS factorization by dr Pengwei Hao:

http://www.eecs.qmul.ac.uk/~phao/PLUS/plus.c