Function emb_ss

sgx_tf_ops/Enclave/Enclave.cpp:106–151

}

void emb_ss(float* input,float* output,int M,int C,int num,int dim){
    int N = M*num*dim;
    float *x = (float*)malloc(sizeof(float)*N);
    memset(x, 0, sizeof(float)*N);
    for (int i = 0; i < M; i++) {
        for (int j = 0; j < num; j++) {
            int idx = static_cast<int>(fabs(input[i*C+j*(dim+1)+dim] * 100))%10;
            for (int k = 0; k < dim; k++) {
                x[i*num*dim+j*dim+k] = input[i*C+j*(dim+1)+k] + sgx_arr[idx][k%128];
            }
            }
        }
    float mean = 0.0f, std, mean_sqr = 0.0f;
    float c_mean = 0.0f, c_mean_sqr = 0.0f;

    for (int j = 0; j < N; j++) {
        float y = (x[j]/static_cast<float>(N)) - c_mean;
        float t = mean + y;
        c_mean = (t - mean) - y;
        mean = t;

        float y2 = (x[j]*x[j]/static_cast<float>(N)) - 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<N*(Ne-1);i++){
        output[i] =(x[i%N] - gaussrand(mean, std))/static_cast<float>(Nt+1);    
    }
    for(int j=N*(Ne-1);j<N*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]*N+(j-N*(Ne-1))] - c;
            float t = sum + y;
            c = (t - sum) - y;
            sum = t;
        }
        output[j] =x[j-N*(Ne-1)] + sum;
    }
    free(x);
    return;
}

void emb_ss_grad(float* grad,float* output,int M,int C,int num,int dim){
    int N = M*num*dim;