fas.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   a:  coefficients for equation on grid n, computed with "afromc"
00017   b:  coefficients "   "
00018   ei: coefficients inverse "  "
00019 
00020   scalar parameters:
00021   ------------------
00022   n:      dimension of arrays
00023 
00024  ON OUTPUT:
00025 
00026   everything the same, except
00027 
00028   p:      contains the pressure 
00029 */
00030 
00031 #ifndef MG_H
00032 #define MG_H
00033 
00034 #define NGERROR         -1
00035 #define NXERROR         -2
00036 #define NYERROR         -3
00037 #define MEMERROR        -5
00038 
00039 
00040 extern void relax(real2D u, real2D rhs, 
00041                   real2D a, real2D b, real2D ei, real2D cc,
00042                   int n1, int n2);
00043 
00044 extern int computeng(int n1, int n2);
00045 
00046 extern int initpressure(int n1, int n2, int ng);
00047 
00048 extern void mginit(real2D u, real2D p, real2D a, real2D b, real2D ei, 
00049                    int n1, int n2, int ng); 
00050 
00051 extern void mgvcycle(int n, int nx, int ny, int npre, int npost);
00052 
00053 extern void mgfas(real2D u, real2D p, real2D a, real2D b, real2D ei, 
00054                   int n1, int n2, int ng, int npre, int npost, int ncycle);
00055 #endif
00056