Ignore:
Timestamp:
09/14/12 16:29:23 (13 years ago)
Author:
Mathieu Morlighem
Message:

BUG: Fixed parallel kriging with new iomodel->dependents

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/modules/Krigingx/pKrigingx.cpp

    r13216 r13376  
    5656                /*partition loop across threads: */
    5757                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    58                         _pprintString_("\r      interpolation progress: "<<setw(6)<<setprecision(2)<<double(idx)/double(n_interp)*100<<"%");
     58                        _pprintLine_("      interpolation progress: "<<setw(6)<<setprecision(2)<<double(idx)/double(n_interp)*100<<"%");
    5959                        observations->InterpolationKriging(&predictions[idx],&error[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata,variogram);
    6060                }
    61                 _pprintLine_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%");
     61                _pprintLine_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%");
    6262
    6363#ifdef _HAVE_MPI_
     
    7474                /*partition loop across threads: */
    7575                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    76                         _pprintString_("\r      interpolation progress: "<<setw(6)<<setprecision(2)<<double(idx)/double(n_interp)*100<<"%");
     76                        _pprintLine_("      interpolation progress: "<<setw(6)<<setprecision(2)<<double(idx)/double(n_interp)*100<<"%");
    7777                        observations->InterpolationNearestNeighbor(&predictions[idx],x_interp[idx],y_interp[idx],radius);
    7878                }
    79                 _pprintLine_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%");
     79                _pprintLine_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%");
    8080
    8181#ifdef _HAVE_MPI_
     
    9191                /*partition loop across threads: */
    9292                for(int idx=my_rank;idx<n_interp;idx+=num_procs){
    93                         _pprintString_("\r      interpolation progress: "<<setw(6)<<setprecision(2)<<double(idx)/double(n_interp)*100<<"%");
     93                        _pprintLine_("      interpolation progress: "<<setw(6)<<setprecision(2)<<double(idx)/double(n_interp)*100<<"%");
    9494                        observations->InterpolationIDW(&predictions[idx],x_interp[idx],y_interp[idx],radius,mindata,maxdata,power);
    9595                }
    96                 _pprintLine_("\r      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%");
     96                _pprintLine_("      interpolation progress: "<<fixed<<setw(6)<<setprecision(2)<<100.<<"%");
    9797
    9898#ifdef _HAVE_MPI_
Note: See TracChangeset for help on using the changeset viewer.