File: USERD_get_nfaced_nodes_per_face.H

package info (click to toggle)
openfoam 4.1%2Bdfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 163,028 kB
  • ctags: 58,990
  • sloc: cpp: 830,760; sh: 10,227; ansic: 8,215; xml: 745; lex: 437; awk: 194; sed: 91; makefile: 77; python: 18
file content (74 lines) | stat: -rw-r--r-- 1,934 bytes parent folder | download | duplicates (3)
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
72
73
74

int USERD_get_nfaced_nodes_per_face
(
    int part_number,
    int *nfaced_npf_array
)
{
#ifdef ENSIGHTDEBUG
    Info<< "Entering: USERD_get_nfaced_nodes_per_face"
        << ", part_number = " << part_number
        << endl
        << flush;
#endif

    if (part_number == 1)
    {
        const cellShapeList& cellShapes = meshPtr->cellShapes();
        const cellList& cells = meshPtr->cells();
        const faceList& faces = meshPtr->faces();

        label nCells = cellShapes.size();
        label nFaced = 0;
        for (label n=0; n<nCells; n++)
        {
            label nFacesInCell = cells[n].size();
            labelList points = cellShapes[n];
            label nPoints = points.size();

            if ((nFacesInCell == 6) && (nPoints == 8))
            {}
            else if ((nFacesInCell == 4) && (nPoints == 4))
            {}
            else if (nFacesInCell == 5)
            {
                if (nPoints == 6)
                {}
                else if (nPoints == 5)
                {}
                else
                {
                    for (label i=0; i<nFacesInCell; i++)
                    {
                        label facei = cells[n][i];
                        label nFacePoints = faces[facei].size();
                        nfaced_npf_array[nFaced++] = nFacePoints;
                    }
                }
            }
            else
            {
                for (label i=0; i<nFacesInCell; i++)
                {
                    label facei = cells[n][i];
                    label nFacePoints = faces[facei].size();
                    nfaced_npf_array[nFaced++] = nFacePoints;
                }
            }
        }

    }
    else if (part_number < nPatches+2)
    {
        return Z_ERR;
    }
    else
    {
        return Z_ERR;
    }
#ifdef ENSIGHTDEBUG
    Info<< "Exiting: USERD_get_nfaced_nodes_per_face" << endl
        << flush;
#endif
    return Z_OK;
}