/* aux.h: Header for aux.c
 * 
 * Madhu Nayakkankuppam, 10/09/2002
 *
 */

#ifndef _AUX_H
#define _AUX_H
#endif

#define PI (3.14159265358979323846)

/*--------------------------------------------------------------------------*/

void allocate_matrix(int m,       /* # rows */
                     int n,       /* # cols */
                     double ***a  /* a is ptr to array of row ptrs */
                    );

void free_matrix(double **a    /* pointer to array of row pointers */
                );

void rhs(double **f,   /* rhs function */
         int nx,       /* # grid points in x-direction */
         int local_nx, /* local # grid points in y-direction */
         double h,     /* grid spacing */
         double ys     /* starting value of y */
        );

void bdry(double **a,    /* solution matrix */
          int nx,        /* # grid points along x */
          int local_nx,  /* # grid points along y */
          int my_rank,   /* process rank */
          int p          /* total number of processes */
         );

void init(double **a,   /* solution matrix */
          int nx,       /* # grid points along x */
          int local_nx  /* # grid points along y */
         );

double diff(double **a,    /* solution matrix */
            double **b,    /* solution matrix */
            int nx,        /* # points along x */
            int local_nx   /* # points along y */
           );

void dump_matrix(double **a,   /* matrix to be printed */
                 int m,        /* number of rows */
                 int n         /* number of columns */
                );