Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save wrathematics/1c0dedae491458cb816372a099b6a881 to your computer and use it in GitHub Desktop.
Save wrathematics/1c0dedae491458cb816372a099b6a881 to your computer and use it in GitHub Desktop.
static inline int get_rank(int m, int n, double *qr, double tol)
{
int i;
const double minval = fabs(tol*qr[0]);
const int minmn = m > n ? n : m;
for (i=1; i<minmn; i++)
{
if (fabs(qr[i + m*i]) < minval)
return i;
}
return minmn;
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment