--- /dev/null
+#ifndef select_EMGrank_H
+#define select_EMGrank_H
+
+#include "ioutils.h"
+
+void EMGrank(
+ // IN parameters
+ const Real* Pi,
+ const Real* Rho,
+ Int mini,
+ Int maxi,
+ const Real* X,
+ const Real* Y,
+ Real tau,
+ const Int* rank,
+ // OUT parameters
+ Real* phi,
+ Real* LLF,
+ // additional size parameters
+ mwSize n,
+ mwSize p,
+ mwSize m,
+ mwSize k);
+
+#endif