objectivefunctionC

PURPOSE ^

figure out number of parameters to search for

SYNOPSIS ^

function J =objectivefunctionC(search_vector,m,u_g_obs,parameters, optscaln, fit, gradients,control_types,analysis);

DESCRIPTION ^

figure out number of parameters to search for

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function J =objectivefunctionC(search_vector,m,u_g_obs,parameters, optscaln, fit, gradients,control_types,analysis);
0002 
0003 %figure out number of parameters to search for
0004 numparams=length(control_types);
0005 
0006 %Go through parameters, and update along gradients, multiplying by search_vector.
0007 for i=1:numparams,
0008     parameters{i}=parameters{i}+search_vector(i)*optscaln(i)*gradients{i};
0009 end
0010 
0011 %From parameters, build inputs for icediagnostic_core, using those new scaled parameters.
0012 for i=1:numparams,
0013     eval(['inputs.' control_types{i} '=parameters{' num2str(i) '};']);
0014 end
0015 
0016 %Run diagnostic with updated parameters.
0017 u_g =cielodiagnostic_core_nonlinear(m,m.params,inputs,analysis); u_g=full(u_g);
0018 
0019 %Compute misfit for this velocity field.
0020 inputs.fit=fit;
0021 J=Misfit(m.bgpdt,m.bgpdtb, m.est,m.ept, m.mpt, m.geom3, u_g,u_g_obs,inputs,analysis);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003