int main() { double Delta_i, Delta_f, U_i; double* y; Delta_i=0.8; Delta_f=1; N=2*N_half+1; Ns=N*N; y=(double*)malloc(sizeof(double)*4*Ns); Delta=Delta_f; U_f=inity(y); Delta=Delta_i; U_i=inity(y); printf("Delta_i=%lf U_i=%lf Delta_f=%lf U_f=%lf\n",Delta_i,U_i,Delta_f,U_f); evolution(y); free(y); return 0; }
int main() { double* y; long int N_half,N,Ns; double U; N_half=400; N=2*N_half+1; Ns=N*N; y=(double*)malloc(sizeof(double)*4*Ns); U=inity(y); free(y); return 0; }
int main() { double* y; double t; N_half=N/2; Ns=N*N; t=0; y=(double*) malloc(sizeof(double)*16*Ns); //provide y at t=0,and U. comment off if the values read from data. U=inity(y); printf("U=%f\n",U); omp_set_num_threads(3); evolution(y,t); free(y); return 0; }
int main(int argc, char **argv){ int N = 64; int M = 28815; int O = 28768; int N2 = N/2; int L = (M/N2); int P = 4; clock_t start; clock_t stop; int i; int j = 0; int k = 0; int n; int idx; double x1[N]; double y1[N]; double y[O]; double *temparray = (double *) malloc(sizeof(double) * N); double *G = (double *) malloc(sizeof(double) * N); double *D = (double *) calloc(N,sizeof(double)); double *past_tail = (double *) calloc(N2,sizeof(double)); double X1[2*N+1]; double S[2*N+1]; double tempy1[2*N+1]; for (k = 0; k < P; k++) { n = N2*k; memcpy(temparray,x+(n),sizeof(double)*N); windowingFilter(h,temparray,x1,N); preFFT(x1,X1,N); modelfft(X1, N, 1); NoiseSpecterEstimate(D,X1,P,N); } inity(y,P*N2); idx = P*N2; for (k = P; k < L-1; k++) { n = N2*k; memcpy(temparray,x+(n),sizeof(double)*N); windowingFilter(h,temparray,x1,N); preFFT(x1,X1,N); modelfft(X1, N, 1); start = clock(); SSMethod(S,G,D,X1,N); stop = clock(); printf("Time start: %f\n", (double)start); printf("Time stop: %f\n", (double)stop); printf("Time elapsed: %f\n", ((double)stop - start) / CLOCKS_PER_SEC); memcpy(tempy1,S,sizeof(double)*2*N+1); modelfft(tempy1, N, -1); postIFFT(tempy1,y1,N); cascadeOUT(y,past_tail,y1,idx,N); idx+=(N/2); j = 0; for (i = N2;i < N;i++) { past_tail[j] = y1[i]; j++; } } // FILE *fp; // fp = fopen("result.dat","w"); // for(i=0; i<O; i++) { // fprintf(fp,"%.8f\n",y[i]); // } // fclose(fp); // for(i=O-100; i<O; i++) { // printf("%.8f\n",y[i]); // } free(temparray); free(G); free(D); free(past_tail); return 0; }