SphericalHarmonics(unsigned n1, unsigned n2)
 {
     m_n1 = n1;
     m_n2 = n2;
     m_Cell.setBasisVectors(1.0,1.0,1.0,90.0,90.0,90.0);
     m_Energies = Energies(0.0,0.0,1);
 };
void PlotVlasovPoisson(void *vf, real *w) {
  real k_energy = 0, e_energy = 0, t_energy = 0;
  
  field *f = vf;
  
  Energies(f, w, k_energy, e_energy, t_energy,1);
  vf = f;
}
 ConstantFunction()
 {
     m_Cell.setBasisVectors(1.0,1.0,1.0,90.0,90.0,90.0);
     m_Energies = Energies(0.0,0.0,1);
 };
void TwoDimensionalCut::setEnergyPoints(double min, double max, std::size_t points)
{
  this->cut->setEnergies(Energies(min, max, points));
}