/* Wrapper for the HSTCYL function. 
   Alejandro Luque Estepa - 2006
*/

#include "stdlib.h"
#include "stdio.h"
#include "../include/fishpack.h"
#include "math.h"
#include "assert.h"

const char *hstcyl_error_str[FISH_ERROR_MAX] = 
  {"No error",
   "A .LT. 0",
   "A .GE. B",
   "MBDCND .LT. 1 OR MBDCND .GT. 6",
   "C .GE. D",
   "N .LE. 2",
   "NBDCND .LT. 0 OR NBDCND .GT. 4",
   "A = 0 AND MBDCND = 1,2,3, OR 4",
   "A .GT. 0 AND MBDCND .GE. 5",
   "M .LE. 2",
   "IDIMF .LT. M",
   "LAMBDA .GT. 0",
   "A=0, MBDCND .GE. 5, ELMBDA .NE. 0"};

int number=0;

void mhstcyl_(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, double *pertrb, int *ierror);

/* Calls hstcyl and checks for errors. If an error occurs, prints the
 corresponding message and exists. */
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)
{
  double pertrb;
  int ierror;

  mhstcyl_(&r0, &r1, &nr, &rbndcnd, bcr0, bcr1, 
	   &z0, &z1, &nz, &zbndcnd, bcz0, bcz1,
	   &lambda, &s, f, &idimf, &pertrb, &ierror);
  number+=1;

/*  fprintf(stderr, "# calls mhstcyl       =%d\n",number); */
  if (pertrb != 0.0) {
    fprintf (stderr, "%s: ERROR: Undefined solution to the Poisson equation.",
	     __func__);
    fprintf (stderr, "\nTry to change your boundary conditions.\n");
    exit(-1);
  }

  if (0 != ierror) {
    if (ierror <= FISH_ERROR_MAX) {
      fprintf (stderr, "%s: ERROR: %s\n", __func__, hstcyl_error_str[ierror]);
      fprintf (stderr, "The call to fishpack was\n");
      fprintf (stderr, "mhstcyl_ (A = %f, B = %f, M = %d,\n",  r0, r1, nr);
      fprintf (stderr, "         MBDCND = %d, ...\n", rbndcnd);
      fprintf (stderr, "         C = %f, D = %f, N = %d,\n",  z0, z1, nz);
      fprintf (stderr, "         NBDCND = %d, ...\n", zbndcnd);
      fprintf (stderr, "         ELMBDA = %f, ...\n", lambda);
      fprintf (stderr, "         ES = %f, ...\n", s);
      fprintf (stderr, "         IDIMF = %d, ...)\n", idimf);
    } else {
      fprintf (stderr, "%s: Unknown error #%d\n", __func__, ierror);
    }

    exit(-1);
  }
}