void Bistrya(int m, int n)
{
i=m;
k=n;
dopoln=matriza[(i+k)/ 2 ][j];
do 
{
  while (fabs(matriza[i][j])<fabs(dopoln) && i<=n)i++;
  while (fabs(matriza[k][j])>fabs(dopoln) && k>=m)k--;
       if (i<=k)
		{
		prom=matriza[i][j];
		matriza[i][j]=matriza[k][j];
		matriza[k][j]=prom;
		i++;k--;srav++;per++;
		}

}
while (i<=k);
if (m < k) 
Bistrya(m, k); 
   if (i < n)
     Bistrya(i, n); 

typedef int (*compare)(void *, void* );   // указатель на функцию сравнения

// в нашем случае функция сравнения такая
int my_cpm(void* left, void* right)
{
  double l = fabs(*((double*) left));
  double r = fabs(*((double*) right));
  return l > r ? - 1  : l < r;
}

void swap(char* lptr, char* rptr, size_t sz)
{
  char tmp;
  for (int i =  0 ; i < sz; ++i)
  {
    tmp = lptr[i];
    lptr[i] = rptr[i];
    rptr[i] = tmp;
  }
}

int partition(void* ptr, size_t rows, size_t cols, size_t col, size_t sz, compare cmp)
{
  char* pptr = (char*) ptr + col * sz;
  int i = - 1 ;
  int j = rows;
  while (true)
  {
     char* rptr, lptr =  0 ;
     for(--j; j > 0 ; --j)
     {
        rptr = pptr + j * cols * sz;
        if (cmp(pptr, rptr) >= 0 )
          break;
     }
     for(++i; i != j; ++i)
     {
        lptr = pptr + i * cols * sz;
        if (cmp(pptr, lptr) <= 0 )
          break;
     }
     if (i < j)
       swap(i, j, sz);
     else 
       return j;
  }
}

void column_qsort(void* ptr, size_t rows, size_t cols, size_t col, size_t sz, compare cmp)
{
  if (rows <=  1 )
    return;
  int qrow = partition(pptr, rows, cols, col, sz, cmp);
  column_qsort(pptr, qrow, cols, col, sz, cmp);
  column_qsort(((char*)pptr + (qrow +  1 ) * cols * sz), (rows - qrow), cols, col, sz, cmp)
}

int main()
{
  double matrix[ 10 ][ 5 ];
  for ( int i =  0 ; i <  10 ; ++i)
    column_qsort(matrix,  10 ,  5 , i, sizeof(double), my_cmp); 
}