source: issm/trunk/src/c/objects/FemModel.cpp@ 6305

Last change on this file since 6305 was 6305, checked in by Mathieu Morlighem, 14 years ago

Added verbose class

File size: 5.7 KB
Line 
1/*!\file FemModel.c
2 * \brief: implementation of the FemModel object
3 */
4
5
6#ifdef HAVE_CONFIG_H
7 #include "config.h"
8#else
9#error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
10#endif
11
12#include "stdio.h"
13#include "../Container/Container.h"
14#include "../modules/ModelProcessorx/ModelProcessorx.h"
15#include "./objects.h"
16#include "../include/include.h"
17#include "../EnumDefinitions/EnumDefinitions.h"
18#include "../modules/modules.h"
19
20/*Object constructors and destructor*/
21/*FUNCTION FemModel::constructor {{{1*/
22FemModel::FemModel(ConstDataHandle IOMODEL,const int in_solution_type,const int* analyses,const int nummodels){
23
24 /*intermediary*/
25 int i;
26 int analysis_type;
27
28 /*Initialize internal data: */
29 this->nummodels=nummodels;
30 this->solution_type=in_solution_type;
31 this->analysis_counter=nummodels-1; //point to last analysis_type carried out.
32 this->results=new DataSet(); //not initialized by CreateDataSets
33
34 /*Dynamically allocate whatever is a list of length nummodels: */
35 analysis_type_list=(int*)xmalloc(nummodels*sizeof(int));
36 m_nodesets=(NodeSets**)xmalloc(nummodels*sizeof(NodeSets*));
37 m_ys=(Vec*)xmalloc(nummodels*sizeof(Vec));
38
39 /*Initialize: */
40 for(i=0;i<nummodels;i++)analysis_type_list[i]=analyses[i];
41 for(i=0;i<nummodels;i++)m_nodesets[i]=NULL;
42 for(i=0;i<nummodels;i++)m_ys[i]=NULL;
43
44 /*create datasets for all analyses*/
45 ModelProcessorx(&this->elements,&this->nodes,&this->vertices,&this->materials,&this->constraints,&this->loads,&this->parameters,IOMODEL,this->solution_type,nummodels,analyses);
46
47 /*do the post-processing of the datasets to get an FemModel that can actually run analyses: */
48 for(i=0;i<nummodels;i++){
49
50 if(VerbModProc()) _printf_(" Processing finite element model of analysis %s:\n",EnumToString(analysis_type_list[i]));
51 analysis_type=analysis_type_list[i];
52 this->SetCurrentConfiguration(analysis_type);
53
54 if(i==0){
55 if(VerbModProc()) _printf_(" create vertex degrees of freedom\n");
56 VerticesDofx(vertices,parameters); //only call once, we only have one set of vertices
57 }
58
59 if(VerbModProc()) _printf_(" resolve node constraints\n");
60 SpcNodesx(nodes,constraints,analysis_type);
61
62 if(VerbModProc()) _printf_(" create nodal degrees of freedom\n");
63 NodesDofx(nodes,parameters,analysis_type);
64
65 if(VerbModProc()) _printf_(" create nodal constraints vector\n");
66 CreateNodalConstraintsx(&m_ys[i],nodes,analysis_type);
67
68 if(VerbModProc()) _printf_(" create node sets\n");
69 BuildNodeSetsx(&m_nodesets[i], nodes,analysis_type);
70
71 if(VerbModProc()) _printf_(" configuring element and loads\n");
72 ConfigureObjectsx(elements, loads, nodes, vertices, materials,parameters);
73 }
74}
75/*}}}1*/
76/*FUNCTION FemModel::destructor {{{1*/
77FemModel::~FemModel(){
78
79 /*Intermediary*/
80 int i;
81
82 /*Delete all the datasets: */
83 xfree((void**)&analysis_type_list);
84 delete elements;
85 delete nodes;
86 delete vertices;
87 delete constraints;
88 delete loads;
89 delete materials;
90 delete parameters;
91 delete results;
92
93 for(i=0;i<nummodels;i++){
94 NodeSets* temp_nodesets=m_nodesets[i];
95 delete temp_nodesets;
96 Vec temp_ys=m_ys[i];
97 VecFree(&temp_ys);
98 }
99
100 /*Delete dynamically allocated arrays: */
101 xfree((void**)&m_nodesets);
102 xfree((void**)&m_ys);
103
104}
105/*}}}1*/
106
107/*Object management*/
108/*FUNCTION FemModel::Echo {{{1*/
109void FemModel::Echo(void){
110
111 printf("FemModel echo: \n");
112 printf(" number of fem models: %i\n",nummodels);
113 printf(" analysis_type_list: \n");
114 for(int i=0;i<nummodels;i++)printf(" %i: %s\n",i,EnumToString(analysis_type_list[i]));
115 printf(" current analysis_type: \n");
116 printf(" %i: %s\n",analysis_counter,EnumToString(analysis_type_list[analysis_counter]));
117
118
119}
120/*}}}*/
121
122/*Numerics: */
123/*FUNCTION FemModel::SetCurrentConfiguration(int configuration_type,int analysis_type){{{1*/
124void FemModel::SetCurrentConfiguration(int configuration_type,int analysis_type){
125
126 int verbose=0;
127
128 /*retrieve parameters: */
129 this->parameters->FindParam(&verbose,VerboseEnum);
130
131 /*Use configuration_type to setup the analysis counter, the configurations of objects etc ... but use
132 * analysis_type to drive the element numerics. This allows for use of 1 configuration_type for several
133 * analyses. For example: do a SurfaceSlopeX, SurfaceSlopeY, BedSlopeX and BedSlopeY analysis using the
134 * Slope configuration.*/
135
136 int found=-1;
137 for(int i=0;i<nummodels;i++){
138 if (analysis_type_list[i]==configuration_type){
139 found=i;
140 break;
141 }
142 }
143 if(found!=-1) analysis_counter=found;
144 else ISSMERROR("Could not find alias for analysis_type %s in list of FemModel analyses",EnumToString(configuration_type));
145
146 /*activate matrices/vectors: */
147 nodesets=m_nodesets[analysis_counter];
148 ys=m_ys[analysis_counter];
149
150 /*Now, plug analysis_counter and analysis_type inside the parameters: */
151 this->parameters->SetParam(analysis_counter,AnalysisCounterEnum);
152 this->parameters->SetParam(analysis_type,AnalysisTypeEnum);
153 this->parameters->SetParam(configuration_type,ConfigurationTypeEnum);
154
155 /*configure elements, loads and nodes, for this new analysis: */
156 this->elements->SetCurrentConfiguration(elements,loads, nodes,vertices, materials,parameters);
157 this->nodes->SetCurrentConfiguration(elements,loads, nodes,vertices, materials,parameters);
158 this->loads->SetCurrentConfiguration(elements, loads, nodes,vertices, materials,parameters);
159
160 /*take care of petsc options, that depend on this analysis type: */
161 PetscOptionsFromAnalysis(this->parameters,analysis_type);
162
163 if(verbose){
164 _printf_(" petsc Options set for analysis type: %s\n",EnumToString(analysis_type));
165 }
166
167
168}
169/*}}}1*/
170/*FUNCTION FemModel::SetCurrentConfiguration(int configuration_type){{{1*/
171void FemModel::SetCurrentConfiguration(int configuration_type){
172
173 /*overload: analysis_type = configuration_type: */
174 this->SetCurrentConfiguration(configuration_type,configuration_type);
175}
176/*}}}1*/
Note: See TracBrowser for help on using the repository browser.