Changeset 15229


Ignore:
Timestamp:
06/10/13 13:48:39 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: improved parallel kriging

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/contrib/gslib/pkriging.m

    r13730 r15229  
    66
    77options=pairoptions(varargin{:});
    8 cluster=getfieldvalue(options,'cluster',generic('np',10));
     8cluster=getfieldvalue(options,'cluster',generic('np',1));
    99options=removefield(options,'cluster',0);
    1010name   = ['krig' num2str(feature('GetPid'))];
     
    3535%Launch job on remote cluster
    3636BuildKrigingQueueScript(cluster,name,'',1,0,0); %gather, valgrind, gprof
    37 tic
    3837LaunchQueueJob(cluster,name,name,{[name '.bin'] [name '.queue']});
    39 toc
    40 choice=input('Is the job successfully completed? (y/n)','s');
     38
     39%Call waitonlock
     40md=model; md.cluster=cluster; md.settings.waitonlock=Inf;
     41waitonlock(md)
     42
     43%Download
    4144Download(cluster,name,{[name '.outbin']});
    4245structure=parseresultsfromdisk([name '.outbin'],0);
     
    5053
    5154%Process results
    52 B=structure.AutodiffForward;
     55B=structure.(EnumToString(0));
    5356B=reshape(B,size(x_interp,2),size(x_interp,1))';
    54 E=structure.AutodiffIsautodiff;
     57E=structure.(EnumToString(1));
    5558E=reshape(E,size(x_interp,2),size(x_interp,1))';
Note: See TracChangeset for help on using the changeset viewer.