cieloprognostic_core

PURPOSE ^

INPUT function rh_g =cieloprognostic_core_linear(m,params,inputs,analysis)

SYNOPSIS ^

function h_g =cieloprognostic_core_linear(m,params,inputs,analysis)

DESCRIPTION ^

INPUT function rh_g =cieloprognostic_core_linear(m,params,inputs,analysis)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function h_g =cieloprognostic_core_linear(m,params,inputs,analysis)
0002 %INPUT function rh_g =cieloprognostic_core_linear(m,params,inputs,analysis)
0003 
0004     %recover flag_y_s from m.y_s
0005     if size(m.y_s,1)==0,
0006         flag_y_s= 0;
0007     else
0008         flag_y_s= 1;
0009     end
0010 
0011     %stiffness and load generation only:
0012     params.kflag=1; params.pflag=1;
0013     params.ktflag=0;
0014 
0015     %system matrices
0016     [K_gg , p_g , dK_gg]=Emg(m.bgpdt,m.bgpdtb, m.est,m.lst,m.ept,m.mpt,m.geom3,params,inputs,analysis);
0017 
0018     %Reduce tangent matrix from g size to f size
0019     [K_ff, K_fs] = Reducematrixfromgtof( K_gg, m.G_mn, flag_y_s ); 
0020     clear K_gg dK_gg
0021     
0022     %Reduce load from g size to f size
0023     [rp_f] = Reducerightside( rp_g, m.G_mn, rK_fs, m.y_s, flag_y_s );
0024     [clear p_g K_fs
0025 
0026     %Solve
0027     [h_f]=Solver(K_ff,p_f,params.solverstring);
0028     clear K_ff p_f
0029     
0030     %Merge back to g set
0031     h_g= Mergesolvec( h_f, m.G_mn, m.y_s, m.uset); 
0032     
0033     if count>2,
0034         clear h_f
0035     end

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