0001 function J =objectivefunctionC(search_vector,m,u_g_obs,parameters, optscaln, fit, gradients,control_types,analysis);
0002
0003
0004 numparams=length(control_types);
0005
0006
0007 for i=1:numparams,
0008 parameters{i}=parameters{i}+search_vector(i)*optscaln(i)*gradients{i};
0009 end
0010
0011
0012 for i=1:numparams,
0013 eval(['inputs.' control_types{i} '=parameters{' num2str(i) '};']);
0014 end
0015
0016
0017 u_g =cielodiagnostic_core_nonlinear(m,m.params,inputs,analysis); u_g=full(u_g);
0018
0019
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);