Heed::PointCoorMesh< T, D > Class Template Reference

#include <tline.h>

Public Member Functions
long	get_qi (void) const
T	get_xmin (void) const
T	get_xmax (void) const
void	get_scoor (long n, T &b) const
virtual int	get_interval (long n, T &b1, T &b2) const
int	get_interval (T x, long &n1) const
int	get_interval (T x, long &n1, T &b1, long &n2, T &b2) const
int	get_interval_extrap (T x, long &n1, T &b1, long &n2, T &b2) const
int	get_step (long n, T &fstep) const
	PointCoorMesh (void)
	PointCoorMesh (long fq, D *famesh)
virtual	~PointCoorMesh ()
void	check (void)
virtual void	print (std::ostream &file) const

Detailed Description

template<class T, class D>
class Heed::PointCoorMesh< T, D >

Generic mesh with arbitrary steps. The array determining the step edges is located somewhere outside. In object of this class only the raw pointer is contained with consequences:

Definition at line 389 of file tline.h.

Constructor & Destructor Documentation

PointCoorMesh() [1/2]

template<class T, class D>

Heed::PointCoorMesh< T, D >::PointCoorMesh ( void )

inline

Definition at line 427 of file tline.h.

      : q(0), xmin(0), xmax(0), x_old(0), n_old(-1), amesh(NULL) {
    ;
  }

PointCoorMesh() [2/2]

template<class T, class D>

Heed::PointCoorMesh< T, D >::PointCoorMesh	(	long	fq,
		D *	famesh )

Definition at line 459 of file tline.h.

    : q(fq), x_old(0), n_old(-1) {
  if (q <= 1) {
    mcerr << "ERROR in PointCoorMesh<T,D>::PointCoorMesh<T,D>:\n"
          << "q <= 1\n";
    Iprintn(mcerr, q);
    spexit(mcerr);
  }
#ifndef TLINE_REDUCE_TO_RAW_ARR
  // amesh.put( famesh );
  amesh = famesh;
  xmin = (*amesh)[0];
  xmax = (*amesh)[q - 1];
#else
  amesh = &((*famesh)[0]);
  xmin = amesh[0];
  xmax = amesh[q - 1];
#endif
  // check consistence
  if (xmin > xmax) {
    mcerr << "ERROR in PointCoorMesh<T,D>::PointCoorMesh<T,D>:\n"
          << "xmin > xmax\n";
    Iprint2n(mcerr, xmin, xmax);
    spexit(mcerr);
  }
#ifdef CHECK_POINT_MESH
  long n;
  for (n = 0; n < q - 1; n++) {
#ifndef TLINE_REDUCE_TO_RAW_ARR
    if ((*amesh)[n] >= (*amesh)[n + 1])
#else
    if (amesh[n] >= amesh[n + 1])
#endif
    {
      mcerr << "ERROR in PointCoorMesh<T,D>::PointCoorMesh<T,D>:\n"
            << "amesh[n] >= amesh[n+1]\n";
#ifndef TLINE_REDUCE_TO_RAW_ARR
      Iprint3n(mcerr, n, (*amesh)[n], (*amesh)[n + 1]);
#else
      Iprint3n(mcerr, n, amesh[n], amesh[n + 1]);
#endif
      spexit(mcerr);
    }
  }
#endif
}

~PointCoorMesh()

template<class T, class D>

virtual Heed::PointCoorMesh< T, D >::~PointCoorMesh ( )

inlinevirtual

Definition at line 435 of file tline.h.

{}

Member Function Documentation

check()

template<class T, class D>

void Heed::PointCoorMesh< T, D >::check

(

void

)

Definition at line 507 of file tline.h.

                                    {
  long n;
  for (n = 0; n < q - 1; n++) {
#ifndef TLINE_REDUCE_TO_RAW_ARR
    if ((*amesh)[n] >= (*amesh)[n + 1])
#else
    if (amesh[n] >= amesh[n + 1])
#endif
    {
      mcerr << "ERROR in PointCoorMesh<T,D>::check(void):\n"
            << "amesh[n] >= amesh[n+1]\n";
#ifndef TLINE_REDUCE_TO_RAW_ARR
      Iprint3n(mcerr, n, (*amesh)[n], (*amesh)[n + 1]);
#else
      Iprint3n(mcerr, n, amesh[n], amesh[n + 1]);
#endif
      spexit(mcerr);
    }
  }
}

get_interval() [1/3]

template<class T, class D>

virtual int Heed::PointCoorMesh< T, D >::get_interval ( long n, T & b1, T & b2 ) const

