[routines][generic] Add functionality to allow functions on 1-d array to work on any axis

[forfudan]

1. Adds functions that are able to apply any functions working on 1-d array to any axis, with or without dimension reduction.

   • Add apply_func_on_array_with_dim_reduction() and apply_func_on_array_without_dim_reduction(). They try to utilize parallelization as much as possible.
   • In future, we only need to focus on writing (and optimizing) functions for 1-d arrays. Operating along certain axis can be easily achieved by applying the function, e.g., `apply_func_on_array_with_dim_reduction[max_1d](array, axis=axis).

2. Implements this approach on functions in statistics.averages module and on the sort function. The sort function gain speed increase compared to the old method and is quicker than numpy.sort for large arrays.

[forfudan]

@shivasankarka Let me resolve conflicts first.

[forfudan]

@shivasankarka Conflicts resolved!