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

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

merged trunk-jpl and trunk for revision 11236
File size: 5.4 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 "../io/io.h"
16#include "./objects.h"
17#include "../include/include.h"
18#include "../EnumDefinitions/EnumDefinitions.h"
19#include "../modules/modules.h"
20
21/*Object constructors and destructor*/
22/*FUNCTION FemModel::constructor {{{1*/
23FemModel::FemModel(char* inputfilename, char* outputfilename, const int in_solution_type,const int* analyses,const int nummodels){
24#ifdef _PARALLEL_
25
26 /*intermediary*/
27 int i;
28 int analysis_type;
29 FILE* IOMODEL;
30 extern int my_rank;
31
32 /*Open input file on cpu 0: */
33 if(my_rank==0) IOMODEL= pfopen(inputfilename ,"rb");
34
35 /*Initialize internal data: */
36 this->nummodels=nummodels;
37 this->solution_type=in_solution_type;
38 this->analysis_counter=nummodels-1; //point to last analysis_type carried out.
39 this->results=new Results(); //not initialized by CreateDataSets
40
41 /*Dynamically allocate whatever is a list of length nummodels: */
42 analysis_type_list=(int*)xmalloc(nummodels*sizeof(int));
43
44 /*Initialize: */
45 for(i=0;i<nummodels;i++)analysis_type_list[i]=analyses[i];
46
47 /*create datasets for all analyses*/
48 ModelProcessorx(&this->elements,&this->nodes,&this->vertices,&this->materials,&this->constraints,&this->loads,&this->parameters,IOMODEL,this->solution_type,nummodels,analyses);
49
50 /*do the post-processing of the datasets to get an FemModel that can actually run analyses: */
51 for(i=0;i<nummodels;i++){
52
53 _printf_(VerboseMProcessor()," Processing finite element model of analysis %s:\n",EnumToStringx(analysis_type_list[i]));
54 analysis_type=analysis_type_list[i];
55 this->SetCurrentConfiguration(analysis_type);
56
57 if(i==0){
58 _printf_(VerboseMProcessor()," creating vertex degrees of freedom\n");
59 VerticesDofx(vertices,parameters); //only call once, we only have one set of vertices
60 }
61
62 _printf_(VerboseMProcessor()," resolving node constraints\n");
63 SpcNodesx(nodes,constraints,parameters,analysis_type);
64
65 _printf_(VerboseMProcessor()," creating nodal degrees of freedom\n");
66 NodesDofx(nodes,parameters,analysis_type);
67
68 _printf_(VerboseMProcessor()," configuring element and loads\n");
69 ConfigureObjectsx(elements, loads, nodes, vertices, materials,parameters);
70 }
71
72 /*Close input file descriptors: */
73 if(my_rank==0) pfclose(IOMODEL,inputfilename);
74
75 /*Add output file name to parameters: */
76 this->parameters->AddObject(new StringParam(OutputfilenameEnum,outputfilename));
77
78#endif
79
80}
81
82/*}}}1*/
83/*FUNCTION FemModel::destructor {{{1*/
84FemModel::~FemModel(){
85
86 /*Intermediary*/
87 int i;
88
89 /*Delete all the datasets: */
90 xfree((void**)&analysis_type_list);
91 delete elements;
92 delete nodes;
93 delete vertices;
94 delete constraints;
95 delete loads;
96 delete materials;
97 delete parameters;
98 delete results;
99
100}
101/*}}}1*/
102
103/*Object management*/
104/*FUNCTION FemModel::Echo {{{1*/
105void FemModel::Echo(void){
106
107 printf("FemModel echo: \n");
108 printf(" number of fem models: %i\n",nummodels);
109 printf(" analysis_type_list: \n");
110 for(int i=0;i<nummodels;i++)printf(" %i: %s\n",i,EnumToStringx(analysis_type_list[i]));
111 printf(" current analysis_type: \n");
112 printf(" %i: %s\n",analysis_counter,EnumToStringx(analysis_type_list[analysis_counter]));
113
114}
115/*}}}*/
116
117/*Numerics: */
118/*FUNCTION FemModel::SetCurrentConfiguration(int configuration_type,int analysis_type){{{1*/
119void FemModel::SetCurrentConfiguration(int configuration_type,int analysis_type){
120
121 /*Use configuration_type to setup the analysis counter, the configurations of objects etc ... but use
122 * analysis_type to drive the element numerics. This allows for use of 1 configuration_type for several
123 * analyses. For example: do a SurfaceSlopeX, SurfaceSlopeY, BedSlopeX and BedSlopeY analysis using the
124 * Slope configuration.*/
125
126 int found=-1;
127 for(int i=0;i<nummodels;i++){
128 if (analysis_type_list[i]==configuration_type){
129 found=i;
130 break;
131 }
132 }
133 if(found!=-1) analysis_counter=found;
134 else _error_("Could not find alias for analysis_type %s in list of FemModel analyses",EnumToStringx(configuration_type));
135
136 /*Now, plug analysis_counter and analysis_type inside the parameters: */
137 this->parameters->SetParam(analysis_counter,AnalysisCounterEnum);
138 this->parameters->SetParam(analysis_type,AnalysisTypeEnum);
139 this->parameters->SetParam(configuration_type,ConfigurationTypeEnum);
140
141 /*configure elements, loads and nodes, for this new analysis: */
142 this->elements->SetCurrentConfiguration(elements,loads, nodes,vertices, materials,parameters);
143 this->nodes->SetCurrentConfiguration(elements,loads, nodes,vertices, materials,parameters);
144 this->loads->SetCurrentConfiguration(elements, loads, nodes,vertices, materials,parameters);
145
146 /*take care of petsc options, that depend on this analysis type (present only after model processor)*/
147 if(this->parameters->Exist(PetscOptionsStringsEnum)){
148 PetscOptionsFromAnalysis(this->parameters,analysis_type);
149 _printf_(VerboseSolver()," petsc Options set for analysis type: %s\n",EnumToStringx(analysis_type));
150 }
151
152}
153/*}}}1*/
154/*FUNCTION FemModel::SetCurrentConfiguration(int configuration_type){{{1*/
155void FemModel::SetCurrentConfiguration(int configuration_type){
156
157 /*overload: analysis_type = configuration_type: */
158 this->SetCurrentConfiguration(configuration_type,configuration_type);
159}
160/*}}}1*/
Note: See TracBrowser for help on using the repository browser.