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

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

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