1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
/*
* MLACTS.C - test the mapping from LR to AC meshes
*
* Source Version: 2.0
* Software Release #92-0043
*
*/
#include "cpyright.h"
#include "pml.h"
/*--------------------------------------------------------------------------*/
/*--------------------------------------------------------------------------*/
/* MAIN - a sample program */
main()
{PM_mesh_topology *mt;
REAL *x, *y, *px, *py;
int i, j, k, l;
int nd, nc, np;
int kmax, lmax, k1, k2, l1, l2, ord;
long *ncell, *pcell;
kmax = 6;
lmax = 8;
k1 = 2;
k2 = 4;
l1 = 1;
l2 = 4;
x = FMAKE_N(REAL, kmax*lmax, "MLACTS.C:x");
y = FMAKE_N(REAL, kmax*lmax, "MLACTS.C:y");
for (l = 0; l < lmax; l++)
for (k = 0; k < kmax; k++)
{i = l*kmax + k;
x[i] = k + 1.0;
y[i] = l + 1.0;};
px = x;
py = y;
ord = NGB_CELL;
mt = PM_lr_ac_mesh_2d(&px, &py, kmax, lmax, k1, k2, l1, l2, ord);
/* print cell info */
nd = mt->n_dimensions;
for (k = nd; 0 < k; k--)
{ncell = mt->boundaries[k];
nc = mt->n_cells[k];
np = mt->n_bound_params[k];
for (i = 0; i < nc; i++)
{pcell = ncell + i*np;
for (j = 0; j < np; j++)
PRINT(STDOUT, " %4d", pcell[j]);
PRINT(STDOUT, "\n");};
PRINT(STDOUT, "\n");};
SFREE(x);
SFREE(y);
SFREE(px);
SFREE(py);
return(0);}
/*--------------------------------------------------------------------------*/
/*--------------------------------------------------------------------------*/
|
