Changeset 15104 for issm/trunk-jpl/src/c/modules/Krigingx/pKrigingx.cpp
- Timestamp:
- 05/24/13 19:11:47 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/modules/Krigingx/pKrigingx.cpp ¶
r15100 r15104 68 68 /*partition loop across threads: */ 69 69 for(int idx=my_rank;idx<n_interp;idx+=num_procs){ 70 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"% " << "\n");70 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n"); 71 71 observations->InterpolationKriging(&predictions[idx],&error[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata,variogram); 72 72 } 73 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"% " << "\n");73 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n"); 74 74 75 75 #ifdef _HAVE_MPI_ … … 86 86 /*partition loop across threads: */ 87 87 for(int idx=my_rank;idx<n_interp;idx+=num_procs){ 88 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"% " << "\n");88 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n"); 89 89 observations->InterpolationV4(&predictions[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata); 90 90 } 91 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"% " << "\n");91 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n"); 92 92 93 93 #ifdef _HAVE_MPI_ … … 101 101 /*partition loop across threads: */ 102 102 for(int idx=my_rank;idx<n_interp;idx+=num_procs){ 103 _printf0_(" interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"% " << "\n");103 _printf0_(" interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n"); 104 104 observations->InterpolationNearestNeighbor(&predictions[idx],x_interp[idx],y_interp[idx],radius); 105 105 } 106 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"% " << "\n");106 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n"); 107 107 108 108 #ifdef _HAVE_MPI_ … … 118 118 /*partition loop across threads: */ 119 119 for(int idx=my_rank;idx<n_interp;idx+=num_procs){ 120 _printf0_(" interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"% " << "\n");120 _printf0_(" interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%\n"); 121 121 observations->InterpolationIDW(&predictions[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata,power); 122 122 } 123 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"% " << "\n");123 _printf0_(" interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%\n"); 124 124 125 125 #ifdef _HAVE_MPI_ … … 142 142 143 143 MPI_Barrier(MPI_COMM_WORLD); finish=MPI_Wtime( ); 144 _printf_("\n " << setw(34) << left << "Observation fitering elapsed time: " << finish_init-start_init << " seconds \n " << "\n");145 _printf_(" " << setw(34) << left << "Kriging prediction elapsed time: " << finish_core-start_core << " seconds \n " << "\n");146 _printf_("\n " << "Total elapsed time " << int((finish-start)/3600) << " hrs " << int(int(finish-start)%3600/60) << " min " << int(finish-start)%60 << " sec\n\n " << "\n");144 _printf_("\n " << setw(34) << left << "Observation fitering elapsed time: " << finish_init-start_init << " seconds \n\n"); 145 _printf_(" " << setw(34) << left << "Kriging prediction elapsed time: " << finish_core-start_core << " seconds \n\n"); 146 _printf_("\n " << "Total elapsed time " << int((finish-start)/3600) << " hrs " << int(int(finish-start)%3600/60) << " min " << int(finish-start)%60 << " sec\n\n\n"); 147 147 return 1; 148 148 #else
Note:
See TracChangeset
for help on using the changeset viewer.