Changeset 4193 for issm/trunk/src/m/solutions/SpawnCore.m
- Timestamp:
- 06/24/10 13:46:42 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/SpawnCore.m
r4131 r4193 1 function responses=SpawnCore( models,variables,variabledescriptors,counter);1 function responses=SpawnCore(femmodel,variables,variabledescriptors,counter); 2 2 %SPAWNCORE - for Qmu analysis, using Dakota. Spawn the core solution. 3 3 % 4 4 % Usage: 5 % responses=SpawnCore( models,variables,variabledescriptors)5 % responses=SpawnCore(femmodel,variables,variabledescriptors) 6 6 % 7 7 8 8 %recover parameters 9 verbose=models.dh.parameters.Verbose; 10 parameters=models.dh.parameters; 11 responsedescriptors=models.dh.parameters.responsedescriptors; 12 npart=models.dh.parameters.qmu_npart; 13 part=models.dh.parameters.qmu_part+1; %C indexing 14 numberofnodes=models.dh.parameters.NumberOfNodes; 9 verbose=femmodel.parameters.Verbose; 10 responsedescriptors=femmodel.parameters.responsedescriptors; 11 npart=femmodel.parameters.qmu_npart; 12 part=femmodel.parameters.qmu_part+1; %C indexing 13 numberofnodes=femmodel.parameters.NumberOfNodes; 15 14 16 disp ([' qmu iteration: ',num2str(counter)]);15 displaystring(verbose,'%s%i',[' qmu iteration:'],counter); 17 16 18 %first update the inputs to the modelsusing the variables provided to us by dakota.17 %first update the inputs to the femmodel using the variables provided to us by dakota. 19 18 count=1; 20 19 while count<=numel(variables) … … 27 26 else 28 27 root=qmuroot(descriptor); 29 param= parameters.(root);28 param=femmodel.parameters.(root); 30 29 31 30 %next npart values in variables are partition values for this param, collect them. … … 47 46 if analysis_type==DiagnosticAnalysisEnum(), 48 47 49 results=diagnostic_core( models);48 results=diagnostic_core(femmodel); 50 49 51 50 else … … 54 53 55 54 %process results 56 processedresults=processresults( models,results);55 processedresults=processresults(femmodel,results); 57 56 58 57 %now process the results to get response function values … … 60 59 for i=1:numel(responsedescriptors), 61 60 descriptor=responsedescriptors{i}; 62 responses(i)=qmuresponse( models,results,processedresults,descriptor);61 responses(i)=qmuresponse(femmodel,results,processedresults,descriptor); 63 62 end
Note:
See TracChangeset
for help on using the changeset viewer.