source: issm/trunk/src/m/solutions/surfaceslope.m@ 4294

Last change on this file since 4294 was 4294, checked in by Mathieu Morlighem, 15 years ago

Homogenized with parallel, added missing timing and fixed diagnostic

File size: 1.0 KB
RevLine 
[4113]1function md=surfaceslopecompute(md);
2%SLOPECOMPUTE - compute the surface slope of a model
[4102]3%
4% Usage:
[4113]5% md=surfaceslope(md)
[4102]6%
7 %timing
8 t1=clock;
9
[4294]10 numanalyses=1;
11 analyses=[SlopeAnalysisEnum];
12 solution_type=SlopeSolutionEnum;
[4102]13
14 displaystring(md.verbose,'%s',['create finite element model']);
[4294]15 femmodel=NewFemModel(md,solution_type,analyses,1);
[4102]16
17 %retrieve parameters
18 verbose=femmodel.parameters.Verbose;
19 qmu_analysis=femmodel.parameters.QmuAnalysis;
20
21 %compute solution
22 if ~qmu_analysis,
23 displaystring(verbose,'%s',['call computational core']);
[4113]24 femmodel=surfaceslope_core(femmodel);
25
26 displaystring(verbose,'%s',['write results']);
[4184]27 md.results.surfaceslopecompute=OutputResults(femmodel.elements, femmodel.nodes , femmodel.vertices , femmodel.loads , femmodel.materials, femmodel.parameters, femmodel.results);
[4113]28
[4102]29 else
30 %launch dakota driver for diagnostic core solution
[4193]31 Qmu(femmodel,femmodel.parameters);
[4102]32 end
33
34 %stop timing
35 t2=clock;
36 displaystring(md.verbose,'\n%s\n',['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
Note: See TracBrowser for help on using the repository browser.