CopiedPointCoorMesh< T > Class Template Reference

#include <tline.h>

Inheritance diagram for CopiedPointCoorMesh< T >:

Public Member Functions
long	get_qi (void) const
T	get_xmin (void) const
T	get_xmax (void) const
const DynLinArr< T > &	get_mesh (void) const
void	get_scoor (long n, T &b) const
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
	CopiedPointCoorMesh (void)
	CopiedPointCoorMesh (long fq, T *famesh)
	CopiedPointCoorMesh (DynLinArr< T > &famesh)
	CopiedPointCoorMesh (const CopiedPointCoorMesh< T > &f)
CopiedPointCoorMesh< T > &	operator= (const CopiedPointCoorMesh< T > &f)
void	print (std::ostream &file) const

Detailed Description

template<class T>
class CopiedPointCoorMesh< T >

Definition at line 760 of file tline.h.

Constructor & Destructor Documentation

CopiedPointCoorMesh() [1/4]

template<class T >

CopiedPointCoorMesh< T >::CopiedPointCoorMesh ( void  )

inline

Definition at line 795 of file tline.h.

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

CopiedPointCoorMesh() [2/4]

template<class T >

CopiedPointCoorMesh< T >::CopiedPointCoorMesh	(	long	fq,
		T *	famesh
	)

Definition at line 861 of file tline.h.

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

CopiedPointCoorMesh() [3/4]

template<class T >

CopiedPointCoorMesh< T >::CopiedPointCoorMesh

(

DynLinArr< T > &

famesh

)

Definition at line 905 of file tline.h.

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

CopiedPointCoorMesh() [4/4]

template<class T >

CopiedPointCoorMesh< T >::CopiedPointCoorMesh ( const CopiedPointCoorMesh< T > &  f )

inline

Definition at line 795 of file tline.h.

      : q(f.q),
        mesh(f.mesh),
        xmin(f.xmin),
        xmax(f.xmax),
        x_old(f.x_old),
        n_old(f.n_old)
#ifdef TLINE_REDUCE_TO_RAW_ARR
        amesh(&mesh)
#endif
        {
  }

Member Function Documentation

get_interval() [1/3]

template<class T >

int CopiedPointCoorMesh< T >::get_interval ( long  n, T &  b1, T &  b2  ) const

inline

Definition at line 775 of file tline.h.

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

Referenced by CopiedPointCoorMesh< T >::CopiedPointCoorMesh(), and VirtCopiedPointCoorMesh< T >::get_interval().

get_interval() [2/3]

template<class T >

int CopiedPointCoorMesh< T >::get_interval	(	T	x,
		long &	n1
	)		const