inlinevirtual

Definition at line 401 of file tline.h.

                                                       {
    if (n < 0 || n >= q - 1) return 0;
#ifndef TLINE_REDUCE_TO_RAW_ARR
    b1 = (*amesh)[n];
    b2 = (*amesh)[n + 1];
#else
    b1 = amesh[n];
    b2 = amesh[n + 1];
#endif
    return 1;
  }

Referenced by get_step(), and Heed::HeedParticle_BGM::physics().

get_interval() [2/3]

template<class T, class D>

int Heed::PointCoorMesh< T, D >::get_interval	(	T	x,
		long &	n1 ) const

Definition at line 529 of file tline.h.

                                                         {
  if (x < xmin || x >= xmax) {
    n1 = 0;
    return 0;
  }
#ifndef TLINE_REDUCE_TO_RAW_ARR
  if (n_old >= 0 && x_old <= x) {
    n1 = t_find_interval_end<T, D>(x, q, *amesh, n_old);
  } else {
    n1 = t_find_interval<T, D>(x, q, *amesh);
  }
// n1 = t_find_interval< T , D >(x, q, amesh);
#else
  if (n_old >= 0 && x_old <= x) {
    n1 = t_find_interval_end<T, T*>(x, q, amesh, n_old);
  } else {
    n1 = t_find_interval<T, T*>(x, q, amesh);
  }
// n1 = t_find_interval< T , T* >(x, q, &amesh);
#endif
 
  if (n1 < 0) {
    mcerr << "ERROR in PointCoorMesh<T,D>::get_interval:\n"
          << "n1 < 0\n";
    print(mcerr);
    Iprintn(mcerr, n1);
    spexit(mcerr);
  }
  n_old = n1;
  x_old = x;
  return 1;
}

get_interval() [3/3]

template<class T, class D>

int Heed::PointCoorMesh< T, D >::get_interval	(	T	x,
		long &	n1,
		T &	b1,
		long &	n2,
		T &	b2 ) const

Definition at line 563 of file tline.h.

                                                   {
  if (x < xmin || x >= xmax) {
    n1 = 0;
    n2 = 0;
    b1 = 0;
    b2 = 0;
    return 0;
  }
#ifndef TLINE_REDUCE_TO_RAW_ARR
  if (n_old >= 0 && x_old <= x) {
    n1 = t_find_interval_end<T, D>(x, q, *amesh, n_old);
  } else {
    n1 = t_find_interval<T, D>(x, q, *amesh);
  }
// n1 = t_find_interval< T , D >(x, q, amesh);
#else
  if (n_old >= 0 && x_old <= x) {
    n1 = t_find_interval_end<T, T*>(x, q, amesh, n_old);
  } else {
    n1 = t_find_interval<T, T*>(x, q, amesh);
  }
// n1 = t_find_interval< T , T* >(x, q, &amesh);
#endif
  n2 = n1 + 1;
  if (n1 < 0 || n1 >= q || n2 < 0 || n2 >= q) {
    mcerr << "ERROR in PointCoorMesh<T,D>::get_interval:\n"
          << "n1 < 0 || n1 >= q || n2 < 0 || n2 >= q\n";
    print(mcerr);
    Iprint2n(mcerr, n1, n2);
    spexit(mcerr);
  }
#ifndef TLINE_REDUCE_TO_RAW_ARR
  b1 = (*amesh)[n1];
  b2 = (*amesh)[n2];
#else
  b1 = amesh[n1];
  b2 = amesh[n2];
#endif
  if (b1 < xmin || b1 > xmax || b2 < xmin || b2 > xmax) {
    mcerr << "ERROR in PointCoorMesh<T,D>::get_interval:\n"
          << "b1 < xmin || b1 > xmax || b2 < xmin || b2 > xmax\n";
    print(mcerr);
    Iprint2n(mcerr, b1, b2);
    spexit(mcerr);
  }
  n_old = n1;
  x_old = x;
  return 1;
}

get_interval_extrap()

template<class T, class D>

int Heed::PointCoorMesh< T, D >::get_interval_extrap	(	T	x,
		long &	n1,
		T &	b1,
		long &	n2,
		T &	b2 ) const

