Changeset 4453
- Timestamp:
- 07/07/10 17:35:32 (15 years ago)
- Location:
- issm/trunk/src
- Files:
-
- 58 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/modules/InputUpdateFromDakotax/InputUpdateFromDakotax.cpp
r4439 r4453 9 9 #include "../../EnumDefinitions/EnumDefinitions.h" 10 10 11 void InputUpdateFromDakotax(FemModel* femmodel, double* variables,char**variables_descriptors,int numvariables){11 void InputUpdateFromDakotax(Elements* elements,Nodes* nodes,Vertices* vertices,Loads* loads,Materials* materials,Parameters* parameters,double* variables,char* *variables_descriptors,int numvariables){ 12 12 13 13 int verbose; … … 17 17 18 18 /*retrieve parameters: */ 19 femmodel->parameters->FindParam(&verbose,VerboseEnum);20 femmodel->parameters->FindParam(&qmu_npart,QmuNPartEnum);21 femmodel->parameters->FindParam(&qmu_part,&dummy,QmuPartEnum);19 parameters->FindParam(&verbose,VerboseEnum); 20 parameters->FindParam(&qmu_npart,QmuNPartEnum); 21 parameters->FindParam(&qmu_part,&dummy,QmuPartEnum); 22 22 23 23 ISSMERROR(" not supported yet!"); -
issm/trunk/src/c/modules/InputUpdateFromDakotax/InputUpdateFromDakotax.h
r4236 r4453 9 9 #include "../../Container/Container.h" 10 10 11 void InputUpdateFromDakotax(FemModel* femmodel,double* variables,char**variables_descriptors,int numvariables);11 void InputUpdateFromDakotax(Elements* elements,Nodes* nodes,Vertices* vertices,Loads* loads,Materials* materials,Parameters* parameters,double* variables,char* *variables_descriptors,int numvariables); 12 12 13 13 #endif /* _INPUTUPDATEFROMDAKOTAXX_H */ -
issm/trunk/src/c/modules/Qmux/SpawnCoreParallel.cpp
r4439 r4453 61 61 62 62 /*Modify core inputs in objects contained in femmodel, to reflect the dakota variables inputs: */ 63 InputUpdateFromDakotax(femmodel->elements,femmodel->nodes,femmodel->vertices, loads,femmodel->materials,femmodel->parameters,variables,variables_descriptors,numvariables);63 InputUpdateFromDakotax(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,variables,variables_descriptors,numvariables); 64 64 65 65 /*Run the core solution sequence: */ 66 66 if(verbose)_printf_("%s%s%s\n","Starting ",EnumAsString(solution_type)," core:"); 67 SolutionConfiguration(NULL,NULL,&solutioncore ); solution_core(femmodel);67 SolutionConfiguration(NULL,NULL,&solutioncore,solution_type); solutioncore(femmodel); 68 68 69 69 /*compute responses: */ 70 70 if(verbose)_printf_("compute dakota responses:\n"); 71 DakotaResponsesx(responses,femmodel->elements,femmodel->nodes,femmodel->vertices, loads,femmodel->materials,femmodel->parameters,responses_descriptors,numresponses,);71 DakotaResponsesx(responses,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,responses_descriptors,numresponses); 72 72 73 73 /*Free ressources:*/ -
issm/trunk/src/c/modules/modules.h
r4168 r4453 82 82 #include "./GetVectorFromInputsx/GetVectorFromInputsx.h" 83 83 #include "./InputConvergencex/InputConvergencex.h" 84 #include "./DakotaResponsesx/DakotaResponsesx.h" 84 85 85 86 #endif -
issm/trunk/src/mex/AverageFilter/AverageFilter.cpp
r2775 r4453 17 17 { 18 18 19 int i,j ,noerr;19 int i,j; 20 20 21 21 /* required input: */ -
issm/trunk/src/mex/Bamg/Bamg.cpp
r3579 r4453 11 11 12 12 /*diverse: */ 13 int noerr=1;14 13 int i; 15 14 int lines,cols; -
issm/trunk/src/mex/BuildNodeSets/BuildNodeSets.cpp
r4211 r4453 7 7 8 8 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 9 10 /*diverse: */11 int noerr=1;12 9 13 10 /*input datasets: */ -
issm/trunk/src/mex/ComputeBasalStress/ComputeBasalStress.cpp
r4404 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/ComputePressure/ComputePressure.cpp
r4404 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/ConfigureObjects/ConfigureObjects.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/ControlOptimization/ControlOptimization.cpp
r3841 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/CostFunction/CostFunction.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/DakotaResponses/DakotaResponses.cpp
r4449 r4453 1 /*\file InputUpdateFromDakota.c2 *\brief: update inputs using dakota input variables1 /*\file DakotaResponses.c 2 *\brief: compute dakota responses after a run 3 3 */ 4 4 5 #include "./ InputUpdateFromDakota.h"5 #include "./DakotaResponses.h" 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr=1; 9 int i; 11 10 12 11 /*input datasets: */ … … 22 21 int numresponses; 23 22 mxArray* pfield=NULL; 23 int stringlength; 24 24 25 25 /*Boot module: */ … … 27 27 28 28 /*checks on arguments on the matlab side: */ 29 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,& InputUpdateFromDakotaUsage);29 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&DakotaResponsesUsage); 30 30 31 31 /*Input datasets: */ … … 75 75 } 76 76 77 void InputUpdateFromDakotaUsage(void)77 void DakotaResponsesUsage(void) 78 78 { 79 79 _printf_("\n"); 80 _printf_(" usage: [ elements loads] = %s(elements,nodes,vertices,loads,materials,parameters,variables,variables_descriptors);\n",__FUNCT__);80 _printf_(" usage: [responses] = %s(elements,nodes,vertices,loads,materials,parameters,responses_descritpors,numresponses);\n",__FUNCT__); 81 81 _printf_("\n"); 82 82 } -
issm/trunk/src/mex/Du/Du.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Echo/Echo.cpp
r2333 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/FieldAverageOntoVertices/FieldAverageOntoVertices.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/FieldDepthAverage/FieldDepthAverage.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/GetSolutionFromInputs/GetSolutionFromInputs.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Gradj/Gradj.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr = 1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/HoleFiller/HoleFiller.cpp
r2333 r4453 17 17 { 18 18 19 int i,j ,noerr;19 int i,j; 20 20 21 21 /* required input: */ -
issm/trunk/src/mex/InputControlConstrain/InputControlConstrain.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr = 1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/InputDepthAverage/InputDepthAverage.cpp
r4218 r4453 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr = 1; 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 enum_type; 17 int average_enum_type; 11 18 12 /*input datasets: */ 13 Elements* elements = NULL; 14 Nodes *nodes = NULL; 15 Vertices *vertices = NULL; 16 Loads *loads = NULL; 17 Materials *materials = NULL; 18 Parameters *parameters = NULL; 19 int enum_type; 20 int average_enum_type; 19 /*Boot module: */ 20 MODULEBOOT(); 21 21 22 /*Boot module: */23 MODULEBOOT();22 /*checks on arguments on the matlab side: */ 23 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputDepthAverageUsage); 24 24 25 /*checks on arguments on the matlab side: */ 26 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputDepthAverageUsage); 25 /*Input datasets: */ 26 FetchData((DataSet**)&elements,ELEMENTSIN); 27 FetchData((DataSet**)&nodes,NODESIN); 28 FetchData((DataSet**)&vertices,VERTICESIN); 29 FetchData((DataSet**)&loads,LOADSIN); 30 FetchData((DataSet**)&materials,MATERIALSIN); 31 FetchParams(¶meters,PARAMETERSIN); 32 FetchData(&enum_type,ENUMTYPE); 33 FetchData(&average_enum_type,AVERAGEENUMTYPE); 27 34 28 /*Input datasets: */ 29 FetchData((DataSet**)&elements,ELEMENTSIN); 30 FetchData((DataSet**)&nodes,NODESIN); 31 FetchData((DataSet**)&vertices,VERTICESIN); 32 FetchData((DataSet**)&loads,LOADSIN); 33 FetchData((DataSet**)&materials,MATERIALSIN); 34 FetchParams(¶meters,PARAMETERSIN); 35 FetchData(&enum_type,ENUMTYPE); 36 FetchData(&average_enum_type,AVERAGEENUMTYPE); 35 /*call "x" code layer*/ 36 InputDepthAveragex(elements,nodes,vertices,loads, materials,parameters,enum_type,average_enum_type); 37 37 38 /*call "x" code layer*/ 39 InputDepthAveragex(elements,nodes,vertices,loads, materials,parameters,enum_type,average_enum_type); 38 /*write output datasets: */ 39 WriteData(ELEMENTS,elements); 40 WriteData(NODES,nodes); 41 WriteData(VERTICES,vertices); 42 WriteData(LOADS,loads); 43 WriteData(MATERIALS,materials); 44 WriteParams(PARAMETERS,parameters); 40 45 41 /*write output datasets: */42 WriteData(ELEMENTS,elements);43 WriteData(NODES,nodes);44 WriteData(VERTICES,vertices);45 WriteData(LOADS,loads);46 WriteData(MATERIALS,materials);47 WriteParams(PARAMETERS,parameters);46 /*Free ressources: */ 47 delete elements; 48 delete nodes; 49 delete vertices; 50 delete loads; 51 delete materials; 52 delete parameters; 48 53 49 /*Free ressources: */ 50 delete elements; 51 delete nodes; 52 delete vertices; 53 delete loads; 54 delete materials; 55 delete parameters; 56 57 /*end module: */ 58 MODULEEND(); 54 /*end module: */ 55 MODULEEND(); 59 56 } 60 57 -
issm/trunk/src/mex/InputDuplicate/InputDuplicate.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/InputExtrude/InputExtrude.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/InputScale/InputScale.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/InputToResult/InputToResult.cpp
r4228 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/InputUpdateFromConstant/InputUpdateFromConstant.cpp
r4218 r4453 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr=1; 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 double constant; 17 int name; 11 18 12 /*input datasets: */ 13 Elements* elements=NULL; 14 Nodes* nodes=NULL; 15 Vertices* vertices=NULL; 16 Loads* loads=NULL; 17 Materials* materials=NULL; 18 Parameters* parameters=NULL; 19 double constant; 20 int name; 19 /*Boot module: */ 20 MODULEBOOT(); 21 21 22 /*Boot module: */23 MODULEBOOT();22 /*checks on arguments on the matlab side: */ 23 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputUpdateFromConstantUsage); 24 24 25 /*checks on arguments on the matlab side: */ 26 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputUpdateFromConstantUsage); 25 /*Input datasets: */ 26 FetchData((DataSet**)&elements,ELEMENTSIN); 27 FetchData((DataSet**)&nodes,NODESIN); 28 FetchData((DataSet**)&vertices,VERTICESIN); 29 FetchData((DataSet**)&loads,LOADSIN); 30 FetchData((DataSet**)&materials,MATERIALSIN); 31 FetchParams(¶meters,PARAMETERSIN); 32 FetchData(&constant,CONSTANT); 33 FetchData(&name,NAME); 27 34 28 /*Input datasets: */ 29 FetchData((DataSet**)&elements,ELEMENTSIN); 30 FetchData((DataSet**)&nodes,NODESIN); 31 FetchData((DataSet**)&vertices,VERTICESIN); 32 FetchData((DataSet**)&loads,LOADSIN); 33 FetchData((DataSet**)&materials,MATERIALSIN); 34 FetchParams(¶meters,PARAMETERSIN); 35 FetchData(&constant,CONSTANT); 36 FetchData(&name,NAME); 35 /*!Generate internal degree of freedom numbers: */ 36 InputUpdateFromConstantx(elements,nodes,vertices,loads, materials,parameters,constant,name); 37 37 38 /*!Generate internal degree of freedom numbers: */ 39 InputUpdateFromConstantx(elements,nodes,vertices,loads, materials,parameters,constant,name); 38 /*write output datasets: */ 39 WriteData(ELEMENTS,elements); 40 WriteData(LOADS,loads); 40 41 41 /*write output datasets: */ 42 WriteData(ELEMENTS,elements); 43 WriteData(LOADS,loads); 42 /*Free ressources: */ 43 delete elements; 44 delete nodes; 45 delete vertices; 46 delete loads; 47 delete materials; 48 delete parameters; 44 49 45 /*Free ressources: */ 46 delete elements; 47 delete nodes; 48 delete vertices; 49 delete loads; 50 delete materials; 51 delete parameters; 52 53 /*end module: */ 54 MODULEEND(); 50 /*end module: */ 51 MODULEEND(); 55 52 } 56 53 -
issm/trunk/src/mex/InputUpdateFromDakota/InputUpdateFromDakota.cpp
r4449 r4453 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr=1; 9 int i; 11 10 12 11 /*input datasets: */ … … 20 19 char** variables_descriptors=NULL; 21 20 char* string=NULL; 21 int stringlength; 22 22 int numvariables; 23 23 mxArray* pfield=NULL; … … 38 38 39 39 /*dakota input: */ 40 FetchData(&variables,&numvariables,VARIABLES IN);40 FetchData(&variables,&numvariables,VARIABLES); 41 41 42 42 variables_descriptors=(char**)xmalloc(numvariables*sizeof(char*)); -
issm/trunk/src/mex/InputUpdateFromSolution/InputUpdateFromSolution.cpp
r4218 r4453 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr=1; 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 Vec solution=NULL; 11 17 12 /*input datasets: */ 13 Elements* elements=NULL; 14 Nodes* nodes=NULL; 15 Vertices* vertices=NULL; 16 Loads* loads=NULL; 17 Materials* materials=NULL; 18 Parameters* parameters=NULL; 19 Vec solution=NULL; 18 /*Boot module: */ 19 MODULEBOOT(); 20 20 21 /*Boot module: */22 MODULEBOOT();21 /*checks on arguments on the matlab side: */ 22 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputUpdateFromSolutionUsage); 23 23 24 /*checks on arguments on the matlab side: */ 25 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputUpdateFromSolutionUsage); 24 /*Input datasets: */ 25 FetchData((DataSet**)&elements,ELEMENTSIN); 26 FetchData((DataSet**)&nodes,NODESIN); 27 FetchData((DataSet**)&vertices,VERTICESIN); 28 FetchData((DataSet**)&loads,LOADSIN); 29 FetchData((DataSet**)&materials,MATERIALSIN); 30 FetchParams(¶meters,PARAMETERSIN); 31 FetchData(&solution,SOLUTION); 26 32 27 /*Input datasets: */ 28 FetchData((DataSet**)&elements,ELEMENTSIN); 29 FetchData((DataSet**)&nodes,NODESIN); 30 FetchData((DataSet**)&vertices,VERTICESIN); 31 FetchData((DataSet**)&loads,LOADSIN); 32 FetchData((DataSet**)&materials,MATERIALSIN); 33 FetchParams(¶meters,PARAMETERSIN); 34 FetchData(&solution,SOLUTION); 33 /*!Generate internal degree of freedom numbers: */ 34 InputUpdateFromSolutionx(elements,nodes,vertices,loads, materials,parameters,solution); 35 35 36 /*!Generate internal degree of freedom numbers: */37 InputUpdateFromSolutionx(elements,nodes,vertices,loads, materials,parameters,solution);36 /*write output datasets: */ 37 WriteData(ELEMENTS,elements); 38 38 39 /*write output datasets: */ 40 WriteData(ELEMENTS,elements); 39 /*Free ressources: */ 40 delete elements; 41 delete nodes; 42 delete vertices; 43 delete loads; 44 delete materials; 45 delete parameters; 46 VecFree(&solution); 41 47 42 /*Free ressources: */ 43 delete elements; 44 delete nodes; 45 delete vertices; 46 delete loads; 47 delete materials; 48 delete parameters; 49 VecFree(&solution); 50 51 /*end module: */ 52 MODULEEND(); 48 /*end module: */ 49 MODULEEND(); 53 50 } 54 51 -
issm/trunk/src/mex/InputUpdateFromVector/InputUpdateFromVector.cpp
r4218 r4453 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr=1; 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 double* vector=NULL; 17 int dummy; 18 int NameEnum; 19 int TypeEnum; 11 20 12 /*input datasets: */ 13 Elements* elements=NULL; 14 Nodes* nodes=NULL; 15 Vertices* vertices=NULL; 16 Loads* loads=NULL; 17 Materials* materials=NULL; 18 Parameters* parameters=NULL; 19 double* vector=NULL; 20 int dummy; 21 int NameEnum; 22 int TypeEnum; 21 /*Boot module: */ 22 MODULEBOOT(); 23 23 24 /*Boot module: */25 MODULEBOOT();24 /*checks on arguments on the matlab side: */ 25 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputUpdateFromVectorUsage); 26 26 27 /*checks on arguments on the matlab side: */ 28 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&InputUpdateFromVectorUsage); 27 /*Input datasets: */ 28 FetchData((DataSet**)&elements,ELEMENTSIN); 29 FetchData((DataSet**)&nodes,NODESIN); 30 FetchData((DataSet**)&vertices,VERTICESIN); 31 FetchData((DataSet**)&loads,LOADSIN); 32 FetchData((DataSet**)&materials,MATERIALSIN); 33 FetchParams(¶meters,PARAMETERSIN); 34 FetchData(&vector,&dummy,VECTOR); 35 FetchData(&NameEnum,NAME); 36 FetchData(&TypeEnum,TYPE); 29 37 30 /*Input datasets: */ 31 FetchData((DataSet**)&elements,ELEMENTSIN); 32 FetchData((DataSet**)&nodes,NODESIN); 33 FetchData((DataSet**)&vertices,VERTICESIN); 34 FetchData((DataSet**)&loads,LOADSIN); 35 FetchData((DataSet**)&materials,MATERIALSIN); 36 FetchParams(¶meters,PARAMETERSIN); 37 FetchData(&vector,&dummy,VECTOR); 38 FetchData(&NameEnum,NAME); 39 FetchData(&TypeEnum,TYPE); 38 /*Check that type is one of Constant, Vertex or Element: */ 39 if ((TypeEnum!=ConstantEnum) && (TypeEnum!=VertexEnum) && (TypeEnum!=ElementEnum)){ 40 ISSMERROR("%s%s%s\n","Type of input can only be a constant, a vertex or an element vector. Right now, you input a ",EnumAsString(TypeEnum)," type data for input update"); 41 } 40 42 41 /*Check that type is one of Constant, Vertex or Element: */ 42 if ((TypeEnum!=ConstantEnum) && (TypeEnum!=VertexEnum) && (TypeEnum!=ElementEnum)){ 43 ISSMERROR("%s%s%s\n","Type of input can only be a constant, a vertex or an element vector. Right now, you input a ",EnumAsString(TypeEnum)," type data for input update"); 44 } 43 /*!Generate internal degree of freedom numbers: */ 44 InputUpdateFromVectorx(elements,nodes,vertices,loads, materials,parameters,vector,NameEnum,TypeEnum); 45 45 46 /*!Generate internal degree of freedom numbers: */ 47 InputUpdateFromVectorx(elements,nodes,vertices,loads, materials,parameters,vector,NameEnum,TypeEnum); 46 /*write output datasets: */ 47 WriteData(ELEMENTS,elements); 48 WriteData(NODES,nodes); 49 WriteData(VERTICES,vertices); 50 WriteData(LOADS,loads); 51 WriteData(MATERIALS,materials); 52 WriteParams(PARAMETERS,parameters); 48 53 49 /*write output datasets: */ 50 WriteData(ELEMENTS,elements); 51 WriteData(NODES,nodes); 52 WriteData(VERTICES,vertices); 53 WriteData(LOADS,loads); 54 WriteData(MATERIALS,materials); 55 WriteParams(PARAMETERS,parameters); 54 /*Free ressources: */ 55 delete elements; 56 delete nodes; 57 delete vertices; 58 delete loads; 59 delete materials; 60 delete parameters; 61 xfree((void**)&vector); 56 62 57 /*Free ressources: */ 58 delete elements; 59 delete nodes; 60 delete vertices; 61 delete loads; 62 delete materials; 63 delete parameters; 64 xfree((void**)&vector); 65 66 /*end module: */ 67 MODULEEND(); 63 /*end module: */ 64 MODULEEND(); 68 65 } 69 66 -
issm/trunk/src/mex/MassFlux/MassFlux.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Mergesolutionfromftog/Mergesolutionfromftog.cpp
r3960 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Misfit/Misfit.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/MpcNodes/MpcNodes.cpp
r4217 r4453 7 7 8 8 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 9 10 /*diverse: */11 int noerr=1;12 9 13 10 /*input datasets: */ -
issm/trunk/src/mex/NormalizeConstraints/NormalizeConstraints.cpp
r2333 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Orth/Orth.cpp
r2333 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/OutputResults/OutputResults.cpp
r4306 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/OutputRifts/OutputRifts.cpp
r4214 r4453 8 8 9 9 /*diverse: */ 10 int noerr=1;11 10 int numrifts; 12 11 -
issm/trunk/src/mex/ParameterOutput/ParameterOutput.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/PenaltyConstraints/PenaltyConstraints.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/PenaltySystemMatrices/PenaltySystemMatrices.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*output: */ -
issm/trunk/src/mex/PropagateFlagsFromConnectivity/PropagateFlagsFromConnectivity.cpp
r2888 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 12 8 13 9 /*input/output datasets: */ -
issm/trunk/src/mex/Qmu/Qmu.cpp
r4211 r4453 12 12 13 13 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 14 15 /*diverse: */16 int noerr = 1;17 14 18 15 /*input datasets: */ -
issm/trunk/src/mex/Reduceloadfromgtof/Reduceloadfromgtof.cpp
r3960 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Reducematrixfromgtof/Reducematrixfromgtof.cpp
r2333 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Reducevectorgtof/Reducevectorgtof.cpp
r2333 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Reducevectorgtos/Reducevectorgtos.cpp
r3960 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Solver/Solver.cpp
r3715 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/SpcNodes/SpcNodes.cpp
r4217 r4453 7 7 8 8 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 9 10 /*diverse: */11 int noerr=1;12 9 13 10 /*input datasets: */ -
issm/trunk/src/mex/SystemMatrices/SystemMatrices.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/Test/Test.cpp.grid
r2994 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*Boot module: */ -
issm/trunk/src/mex/TriMesh/TriMesh.cpp
r3755 r4453 48 48 double** riftsgridsy=NULL; 49 49 50 /* error handling: */51 int noerr=1;52 53 50 /* Triangle structures: */ 54 51 struct triangulateio in,out; … … 99 96 100 97 /*Start reading the domain outline file: */ 101 noerr=DomainOutlineRead(&nprof,&profngrids,&pprofx,&pprofy,domainname); 102 if(!noerr){ 98 if(!DomainOutlineRead(&nprof,&profngrids,&pprofx,&pprofy,domainname)){ 103 99 printf("%s%s%s\n",__FUNCT__," error message reading domain outline ",domainname); 104 100 mexErrMsgTxt(" "); … … 107 103 /*Read rifts file if present: */ 108 104 if(riftname){ 109 noerr=DomainOutlineRead(&numrifts,&riftsnumgrids,&riftsgridsx,&riftsgridsy,riftname); 110 if(!noerr){ 105 if(!DomainOutlineRead(&numrifts,&riftsnumgrids,&riftsgridsx,&riftsgridsy,riftname)){ 111 106 printf("%s%s%s\n",__FUNCT__," error message reading rifts outline ",riftname); 112 107 mexErrMsgTxt(" "); … … 268 263 269 264 /*Associate elements with segments: */ 270 noerr=AssociateSegmentToElement(&segments,out.numberofsegments,index,out.numberoftriangles);265 AssociateSegmentToElement(&segments,out.numberofsegments,index,out.numberoftriangles); 271 266 272 267 /*Order segments so that their normals point outside the domain: */ -
issm/trunk/src/mex/TriMeshNoDensity/TriMeshNoDensity.cpp
r3755 r4453 46 46 double** riftsgridsy=NULL; 47 47 48 /* error handling: */49 int noerr=1;50 51 48 /* Triangle structures: */ 52 49 struct triangulateio in,out; … … 89 86 90 87 /*Start reading the domain outline file: */ 91 noerr=DomainOutlineRead(&nprof,&profngrids,&pprofx,&pprofy,domainname); 92 if(!noerr){ 88 if(!DomainOutlineRead(&nprof,&profngrids,&pprofx,&pprofy,domainname)){ 93 89 printf("%s%s%s\n",__FUNCT__," error message reading domain outline ",domainname); 94 90 mexErrMsgTxt(" "); … … 97 93 /*Read rifts file if present: */ 98 94 if(riftname){ 99 noerr=DomainOutlineRead(&numrifts,&riftsnumgrids,&riftsgridsx,&riftsgridsy,riftname); 100 if(!noerr){ 95 if(!DomainOutlineRead(&numrifts,&riftsnumgrids,&riftsgridsx,&riftsgridsy,riftname)){ 101 96 printf("%s%s%s\n",__FUNCT__," error message reading rifts outline ",riftname); 102 97 mexErrMsgTxt(" "); … … 259 254 260 255 /*Associate elements with segments: */ 261 noerr=AssociateSegmentToElement(&segments,out.numberofsegments,index,out.numberoftriangles);256 AssociateSegmentToElement(&segments,out.numberofsegments,index,out.numberoftriangles); 262 257 263 258 /*Order segments so that their normals point outside the domain: */ -
issm/trunk/src/mex/TriMeshProcessRifts/TriMeshProcessRifts.cpp
r3755 r4453 53 53 int numrifts; 54 54 55 /* error handling: */56 int noerr=1;57 58 55 /* verify correct usage: */ 59 56 if (nlhs==0 && nrhs==0) { … … 176 173 177 174 /*Order segments so that their normals point outside the domain: */ 178 noerr=OrderSegments(&segments_in,num_seg, index_in,nel);175 OrderSegments(&segments_in,num_seg, index_in,nel); 179 176 180 177 -
issm/trunk/src/mex/TriMeshRefine/TriMeshRefine.cpp
r3755 r4453 43 43 char* order=NULL; 44 44 45 /* error handling: */46 int noerr=1;47 48 45 /* Triangle structures: */ 49 46 struct triangulateio in,out; … … 276 273 277 274 /*Associate elements with segments: */ 278 noerr=AssociateSegmentToElement(&segments,out.numberofsegments,index,out.numberoftriangles);275 AssociateSegmentToElement(&segments,out.numberofsegments,index,out.numberoftriangles); 279 276 280 277 /*Order segments so that their normals point outside the domain: */ 281 278 if(!strcmp(order,"yes")){ 282 noerr=OrderSegments(&segments,out.numberofsegments, index,out.numberoftriangles);279 OrderSegments(&segments,out.numberofsegments, index,out.numberoftriangles); 283 280 } 284 281 -
issm/trunk/src/mex/UpdateGeometry/UpdateGeometry.cpp
r4218 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr = 1;11 8 12 9 /*input datasets: */ -
issm/trunk/src/mex/UpdateVertexPositions/UpdateVertexPositions.cpp
r4218 r4453 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 8 9 /*diverse: */ 10 int noerr = 1; 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; 11 16 12 /*input datasets: */ 13 Elements* elements = NULL; 14 Nodes *nodes = NULL; 15 Vertices *vertices = NULL; 16 Loads *loads = NULL; 17 Materials *materials = NULL; 18 Parameters *parameters = NULL; 17 /*Boot module: */ 18 MODULEBOOT(); 19 19 20 /*Boot module: */21 MODULEBOOT();20 /*checks on arguments on the matlab side: */ 21 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&UpdateVertexPositionsUsage); 22 22 23 /*checks on arguments on the matlab side: */ 24 CheckNumMatlabArguments(nlhs,NLHS,nrhs,NRHS,__FUNCT__,&UpdateVertexPositionsUsage); 23 /*Input datasets: */ 24 FetchData((DataSet**)&elements,ELEMENTSIN); 25 FetchData((DataSet**)&nodes,NODESIN); 26 FetchData((DataSet**)&vertices,VERTICESIN); 27 FetchData((DataSet**)&loads,LOADSIN); 28 FetchData((DataSet**)&materials,MATERIALSIN); 29 FetchParams(¶meters,PARAMETERSIN); 25 30 26 /*Input datasets: */ 27 FetchData((DataSet**)&elements,ELEMENTSIN); 28 FetchData((DataSet**)&nodes,NODESIN); 29 FetchData((DataSet**)&vertices,VERTICESIN); 30 FetchData((DataSet**)&loads,LOADSIN); 31 FetchData((DataSet**)&materials,MATERIALSIN); 32 FetchParams(¶meters,PARAMETERSIN); 31 /*call "x" code layer*/ 32 UpdateVertexPositionsx(elements,nodes,vertices,loads, materials,parameters); 33 33 34 /*call "x" code layer*/ 35 UpdateVertexPositionsx(elements,nodes,vertices,loads, materials,parameters); 34 /*write output datasets: */ 35 WriteData(ELEMENTS,elements); 36 WriteData(NODES,nodes); 37 WriteData(VERTICES,vertices); 38 WriteData(LOADS,loads); 39 WriteData(MATERIALS,materials); 40 WriteParams(PARAMETERS,parameters); 36 41 37 /*write output datasets: */38 WriteData(ELEMENTS,elements);39 WriteData(NODES,nodes);40 WriteData(VERTICES,vertices);41 WriteData(LOADS,loads);42 WriteData(MATERIALS,materials);43 WriteParams(PARAMETERS,parameters);42 /*Free ressources: */ 43 delete elements; 44 delete nodes; 45 delete vertices; 46 delete loads; 47 delete materials; 48 delete parameters; 44 49 45 /*Free ressources: */ 46 delete elements; 47 delete nodes; 48 delete vertices; 49 delete loads; 50 delete materials; 51 delete parameters; 52 53 /*end module: */ 54 MODULEEND(); 50 /*end module: */ 51 MODULEEND(); 55 52 } 56 53 -
issm/trunk/src/mex/VerticesDof/VerticesDof.cpp
r4213 r4453 6 6 7 7 void mexFunction( int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]){ 8 9 /*diverse: */10 int noerr=1;11 8 12 9 /*input datasets: */
Note:
See TracChangeset
for help on using the changeset viewer.