Definition at line 950 of file tline.h.

                                                            {
  if (x < xmin || x >= xmax) {
    n1 = 0;
    return 0;
  }
#ifdef TLINE_REDUCE_TO_RAW_ARR
  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 , D >(x, q, amesh);
#else
  if (n_old >= 0 && x_old <= x) {
    n1 = t_find_interval_end<T, DynLinArr<T> >(x, q, mesh, n_old);
  } else {
    n1 = t_find_interval<T, DynLinArr<T> >(x, q, mesh);
  }
//n1 = t_find_interval< T , T* >(x, q, &amesh);
#endif
 
  if (n1 < 0) {
    mcerr << "ERROR in CopiedPointCoorMesh<T>::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 >

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

Definition at line 984 of file tline.h.

                                                      {
  if (x < xmin || x >= xmax) {
    n1 = 0;
    n2 = 0;
    b1 = 0;
    b2 = 0;
    return 0;
  }
#ifdef TLINE_REDUCE_TO_RAW_ARR
  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 , D >(x, q, amesh);
#else
  if (n_old >= 0 && x_old <= x) {
    n1 = t_find_interval_end<T, DynLinArr<T> >(x, q, mesh, n_old);
  } else {
    n1 = t_find_interval<T, DynLinArr<T> >(x, q, mesh);
  }
//n1 = t_find_interval< T , T* >(x, q, &amesh);
#endif
  n2 = n1 + 1;
#ifdef TLINE_REDUCE_TO_RAW_ARR
  b1 = amesh[n1];
  b2 = amesh[n2];
#else
  b1 = mesh[n1];
  b2 = mesh[n2];
#endif
  if (n1 < 0 || n1 >= q || n2 < 0 || n2 >= q || b1 < xmin || b1 > xmax ||
      b2 < xmin || b2 > xmax) {
    mcerr << "ERROR in CopiedPointCoorMesh<T>::get_interval:\n"
          << "n1 < 0 || n1 >= q || n2 < 0 || n2 >= q || b1 < xmin || b1 > xmax "
             "|| b2 < xmin || b2 > xmax\n";
    print(mcerr);
    Iprint4n(mcerr, n1, n2, b1, b2);
    spexit(mcerr);
  }
  n_old = n1;
  x_old = x;
  return 1;
}

get_interval_extrap()

template<class T >

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

Definition at line 1031 of file tline.h.

                                                             {
  int i_ret = 1;
 
  if (x < xmin) {
    i_ret = 0;
    n1 = 0;
    n2 = 1;
    b1 = xmin;
#ifdef TLINE_REDUCE_TO_RAW_ARR
    b2 = (*amesh)[1];
#else
    b2 = mesh[1];
#endif
  } else if (x >= xmax) {
    i_ret = 2;
    n1 = q - 2;
    n2 = q - 1;
#ifdef TLINE_REDUCE_TO_RAW_ARR
    b1 = (*amesh)[q - 2];
#else
    b1 = mesh[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, DynLinArr<T> >(x, q, mesh, n_old);
    } else {
      n1 = t_find_interval<T, DynLinArr<T> >(x, q, mesh);
    }
//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;
#ifndef TLINE_REDUCE_TO_RAW_ARR
    b1 = mesh[n1];
    b2 = mesh[n2];
#else
    b1 = amesh[n1];
    b2 = amesh[n2];
#endif
    if (n1 < 0 || n1 >= q || n2 < 0 || n2 >= q || b1 < xmin || b1 > xmax ||
        b2 < xmin || b2 > xmax) {
      mcerr << "ERROR in CopiedPointCoorMesh<T>::get_interval:\n"
            << "n1 < 0 || n1 >= q || n2 < 0 || n2 >= q || b1 < xmin || b1 > "
               "xmax || b2 < xmin || b2 > xmax\n";
      print(mcerr);
      Iprint4n(mcerr, n1, n2, b1, b2);
      spexit(mcerr);
    }
    n_old = n1;
    x_old = x;
  }
  return i_ret;
 
}

Referenced by VirtCopiedPointCoorMesh< T >::get_interval_extrap().

get_mesh()

template<class T >

const DynLinArr< T > & CopiedPointCoorMesh< T >::get_mesh ( void  ) const

inline

Definition at line 766 of file tline.h.

{ return mesh; }

Referenced by apeq_mant(), operator!=(), and operator==().

get_qi()

template<class T >

long CopiedPointCoorMesh< T >::get_qi ( void  ) const

inline

Definition at line 763 of file tline.h.

{ return q - 1; }

Referenced by apeq_mant(), VirtCopiedPointCoorMesh< T >::get_qi(), operator!=(), and operator==().

get_scoor()

template<class T >

void CopiedPointCoorMesh< T >::get_scoor ( long  n, T &  b  ) const

inline

Definition at line 767 of file tline.h.

      {
    b = (*amesh)[n];
  }
#else
  { b = mesh[n]; }

Referenced by VirtCopiedPointCoorMesh< T >::get_scoor().

get_step()

template<class T >

int CopiedPointCoorMesh< T >::get_step ( long  n, T &  fstep  ) const

inline

Definition at line 795 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;
  }

Referenced by VirtCopiedPointCoorMesh< T >::get_step().

get_xmax()

template<class T >

T CopiedPointCoorMesh< T >::get_xmax ( void  ) const

inline

Definition at line 765 of file tline.h.

{ return xmax; }

Referenced by apeq_mant(), VirtCopiedPointCoorMesh< T >::get_xmax(), operator!=(), and operator==().

get_xmin()

template<class T >

T CopiedPointCoorMesh< T >::get_xmin ( void  ) const

inline

Definition at line 764 of file tline.h.

{ return xmin; }

Referenced by apeq_mant(), VirtCopiedPointCoorMesh< T >::get_xmin(), operator!=(), and operator==().

operator=()

template<class T >

CopiedPointCoorMesh< T > & CopiedPointCoorMesh< T >::operator= ( const CopiedPointCoorMesh< T > &  f )

inline

Definition at line 827 of file tline.h.

                                                                     {
    if (this != &f) {
      q = f.q;
      mesh = f.mesh;
      xmin = f.xmin;
      xmax = f.xmax;
      x_old = f.x_old;
      n_old = f.n_old;
#ifdef TLINE_REDUCE_TO_RAW_ARR
      amesh = &mesh;
#endif
    }
    return *this;
  }

print()

template<class T >

void CopiedPointCoorMesh< T >::print

(

std::ostream &

file

)

const

Definition at line 1096 of file tline.h.

                                                         {
  Ifile << "CopiedPointCoorMesh<T>:\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);
 
  Ifile << "mesh= " << mesh << '\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 operator<<(), and VirtCopiedPointCoorMesh< T >::print().

---

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

• tline.h