#include <MsmMap.h>

Inheritance diagram for msm::GridFixed< T, N >:

Public Member Functions
	GridFixed ()

void	init (const IndexRange &n)

void	set (int pia, int pni, int pja, int pnj, int pka, int pnk)

void	setbounds (int pia, int pib, int pja, int pjb, int pka, int pkb)

const T &	operator() (int i, int j, int k) const

const T &	operator() (const Ivec &n) const

const T &	elem (int i, int j, int k) const

T &	operator() (int i, int j, int k)

T &	operator() (const Ivec &n)

T &	elem (int i, int j, int k)

int	flatindex (int i, int j, int k) const

const T *	buffer () const

T *	buffer ()

void	reset (const T &t)

void	updateLower (const Ivec &n)

GridFixed< T, N > &	operator+= (const GridFixed< T, N > &g)

void	extract (GridFixed< T, N > &g)

Public Member Functions inherited from msm::IndexRange
	IndexRange ()

void	set (int pia, int pni, int pja, int pnj, int pka, int pnk)

void	setbounds (int pia, int pib, int pja, int pjb, int pka, int pkb)

int	ia () const

int	ib () const

int	ja () const

int	jb () const

int	ka () const

int	kb () const

int	ni () const

int	nj () const

int	nk () const

int	nn () const

Ivec	lower () const

Ivec	extent () const

int	operator<= (const IndexRange &n)

Friends
class	Grid< T >

Additional Inherited Members
Protected Attributes inherited from msm::IndexRange
Ivec	nlower

Ivec	nextent

Detailed Description

template<class T, int N>
class msm::GridFixed< T, N >

Definition at line 467 of file MsmMap.h.

Constructor & Destructor Documentation

template<class T, int N>

msm::GridFixed< T, N >::GridFixed ( )

inline

Definition at line 470 of file MsmMap.h.

470 { }

Member Function Documentation

template<class T, int N>

const T* msm::GridFixed< T, N >::buffer ( ) const

inline

Definition at line 531 of file MsmMap.h.

Referenced by msm::Grid< Vtype >::extract(), and msm::Grid< Vtype >::operator+=().

531 { return gdata; }

template<class T, int N>

T* msm::GridFixed< T, N >::buffer ( )

inline

Definition at line 532 of file MsmMap.h.

532 { return gdata; }

template<class T, int N>

const T& msm::GridFixed< T, N >::elem	(	int	i,
		int	j,
		int	k
	)		const

inline

Definition at line 494 of file MsmMap.h.

References msm::GridFixed< T, N >::flatindex(), msm::IndexRange::ia(), msm::IndexRange::ib(), msm::IndexRange::ja(), msm::IndexRange::jb(), msm::IndexRange::ka(), msm::IndexRange::kb(), and NAMD_die().

Referenced by msm::GridFixed< T, N >::operator()().

                                                {
         if (i<ia() || i>ib() || j<ja() || j>jb() || k<ka() || k>kb()) {
           char msg[200];
           snprintf(msg, sizeof(msg), "Grid indexing:\n"
               "ia=%d, ib=%d, i=%d\n"
               "ja=%d, jb=%d, j=%d\n"
               "ka=%d, kb=%d, k=%d\n",
               ia(), ib(), i, ja(), jb(), j, ka(), kb(), k);
           NAMD_die(msg);
         }
         return gdata[flatindex(i,j,k)];
       }

template<class T, int N>

T& msm::GridFixed< T, N >::elem	(	int	i,
		int	j,
		int	k
	)

inline

Definition at line 516 of file MsmMap.h.

References msm::GridFixed< T, N >::flatindex(), msm::IndexRange::ia(), msm::IndexRange::ib(), msm::IndexRange::ja(), msm::IndexRange::jb(), msm::IndexRange::ka(), msm::IndexRange::kb(), and NAMD_die().

                                    {
         if (i<ia() || i>ib() || j<ja() || j>jb() || k<ka() || k>kb()) {
           char msg[200];
           snprintf(msg, sizeof(msg), "Grid indexing:\n"
               "ia=%d, ib=%d, i=%d\n"
               "ja=%d, jb=%d, j=%d\n"
               "ka=%d, kb=%d, k=%d\n",
               ia(), ib(), i, ja(), jb(), j, ka(), kb(), k);
           NAMD_die(msg);
         }
         return gdata[flatindex(i,j,k)];
       }

template<class T, int N>

void msm::GridFixed< T, N >::extract ( GridFixed< T, N > & g )

inline

Definition at line 571 of file MsmMap.h.

References ASSERT, msm::Ivec::i, msm::IndexRange::IndexRange(), msm::Ivec::j, msm::Ivec::k, msm::IndexRange::nextent, msm::IndexRange::ni(), and msm::IndexRange::nlower.

                                       {
         ASSERT(IndexRange(g) <= IndexRange(*this));
         int gni = g.nextent.i;
         int gnj = g.nextent.j;
         int gnk = g.nextent.k;
         int index = 0;
         int ni = nextent.i;
         int nij = nextent.i * nextent.j;
         int koff = (g.nlower.k - nlower.k) * nij
           + (g.nlower.j - nlower.j) * ni + (g.nlower.i - nlower.i);
         T *gbuf = g.gdata.buffer();
         const T *buf = gdata.buffer();
         for (int k = 0;  k < gnk;  k++) {
           int jkoff = k * nij + koff;
           for (int j = 0;  j < gnj;  j++) {
             int ijkoff = j * ni + jkoff;
             for (int i = 0;  i < gni;  i++, index++) {
               gbuf[index] = buf[i + ijkoff];
             }
           }
         }
       }

