readSample.m

Go to the documentation of this file.

%>  \brief
%>  Return a list of objects of superclass [pm.sampling.FileContentsSample](@ref FileContentsSample)
%>  containing the contents of a (set of) ParaMonte simulation output sample file(s) whose paths match
%>  the specified input ``pattern`` or the simulation specification ``sampler.spec.outputFileName``.<br>
%>
%>  \warning
%>  This function is to be only used for post-processing of the output sample file(s) of an already finished simulation.<br>
%>  Although possible, this method is **not** meant to be called by all processes in MPI-parallel simulations.<br>
%>
%>  \param[in]  sampler :   The input object of superclass [pm.sampling.Sampler](@ref Sampler)
%>                          whose type and properties determine the type of the output object(s).<br>
%>                          (**optional**. If empty, it is set to [pm.sampling.Sampler()](@ref Sampler).)
%>  \param[in]  pattern :   The input scalar MATLAB string containing the pattern matching
%>                          the desired sample file(s) whose contents is to be read.<br>
%>                          The specified ``pattern`` only needs to partially identify
%>                          the name of the simulation to which the sample file belongs.<br>
%>                          For example, specifying ``"./mydir/mysim"`` as input will
%>                          lead to a search for file(s) beginning with ``"mysim"`` and
%>                          ending with ``"_sample.txt"`` inside the directory ``"./mydir/"``.<br>
%>                          If there are multiple files matching in the input ``pattern``,
%>                          then all such files will be read and returned as elements of a list.<br>
%>                          If the specified pattern is a valid existing URL, the file will be
%>                          downloaded as a temporary file to the local system, its contents
%>                          shall be parsed and the file will be subsequently removed.<br>
%>                          If the input ``pattern`` is empty, then the method will search
%>                          for any possible candidate files with the appropriate suffix
%>                          in the current working directory.<br>
%>                          (**optional**. If empty, it is set to ``sampler.spec.outputFileName`` or if empty, it is set to ``"./"``.)
%>  \param[in]  sep     :   The input MATLAB string containing the field separator used in the file(s).<br>
%>                          (**optional**. If empty, it is set to ``sampler.spec.outputSeparator`` or if empty, it is automatically inferred.)
%>
%>  \return
%>  ``sampleList``      :   The output MATLAB cell array of objects of superclass [pm.sampling.FileContentsSample](@ref FileContentsSample),
%>                          each of which corresponds to the contents of a unique ParaMonte sampler sample file.<br>
%>                          <ol>
%>                              <li>    If the input argument ``sampler`` is of type [pm.sampling.Sampler](@ref Sampler),
%>                                      then the output array elements are of type [pm.sampling.FileContentsSample](@ref FileContentsSample).<br>
%>                              <li>    If the input argument ``sampler`` is of type [pm.sampling.Paradram](@ref Paradram),
%>                                      then the output array elements are of type [pm.sampling.FileContentsSample](@ref FileContentsSample).<br>
%>                          </ol>
%>
%>  \interface{readSample}
%>  \code{.m}
%>
%>      sampleList = pm.sampling.readSample();
%>      sampleList = pm.sampling.readSample([]);
%>      sampleList = pm.sampling.readSample([], []);
%>      sampleList = pm.sampling.readSample([], [], []);
%>      sampleList = pm.sampling.readSample([], pattern);
%>      sampleList = pm.sampling.readSample([], [], sep);
%>      sampleList = pm.sampling.readSample([], pattern, sep);
%>      sampleList = pm.sampling.readSample(sampler, pattern, sep);
%>
%>  \endcode
%>
%>  \note
%>  See the documentation of the subclasses of [pm.sampling.Sampler](@ref Sampler)
%>  (e.g., [pm.sampling.Paradram](@ref Paradram)) for example usage in action.<br>
%>
%>  \example{readSample}
%>  \code{.m}
%>
%>      sampler = pm.sampling.Sampler();
%>
%>      sampleList = pm.sampling.readSample([], "./out/test_run_");
%>      sampleList = pm.sampling.readSample(sampler, "./out/test_run_");
%>
%>      sampler.spec.outputFileName = "./out/test_run_";
%>
%>      sampleList = pm.sampling.readSample(sampler, []);
%>      sampleList = pm.sampling.readSample(sampler);
%>
%>      sampler.spec.outputSeparator = ",";
%>      pm.sampling.readSample(sampler, "./out/test_run_");
%>      pm.sampling.readSample(sampler, "./out/test_run_", []);
%>      sampleList = pm.sampling.readSample(sampler, "./out/test_run_", ",");
%>
%>  \endcode
%>
%>  \final{readSample}
%>
%>  \author
%>  \AmirShahmoradi, Saturday Nov 2 2024, 11:11 PM, Dallas, TX.<br>
%>  \FatemehBagheri, May 20 2024, 1:25 PM, NASA Goddard Space Flight Center (GSFC), Washington, D.C.<br>
function sampleList = readSample(sampler, pattern, sep)
 
    if  nargin < 1
        sampler = [];
    end
    if  isempty(sampler)
        sampler = pm.sampling.Sampler();
    end
 
    if  nargin < 2
        if  0 < pm.array.len(sampler.spec.outputFileName)
            pattern = string(sampler.spec.outputFileName);
        else
            pattern = [];
        end
    end
 
    if  nargin < 3
        if  0 < pm.array.len(sampler.spec.outputSeparator)
            sep = string(sampler.spec.outputSeparator);
        else
            sep = [];
        end
    end
 
    pathList = sampler.findfile("sample", pattern);
    sampleList = cell(numel(pathList), 1);
 
    for ifile = numel(pathList) : -1 : 1
 
        if ~sampler.silent
            disp("processing file: """ + pathList(ifile) + """");
        end
 
        if  isequal(class(sampler), "pm.sampling.Sampler")
            sampleList{ifile} = pm.sampling.FileContentsSample(pathList(ifile), sampler.silent, sep);
        elseif isequal(class(sampler), "pm.sampling.Paradram")
            sampleList{ifile} = pm.sampling.FileContentsSampleDRAM(pathList(ifile), sampler.silent, sep);
        else
            help("pm.sampling.readSample")
            disp("class(sampler)");
            disp( class(sampler) );
            error   ( newline ...
                    + "Invalid input ``sampler`` object type." + newline ...
                    + "See the documentation displayed above for more " + newline ...
                    + "information on possible values for ``sampler`` argument." + newline ...
                    + newline ...
                    );
        end
 
    end
 
end