#ifndef _FISPHACK_H_

/* Arguments for the boundary conditions in fishpack:

   The first name is the b.c. at the lowest value of the variable,

   the second one is the b.c. at the highest value.

     DIR = Dirichlet, 

     NEU = Neumann, 

     UNS = Unspecified (when r = 0 for cylindrical calculations).

*/

#define FISH_PERIODIC 0

#define FISH_DIR_DIR  1

#define FISH_DIR_NEU  2

#define FISH_NEU_NEU  3

#define FISH_NEU_DIR  4

#define FISH_UNS_DIR  5

#define FISH_UNS_NEU  6

/* This is the maximum number of grid cells that fishpack can handle

   without accumulating larger and larger roundoff errors

   (it's only an approximation: not checked throughly so if you notice

   that you are getting strange things such as artificial lines in

   the electric field, this is the first thing you should check). */

#define FISH_MAX_GRIDPOINTS   1700

#define LOG2 0.69314718055994530941723212145818

#define FISH_WORK(M_, N_)  (13 * (M_) + 4 * (N_)		\

			    + (M_) * (int) (log(N_) / LOG2))

#define FISH_ERROR_MAX 13

const char *hstcyl_error_str[FISH_ERROR_MAX];

void 

fish_hstcyl (double r0, double r1, int nr, 

	     int rbndcnd, double *bcr0, double *bcr1,

	     double z0, double z1, int nz, 

	     int zbndcnd, double *bcz0, double *bcz1,

	     double lambda, double s, double *f, int idimf);

void 

fish_hstcrt (double r0, double r1, int nr, 

	     int rbndcnd, double *bcr0, double *bcr1,

	     double z0, double z1, int nz, 

	     int zbndcnd, double *bcz0, double *bcz1,

	     double lambda, double *f, int idimf);

#define _FISHPACK_H_

#endif /* _FISHPACK_H_ */