app-bholger/gfem/gfemBasisOne.h

Go to the documentation of this file.

#pragma once

#include <string>

#include "basics.hh"
#include "geometry.hh"
#include "space.hh"
#include "graphics.hh"
#include "toolbox.hh"
#include "function.hh"
#include "operator.hh"
#include "linearFEM.hh"
#include "linDG1D.hh"
#include "hp2D.hh"
#include "gfemIdentity.h"
#include "gfemLaplace.h"
#include "hp2D/formula.hh"
#include "PC_UC_Solver.h"

namespace concepts {
namespace gfem {

class GfemBasisOne : public PC_UC_Solver {
  public:

  GfemBasisOne(const concepts::Space<concepts::Real>& space, 
          EquationType type, 
          const MID& attrToEpsilon, 
          double omega,
          const concepts::BoundaryConditions& bc) 
    : PC_UC_Solver(space, type, attrToEpsilon, omega, bc)
  {

  }
 
  ~GfemBasisOne() {
    for(int i=0; i < (int)uc_sols.size(); ++i)
      delete uc_sols[i];
    uc_sols.clear();
  }

  void releaseContents() {
    uc_sols.clear();
  }

  void solve(concepts::Cmplx k1, concepts::Cmplx k2) {
    UC_Sol* ucSol = new UC_Sol(space, k1, k2);
    PC_UC_Solver::solve(k1, k2, ucSol->vec);

    uc_sols.push_back(ucSol);
  }
 
  void addSolution(UC_Sol* sol) {
    uc_sols.push_back(sol);
  }

  void orthogonalize();

  const UC_Sol& getSolution(int i) const {
    return *(uc_sols[i]);
  }

  UC_Sol& getSolution(int i) {
    return *(uc_sols[i]);
  }

  int size() const {
    return uc_sols.size();
  } 

  void cutoff(double sigma_thresh);

  void printSigmas() const;

  void setAllKs(Cmplx k1, Cmplx k2);

private:
  void transformReal();
  void sort();

  std::vector< UC_Sol* > uc_sols;
};

} // namespace
} // namespace