poisson.h

Go to the documentation of this file.

#ifndef _POISSON_H_

#include "cstream.h"
#include "rz_array.h"

#ifndef _GRID_H_
#include "grid.h"
#endif

typedef struct pois_grid_t pois_grid_t;

struct pois_grid_t {
  RECT_COORDS;
  LEAF_FIELDS(pois_grid_t);
  int ext_bound;

  rz_array_t *phi;
  rz_array_t *charge;
  rz_array_t *error;
};

typedef struct pois_problem_t pois_problem_t;

struct pois_problem_t {
  int max_level;
  int extra_levels;
  double max_error;

  int bnd_right;
  int bnd_top;
  int bnd_bottom;
};

typedef struct pois_boundaries_t {
  REAL *left, *right, *top, *bottom;
  int r, z;
} pois_boundaries_t;

#define UNCHECK_ER_RZ(grid_, ir_, iz_)             \
  (((RZ(grid_->phi, (ir_) + 1, iz_) - RZ(grid_->phi, ir_, iz_))      \
    / dr[grid_->level]))

#define UNCHECK_EZ_RZ(grid_, ir_, iz_)             \
  (((RZ(grid_->phi, ir_, (iz_) + 1) - RZ(grid_->phi, ir_, iz_))   \
    / dz[grid_->level]))

#define ER_RZ(grid_, ir_, iz_)         \
  (((ir_) < ((grid_)->r1 + 1) && ((ir_) > ((grid_)->r0 - 2)))? \
   UNCHECK_ER_RZ(grid_, ir_, iz_): 0)

#define EZ_RZ(grid_, ir_, iz_)         \
  (((iz_) < ((grid_)->z1 + 1) && ((iz_) > ((grid_)->z0 - 2)))? \
   UNCHECK_EZ_RZ(grid_, ir_, iz_): 0)

#define _POISSON_H_
#endif