[v0.6.3.9000] Singular Matrices in c(1, gen)

[Chris1221]

Dealing with this by Moore-Penrose pseudo-inverse of matrix.

arma::pinv()

[Chris1221]

Fixed by reporting 0 if vector is all 0, 1, or 2. See here.

bool zero = arma::all(X == 0);
bool one = arma::all(X == 1);
bool two = arma::all(X == 2);

bool singular = (zero || one || two);