Ignore:
Timestamp:
07/04/13 15:35:35 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: no need to marshall edges, they are constructed on the fly

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/classes/mesh.m

    r15131 r15432  
    116116                                        if md.prognostic.stabilization==3,
    117117                                                md = checkfield(md,'mesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d meshes');
    118                                                 md = checkfield(md,'mesh.edges','size',[NaN 4]);
    119                                                 md = checkfield(md,'mesh.edges(:,1:3)','>',0);
    120118                                        end
    121119                                case BalancethicknessSolutionEnum(),
    122120                                        if md.balancethickness.stabilization==3,
    123121                                                md = checkfield(md,'mesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d meshes');
    124                                                 md = checkfield(md,'mesh.edges','size',[NaN 4]);
    125                                                 md = checkfield(md,'mesh.edges(:,1:3)','>',0);
    126122                                        end
    127123                                case TransientSolutionEnum(),
    128124                                        if md.transient.isprognostic & md.prognostic.stabilization==3,
    129125                                                md = checkfield(md,'mesh.dimension','values',2,'message','Discontinuous Galerkin only supported for 2d meshes');
    130                                                 md = checkfield(md,'mesh.edges','size',[NaN 4]);
    131                                                 md = checkfield(md,'mesh.edges(:,1:3)','>',0);
    132126                                        end
    133127                                case ThermalSolutionEnum(),
     
    196190                        WriteData(fid,'object',obj,'fieldname','numberofelements','format','Integer');
    197191                        WriteData(fid,'object',obj,'fieldname','numberofvertices','format','Integer');
    198                         WriteData(fid,'object',obj,'fieldname','numberofedges','format','Integer');
    199192                        WriteData(fid,'object',obj,'fieldname','elementonbed','format','BooleanMat','mattype',2);
    200193                        WriteData(fid,'object',obj,'fieldname','elementonsurface','format','BooleanMat','mattype',2);
     
    203196                        WriteData(fid,'object',obj,'fieldname','lowerelements','format','DoubleMat','mattype',2);
    204197                        WriteData(fid,'object',obj,'fieldname','upperelements','format','DoubleMat','mattype',2);
    205                         WriteData(fid,'object',obj,'fieldname','edges','format','DoubleMat','mattype',3);
    206198                        WriteData(fid,'object',obj,'fieldname','elementconnectivity','format','DoubleMat','mattype',3);
    207199                        WriteData(fid,'object',obj,'fieldname','average_vertex_connectivity','format','Integer');
Note: See TracChangeset for help on using the changeset viewer.