Changeset 18197
- Timestamp:
- 06/30/14 14:38:11 (11 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/flowequation.m
r18185 r18197 6 6 classdef flowequation 7 7 properties (SetAccess=public) 8 isSIA = 0;9 isSSA = 0;10 isL1L2 = 0;11 isHO = 0;12 isFS = 0;13 fe_SSA = '';14 fe_HO = '';15 fe_FS = '';16 XTH_r= 1.;17 XTH_theta = 0.;18 vertex_equation = NaN;19 element_equation = NaN;20 borderSSA = NaN;21 borderHO = NaN;22 borderFS = NaN;8 isSIA = 0; 9 isSSA = 0; 10 isL1L2 = 0; 11 isHO = 0; 12 isFS = 0; 13 fe_SSA = ''; 14 fe_HO = ''; 15 fe_FS = ''; 16 augmented_lagrangian_r = 1.; 17 XTH_theta = 0.; 18 vertex_equation = NaN; 19 element_equation = NaN; 20 borderSSA = NaN; 21 borderHO = NaN; 22 borderFS = NaN; 23 23 end 24 24 methods (Static) … … 135 135 md = checkfield(md,'fieldname','flowequation.fe_HO' ,'values',{'P1','P1bubble','P1bubblecondensed','P1xP2','P2xP1','P2','P2bubble','P1xP3','P2xP4'}); 136 136 md = checkfield(md,'fieldname','flowequation.fe_FS' ,'values',{'P1P1','P1P1GLS','MINIcondensed','MINI','TaylorHood','LATaylorHood','XTaylorHood','OneLayerP4z','CrouzeixRaviart'}); 137 md = checkfield(md,'fieldname','flowequation. XTH_r','numel',[1],'>',0.);137 md = checkfield(md,'fieldname','flowequation.augmented_lagrangian_r','numel',[1],'>=',0.); 138 138 md = checkfield(md,'fieldname','flowequation.XTH_theta','numel',[1],'>=',0.,'<',0.5); 139 139 md = checkfield(md,'fieldname','flowequation.borderSSA','size',[md.mesh.numberofvertices 1],'values',[0 1]); … … 192 192 WriteData(fid,'enum',FlowequationFeHOEnum() ,'data',StringToEnum(obj.fe_HO) ,'format','Integer'); 193 193 WriteData(fid,'enum',FlowequationFeFSEnum() ,'data',StringToEnum(obj.fe_FS) ,'format','Integer'); 194 WriteData(fid,'enum',AugmentedLagrangianREnum() ,'data',obj.XTH_r ,'format','Double');194 WriteData(fid,'enum',AugmentedLagrangianREnum(),'data',obj.augmented_lagrangian_r ,'format','Double'); 195 195 WriteData(fid,'enum',AugmentedLagrangianThetaEnum() ,'data',obj.XTH_theta ,'format','Double'); 196 196 WriteData(fid,'object',obj,'fieldname','borderSSA','format','DoubleMat','mattype',1); -
issm/trunk-jpl/src/m/classes/flowequation.py
r18179 r18197 17 17 def __init__(self): # {{{ 18 18 19 self.isSIA = 020 self.isSSA = 021 self.isL1L2 = 022 self.isHO = 023 self.isFS = 024 self.fe_SSA = ''25 self.fe_HO = ''26 self.fe_FS = ''27 self. XTH_r= 1.28 self.XTH_theta = 0.29 self.vertex_equation = float('NaN')30 self.element_equation = float('NaN')31 self.borderSSA = float('NaN')32 self.borderHO = float('NaN')33 self.borderFS = float('NaN')19 self.isSIA = 0 20 self.isSSA = 0 21 self.isL1L2 = 0 22 self.isHO = 0 23 self.isFS = 0 24 self.fe_SSA = '' 25 self.fe_HO = '' 26 self.fe_FS = '' 27 self.augmented_lagrangian_r = 1. 28 self.XTH_theta = 0. 29 self.vertex_equation = float('NaN') 30 self.element_equation = float('NaN') 31 self.borderSSA = float('NaN') 32 self.borderHO = float('NaN') 33 self.borderFS = float('NaN') 34 34 35 35 #set defaults … … 82 82 md = checkfield(md,'fieldname','flowequation.borderHO','size',[md.mesh.numberofvertices],'values',[0,1]) 83 83 md = checkfield(md,'fieldname','flowequation.borderFS','size',[md.mesh.numberofvertices],'values',[0,1]) 84 md = checkfield(md,'fieldname','flowequation. XTH_r','numel',[1],'>',0.)84 md = checkfield(md,'fieldname','flowequation.augmented_lagrangian_r','numel',[1],'>',0.) 85 85 md = checkfield(md,'fieldname','flowequation.XTH_theta','numel',[1],'>=',0.,'<',.5) 86 86 if m.strcmp(md.mesh.domaintype(),'2Dhorizontal'): … … 111 111 WriteData(fid,'enum',FlowequationFeHOEnum() ,'data',StringToEnum(self.fe_HO)[0] ,'format','Integer') 112 112 WriteData(fid,'enum',FlowequationFeFSEnum() ,'data',StringToEnum(self.fe_FS)[0] ,'format','Integer') 113 WriteData(fid,'enum',AugmentedLagrangianREnum() ,'data',self.XTH_r ,'format','Double')113 WriteData(fid,'enum',AugmentedLagrangianREnum(),'data',self.augmented_lagrangian_r ,'format','Double') 114 114 WriteData(fid,'enum',AugmentedLagrangianThetaEnum() ,'data',self.XTH_theta ,'format','Double') 115 115 WriteData(fid,'object',self,'fieldname','borderSSA','format','DoubleMat','mattype',1)
Note:
See TracChangeset
for help on using the changeset viewer.