wzcelltype.cxx

有限元学习研究用源代码(老外的),供科研人员参考

#include "wzcelltype.hxx"

void wzCellTypeInit(); 

int	wzCellTypepoints[IBG2CTMAX];
int	wzCellTypelines[IBG2CTMAX];
int	wzCellTypesides[IBG2CTMAX];
int	wzCellTypeudim[IBG2CTMAX];
int	wzCellTypegdim[IBG2CTMAX];
int	wzCellTypecdim[IBG2CTMAX];
int	wzCellTypesize[IBG2CTMAX];
char	*wzCellTypename[IBG2CTMAX];
int	wzCellTypeepoint1[IBG2CTMAX][IBG2CEMAX];
int	wzCellTypeepoint2[IBG2CTMAX][IBG2CEMAX];
int	wzCellTypeesidel[IBG2CTMAX][IBG2CEMAX];
int	wzCellTypeesider[IBG2CTMAX][IBG2CEMAX];
int	wzCellTypessize[IBG2CTMAX][IBG2CSMAX];
int	wzCellTypesside[IBG2CTMAX][IBG2CSMAX][4];
int	wzCellTypesline[IBG2CTMAX][IBG2CSMAX][4];
int	wzCellTypespoint[IBG2CTMAX][IBG2CSMAX][5];

