source: issm/trunk/src/c/modules/ModelProcessorx/Balancedthickness2/CreateNodesBalancedthickness2.cpp@ 4211

Last change on this file since 4211 was 4211, checked in by Eric.Larour, 15 years ago

New Nodes class
File size: 2.7 KB
Line 
1/*
2 * CreateNodesBalancedthickness2.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 "../../../include/include.h"
12#include "../../MeshPartitionx/MeshPartitionx.h"
13#include "../ModelProcessorx.h"
14
15void CreateNodesBalancedthickness2(Nodes** pnodes, IoModel* iomodel,ConstDataHandle iomodel_handle){
16
17 /*Intermediary*/
18 int i,j;
19 bool continuous_galerkin=false;
20 int node_id;
21 int vertex_id;
22 int io_index;
23
24 /*DataSets: */
25 Nodes* nodes = NULL;
26
27 /*Recover pointer: */
28 nodes=*pnodes;
29
30 /*Create nodes if they do not exist yet*/
31 if(!nodes) nodes = new Nodes(NodesEnum);
32
33 /*Continuous Galerkin partition of nodes: */
34 NodesPartitioning(&iomodel->my_nodes,iomodel->my_elements, iomodel->my_vertices, iomodel->my_bordervertices, iomodel, iomodel_handle,continuous_galerkin);
35
36 /*Create nodes: */
37 IoModelFetchData(&iomodel->gridonbed,NULL,NULL,iomodel_handle,"gridonbed");
38 IoModelFetchData(&iomodel->gridonsurface,NULL,NULL,iomodel_handle,"gridonsurface");
39 IoModelFetchData(&iomodel->gridonhutter,NULL,NULL,iomodel_handle,"gridonhutter");
40 IoModelFetchData(&iomodel->gridonicesheet,NULL,NULL,iomodel_handle,"gridonicesheet");
41 IoModelFetchData(&iomodel->gridoniceshelf,NULL,NULL,iomodel_handle,"gridoniceshelf");
42 IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
43 if (iomodel->dim==3){
44 IoModelFetchData(&iomodel->deadgrids,NULL,NULL,iomodel_handle,"deadgrids");
45 IoModelFetchData(&iomodel->uppernodes,NULL,NULL,iomodel_handle,"uppergrids");
46 ISSMERROR("Not implemented yet");
47 }
48
49 /*Build Nodes dataset -> 3 for each element (Discontinuous Galerkin)*/
50 for (i=0;i<iomodel->numberofelements;i++){
51 for (j=0;j<3;j++){
52
53 if(iomodel->my_nodes[3*i+j]){
54
55 //Get index of the vertex on which the current node is located
56 vertex_id=(int)*(iomodel->elements+3*i+j); //(Matlab indexing)
57 io_index=vertex_id-1; //(C indexing)
58 ISSMASSERT(vertex_id>0 && vertex_id<=iomodel->numberofvertices);
59
60 //Compute Node id
61 node_id=iomodel->nodecounter+3*i+j+1;
62
63 /*Add node to nodes dataset: */
64 nodes->AddObject(new Node(node_id,node_id-1,vertex_id,io_index,iomodel,Balancedthickness2AnalysisEnum));
65
66 }
67 }
68 }
69
70 /*Clean fetched data: */
71 xfree((void**)&iomodel->deadgrids);
72 xfree((void**)&iomodel->gridonbed);
73 xfree((void**)&iomodel->gridonsurface);
74 xfree((void**)&iomodel->gridonhutter);
75 xfree((void**)&iomodel->uppernodes);
76 xfree((void**)&iomodel->gridonicesheet);
77 xfree((void**)&iomodel->gridoniceshelf);
78 xfree((void**)&iomodel->elements);
79
80 /*Assign output pointer: */
81 *pnodes=nodes;
82}
Note: See TracBrowser for help on using the repository browser.