Changeset 642

Timestamp:

05/29/09 11:05:33 (15 years ago)

Author:

jschierm

Message:

Addition of parameter study output file parsing.

File:

• issm/trunk/src/m/solutions/dakota/dakota_out_parse.m (modified) (2 diffs)

issm/trunk/src/m/solutions/dakota/dakota_out_parse.m

@@ -71 +71 @@
     case {'nond_local_reliability'}
         [dresp]=nond_locrel(fidi);
-%     case {'dace','fsu_quasi_mc','fsu_cvt'}
-%     case {'vector_parameter_study','list_parameter_study',...
-%             'centered parameter_study','multidim_parameter_study'}
+     case {'dace','fsu_quasi_mc','fsu_cvt'}
+        [dresp]=param_any(fidi);
+     case {'vector_parameter_study','list_parameter_study',...
+             'centered_parameter_study','multidim_parameter_study'}
+        [dresp]=param_any(fidi);
     otherwise
         error('Unrecognized  method ''%s'' in file %s.',...
@@ -735 +737 @@
 end
 
+%%  function to parse any parameter study output file
+
+function [dresp]=param_any(fidi)
+
+%%  search through the file to find the Dakota output data
+
+[fline]=findline(fidi,'<<<<< Iterator');
+if ~ischar(fline)
+    return
+end
+display(['  ' deblank(fline)]);
+
+%  loop through the file to find the number of function evaluations
+
+[nfeval]=nfeval_read(fidi);
+
+%  loop through the file to find the volumetric uniformity measures
+
+[dresp]=vum_read(fidi,struct([]));
+
+end
+
+%%  function to find and read the volumetric uniformity measures
+
+function [dresp]=vum_read(fidi,dresp)
+
+[fline]=findline(fidi,'The following lists volumetric uniformity measures');
+if ~ischar(fline)
+    return
+end
+dresp(end+1).vum=[];
+
+display('Reading measures for volumetric uniformity.');
+
+fline=fgetl(fidi);
+fline=fgetl(fidi);
+
+while ischar(fline) && ~isempty(fline)
+        [ntokens,tokens]=fltokens(fline);
+    switch lower(tokens{1}{1})
+        case 'chi'
+            dresp.vum.chi=tokens{1}{4};
+        case 'd'
+            dresp.vum.d  =tokens{1}{4};
+        case 'h'
+            dresp.vum.h  =tokens{1}{4};
+        case 'tau'
+            dresp.vum.tau=tokens{1}{4};
+    end
+    fline=fgetl(fidi);
+end
+
+end
+
 %%  function to find a file line starting with a specified string