1 | /*
|
---|
2 | * CreateNodesDiagnosticHutter.c:
|
---|
3 | */
|
---|
4 |
|
---|
5 | #include "../../../DataSet/DataSet.h"
|
---|
6 | #include "../../../toolkits/toolkits.h"
|
---|
7 | #include "../../../io/io.h"
|
---|
8 | #include "../../../EnumDefinitions/EnumDefinitions.h"
|
---|
9 | #include "../../../objects/objects.h"
|
---|
10 | #include "../../../shared/shared.h"
|
---|
11 | #include "../../MeshPartitionx/MeshPartitionx.h"
|
---|
12 | #include "../../../include/include.h"
|
---|
13 | #include "../ModelProcessorx.h"
|
---|
14 |
|
---|
15 | void CreateNodesDiagnosticHutter(DataSet** pelements,DataSet** pnodes, DataSet** pvertices,DataSet** pmaterials, IoModel* iomodel,ConstDataHandle iomodel_handle){
|
---|
16 |
|
---|
17 | /*Intermediary*/
|
---|
18 | int i;
|
---|
19 | int totalnodes;
|
---|
20 | bool continuous_galerkin=true;
|
---|
21 |
|
---|
22 | /*DataSets: */
|
---|
23 | DataSet* nodes = NULL;
|
---|
24 |
|
---|
25 | /*First create nodes*/
|
---|
26 | nodes = new DataSet(NodesEnum);
|
---|
27 |
|
---|
28 | /*Now, is the flag ishutter on? otherwise, do nothing: */
|
---|
29 | if (!iomodel->ishutter)goto cleanup_and_return;
|
---|
30 |
|
---|
31 | /*Recover number of nodes already created in other analyses: */
|
---|
32 | totalnodes=iomodel->nodecounter;
|
---|
33 |
|
---|
34 | /*Continuous Galerkin partition of nodes: */
|
---|
35 | NodesPartitioning(&iomodel->my_nodes,iomodel->my_elements, iomodel->my_vertices, iomodel->my_bordervertices, iomodel, iomodel_handle,continuous_galerkin);
|
---|
36 |
|
---|
37 | /*First fetch data: */
|
---|
38 | if (iomodel->dim==3){
|
---|
39 | IoModelFetchData(&iomodel->deadgrids,NULL,NULL,iomodel_handle,"deadgrids");
|
---|
40 | IoModelFetchData(&iomodel->uppernodes,NULL,NULL,iomodel_handle,"uppergrids");
|
---|
41 | }
|
---|
42 | IoModelFetchData(&iomodel->gridonbed,NULL,NULL,iomodel_handle,"gridonbed");
|
---|
43 | IoModelFetchData(&iomodel->gridonsurface,NULL,NULL,iomodel_handle,"gridonsurface");
|
---|
44 | IoModelFetchData(&iomodel->gridonhutter,NULL,NULL,iomodel_handle,"gridonhutter");
|
---|
45 | IoModelFetchData(&iomodel->gridonicesheet,NULL,NULL,iomodel_handle,"gridonicesheet");
|
---|
46 | IoModelFetchData(&iomodel->gridoniceshelf,NULL,NULL,iomodel_handle,"gridoniceshelf");
|
---|
47 | IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
|
---|
48 | CreateNumberNodeToElementConnectivity(iomodel);
|
---|
49 |
|
---|
50 | for (i=0;i<iomodel->numberofvertices;i++){
|
---|
51 |
|
---|
52 | if(iomodel->my_vertices[i]){
|
---|
53 |
|
---|
54 | /*Add node to nodes dataset: */
|
---|
55 | nodes->AddObject(new Node(i+1,i,iomodel));
|
---|
56 |
|
---|
57 | }
|
---|
58 | }
|
---|
59 |
|
---|
60 | /*Clean fetched data: */
|
---|
61 | xfree((void**)&iomodel->deadgrids);
|
---|
62 | xfree((void**)&iomodel->gridonbed);
|
---|
63 | xfree((void**)&iomodel->gridonsurface);
|
---|
64 | xfree((void**)&iomodel->gridonhutter);
|
---|
65 | xfree((void**)&iomodel->uppernodes);
|
---|
66 | xfree((void**)&iomodel->gridonicesheet);
|
---|
67 | xfree((void**)&iomodel->gridoniceshelf);
|
---|
68 | xfree((void**)&iomodel->elements);
|
---|
69 | xfree((void**)&iomodel->numbernodetoelementconnectivity);
|
---|
70 |
|
---|
71 | /*All our datasets are already order by ids. Set presort flag so that later on, when sorting is requested on these
|
---|
72 | * datasets, it will not be redone: */
|
---|
73 | elements->Presort();
|
---|
74 |
|
---|
75 | cleanup_and_return:
|
---|
76 |
|
---|
77 | /*Assign output pointer: */
|
---|
78 | *pnodes=nodes;
|
---|
79 | }
|
---|