IntNGP.h

Go to the documentation of this file.

// -*- C++ -*-
/***************************************************************************
 *
 * The IPPL Framework
 * 
 *
 * Visit http://people.web.psi.ch/adelmann/ for more details
 *
 ***************************************************************************/

#ifndef INT_NGP_H
#define INT_NGP_H

/* IntNGP.h -- Definition of simple class to perform nearest-grid-point
   interpolation of data for a single particle to or from IPPL Field. */

// include files
#include "Particle/Interpolator.h"
#include "Field/Field.h"


// forward declaration
class IntNGP;

// specialization of InterpolatorTraits

template <class T, unsigned Dim>
struct InterpolatorTraits<T,Dim,IntNGP> {
  typedef NDIndex<Dim> Cache_t;
};


// IntNGP class definition
class IntNGP : public Interpolator {

public:
  // constructor/destructor
  IntNGP() {}
  ~IntNGP() {}

  // gather/scatter functions

  // scatter particle data into Field using particle position and mesh
  template <class FT, unsigned Dim, class M, class C, class PT>
  static
  void scatter(const FT& pdata, Field<FT,Dim,M,C>& f,
               const Vektor<PT,Dim>& ppos, const M& mesh) {
    // find nearest-grid-point for particle position, store in NDIndex obj
    NDIndex<Dim> ngp = FindNGP(mesh, ppos, CenteringTag<C>());
    // scatter data value to Field ... this assumes that the Field
    // data point is local to this processor, if not an error will be printed.

    CompressedBrickIterator<FT,Dim> fiter = getFieldIter(f,ngp);
    *fiter += pdata;

    return;
  }

  // scatter particle data into Field using particle position and mesh
  // and cache mesh information for reuse
  template <class FT, unsigned Dim, class M, class C, class PT>
  static
  void scatter(const FT& pdata, Field<FT,Dim,M,C>& f,
               const Vektor<PT,Dim>& ppos, const M& mesh,
               NDIndex<Dim>& ngp) {
    // find nearest-grid-point for particle position, store in NDIndex obj
    ngp = FindNGP(mesh, ppos, CenteringTag<C>());
    // scatter data value to Field ... this assumes that the Field
    // data point is local to this processor, if not an error will be printed.

    CompressedBrickIterator<FT,Dim> fiter = getFieldIter(f,ngp);
    *fiter += pdata;

    return;
  }

  // scatter particle data into Field using cached mesh information
  template <class FT, unsigned Dim, class M, class C>
  static
  void scatter(const FT& pdata, Field<FT,Dim,M,C>& f,
               const NDIndex<Dim>& ngp) {
    // scatter data value to Field ... this assumes that the Field
    // data point is local to this processor, if not an error will be printed.

    CompressedBrickIterator<FT,Dim> fiter = getFieldIter(f,ngp);
    *fiter += pdata;

    return;
  }

  // gather particle data from Field using particle position and mesh
  template <class FT, unsigned Dim, class M, class C, class PT>
  static
  void gather(FT& pdata, const Field<FT,Dim,M,C>& f,
              const Vektor<PT,Dim>& ppos, const M& mesh) {
    // find nearest-grid-point for particle position, store in NDIndex obj
    NDIndex<Dim> ngp = FindNGP(mesh, ppos, CenteringTag<C>());
    // gather Field value to particle data ... this assumes that the Field
    // data point is local to this processor, if not an error will be printed.

    CompressedBrickIterator<FT,Dim> fiter = getFieldIter(f,ngp);
    pdata = *fiter;

    return;
  }

  // gather particle data from Field using particle position and mesh
  // and cache mesh information for reuse
  template <class FT, unsigned Dim, class M, class C, class PT>
  static
  void gather(FT& pdata, const Field<FT,Dim,M,C>& f,
              const Vektor<PT,Dim>& ppos, const M& mesh,
              NDIndex<Dim>& ngp) {
    // find nearest-grid-point for particle position, store in NDIndex obj
    ngp = FindNGP(mesh, ppos, CenteringTag<C>());
    // gather Field value to particle data ... this assumes that the Field
    // data point is local to this processor, if not an error will be printed.

    CompressedBrickIterator<FT,Dim> fiter = getFieldIter(f,ngp);
    pdata = *fiter;

    return;
  }

  // gather particle data from Field using cached mesh information
  template <class FT, unsigned Dim, class M, class C>
  static
  void gather(FT& pdata, const Field<FT,Dim,M,C>& f,
              const NDIndex<Dim>& ngp) {
    // gather Field value to particle data ... this assumes that the Field
    // data point is local to this processor, if not an error will be printed.

    CompressedBrickIterator<FT,Dim> fiter = getFieldIter(f,ngp);
    pdata = *fiter;

    return;
  }

};

#endif // INT_NGP_H

/***************************************************************************
 * $RCSfile: IntNGP.h,v $   $Author: adelmann $
 * $Revision: 1.1.1.1 $   $Date: 2003/01/23 07:40:28 $
 * IPPL_VERSION_ID: $Id: IntNGP.h,v 1.1.1.1 2003/01/23 07:40:28 adelmann Exp $ 
 ***************************************************************************/