source: issm/trunk/src/c/solutions/prognostic2.cpp@ 4211

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

New Nodes class
File size: 3.5 KB
Line 
1/*!\file: prognostic2.cpp
2 * \brief: prognostic2 solution
3 */
4
5#ifdef HAVE_CONFIG_H
6 #include "config.h"
7#else
8#error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
9#endif
10
11#include "../objects/objects.h"
12#include "../shared/shared.h"
13#include "../DataSet/DataSet.h"
14#include "../EnumDefinitions/EnumDefinitions.h"
15#include "../include/include.h"
16#include "../modules/modules.h"
17#include "./solutions.h"
18
19int main(int argc,char* *argv){
20
21 /*I/O: */
22 FILE* fid=NULL;
23 char* inputfilename=NULL;
24 char* outputfilename=NULL;
25 char* lockname=NULL;
26 bool qmu_analysis=false;
27 bool waitonlock=false;
28
29 /*FemModel: */
30 FemModel* femmodel=NULL;
31
32 /*time*/
33 double start, finish;
34 double start_core, finish_core;
35 double start_init, finish_init;
36
37 int analyses[1]={Prognostic2AnalysisEnum};
38 int solution_type=Prognostic2AnalysisEnum;
39
40 MODULEBOOT();
41
42 #if !defined(_PARALLEL_) || (defined(_PARALLEL_) && !defined(_HAVE_PETSC_))
43 ISSMERROR(" parallel executable was compiled without support of parallel libraries!");
44 #endif
45
46 /*Initialize Petsc and get start time*/
47 PetscInitialize(&argc,&argv,(char *)0,"");
48 MPI_Barrier(MPI_COMM_WORLD); start=MPI_Wtime();
49
50 /*Size and rank: */
51 MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);
52 MPI_Comm_size(MPI_COMM_WORLD,&num_procs);
53
54 _printf_("recover input file name and output file name:\n");
55 inputfilename=argv[2];
56 outputfilename=argv[3];
57 lockname=argv[4];
58
59 /*Initialize model structure: */
60 MPI_Barrier(MPI_COMM_WORLD); start_init=MPI_Wtime();
61
62 /*Open handle to data on disk: */
63 fid=pfopen(inputfilename,"rb");
64
65 _printf_("create finite element model:\n");
66 femmodel=new FemModel(fid,solution_type,analyses,1);
67
68 /*add outputfilename in parameters: */
69 femmodel->parameters->AddObject(new StringParam(OutputFileNameEnum,outputfilename));
70
71 /*get parameters: */
72 femmodel->parameters->FindParam(&qmu_analysis,QmuAnalysisEnum);
73 femmodel->parameters->FindParam(&waitonlock,WaitOnLockEnum);
74
75 MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime();
76
77 /*are we running the solutoin sequence, or a qmu wrapper around it? : */
78 if(!qmu_analysis){
79
80 /*run prognostic2 analysis: */
81 _printf_("call computational core:\n");
82 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
83 prognostic2_core(femmodel);
84 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
85
86 _printf_("write results to disk:\n");
87 OutputResultsx(femmodel->elements, femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,femmodel->results);
88
89 }
90 else{
91
92 /*run qmu analysis: */
93 _printf_("calling qmu analysis on prognostic2 core:\n");
94
95 #ifdef _HAVE_DAKOTA_
96 MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( );
97 Qmux(femmodel,Prognostic2AnalysisEnum,NoneAnalysisEnum);
98 MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( );
99 #else
100 ISSMERROR(" Dakota not present, cannot do qmu!");
101 #endif
102 }
103
104 if (waitonlock>0){
105 _printf_("write lock file:\n");
106 WriteLockFile(lockname);
107 }
108
109 /*Free ressources:*/
110 delete femmodel;
111
112 /*Get finish time and close*/
113 MPI_Barrier(MPI_COMM_WORLD); finish = MPI_Wtime( );
114 _printf_("\n %-34s %f seconds \n","Model initialization elapsed time:",finish_init-start_init);
115 _printf_(" %-34s %f seconds \n","Core solution elapsed time:",finish_core-start_core);
116 _printf_("\n %s %i hrs %i min %i sec\n\n","Total elapsed time:",int((finish-start)/3600),int(int(finish-start)%3600/60),int(finish-start)%60);
117 _printf_("closing MPI and Petsc\n");
118 PetscFinalize();
119
120 /*end module: */
121 MODULEEND();
122
123 return 0; //unix success return;
124}
Note: See TracBrowser for help on using the repository browser.