Changeset 8075
- Timestamp:
- 04/29/11 11:21:20 (14 years ago)
- Location:
- issm/trunk/src/m/solutions
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/steadystate_core.m
r6323 r8075 28 28 29 29 issmprintf(VerboseSolution,'\n%s',[' saving velocities, temperature and pressure for convergence']); 30 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,Vx OldEnum);31 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,Vy OldEnum);32 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,Vz OldEnum);33 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,Pressure OldEnum);30 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,VxPicardEnum); 31 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,VyPicardEnum); 32 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,VzPicardEnum); 33 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,PressurePicardEnum); 34 34 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,TemperatureOldEnum); 35 35 -
issm/trunk/src/m/solutions/steadystateconvergence.m
r4889 r8075 11 11 %Retrieve parameters 12 12 convergencecriterionvalue=femmodel.parameters.EpsRel; 13 velocityenums=[VxEnum;Vx OldEnum;VyEnum;VyOldEnum;VzEnum;VzOldEnum;PressureEnum;PressureOldEnum]; %pairs of enums (new and old) on which to carry out the converence tests13 velocityenums=[VxEnum;VxPicardEnum;VyEnum;VyPicardEnum;VzEnum;VzPicardEnum;PressureEnum;PressurePicardEnum]; %pairs of enums (new and old) on which to carry out the converence tests 14 14 temperatureenums=[TemperatureEnum;TemperatureOldEnum]; 15 15 convergencecriterion=[RelativeEnum]; %criterions for convergence, RelativeEnum or AbsoluteEnum
Note:
See TracChangeset
for help on using the changeset viewer.