Changeset 5776
- Timestamp:
- 09/13/10 14:19:35 (15 years ago)
- Location:
- issm/trunk/src/m
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/NewFemModel.m
r5103 r5776 19 19 20 20 %Initialize some fiels with empty celils 21 femmodel.m_yg=cell(nummodels,1);22 21 femmodel.m_nodesets=cell(nummodels,1); 23 22 femmodel.m_ys=cell(nummodels,1); … … 32 31 femmodel=SetCurrentConfiguration(femmodel,analysis_type); 33 32 34 displaystring(md.verbose,'%s',' generating degrees of freedom...');33 displaystring(md.verbose,'%s',' generating vertices degrees of freedom'); 35 34 if ~isfield(femmodel,'part'), 36 35 [femmodel.vertices,femmodel.part,femmodel.tpart]=VerticesDof(femmodel.vertices, femmodel.parameters); %do not create partition vector twice! we only have one set of vertices! 37 36 end 37 38 displaystring(md.verbose,'%s',' resolve node constraints'); 39 [femmodel.nodes]=SpcNodes(femmodel.nodes,femmodel.constraints,analysis_type); 40 41 displaystring(md.verbose,'%s',' create nodal degrees of freedom'); 38 42 [femmodel.nodes]=NodesDof(femmodel.nodes,femmodel.parameters); 39 43 40 displaystring(md.verbose,'%s',' generating single point constraints...');41 [femmodel.nodes,femmodel.m_yg{i}]=SpcNodes(femmodel.nodes,femmodel.constraints,analysis_type);44 displaystring(md.verbose,'%s',' create nodal constraints vector'); 45 femmodel.m_ys{i}=CreateNodalConstraints(femmodel.nodes,analysis_type); 42 46 43 displaystring(md.verbose,'%s',' generating node sets...');47 displaystring(md.verbose,'%s',' create node sets'); 44 48 femmodel.m_nodesets{i}=BuildNodeSets(femmodel.nodes,analysis_type); 45 49 46 displaystring(md.verbose,'%s',' reducing single point constraints vector...'); 47 femmodel.m_ys{i}=Reducevectorgtos(femmodel.m_yg{i},femmodel.m_nodesets{i}); 48 49 displaystring(md.verbose,'%s',' configuring elements and loads...'); 50 displaystring(md.verbose,'%s',' configuring elements and loads'); 50 51 [femmodel.elements,femmodel.loads,femmodel.nodes,femmodel.parameters] = ConfigureObjects( femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices,femmodel.materials,femmodel.parameters); 51 52 end -
issm/trunk/src/m/solutions/ResetBoundaryConditions.m
r4520 r5776 20 20 %For this analysis_type, free existing boundary condition vectors: 21 21 analysis_counter=femmodel.parameters.AnalysisCounter+1; %matlab indexing on counter 22 23 %global dof set24 femmodel.m_yg{analysis_counter}=[];25 %in the s-set26 22 femmodel.m_ys{analysis_counter}=[]; 27 23 28 % Now, duplicate ug (the solution vector) into the boundary conditions vector on the g-set29 femmodel.m_yg{analysis_counter}=ug;24 %Reduce from g to s set 25 ys=Reducevectorgtos(ug,femmodel.m_nodesets{analysis_counter}); 30 26 31 %Reduce from g to s set 32 femmodel.m_ys{analysis_counter}=Reducevectorgtos(femmodel.m_yg{analysis_counter},femmodel.m_nodesets{analysis_counter}); 27 %in the s-set 28 femmodel.m_ys{analysis_counter}=ys; 29 -
issm/trunk/src/m/solutions/SetCurrentConfiguration.m
r5377 r5776 39 39 %activate matrices and vectors: 40 40 femmodel.nodesets=femmodel.m_nodesets{found}; 41 femmodel.yg=femmodel.m_yg{found};42 41 femmodel.ys=femmodel.m_ys{found}; 43 42 -
issm/trunk/src/m/solvers/solver_adjoint_linear.m
r5693 r5776 5 5 % femmodel =solver_adjoint_linear(femmodel) 6 6 7 [K_gg, p_g,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);7 [K_gg,K_ff,K_fs,p_g,p_f,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters); 8 8 9 9 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, femmodel.nodesets,femmodel.parameters); -
issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m
r5693 r5776 27 27 old_uf=uf; 28 28 29 [K_gg, p_g,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,loads,femmodel.materials,femmodel.parameters);29 [K_gg,K_ff,K_fs,p_g,p_f,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,loads,femmodel.materials,femmodel.parameters); 30 30 31 31 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, femmodel.nodesets,femmodel.parameters); -
issm/trunk/src/m/solvers/solver_linear.m
r5693 r5776 5 5 % femmodel =solver_linear(femmodel) 6 6 7 [K_gg, p_g,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);7 [K_gg,K_ff,K_fs,p_g,p_f,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters); 8 8 9 9 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, femmodel.nodesets,femmodel.parameters); -
issm/trunk/src/m/solvers/solver_thermal_nonlinear.m
r5707 r5776 20 20 while(~converged), 21 21 22 [K_gg, p_g,melting_offset]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);22 [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); 23 23 24 24 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, femmodel.nodesets,femmodel.parameters);
Note:
See TracChangeset
for help on using the changeset viewer.