DataStore.hpp

// Copyright Lucas Mieg, Martin Moennigmann
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
 
#ifndef __INC_OMDC_DATA_STORE__
#define __INC_OMDC_DATA_STORE__
 
#include <Eigen/Core>
 
namespace OMDc {
 
class DataStore
{
public:
    const Eigen::Index m;
    const Eigen::Index r;
    const Eigen::Index n;
    const Eigen::Index p;
private:
    const Eigen::Ref<const Eigen::MatrixXd> S;
    const Eigen::Ref<const Eigen::MatrixXd> U;
    Eigen::MatrixXd LTS;
    const Eigen::MatrixXd Z;
    Eigen::MatrixXd ZW;
    const double normY;
    Eigen::MatrixXd ZTXTL_qr_rep;
    Eigen::MatrixXd ZTYTL_qr_rep_top;
    Eigen::MatrixXd ZTYTL_qr_rep_bottom;
public:
    typedef Eigen::Block<const Eigen::MatrixXd,Eigen::Dynamic,Eigen::Dynamic,true> Block;
    typedef Eigen::Block<const Eigen::Ref<const Eigen::MatrixXd>,Eigen::Dynamic,Eigen::Dynamic,true> RefBlock;
    inline const RefBlock X() const {
        return S.leftCols(m-1);
    };
    inline const RefBlock Y() const {
        return S.rightCols(m-1);
    };
    inline const Block LTX() const {
        return LTS.leftCols(m-1);
    };
    inline const Block LTY() const {
        return LTS.rightCols(m-1);
    };
    DataStore( const Eigen::Ref<const Eigen::MatrixXd> S, const Eigen::Ref<const Eigen::MatrixXd> U, Eigen::Index r);
    void updateSystemMatrices(const Eigen::Ref<const Eigen::MatrixXd> L);
    const Eigen::MatrixXd getSystemMatrix() const;
    const Eigen::MatrixXd& getYstar() const;
    const Eigen::MatrixXd &getProjector() const;
    const Eigen::MatrixXd getInputMatrix() const;
    const double getResidual() const;
    const Eigen::VectorXd getInitialState() const;
};
 
};
 
#endif