template<class T, int N>

int msm::GridFixed< T, N >::flatindex	(	int	i,
		int	j,
		int	k
	)		const

inline

Definition at line 528 of file MsmMap.h.

References msm::IndexRange::ia(), msm::IndexRange::ja(), msm::IndexRange::ka(), msm::IndexRange::ni(), and msm::IndexRange::nj().

Referenced by msm::GridFixed< T, N >::elem(), and msm::GridFixed< T, N >::operator()().

                                                {
         return ((k-ka())*nj() + (j-ja()))*ni() + (i-ia());
       }

template<class T, int N>

void msm::GridFixed< T, N >::init ( const IndexRange & n )

inline

Definition at line 471 of file MsmMap.h.

References ASSERT, msm::IndexRange::extent(), msm::Ivec::i, msm::Ivec::j, msm::Ivec::k, msm::IndexRange::lower(), msm::IndexRange::nextent, and msm::IndexRange::nlower.

                                      {
         nlower = n.lower();
         nextent = n.extent();
         ASSERT(nextent.i * nextent.j * nextent.k <= N);
       }

template<class T, int N>

const T& msm::GridFixed< T, N >::operator()	(	int	i,
		int	j,
		int	k
	)		const

inline

Definition at line 484 of file MsmMap.h.

References msm::GridFixed< T, N >::elem(), and msm::GridFixed< T, N >::flatindex().

Referenced by msm::GridFixed< T, N >::operator()().

                                                      {
 #ifdef DEBUG_MSM
         return elem(i,j,k);
 #else
         return gdata[flatindex(i,j,k)];
 #endif
       }

template<class T, int N>

const T& msm::GridFixed< T, N >::operator() ( const Ivec & n ) const

inline

Definition at line 491 of file MsmMap.h.

References msm::Ivec::i, msm::Ivec::j, msm::Ivec::k, and msm::GridFixed< T, N >::operator()().

                                                {
         return this->operator()(n.i, n.j, n.k);
       }

template<class T, int N>

T& msm::GridFixed< T, N >::operator()	(	int	i,
		int	j,
		int	k
	)

inline

Definition at line 506 of file MsmMap.h.

References msm::GridFixed< T, N >::elem(), and msm::GridFixed< T, N >::flatindex().

                                          {
 #ifdef DEBUG_MSM
         return elem(i,j,k);
 #else
         return gdata[flatindex(i,j,k)];
 #endif
       }

template<class T, int N>

T& msm::GridFixed< T, N >::operator() ( const Ivec & n )

inline

Definition at line 513 of file MsmMap.h.

References msm::Ivec::i, msm::Ivec::j, msm::Ivec::k, and msm::GridFixed< T, N >::operator()().

                                    {
         return this->operator()(n.i, n.j, n.k);
       }

template<class T, int N>

GridFixed<T,N>& msm::GridFixed< T, N >::operator+= ( const GridFixed< T, N > & g )

inline

Definition at line 545 of file MsmMap.h.

References ASSERT, msm::Ivec::i, msm::IndexRange::IndexRange(), msm::Ivec::j, msm::Ivec::k, msm::IndexRange::nextent, msm::IndexRange::ni(), and msm::IndexRange::nlower.

                                                           {
         ASSERT(IndexRange(g) <= IndexRange(*this));
         int gni = g.nextent.i;
         int gnj = g.nextent.j;
         int gnk = g.nextent.k;
         int index = 0;
         int ni = nextent.i;
         int nij = nextent.i * nextent.j;
         int koff = (g.nlower.k - nlower.k) * nij
           + (g.nlower.j - nlower.j) * ni + (g.nlower.i - nlower.i);
         const T *gbuf = g.gdata.buffer();
         T *buf = gdata.buffer();
         for (int k = 0;  k < gnk;  k++) {
           int jkoff = k * nij + koff;
           for (int j = 0;  j < gnj;  j++) {
             int ijkoff = j * ni + jkoff;
             for (int i = 0;  i < gni;  i++, index++) {
               buf[i + ijkoff] += gbuf[index];
             }
           }
         }
         return(*this);
       }

template<class T, int N>

void msm::GridFixed< T, N >::reset ( const T & t )

inline

Definition at line 535 of file MsmMap.h.

References msm::IndexRange::nn().

                              {
         int len = nn();
         for (int n = 0;  n < len;  n++) { gdata[n] = t; }
       }

template<class T, int N>

void msm::GridFixed< T, N >::set	(	int	pia,
		int	pni,
		int	pja,
		int	pnj,
		int	pka,
		int	pnk
	)

inline

Definition at line 476 of file MsmMap.h.

References ASSERT, msm::Ivec::i, msm::Ivec::j, msm::Ivec::k, msm::IndexRange::nextent, and msm::IndexRange::set().

                                                                      {
         IndexRange::set(pia, pni, pja, pnj, pka, pnk);
         ASSERT(nextent.i * nextent.j * nextent.k <= N);
       }

template<class T, int N>

void msm::GridFixed< T, N >::setbounds	(	int	pia,
		int	pib,
		int	pja,
		int	pjb,
		int	pka,
		int	pkb
	)

inline

Definition at line 480 of file MsmMap.h.

References ASSERT, msm::Ivec::i, msm::Ivec::j, msm::Ivec::k, msm::IndexRange::nextent, and msm::IndexRange::setbounds().

                                                                            {
         IndexRange::setbounds(pia, pib, pja, pjb, pka, pkb);
         ASSERT(nextent.i * nextent.j * nextent.k <= N);
       }