Main Page   Data Structures   File List   Data Fields   Globals  

vmg.h File Reference

This graph shows which files directly or indirectly include this file:

Included by dependency graph

Go to the source code of this file.

Defines

#define NGERROR   1
#define NXERROR   2
#define NYERROR   3
#define MEMERROR   4

Functions

int computeng (int n1, int n2)
int initpressure (int n1, int n2, int ng)
void mglin (real2D u, real2D p, real2D cc, real2D a, real2D c, real2D ei, real epsilon, int n1, int n2, int ng)


Define Documentation

#define MEMERROR   4
 

Definition at line 7 of file vmg.h.

#define NGERROR   1
 

Definition at line 4 of file vmg.h.

#define NXERROR   2
 

Definition at line 5 of file vmg.h.

#define NYERROR   3
 

Definition at line 6 of file vmg.h.


Function Documentation

int computeng int    n1,
int    n2
 

Definition at line 35 of file fas.c.

References MIN.

Referenced by initialize().

00036 {
00037   int nmin, n, nn, j;
00038 
00039   nmin = min(n1, n2);
00040   
00041   j = nmin - 2;  n = 1;
00042   while (j > 2) 
00043   {  j /= 2;   n++;  }
00044   
00045   nn = 1;  
00046   for (j = 1; j <= n; j++) 
00047     nn *= 2;       
00048   if((nn + 2) != nmin)
00049      return NGERROR;
00050 
00051   if ((n1 - 2) % (nmin - 2) != 0)
00052      return NXERROR;
00053   if ((n2 - 2) % (nmin - 2) != 0)
00054      return NYERROR;
00055   
00056   return n;
00057 }

int initpressure int    n1,
int    n2,
int    ng
 

Definition at line 61 of file fas.c.

References mymalloc(), real, and real2D.

Referenced by initialize().

00062 {
00063   int j;
00064   
00065   irhs = (real2D *) malloc((ng + 1) * sizeof(real2D));
00066   irho = (real2D *) malloc((ng + 1) * sizeof(real2D));
00067   itau = (real2D *) malloc((ng + 1) * sizeof(real2D));
00068   ip = (real2D *) malloc((ng + 1) * sizeof(real2D));
00069   ia = (real2D *) malloc((ng + 1) * sizeof(real2D));
00070   ib = (real2D *) malloc((ng + 1) * sizeof(real2D));
00071   iei = (real2D *) malloc((ng + 1) * sizeof(real2D));
00072   itemp = (real2D *) malloc((ng + 1) * sizeof(real2D));
00073   icc = (real2D *) malloc((ng + 1) * sizeof(real2D));
00074   
00075   if (icc EQUALS NULL)
00076       return MEMERROR;
00077   
00078   irhs[ng] = (real2D)mymalloc(sizeof(real), n1, n2);
00079   itau[ng] = (real2D)mymalloc(sizeof(real), n1, n2); 
00080   itemp[ng] = (real2D)mymalloc(sizeof(real), n1, n2);
00081 
00082   n1 = n1 / 2 + 1;
00083   n2 = n2 / 2 + 1;
00084 
00085   for (j = ng - 1; j >= 1; j--)
00086     {
00087       irhs[j] = (real2D)mymalloc(sizeof(real), n1, n2); 
00088       irho[j] = (real2D)mymalloc(sizeof(real), n1, n2); 
00089       itau[j] = (real2D)mymalloc(sizeof(real), n1, n2); 
00090       ip[j] = (real2D)mymalloc(sizeof(real), n1, n2); 
00091       ia[j] = (real2D)mymalloc(sizeof(real), n1, n2);
00092       ib[j] = (real2D)mymalloc(sizeof(real), n1, n2); 
00093       iei[j] = (real2D)mymalloc(sizeof(real), n1, n2); 
00094       icc[j] = (real2D)mymalloc(sizeof(real), n1, n2);
00095       itemp[j] = (real2D)mymalloc(sizeof(real), n1, n2);
00096 
00097       if (itemp[j] EQUALS NULL)
00098         return MEMERROR;
00099 
00100       n1 = n1 / 2 + 1;
00101       n2 = n2 / 2 + 1;
00102     }
00103     
00104   return 0;
00105 } /* end initpressure() */

void mglin real2D    u,
real2D    p,
real2D    cc,
real2D    a,
real2D    c,
real2D    ei,
real    epsilon,
int    n1,
int    n2,
int    ng
 

Referenced by pressure().


Generated on Wed Feb 19 22:29:04 2003 for Markers by doxygen1.2.18