This graph shows which files directly or indirectly include this file:
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) |
|
|
|
|
|
|
|
|
|
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 } |
|
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() */ |
|
Referenced by pressure(). |