Heed::EnergyMesh Class Reference

#include <EnergyMesh.h>

Inheritance diagram for Heed::EnergyMesh:

Public Member Functions
	EnergyMesh (void)
	EnergyMesh (double femin, double femax, long fq)
	EnergyMesh (DynLinArr< double > fec)
long	get_q () const
double	get_emin () const
double	get_emax () const
double	get_e (long n) const
double	get_ec (long n) const
const double *	get_ae (void) const
const double *	get_aec (void) const
long	get_interval_number (double ener)
long	get_interval_number_between_centers (double ener)
virtual void	print (std::ostream &file, int l) const
	macro_copy_total (EnergyMesh)
virtual	~EnergyMesh ()

Friends
std::ostream &	operator<< (std::ostream &file, EnergyMesh &f)

Detailed Description

Definition at line 44 of file EnergyMesh.h.

Constructor & Destructor Documentation

EnergyMesh() [1/3]

Heed::EnergyMesh::EnergyMesh ( void  )

inline

Definition at line 46 of file EnergyMesh.h.

: q(0), emin(0.0), emax(0.0) {}

EnergyMesh() [2/3]

Heed::EnergyMesh::EnergyMesh	(	double	femin,
		double	femax,
		long	fq
	)

Definition at line 8 of file EnergyMesh.cpp.

    : q(fq), emin(femin), emax(femax) {
  mfunname("EnergyMesh::EnergyMesh(double femin, double femax, long fq)");
  check_econd21(q, < 0 ||, > pqener - 1, mcerr);
 
  double rk = pow(emax / emin, (1.0 / double(q)));
  double er = emin;
  e[0] = er;
  for (long n = 1; n < q + 1; n++) {
    e[n] = er * rk;
    ec[n - 1] = (e[n - 1] + e[n]) * 0.5;
    er = e[n];
  }
  //pcm_e = PointCoorMesh< double, double[pqener] >( q, &e);
  //pcm_ec = PointCoorMesh< double, double[pqener-1] >( q-1, &ec);
  //mcout<<"EnergyMesh is done"<<endl;
}

EnergyMesh() [3/3]

Heed::EnergyMesh::EnergyMesh

(

DynLinArr< double >

fec

)

Definition at line 26 of file EnergyMesh.cpp.

                                            : q(fec.get_qel()) {
  mfunname("DynLinArr< double > fec");
  check_econd21(q, < 0 ||, > pqener - 1, mcerr);
  check_econd11(q, != 1, mcerr);  // otherwise problems with emin/emax
  if (q > 0) {
    emin = fec[0] - (fec[1] - fec[0]) / 2.0;
    emax = fec[q - 1] + (fec[q - 1] - fec[q - 2]) / 2.0;
    e[0] = emin;
    e[q] = emax;
 
    for (long n = 0; n < q; n++) {
      ec[n] = fec[n];
    }
    for (long n = 1; n < q; n++) {
      e[n] = 0.5 * (fec[n - 1] + fec[n]);
    }
  } else {
    emin = 0.0;
    emax = 0.0;
  }
  //pcm_e = PointCoorMesh< double, double[pqener] >( q, &e);
  //pcm_ec = PointCoorMesh< double, double[pqener-1] >( q-1, &ec);
}

~EnergyMesh()

virtual Heed::EnergyMesh::~EnergyMesh ( )

inlinevirtual

Definition at line 75 of file EnergyMesh.h.

{}

Member Function Documentation

get_ae()

const double * Heed::EnergyMesh::get_ae ( void  ) const

inline

Definition at line 60 of file EnergyMesh.h.

                                          {
    return e;
  }  // left sides

get_aec()

const double * Heed::EnergyMesh::get_aec ( void  ) const

inline

Definition at line 63 of file EnergyMesh.h.

{ return ec; }

get_e()

double Heed::EnergyMesh::get_e ( long  n ) const

inline

Definition at line 54 of file EnergyMesh.h.

                                    {
    return e[n];
  }  // left side of interval

get_ec()

double Heed::EnergyMesh::get_ec ( long  n ) const

inline

Definition at line 57 of file EnergyMesh.h.

                                     {
    return ec[n];
  }  // center of interval