Definition at line 615 of file tline.h.

                                                          {
  int i_ret = 1;
 
  if (x < xmin) {
    i_ret = 0;
    n1 = 0;
    n2 = 1;
    b1 = xmin;
#ifndef TLINE_REDUCE_TO_RAW_ARR
    b2 = (*amesh)[1];
#else
    b2 = amesh[1];
#endif
  } else if (x >= xmax) {
    i_ret = 2;
    n1 = q - 2;
    n2 = q - 1;
#ifndef TLINE_REDUCE_TO_RAW_ARR
    b1 = (*amesh)[q - 2];
#else
    b1 = amesh[q - 2];
#endif
    b2 = xmax;
  } else {
#ifndef TLINE_REDUCE_TO_RAW_ARR
    if (n_old >= 0 && x_old <= x) {
      n1 = t_find_interval_end<T, D>(x, q, *amesh, n_old);
    } else {
      n1 = t_find_interval<T, D>(x, q, *amesh);
    }
// n1 = t_find_interval< T , D >(x, q, amesh);
#else
    if (n_old >= 0 && x_old <= x) {
      n1 = t_find_interval_end<T, T*>(x, q, amesh, n_old);
    } else {
      n1 = t_find_interval<T, T*>(x, q, amesh);
    }
// n1 = t_find_interval< T , T* >(x, q, &amesh);
#endif
    n2 = n1 + 1;
    if (n1 < 0 || n1 >= q || n2 < 0 || n2 >= q) {
      mcerr << "ERROR in PointCoorMesh<T,D>::get_interval:\n"
            << "n1 < 0 || n1 >= q || n2 < 0 || n2 >= q\n";
      print(mcerr);
      Iprint2n(mcerr, n1, n2);
      spexit(mcerr);
    }
#ifndef TLINE_REDUCE_TO_RAW_ARR
    b1 = (*amesh)[n1];
    b2 = (*amesh)[n2];
#else
    b1 = amesh[n1];
    b2 = amesh[n2];
#endif
    if (b1 < xmin || b1 > xmax || b2 < xmin || b2 > xmax) {
      mcerr << "ERROR in PointCoorMesh<T,D>::get_interval:\n"
            << "b1 < xmin || b1 > xmax || b2 < xmin || b2 > xmax\n";
      print(mcerr);
      Iprint2n(mcerr, b1, b2);
      spexit(mcerr);
    }
    n_old = n1;
    x_old = x;
  }
  return i_ret;
}

get_qi()

template<class T, class D>

long Heed::PointCoorMesh< T, D >::get_qi ( void ) const

inline

Definition at line 391 of file tline.h.

{ return q - 1; }

get_scoor()

template<class T, class D>

void Heed::PointCoorMesh< T, D >::get_scoor ( long n, T & b ) const

inline

Definition at line 394 of file tline.h.

                                            {
#ifndef TLINE_REDUCE_TO_RAW_ARR
    b = (*amesh)[n];
#else
    b = amesh[n];
#endif
  }

get_step()

template<class T, class D>

int Heed::PointCoorMesh< T, D >::get_step ( long n, T & fstep ) const

inline

Definition at line 419 of file tline.h.

                                              {
    if (n < 0 || n >= q - 1) return 0;
    T b1, b2;
    get_interval(n, b1, b2);
    fstep = b2 - b1;
    return 1;
  }

get_xmax()

template<class T, class D>

T Heed::PointCoorMesh< T, D >::get_xmax ( void ) const

inline

Definition at line 393 of file tline.h.

{ return xmax; }

get_xmin()

template<class T, class D>

T Heed::PointCoorMesh< T, D >::get_xmin ( void ) const

inline

Definition at line 392 of file tline.h.

{ return xmin; }

print()

template<class T, class D>

void Heed::PointCoorMesh< T, D >::print ( std::ostream & file ) const

virtual

Definition at line 684 of file tline.h.

                                                      {
  Ifile << "PointCoorMesh<T,D>:\n";
  indn.n += 2;
  Ifile << "Type of T is (in internal notations) " << typeid(T).name() << '\n';
  Ifile << "Type of D is (in internal notations) " << typeid(D).name() << '\n';
  Iprint3n(file, q, xmin, xmax);
  Iprint2n(file, n_old, x_old);
#ifndef TLINE_REDUCE_TO_RAW_ARR
  // Ifile << "(*amesh)=" << (*amesh) << '\n';
  Ifile << "(*amesh)=" << (*amesh)[0] << '\n';
#else
  Ifile << "amesh:" << '\n';
  long n;
  indn.n += 2;
  for (n = 0; n < q; n++) {
    Ifile << "n=" << n << " amesh[n]=" << noindent << amesh[n] << yesindent
          << '\n';
  }
  file << yesindent;
  indn.n -= 2;
#endif
  indn.n -= 2;
}

Referenced by get_interval(), get_interval(), get_interval_extrap(), and Heed::operator<<().

---

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

• tline.h