Function ecall_reducemax_one

sgx_tf_ops/Enclave/Enclave.cpp:2571–2634

}

void ecall_reducemax_one(float *input,int N,int M,int C,int L,float *output){
    float *x = (float*)malloc(sizeof(float)*N);
    memset(x, 0, sizeof(float)*N);
    for(int j=0;j<N;j++){
        x[j] = 0.0f;
        float sum = input[insert*N+j];
        float c = 0.0f;
        for(int i=0;i<Nt;i++){
            float y = input[indexs[i]*N+j] - c;
            float t = sum + y;
            c = (t - sum) - y;
            sum = t;
        }
        x[j] = sum;
        
    }
    float *x1 = (float*)malloc(sizeof(float)*M*L);
    memset(x1, 0, sizeof(float)*M*L);

    for (int i = 0; i < M; i++) {
        for (int j = 0; j < L; j++) {
            float max_val = x[i*C*L+j];
            for (int k = 1; k < C; k++) {
                max_val = fmaxf(max_val, x[i*C*L+j+k*L]);
            }
            x1[i*L+j] = max_val;
        }
    }

    float mean = 0.0f, std, mean_sqr = 0.0f;
    float c_mean = 0.0f, c_mean_sqr = 0.0f;

    for (int j = 0; j < M*L; j++) {
        float y = (x1[j]/static_cast<float>(M*L)) - c_mean;
        float t = mean + y;
        c_mean = (t - mean) - y;
        mean = t;

        float y2 = (x1[j]*x1[j]/static_cast<float>(M*L)) - c_mean_sqr;
        float t2 = mean_sqr + y2;
        c_mean_sqr = (t2 - mean_sqr) - y2;
        mean_sqr = t2;
    } 
    std = sqrtf(mean_sqr - mean * mean)/static_cast<float>(Nt);
    mean = mean/static_cast<float>(Nt);
    for(int i=0;i<M*L*(Ne-1);i++){
        output[i] =(x1[i%(M*L)] - gaussrand(mean, std))/static_cast<float>(Nt+1);    
    }
    for(int j=M*L*(Ne-1);j<M*L*Ne;j++){
        output[j] = 0.0f;
        float sum = 0.0f;
        float c = 0.0f;
        for(int i=0;i<Nt;i++){
            float y = -1.0f * output[indexs[i]*M*L+(j-M*L*(Ne-1))] - c;
            float t = sum + y;
            c = (t - sum) - y;
            sum = t;
        }