source: issm/trunk/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateNodesDiagnosticHutter.cpp@ 3999

Last change on this file since 3999 was 3999, checked in by Mathieu Morlighem, 15 years ago

updated CreateNodes
File size: 2.7 KB
Line 
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
15void 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}
Note: See TracBrowser for help on using the repository browser.