Lineq.h

Go to the documentation of this file.

/* 
 * @BEGIN LICENSE
 *
 * Copyright (C) 2014-2015  Ward Poelmans
 *
 * This file is part of v2DM-DOCI.
 * 
 * v2DM-DOCI is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * Foobar is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with Foobar.  If not, see <http://www.gnu.org/licenses/>.
 *
 * @END LICENSE
 */

#ifndef LINEQ_H
#define LINEQ_H

#include <iostream>
#include <vector>

#include "TPM.h"
#include "SUP.h"

namespace doci2DM
{

class Lineq
{
   friend std::ostream &operator<<(std::ostream &output,Lineq &lineq_p);

   public:

      Lineq(int L,int N, bool=false);

      virtual ~Lineq() = default;

      int gN() const;

      int gnr() const;

      int gL() const;

      const TPM &gE(int) const;

      double ge(int) const;

      const TPM &gE_ortho(int) const;

      double ge_ortho(int) const;

      const SUP &gu_0(int) const;

      const SUP &gu_0_ortho(int) const;

      void check(const TPM &tpm) const;

      void orthogonalize();

   private:

      void constr_u_0();

      void orthogonalize_u_0();

      std::vector<TPM> E;
      std::vector<double> e;
      
      std::vector<TPM> E_ortho;
      std::vector<double> e_ortho;

      std::vector<SUP> u_0;

      std::vector<SUP> u_0_ortho;

      int N;

      int L;
};

}

#endif /* LINEQ_H */

/* vim: set ts=3 sw=3 expandtab :*/