Referenced by Heed::HeedDeltaElectron::physics_after_new_speed(), and Heed::HeedDeltaElectron::physics_mrange().

get_emax()

double Heed::EnergyMesh::get_emax ( ) const

inline

Definition at line 53 of file EnergyMesh.h.

{ return emax; }

get_emin()

double Heed::EnergyMesh::get_emin ( ) const

inline

Definition at line 52 of file EnergyMesh.h.

{ return emin; }

get_interval_number()

long Heed::EnergyMesh::get_interval_number

(

double

ener

)

Definition at line 75 of file EnergyMesh.cpp.

                                                {
  if (ener < emin) return -1;
  if (ener > emax) return q;
 
  long n1 = 0;
  long n2 = q;  // right side of last
  long n3 = n1;
  while (n2 - n1 > 1) {
    n3 = n1 + (n2 - n1) / 2;
    if (ener < e[n3]) {
      n2 = n3;
    } else {
      n1 = n3;
    }
  }
  return n1;
}

get_interval_number_between_centers()

long Heed::EnergyMesh::get_interval_number_between_centers

(

double

ener

)

Definition at line 93 of file EnergyMesh.cpp.

                                                                {
  if (ener < ec[0]) return -1;
  if (ener > ec[q - 1]) return q;
 
  long n1 = 0;
  long n2 = q - 1;  // right side of last
  long n3 = n1;
  while (n2 - n1 > 1) {
    n3 = n1 + (n2 - n1) / 2;
    if (ener < ec[n3]) {
      n2 = n3;
    } else {
      n1 = n3;
    }
  }
  return n1;
}

Referenced by Heed::HeedDeltaElectron::physics_after_new_speed(), and Heed::HeedDeltaElectron::physics_mrange().

get_q()

long Heed::EnergyMesh::get_q ( ) const

inline

Definition at line 51 of file EnergyMesh.h.

{ return q; }

Referenced by Heed::HeedParticle::physics(), Heed::HeedParticle_BGM::physics(), Heed::HeedDeltaElectron::physics_after_new_speed(), and Heed::HeedDeltaElectron::physics_mrange().

macro_copy_total()

Heed::EnergyMesh::macro_copy_total

(

EnergyMesh

)

print()

void Heed::EnergyMesh::print ( std::ostream &  file, int  l  ) const

virtual

Definition at line 129 of file EnergyMesh.cpp.

                                                    {
  if (l <= 0) return;
  Ifile << "EnergyMesh (l=" << l << "): \n";
  indn.n += 2;
  Ifile << "emin=" << emin << " emax=" << emax << " quantity of intervals=" << q
        << '\n' << " maximal possible quantity of intervals=" << pqener << '\n';
  if (l > 1) {
    Ifile << " number  left side        center       right side       widht\n";
    for (int n = 0; n < q; n++) {
      Ifile << std::setw(5) << n << std::setw(15) << e[n] << std::setw(15)
            << ec[n] << std::setw(15) << e[n + 1] << std::setw(15)
            << (e[n + 1] - e[n]) << '\n';
    }
  }
  indn.n -= 2;
}

Friends And Related Function Documentation

operator<<

std::ostream & operator<< ( std::ostream &  file, EnergyMesh &  f  )

friend

Definition at line 111 of file EnergyMesh.cpp.

                                                        {
  Ifile << "EnergyMesh: \n";
  indn.n += 2;
  Ifile << "emin=" << f.emin << " emax=" << f.emax
        << " quantity of intervals=" << f.q << '\n'
        << " maximal possible quantity of intervals=" << pqener << '\n';
  Ifile << " number  left side        center       right side       widht\n";
  for (int n = 0; n < f.q; n++) {
    Ifile << std::setw(5) << n << std::setw(15) << f.e[n] << std::setw(15)
          << f.ec[n] << std::setw(15) << f.e[n + 1] << std::setw(15)
          << (f.e[n + 1] - f.e[n]) << '\n';
  }
  //f.pcm_e.print(mcout);
  //f.pcm_ec.print(mcout);
  indn.n -= 2;
  return file;
}

---

The documentation for this class was generated from the following files:

• EnergyMesh.h
• EnergyMesh.cpp