Ignore:
Timestamp:
06/24/10 13:46:42 (15 years ago)
Author:
Eric.Larour
Message:

Some updates on Qmu module: simplifications

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/SpawnCore.m

    r4131 r4193  
    1 function responses=SpawnCore(models,variables,variabledescriptors,counter);
     1function responses=SpawnCore(femmodel,variables,variabledescriptors,counter);
    22%SPAWNCORE - for Qmu analysis, using Dakota. Spawn the core solution.
    33%
    44%   Usage:
    5 %      responses=SpawnCore(models,variables,variabledescriptors)
     5%      responses=SpawnCore(femmodel,variables,variabledescriptors)
    66%
    77
    88%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;
     9verbose=femmodel.parameters.Verbose;
     10responsedescriptors=femmodel.parameters.responsedescriptors;
     11npart=femmodel.parameters.qmu_npart;
     12part=femmodel.parameters.qmu_part+1; %C indexing
     13numberofnodes=femmodel.parameters.NumberOfNodes;
    1514
    16 disp(['   qmu iteration: ',num2str(counter)]);
     15displaystring(verbose,'%s%i',['   qmu iteration:'],counter);
    1716
    18 %first update the inputs to the models using the variables provided to us by dakota.
     17%first update the inputs to the femmodel using the variables provided to us by dakota.
    1918count=1;
    2019while count<=numel(variables)
     
    2726        else
    2827                root=qmuroot(descriptor);
    29                 param=parameters.(root);
     28                param=femmodel.parameters.(root);
    3029
    3130                %next npart values in variables are partition values for this param, collect them.
     
    4746if analysis_type==DiagnosticAnalysisEnum(),
    4847
    49         results=diagnostic_core(models);
     48        results=diagnostic_core(femmodel);
    5049
    5150else
     
    5453
    5554%process results
    56 processedresults=processresults(models,results);
     55processedresults=processresults(femmodel,results);
    5756
    5857%now process the results to get response function values
     
    6059for i=1:numel(responsedescriptors),
    6160        descriptor=responsedescriptors{i};
    62         responses(i)=qmuresponse(models,results,processedresults,descriptor);
     61        responses(i)=qmuresponse(femmodel,results,processedresults,descriptor);
    6362end
Note: See TracChangeset for help on using the changeset viewer.