PODReducer.m

Go to the documentation of this file.

namespace spacereduction{


/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

class PODReducer
  :public spacereduction.BaseSpaceReducer,
   public general.POD,
   public IReductionSummaryPlotProvider {
    public:

        rowvec<double> SingularValues = "[]";
    public:

        .false IncludeInitialSpace = false;
    public:

         function  plotSummary(pm,context) {
            plotSummary@spacereduction.BaseSpaceReducer(this, pm, context);
            if ~isempty(this.SingularValues)
                str = sprintf(" %s: POD singular value decay ",context);
                h = pm.nextPlot(" podreducer_singvals ",str,...
                    " subspace size "," singular values ");
                semilogy(h,this.SingularValues," LineWidth ",2);
/*                  h = pm.nextPlot('podreducer_projerr',str,...
 *                     'subspace size','projection error');
 *                 semilogy(h,this.SingularValues,'LineWidth',2); */
            else
                warning(" spacereduction:PODReducer ",...
                    " Singular value data empty. Not providing summary. ");
            end
        end
    end

    methods(Access=protected)        
        function [V,W] = generateReducedSpaceImpl(this, model, subset)
            /*  Implements the abstract method from BaseSpaceReducer */

            /*  Compile all trajectories into a large vector */
            md = model.Data;
            td = md.TrajectoryData;
            nt = td.getNumTrajectories;
            if nt == 0
                error(" No training data found in ModelData. ");
            end

            /*  Augment block data with fxi values */
            if this.IncludeTrajectoryFxiData
                if isempty(md.TrajectoryFxiData)
                    error(" No training fxi data found in ModelData. ");
                end
                td = data.JoinedBlockData(td, md.TrajectoryData);
            end
            /*  Wrap in finite difference adder */
            if this.IncludeFiniteDifferences
                td = data.FinDiffBlockData(td);
            end

            Vex = [];
            if this.IncludeBSpan
                Vex = md.InputSpaceSpan(subset,:);
            end
            if this.IncludeInitialSpace
                is = this.getInitialSpace(md.TrajectoryData, this, subset);
                Vex = [Vex is];
            end

            [V, this.SingularValues] = this.computePOD(td, Vex, subset);
            this.ProjectionError= flipud(cumsum(flipud(this.SingularValues)));
/*              if ~isempty(Vex)
 *                 o = general.Orthonormalizer;
 *                 V = o.orthonormalize([Vex V]);
 *             end */
            if ~isempty(Vex)
                o = general.Orthonormalizer;
                V = [Vex V];
                id = speye(size(V,2));
                while true  /*  V isn't really orthogonal after orthonormalizing once */

                    diff = max(max(abs(V^t*V-id)));
                    if diff < 1e-9
                        break
                    else
                        V = o.orthonormalize(V);
                    end
                end
            end            
            /*  Galerkin projection */
            W = [];
        end
    end

end


         }

};