Обработка массивов

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
 
const int A = 441;
const double PI = 3.1415926535897931159979635;
 
inline void dwarf_sort(double* array, int size) {
    size_t i = 1;
    while (i < size) {
        if (i == 0) {
            i = 1;
        }
        if (array[i - 1] <= array[i]) {
            ++i;
        }
        else
        {
            long tmp = array[i];
            array[i] = array[i - 1];
            array[i - 1] = tmp;
            --i;
        }
    }
}
 
inline double reduce(double* array, int size) {
 
    size_t i;
    double min = RAND_MAX, sum = 0;
 
    for (i = 0; i < size; ++i) {
        if (array[i] < min && array[i] != 0) {
            min = array[i];
        }
    }
 
    for (i = 0; i < size; ++i) {
        if ((int)(array[i] / min) % 2 == 0) {
            sum += sin(array[i]);
        }
    }
    return sum;
}
 
int main(int argc, char* argv[])
{
    int i, N, j;
    double* M1 = NULL, * M2 = NULL, * M2_copy = NULL;
    double X;
    unsigned int seed = 0;
 
    N = atoi(argv[1]); /* N равен первому параметру командной строки */
    M1 = malloc(N * sizeof(double));
    M2 = malloc(N / 2 * sizeof(double));
    M2_copy = malloc(N / 2 * sizeof(double));
 
    for (i = 0; i < 100; i++) 
    {
        seed = i;
        srand(i); 
        
        /*generate*/
        for (j = 0; j < N; ++j) {
            M1[j] = (rand_r(&seed) % A) + 1;
        }
 
        for (j = 0; j < N / 2; ++j) {
            M2[j] = (rand_r(&seed) % (10 * A)) + 1;
        }
 
        /*map*/
        for (j = 0; j < N; ++j)
        {
            M1[j] = pow(M1[j] / PI, 3);
        }
 
        for (j = 0; j < N / 2; ++j) {
            M2_copy[j] = M2[j];
        }
 
        M2[0] = fabs(tan(M2_copy[0]));
        for (j = 0; j < N / 2; ++j) {
            M2[j] = fabs(tan(M2[j] + M2_copy[j]));
        }
 
        /*merge*/
        for (j = 0; j < N / 2; ++j) {
            M2[j] = pow(M1[j], M2[j]);
        }
 
        /*sort*/
        dwarf_sort(M2, N / 2);
 
        /*sort*/
        X = reduce(M2, N / 2);
    }
    
    printf("\nN=%d.\n", N); 
    printf("X=%f\n", X);
    return 0;
}