source: issm/trunk/src/m/classes/public/importancefactors.m@ 2049

function factors=importancefactors(md,variablename,responsename)
%IMPORTANCEFACTORS - compute importance factors for a certain variable and response.
%
%   Usage:
%      factors=importancefactors(md,variablename,responsename)
%
%
%   Example: factors=importancefactors(md,'drag','max_vel');
%


variablenamelength=length(variablename);

%go through all response functions and find the one corresponding to the correct responsename
responsefunctions=md.dakotaresults.dresp_out;
found=0;
for i=1:length(responsefunctions),
        if strcmpi(responsefunctions(i).descriptor,responsename),
                found=i;
                break;
        end
end
if ~found,
        error('importancefactors error message: could not find correct response function');
end
responsefunctions=responsefunctions(found);
nfun=size(responsefunctions.var,1);

%Now recover response to the correct design variable
importancefactors=zeros(1,0);
count=0;
for i=1:nfun,
        desvar=responsefunctions.var{i};
        if strncmpi(desvar,variablename,variablenamelength),
                importancefactors(end+1)=responsefunctions.impfac(i);
                count=count+1;
        end
end

if count==0,
        error('importancefactors error message: could not find to response functions with corresponding design variable');
end

if count==1, %we have scalar
        factors=importancefactors;
        return;
else
        %Ok, get partitioning, unless already supplied
        if isempty(md.part),
                [epart npart]=MeshPartition(md,md.npart);
        else
                npart=md.part+1;
        end

        %distribute importance factor
        factors=importancefactors(npart)';
end