MatrixDEIM.m

Go to the documentation of this file.

namespace general{


/* (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 MatrixDEIM
  :public general.DEIM {
    public:

        NumRows;
    public:

        matrix<double> Qk = "[]";
    private:

        effNumRows;

        
    public:

        function M = evaluate(colvec<double> x,double t) {
            fx = evaluate@general.DEIM(this, x, t);
            M = reshape(fx, this.effNumRows, []);
        }


        function M = evaluateMulti(colvec<double> x,double t) {
            fx = evaluateMulti@general.DEIM(this, x, t);
            M = reshape(fx, this.effNumRows, []);
        }


        function  setSimilarityTransform(Qk) {
            this.Qk= Qk;
            U = this.U_nonproj;
            if isempty(Qk) || isequal(Qk,1)
                this.effNumRows= this.NumRows;
                unew = U;
            else
                this.effNumRows= size(Qk,2);
                unew = zeros(size(this.Qk,2)^2, size(U,2));
                for i=1:size(U,2)
                    umat = reshape(U(:,i),size(this.Qk,1),[]);
                    umat = this.Qk^t*(umat*this.Qk);
                    unew(:,i) = umat(:);
                end    
            end
            /*  Overwrite the current U with processed U_nonproj. */
            this.U= unew;
        }


        function target = project(V,unused1) {
            target = project@general.DEIM(this, V, [], this.clone);
        }


        function copy = clone() {
            copy = general.MatrixDEIM;
            copy = clone@general.DEIM(this, copy);
            copy.Qk= this.Qk;
            copy.NumRows= this.NumRows;
            copy.effNumRows= this.effNumRows;
        }

    
    protected:

        function  updateOrderData() {
            updateOrderData@general.DEIM(this);

            /*  "special" for jaccompevalwrappers:
             * they only evaluate to nonzero entries, which is why the DEIM approximation
             * U, U_nonproj is of size length(J(J ~= 0)) instead of length(J(:)) */
            if isa(this.f," general.JacCompEvalWrapper ") && ~isempty(this.f.f.JSparsityPattern)
                fsp = this.f.f.JSparsityPattern;
                this.U_nonproj= MatUtils.toSparse(this.U_nonproj, find(fsp), numel(fsp)); /* #ok<FNDSB> */

            end

            /*  Apply partial similarity transform (AFTER original DEIM order
             * update)
             * Checks for having Qk set are done there. */
            this.setSimilarityTransform(this.Qk);

            /*  For now: unset the auxiliary DEIM error estimation matrices
             * as they are not used (cannot be used?) for matrix DEIMs yet. */
            this.M1= [];
            this.M2= [];
/*              this.Uerr1 = [];
 *             this.Uerr2 = []; */
        }
};
}