source: src/LinearAlgebra/Subspace.hpp@ 1fb318

/*
 * Subspace.hpp
 *
 *  Created on: Nov 22, 2010
 *      Author: heber
 */

#ifndef SUBSPACE_HPP_
#define SUBSPACE_HPP_

#include <map>
#include <set>
#include <vector>
#include "Eigenspace.hpp"
#include "MatrixContent.hpp"
#include "VectorContent.hpp"
#include "unittests/SubspaceFactorizerUnittest.hpp"

/** A subset of eigenvectors from an Eigenspace.
 *
 * In this class we regard a sub set of eigenvectors of an Eigenspace
 * which span a subspace of the eigenspace. This is used for diagonalization
 * of the Eigenspace's matrix in linear-scaling, subspace decomposition
 * schemes.
 *
 * Here, beyond the contents of Eigenspace, we need projection matrices from
 * and to this subspace and also mappings from the global indices to the local
 * indices, to identify local eigenvectors in this Subspace with their
 * counterparts in the full Eigenspace.
 *
 */
class Subspace : public Eigenspace
{
  // TODO: Remove if not needed anymore
  friend void SubspaceFactorizerUnittest::SubspaceTest();
public:
  typedef std::map<size_t, size_t> mapping;
  typedef std::set< boost::shared_ptr<Subspace> > subset;

  Subspace(indexset &_s, Eigenspace &_FullSpace);
  ~Subspace();

  // manipulate subsets
  bool addSubset(boost::shared_ptr<Subspace> &_s);
  bool removeSubset(boost::shared_ptr<Subspace> &_s);

  void calculateEigenSubspace();
  void correctEigenvectorsFromSubIndices();
  const MatrixContent & getEigenvectorMatrixInFullSpace();
  const eigenvectorset & getEigenvectorsInFullSpace();
  const VectorContent getEigenvectorParallelToFullOne(size_t i);
  const double getEigenvalueOfEigenvectorParallelToFullOne(size_t i);
  void createLocalMapping();

private:

  void invertLocalToGlobalMapping();
  void getSubspacematrixFromBigmatrix(const MatrixContent & bigmatrix);
  void createProjectionMatrices();
  void projectFullSpaceMatrixToSubspace();

  mapping LocalToGlobal;
  mapping GlobalToLocal;
  subset SubIndices;
  MatrixContent ProjectToSubspace;
  MatrixContent ProjectFromSubspace;
  Eigenspace &FullSpace;

  const VectorContent ZeroVector;
};


#endif /* SUBSPACE_HPP_ */