0001 function [ru_g m]=cielodiagnostic_core_linear(m,params,inputs,analysis)
0002
0003
0004
0005 if isempty(inputs),
0006 clear inputs;
0007 end
0008
0009
0010 if isempty(m.y_s),
0011 flag_y_s= 0;
0012 else
0013 flag_y_s= 1;
0014 end
0015
0016
0017 params.kflag=1; params.pflag=1;
0018 params.ktflag=0;
0019
0020
0021 [rK_gg , rp_g , rdK_gg]=Emg(m.bgpdt,m.bgpdtb, m.est,m.lst,m.ept,m.mpt,m.geom3,params,inputs,analysis);
0022
0023
0024 [rK_ff, rK_fs] = Reducematrixfromgtof( rK_gg, m.G_mn, flag_y_s );
0025
0026
0027 rK_gg=IMdb('drop rK_gg');
0028 if ~isempty(rdK_gg), rdK_gg=IMdb('drop rdK_gg');end;
0029
0030
0031 [rp_f] = Reducerightside( rp_g, m.G_mn, rK_fs, m.y_s, flag_y_s );
0032
0033
0034 [rp_g rK_fs]=IMdb('drop rp_g rK_fs');
0035
0036
0037 [ru_f]=Solver(rK_ff,rp_f,{},params.solverstring);
0038
0039
0040 [rK_ff rp_f]=IMdb('drop rK_ff rp_f');
0041
0042
0043 ru_g= Mergesolvec( ru_f, m.G_mn, m.y_s );
0044
0045
0046 ru_f=IMdb('drop ru_f ');
0047
0048 end