SubspaceData.h

#ifndef LINEARALGEBRA_SRC_MOLPRO_LINALG_ITSOLV_SUBSPACE_SUBSPACEDATA_H
#define LINEARALGEBRA_SRC_MOLPRO_LINALG_ITSOLV_SUBSPACE_SUBSPACEDATA_H
#include <map>
#include <molpro/linalg/itsolv/subspace/Matrix.h>
#include <vector>
 
namespace molpro::linalg::itsolv::subspace {
enum class EqnData { H, S, rhs, value, M };
 
template <typename T = double>
using SubspaceData = std::map<EqnData, Matrix<T>>;
 
template <typename T>
std::vector<EqnData> square_blocks(const SubspaceData<T>& data) {
  auto blocks = std::vector<EqnData>{EqnData::H, EqnData::S};
  if (data.count(EqnData::M))
    blocks.push_back(EqnData::M);
  return blocks;
}
 
template <typename T = double, EqnData... DataTypes>
auto null_data() {
  return SubspaceData<T>{std::make_pair<EqnData, Matrix<T>>(DataTypes, {})...};
}
} // namespace molpro::linalg::itsolv::subspace
 
#endif // LINEARALGEBRA_SRC_MOLPRO_LINALG_ITSOLV_SUBSPACE_SUBSPACEDATA_H