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