app-bholger/gfem/karniadakis.h

Go to the documentation of this file.

#pragma once
#include "basics/typedefs.hh"
#include "polyCoeff.h"
#include <memory>
#include "stdio.h"
#include "assert.h"
#include "algorithm"
#include <vector>

namespace concepts {

class KarniadakisCoeffs {
private:
  KarniadakisCoeffs();
  ~KarniadakisCoeffs();

  KarniadakisCoeffs& operator=(const KarniadakisCoeffs& other);

public:
  KarniadakisCoeffs(const KarniadakisCoeffs& other);
  KarniadakisCoeffs(const KarniadakisCoeffs& other, int nper);

  static KarniadakisCoeffs& getInstance() {
    return g_instance; 
  }

  void reinit(int n);

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

  const PolyCoeff& getCoeffs(int i) const {
    return *kcoeffs[i];
  }

  const PolyCoeff& getCoeffsD(int i) const {
    return *kcoeffsD[i];
  }

  const PolyCoeff& getCoeffs(int i, int i_diff) const {
    if(i_diff == 1)
      return getCoeffs(i);
    return getCoeffsD(i);
  }

  static const int SELMAC = 0;

  const PolyCoeff& getCoeffs(int i, int i_diff, int i_comp, int i_sel) const {
    assert(i < size()); 
    if(i_diff == SELMAC && i_comp == i_sel) 
      return getCoeffsD(i);
    return getCoeffs(i);

  }

  const PolyCoeff& getJacobiCoeffs(int i) const {
    return *pJacobi[i];
  }


private:
  static KarniadakisCoeffs g_instance;
  std::vector<PolyCoeff*> kcoeffs;
  std::vector<PolyCoeff*> kcoeffsD;
  std::vector<PolyCoeff*> pJacobi;
};

} // namespace