void wzCellTypeInit()
{
  int t,s,e,ee;
  // <<< wzCellType0Nothing >>> 
  wzCellTypeName(wzCellType0Nothing)	= "Null";
  wzCellTypePoints(wzCellType0Nothing)	= 0;
  wzCellTypeLines(wzCellType0Nothing)	= 0;
  wzCellTypeSides(wzCellType0Nothing)	= 0;
  wzCellTypeCODIM(wzCellType0Nothing)	= wzIs_Nothing;
  wzCellTypeGDIM(wzCellType0Nothing)	= wzIs_Nothing;
  wzCellTypeCDIM(wzCellType0Nothing)	= 0;
  // <<< wzCellType0Node >>> 
  wzCellTypeName(wzCellType0Node)	= "Node 0D";
  wzCellTypePoints(wzCellType0Node)	= 1;
  wzCellTypeLines(wzCellType0Node)	= 0;
  wzCellTypeSides(wzCellType0Node)	= 0;
  wzCellTypeCODIM(wzCellType0Node)	= wzIsRegion;
  wzCellTypeGDIM(wzCellType0Node)	= 0;
  wzCellTypeCDIM(wzCellType0Node)	= 0;
  // <<< wzCellType1Node >>> 
  wzCellTypeName(wzCellType1Node)	= "Node 1D";
  wzCellTypePoints(wzCellType1Node)	= 1;
  wzCellTypeLines(wzCellType1Node)	= 0;
  wzCellTypeSides(wzCellType1Node)	= 0;
  wzCellTypeCODIM(wzCellType1Node)	= wzIsFace;
  wzCellTypeGDIM(wzCellType1Node)	= 1;
  wzCellTypeCDIM(wzCellType1Node)	= 0;
  // <<< wzCellType2Node >>> 
  wzCellTypeName(wzCellType2Node)	= "Node 2D";
  wzCellTypePoints(wzCellType2Node)	= 1;
  wzCellTypeLines(wzCellType2Node)	= 0;
  wzCellTypeSides(wzCellType2Node)	= 0;
  wzCellTypeCODIM(wzCellType2Node)	= wzIsEdge;
  wzCellTypeGDIM(wzCellType2Node)	= 2;
  wzCellTypeCDIM(wzCellType2Node)	= 0;
  // <<< wzCellType3Node >>> 
  wzCellTypeName(wzCellType3Node)	= "Node 3D";
  wzCellTypePoints(wzCellType3Node)	= 1;
  wzCellTypeLines(wzCellType3Node)	= 0;
  wzCellTypeSides(wzCellType3Node)	= 0;
  wzCellTypeCODIM(wzCellType3Node)	= wzIsVertex;
  wzCellTypeGDIM(wzCellType3Node)	= 3;
  wzCellTypeCDIM(wzCellType3Node)	= 0;
  // <<< wzCellType1Line >>> 
  wzCellTypeName(wzCellType1Line)	= "Edge 1D";
  wzCellTypePoints(wzCellType1Line)	= 2;
  wzCellTypeLines(wzCellType1Line)	= 1;
  wzCellTypeSides(wzCellType1Line)	= 2;
  wzCellTypeCODIM(wzCellType1Line)	= wzIsRegion;
  wzCellTypeGDIM(wzCellType1Line)	= 1;
  wzCellTypeCDIM(wzCellType1Line)	= 1;
  wzCellTypeEPoint1(wzCellType1Line,0)	= 0;
  wzCellTypeEPoint2(wzCellType1Line,0)	= 1;
  wzCellTypeSSize(wzCellType1Line,0)	= 1;
  wzCellTypeSSize(wzCellType1Line,1)	= 1;
  wzCellTypeSPoint1(wzCellType1Line,0,0)	= 0;
  wzCellTypeSPoint1(wzCellType1Line,1,0)	= 1;
  wzCellTypeSPoint1(wzCellType1Line,0,1)	= 0;
  wzCellTypeSPoint1(wzCellType1Line,1,1)	= 1;
  // <<< wzCellType2Line >>> 
  wzCellTypeName(wzCellType2Line)	= "Edge 2D";
  wzCellTypePoints(wzCellType2Line)	= 2;
  wzCellTypeLines(wzCellType2Line)	= 1;
  wzCellTypeSides(wzCellType2Line) 	= 2;
  wzCellTypeCODIM(wzCellType2Line)	= wzIsFace;
  wzCellTypeGDIM(wzCellType2Line)	= 2;
  wzCellTypeCDIM(wzCellType2Line)	= 1;
  wzCellTypeEPoint1(wzCellType2Line,0)	= 0;
  wzCellTypeEPoint2(wzCellType2Line,0)	= 1;
  wzCellTypeSSize(wzCellType2Line,0)	= 1;
  wzCellTypeSSize(wzCellType2Line,1)	= 1;
  wzCellTypeSPoint1(wzCellType2Line,0,0)	= 0;
  wzCellTypeSPoint1(wzCellType2Line,1,0)	= 1;
  wzCellTypeSPoint1(wzCellType2Line,0,1)	= 0;
  wzCellTypeSPoint1(wzCellType2Line,1,1)	= 1;
  // <<< wzCellType3Line >>> 
  wzCellTypeName(wzCellType3Line)	= "Edge 3D";
  wzCellTypePoints(wzCellType3Line)	= 2;
  wzCellTypeLines(wzCellType3Line)	= 1;
  wzCellTypeSides(wzCellType3Line)	= 2;
  wzCellTypeCODIM(wzCellType3Line)	= wzIsEdge;
  wzCellTypeGDIM(wzCellType3Line)	= 3;
  wzCellTypeCDIM(wzCellType3Line)	= 1;
  wzCellTypeEPoint1(wzCellType3Line,0)	= 0;
  wzCellTypeEPoint2(wzCellType3Line,0)	= 1;
  wzCellTypeSSize(wzCellType3Line,0)	= 1;
  wzCellTypeSSize(wzCellType3Line,1)	= 1;
  wzCellTypeSPoint1(wzCellType3Line,0,0)	= 0;
  wzCellTypeSPoint1(wzCellType3Line,1,0)	= 1;
  wzCellTypeSPoint1(wzCellType3Line,0,1)	= 0;
  wzCellTypeSPoint1(wzCellType3Line,1,1)	= 1;
  // <<< wzCellType2Triangle >>> 
  wzCellTypeName(wzCellType2Triangle)	= "Triangle 2D";
  wzCellTypePoints(wzCellType2Triangle)	= 3;
  wzCellTypeLines(wzCellType2Triangle)	= 3;
  wzCellTypeSides(wzCellType2Triangle)	= 3;
  wzCellTypeCODIM(wzCellType2Triangle)	= wzIsRegion;
  wzCellTypeGDIM(wzCellType2Triangle)	= 2;
  wzCellTypeCDIM(wzCellType2Triangle)	= 2;
  wzCellTypeEPoint1(wzCellType2Triangle,0)	= 1;
  wzCellTypeEPoint2(wzCellType2Triangle,0)	= 2;
  wzCellTypeEPoint1(wzCellType2Triangle,1)	= 2;
  wzCellTypeEPoint2(wzCellType2Triangle,1)	= 0;
  wzCellTypeEPoint1(wzCellType2Triangle,2)	= 0;
  wzCellTypeEPoint2(wzCellType2Triangle,2)	= 1;
  wzCellTypeSSize(wzCellType2Triangle,0)	= 2;
  wzCellTypeSSize(wzCellType2Triangle,1) 	= 2;
  wzCellTypeSSize(wzCellType2Triangle,2)	= 2;
  wzCellTypeSSide(wzCellType2Triangle,0,0)	= 2;
  wzCellTypeSSide(wzCellType2Triangle,0,1)	= 1;
  wzCellTypeSSide(wzCellType2Triangle,1,0)	= 0;
  wzCellTypeSSide(wzCellType2Triangle,1,1)	= 2;
  wzCellTypeSSide(wzCellType2Triangle,2,0)	= 1;
  wzCellTypeSSide(wzCellType2Triangle,2,1)	= 0;
  wzCellTypeSPoint1(wzCellType2Triangle,0,0)	= 1;
  wzCellTypeSPoint1(wzCellType2Triangle,0,1)	= 2;
  wzCellTypeSPoint1(wzCellType2Triangle,1,0)	= 2;
  wzCellTypeSPoint1(wzCellType2Triangle,1,1)	= 0;
  wzCellTypeSPoint1(wzCellType2Triangle,2,0)	= 0;
  wzCellTypeSPoint1(wzCellType2Triangle,2,1)	= 1;
  // <<< wzCellType3Triangle >>> 
  wzCellTypeName(wzCellType3Triangle)	= "Triangle 3D";
  wzCellTypePoints(wzCellType3Triangle)	= 3;
  wzCellTypeLines(wzCellType3Triangle)	= 3;
  wzCellTypeSides(wzCellType3Triangle)	= 3;
  wzCellTypeCODIM(wzCellType3Triangle)	= wzIsFace;
  wzCellTypeGDIM(wzCellType3Triangle)	= 3;
  wzCellTypeCDIM(wzCellType3Triangle)	= 2;
  wzCellTypeEPoint1(wzCellType3Triangle,0)	= 1;
  wzCellTypeEPoint2(wzCellType3Triangle,0)	= 2;
  wzCellTypeEPoint1(wzCellType3Triangle,1)	= 2;
  wzCellTypeEPoint2(wzCellType3Triangle,1)	= 0;
  wzCellTypeEPoint1(wzCellType3Triangle,2)	= 0;
  wzCellTypeEPoint2(wzCellType3Triangle,2)	= 1;
  wzCellTypeSSize(wzCellType3Triangle,0)	= 2;
  wzCellTypeSSize(wzCellType3Triangle,1) 	= 2;
  wzCellTypeSSize(wzCellType3Triangle,2)	= 2;
  wzCellTypeSSide(wzCellType3Triangle,0,0)	= 2;
  wzCellTypeSSide(wzCellType3Triangle,0,1)	= 1;
  wzCellTypeSSide(wzCellType3Triangle,1,0)	= 0;
  wzCellTypeSSide(wzCellType3Triangle,1,1)	= 2;
  wzCellTypeSSide(wzCellType3Triangle,2,0)	= 1;
  wzCellTypeSSide(wzCellType3Triangle,2,1)	= 0;
  wzCellTypeSPoint1(wzCellType3Triangle,0,0)	= 1;
  wzCellTypeSPoint1(wzCellType3Triangle,0,1)	= 2;
  wzCellTypeSPoint1(wzCellType3Triangle,1,0)	= 2;
  wzCellTypeSPoint1(wzCellType3Triangle,1,1)	= 0;
  wzCellTypeSPoint1(wzCellType3Triangle,2,0)	= 0;
  wzCellTypeSPoint1(wzCellType3Triangle,2,1)	= 1;
  // <<< wzCellType2Rectangle >>> 
  wzCellTypeName(wzCellType2Rectangle)	= "Rectangle 2D";
  wzCellTypePoints(wzCellType2Rectangle)	= 4;
  wzCellTypeLines(wzCellType2Rectangle)	= 4;
  wzCellTypeSides(wzCellType2Rectangle)	= 4;
  wzCellTypeCODIM(wzCellType2Rectangle)	= wzIsRegion;
  wzCellTypeGDIM(wzCellType2Rectangle)	= 2;
  wzCellTypeCDIM(wzCellType2Rectangle)	= 2;
  wzCellTypeEPoint1(wzCellType2Rectangle,0)	= 0;
  wzCellTypeEPoint2(wzCellType2Rectangle,0)	= 1;
  wzCellTypeEPoint1(wzCellType2Rectangle,1)	= 1;