imsort.cpp

Go to the documentation of this file.

/*
 * $Id$
 * Author: David Fournier
 *
 * Copyright (c) 2009-2012 ADMB Foundation
 */
#include <fvar.hpp>

imatrix sort(const imatrix& m, int col, int NSTACK)
{
  ivector ind(m.rowmin(),m.rowmax());
  ivector ind1(m.rowmin(),m.rowmax());
  ivector ind2(m.rowmin(),m.rowmax());
  int i;
  for (i=m.rowmin();i<=m.rowmax();i++)
  {
    ind1(i)=m(i).indexmin();
    ind2(i)=m(i).indexmax();
  }
  //const ivector& iv=column(m,col);
  sort(column(m,col),ind);
  imatrix tmp(m.rowmin(),m.rowmax(),ind1,ind2);
  for (i=m.rowmin();i<=m.rowmax();i++)
  {
    tmp(i)=m(ind(i));
  }
  return tmp;
}