mg.h Source File

00001 /* --------------------------------------------------------------------- */
00002 /* 
00003    SCCS information: %W% %G%
00004 */
00005 /* --------------------------------------------------------------------- */
00006 
00007 /*  These subroutines invert the pressure equation
00008 
00009   ON INPUT:
00010  
00011   nx by ny double arrays:
00012   ---------------------
00013   u:  initial right hand side (divergence field)
00014   p:  anything
00015   cc: concentrations or volume fractions of phase 1
00016 
00017   scalar parameters:
00018   ------------------
00019   n:      dimension of arrays
00020 
00021  ON OUTPUT:
00022 
00023   everything the same, except
00024 
00025   p:      contains the pressure 
00026 */
00027 
00028 #ifndef MG_H
00029 #define MG_H
00030 
00031 #define NGERROR         -1
00032 #define NXERROR         -2
00033 #define NYERROR         -3
00034 #define MEMERROR        -5
00035 
00036 
00037 extern void relax(real2D p, real2D rhs, real2D ap,
00038                   real2D c0, real2D c1, real2D c2, real2D c3, real2D c4,
00039                   int nx, int ny);
00040 
00041 extern void resid(real2D res, real2D p, real2D rhs, real2D ap, 
00042                   real2D c0, real2D c1, real2D c2, real2D c3, real2D c4,
00043                   int nx, int ny);
00044 
00045 extern int computeng(int n1, int n2);
00046 
00047 extern int initpressure(int n1, int n2, int ng);
00048 
00049 extern int mgsolve(real2D u, real2D p, real2D ap,
00050                    real2D c0, real2D c1, real2D c2, real2D c3, real2D c4,
00051                    interface in, real lredis, real maxdiv,
00052                    int n1, int n2, int ng) ;
00053 
00054 extern void mgvcycle(int n, int nx, int ny, int npre, int npost);
00055 
00056 #endif
00057