ParticleLayoutFromGrid.h

Go to the documentation of this file.

//
// C++ Interface: ParticleLayoutFromGrid
//
// Description:
//
//
// Author: Roman Geus <geus@maxwell>, (C) 2005
//
// Copyright: See COPYING file that comes with this distribution
//
//
#ifndef ParticleLayoutFromGrid_H
#define ParticleLayoutFromGrid_H

#include "Ippl.h"
#include "extpde.h"

class ParticleLayoutFromGrid : public ParticleLayout<double,3> {
public:
    typedef int pair_t;
    typedef pair_t* pair_iterator;
    typedef ParticleAttrib<SingleParticlePos_t> ParticlePos_t;
    typedef ParticleAttrib<Index_t>             ParticleIndex_t;
    typedef Vektor<double,3>                    Vector_t;
    ParticleLayoutFromGrid(Grid* grid, Domain* domain, D3vector& local_min, D3vector& local_max) :
        grid_(grid),
        geom_domain_(domain),
        my_corner_min_(local_min),
        my_corner_max_(local_max)
    {}
    void update(IpplParticleBase< ParticleLayoutFromGrid >& particles);
    inline bool is_local_pos(const Vector_t& x) const {
      const D3vector y(x[0],x[1],x[2]);
      return is_inside_box(my_corner_min_, my_corner_max_, y);
    }
    inline bool is_local_pos(const D3vector& x) const {
      return is_inside_box(my_corner_min_, my_corner_max_, x);
    }

    bool is_in_domain(const D3vector& x) const {
        return geom_domain_->point_in_domain(x);
    }
    bool is_in_domain(const Vector_t& x) const {
      const D3vector y(x[0],x[1],x[2]);
      return geom_domain_->point_in_domain(y);
    }

    inline static bool is_inside_box(const D3vector& corner_min, const D3vector& corner_max, const D3vector& x) {
        return corner_min[0] <= x[0] && x[0] < corner_max[0] &&
               corner_min[1] <= x[1] && x[1] < corner_max[1] &&
               corner_min[2] <= x[2] && x[2] < corner_max[2];
    }

private:
    void apply_bconds(ParticlePos_t& R);
    size_t redistribute_particles(IpplParticleBase< ParticleLayoutFromGrid >& particles);
    Grid* grid_;
    Domain* geom_domain_;
    D3vector my_corner_min_;
    D3vector my_corner_max_;
};

#endif