1 | function md=diagnostic(md);
|
---|
2 | %DIAGNOSTIC - compute the velocity field of a model
|
---|
3 | %
|
---|
4 | % Usage:
|
---|
5 | % md=diagnostic(md)
|
---|
6 | %
|
---|
7 | %timing
|
---|
8 | t1=clock;
|
---|
9 |
|
---|
10 | analysis_types=[DiagnosticHorizAnalysisEnum,DiagnosticVertAnalysisEnum,DiagnosticStokesAnalysisEnum,DiagnosticHutterAnalysisEnum,SlopeAnalysisEnum];
|
---|
11 | solution_type=DiagnosticAnalysisEnum;
|
---|
12 |
|
---|
13 | displaystring(md.verbose,'%s',['create finite element model']);
|
---|
14 | femmodel=NewFemModel(md,solution_type,analysis_types);
|
---|
15 |
|
---|
16 | %retrieve parameters
|
---|
17 | verbose=femmodel.parameters.Verbose;
|
---|
18 | qmu_analysis=femmodel.parameters.QmuAnalysis;
|
---|
19 | control_analysis=femmodel.parameters.ControlAnalysis;
|
---|
20 |
|
---|
21 | %compute solution
|
---|
22 | if ~qmu_analysis,
|
---|
23 | if ~control_analysis,
|
---|
24 |
|
---|
25 | displaystring(verbose,'%s',['call computational core']);
|
---|
26 | femmodel=diagnostic_core(femmodel);
|
---|
27 |
|
---|
28 | else,
|
---|
29 |
|
---|
30 | displaystring(verbose,'%s',['call computational core']);
|
---|
31 | femmodel=control_core(femmodel);
|
---|
32 |
|
---|
33 | end
|
---|
34 |
|
---|
35 | displaystring(verbose,'%s',['write results']);
|
---|
36 | OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
|
---|
37 |
|
---|
38 | else
|
---|
39 | %launch dakota driver for diagnostic core solution
|
---|
40 | Qmu(femmodel);
|
---|
41 | end
|
---|
42 |
|
---|
43 | %stop timing
|
---|
44 | t2=clock;
|
---|
45 | displaystring(md.verbose,'\n%s\n',['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
|
---|