Changeset 6581 for issm/trunk/src/m/solvers/solver_thermal_nonlinear.m
- Timestamp:
- 11/15/10 13:55:35 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solvers/solver_thermal_nonlinear.m
r6323 r6581 25 25 if kffpartitioning, 26 26 [K_gg,K_ff,K_fs,p_g,p_f,melting_offset]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters); 27 p_f = Reduceload( p_f, K_fs, femmodel.ys ,femmodel.parameters);27 p_f = Reduceload( p_f, K_fs, femmodel.ys); 28 28 else 29 29 [K_gg,K_ff,K_fs,p_g,p_f,melting_offset]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters); 30 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, femmodel.nodesets ,femmodel.parameters);31 p_f = Reduceloadfromgtof( p_g, K_fs, femmodel.ys, femmodel.nodesets ,femmodel.parameters);30 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, femmodel.nodesets); 31 p_f = Reduceloadfromgtof( p_g, K_fs, femmodel.ys, femmodel.nodesets); 32 32 end 33 33 34 34 issmprintf(VerboseSolver(),'%s%g',' condition number of stiffness matrix: ',condest(K_ff)); 35 35 t_f=Solver(K_ff,p_f,[],femmodel.parameters); 36 t_g= Mergesolutionfromftog( t_f, femmodel.ys, femmodel.nodesets ,femmodel.parameters);36 t_g= Mergesolutionfromftog( t_f, femmodel.ys, femmodel.nodesets); 37 37 38 38 [femmodel.elements,femmodel.materials]=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,t_g);
Note:
See TracChangeset
for help on using the changeset viewer.