DataStore.hpp

// Copyright 2026 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_OMD_DATA_STORE__
#define __INC_OMD_DATA_STORE__
 
#include <Eigen/Core>
 
namespace OMD {
 
class DataStore
{
public:
    const Eigen::Index m;
    const Eigen::Index r;
    const Eigen::Index n;
protected:
    const Eigen::Ref<const Eigen::MatrixXd> S;
    Eigen::MatrixXd LTS;
    const double normY;
    Eigen::MatrixXd W;
    Eigen::MatrixXd XTL_qr_rep;
    Eigen::MatrixXd YTL_qr_rep_top;
    Eigen::MatrixXd YTL_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, Eigen::Index r);
    void updateSystemMatrices(const Eigen::Ref<const Eigen::MatrixXd> L);
    const Eigen::MatrixXd getSystemMatrix() const;
    const double getResidual() const;
    inline const Eigen::MatrixXd& getProjector() const {return W;};
    inline const Eigen::MatrixXd& getYstar() const {return YTL_qr_rep_bottom;};
};
 
};
 
#endif