source: issm/trunk/src/mex/ComputeBasalStress/ComputeBasalStress.cpp@ 8910

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

Reorganized io/ directory quite a bit.
Separated Matlab io from the other Disk related io.
Renamed WriteData to WriteMatlabData and FetchData to FetchMatlabData.
Folded WriteParams and FetchParams into FetchMatlabData and WriteMatlabData.

File size: 1.7 KB
Line 
1/*\file ComputeBasalStress.c
2 *\brief: recover pressure from elements
3 */
4
5#include "./ComputeBasalStress.h"
6
7void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){
8
9 /*input datasets: */
10 Elements* elements=NULL;
11 Nodes* nodes=NULL;
12 Vertices* vertices=NULL;
13 Loads* loads=NULL;
14 Materials* materials=NULL;
15 Parameters* parameters=NULL;
16 int numberofnodes;
17
18 /* output datasets: */
19 Vec sigma_g=NULL;
20
21 /*Boot module: */
22 MODULEBOOT();
23
24 /*checks on arguments on the matlab side: */
25 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&ComputeBasalStressUsage);
26
27 /*Input datasets: */
28 FetchMatlabData((DataSet**)&elements,ELEMENTS);
29 FetchMatlabData((DataSet**)&nodes,NODES);
30 FetchMatlabData((DataSet**)&vertices,VERTICES);
31 FetchMatlabData((DataSet**)&loads,LOADS);
32 FetchMatlabData((DataSet**)&materials,MATERIALS);
33 FetchMatlabData(&parameters,PARAMETERS);
34
35 /*configure: */
36 elements-> Configure(elements,loads, nodes,vertices, materials,parameters);
37 nodes-> Configure(elements,loads, nodes,vertices, materials,parameters);
38 loads-> Configure(elements, loads, nodes,vertices, materials,parameters);
39
40 /*!Generate internal degree of freedom numbers: */
41 ComputeBasalStressx(&sigma_g, elements,nodes,vertices,loads,materials,parameters);
42
43 /*write output datasets: */
44 WriteMatlabData(SIGMA,sigma_g);
45
46 /*Free ressources: */
47 delete nodes;
48 delete vertices;
49 delete elements;
50 delete materials;
51 delete loads;
52 delete parameters;
53 VecFree(&sigma_g);
54
55 /*end module: */
56 MODULEEND();
57}
58
59void ComputeBasalStressUsage(void) {
60 printf("\n");
61 printf(" usage: [p_g] = %s(elements, nodes, vertices, loads, materials, params);\n",__FUNCT__);
62 printf("\n");
63}
Note: See TracBrowser for help on using the repository browser.