Changeset 17668

Timestamp:

04/08/14 10:25:18 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added distance output for pKriging

File:

• issm/trunk-jpl/src/c/modules/Krigingx/pKrigingx.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/Krigingx/pKrigingx.cpp

@@ -108 +108 @@
                 xDelete<double>(predictions); predictions=sumpredictions;
         }
+        else if(strcmp(output,"distance")==0){
+
+                /*partition loop across threads: */
+                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
+                        _printf0_("      interpolation progress: "<<setw(6)<<setprecision(4)<<double(idx)/double(n_interp)*100.<<"%  \n");
+                        observations->Distances(&predictions[idx],&x_interp[idx],&y_interp[idx],1,radius);
+                }
+                _printf0_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(4)<<100.<<"%  \n");
+
+                double *sumpredictions =xNew<double>(n_interp);
+                ISSM_MPI_Allreduce(predictions,sumpredictions,n_interp,ISSM_MPI_PDOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
+                xDelete<double>(predictions); predictions=sumpredictions;
+        }
         else if(strcmp(output,"idw")==0){
                 double power;