wzcelltype.cxx

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

  wzCellTypeEPoint2(wzCellType2Rectangle,1)	= 2;
  wzCellTypeEPoint1(wzCellType2Rectangle,2)	= 2;
  wzCellTypeEPoint2(wzCellType2Rectangle,2)	= 3;
  wzCellTypeEPoint1(wzCellType2Rectangle,3)	= 3;
  wzCellTypeEPoint2(wzCellType2Rectangle,3)	= 0;
  wzCellTypeSSize(wzCellType2Rectangle,0)	= 2;
  wzCellTypeSSize(wzCellType2Rectangle,1) 	= 2;
  wzCellTypeSSize(wzCellType2Rectangle,2)	= 2;
  wzCellTypeSSize(wzCellType2Rectangle,3)	= 2;
  wzCellTypeSSide(wzCellType2Rectangle,0,0)	= 3;
  wzCellTypeSSide(wzCellType2Rectangle,0,1)	= 1;
  wzCellTypeSSide(wzCellType2Rectangle,1,0)	= 0;
  wzCellTypeSSide(wzCellType2Rectangle,1,1)	= 2;
  wzCellTypeSSide(wzCellType2Rectangle,2,0)	= 1;
  wzCellTypeSSide(wzCellType2Rectangle,2,1)	= 3;
  wzCellTypeSSide(wzCellType2Rectangle,3,0)	= 2;
  wzCellTypeSSide(wzCellType2Rectangle,3,1)	= 0;
  wzCellTypeSPoint1(wzCellType2Rectangle,0,0)	= 0;
  wzCellTypeSPoint1(wzCellType2Rectangle,0,1)	= 1;
  wzCellTypeSPoint1(wzCellType2Rectangle,1,0)	= 1;
  wzCellTypeSPoint1(wzCellType2Rectangle,1,1)	= 2;
  wzCellTypeSPoint1(wzCellType2Rectangle,2,0)	= 2;
  wzCellTypeSPoint1(wzCellType2Rectangle,2,1)	= 3;
  wzCellTypeSPoint1(wzCellType2Rectangle,3,0)	= 3;
  wzCellTypeSPoint1(wzCellType2Rectangle,3,1)	= 0;
  // <<< wzCellType3Rectangle >>> 
  wzCellTypeName(wzCellType3Rectangle)		= "Rectangle 3D";
  wzCellTypePoints(wzCellType3Rectangle)	= 4;
  wzCellTypeLines(wzCellType3Rectangle)	= 4;
  wzCellTypeSides(wzCellType3Rectangle)	= 4;
  wzCellTypeCODIM(wzCellType3Rectangle)	= wzIsFace;
  wzCellTypeGDIM(wzCellType3Rectangle)	= 3;
  wzCellTypeCDIM(wzCellType3Rectangle)	= 2;
  wzCellTypeEPoint1(wzCellType3Rectangle,0)	= 0;
  wzCellTypeEPoint2(wzCellType3Rectangle,0)	= 1;
  wzCellTypeEPoint1(wzCellType3Rectangle,1)	= 1;
  wzCellTypeEPoint2(wzCellType3Rectangle,1)	= 2;
  wzCellTypeEPoint1(wzCellType3Rectangle,2)	= 2;
  wzCellTypeEPoint2(wzCellType3Rectangle,2)	= 3;
  wzCellTypeEPoint1(wzCellType3Rectangle,3)	= 3;
  wzCellTypeEPoint2(wzCellType3Rectangle,3)	= 0;
  wzCellTypeSSize(wzCellType3Rectangle,0)	= 2;
  wzCellTypeSSize(wzCellType3Rectangle,1) 	= 2;
  wzCellTypeSSize(wzCellType3Rectangle,2)	= 2;
  wzCellTypeSSize(wzCellType3Rectangle,3)	= 2;
  wzCellTypeSSide(wzCellType3Rectangle,0,0)	= 3;
  wzCellTypeSSide(wzCellType3Rectangle,0,1)	= 1;
  wzCellTypeSSide(wzCellType3Rectangle,1,0)	= 0;
  wzCellTypeSSide(wzCellType3Rectangle,1,1)	= 2;
  wzCellTypeSSide(wzCellType3Rectangle,2,0)	= 1;
  wzCellTypeSSide(wzCellType3Rectangle,2,1)	= 3;
  wzCellTypeSSide(wzCellType3Rectangle,3,0)	= 2;
  wzCellTypeSSide(wzCellType3Rectangle,3,1)	= 0;
  wzCellTypeSPoint1(wzCellType3Rectangle,0,0)	= 0;
  wzCellTypeSPoint1(wzCellType3Rectangle,0,1)	= 1;
  wzCellTypeSPoint1(wzCellType3Rectangle,1,0)	= 1;
  wzCellTypeSPoint1(wzCellType3Rectangle,1,1)	= 2;
  wzCellTypeSPoint1(wzCellType3Rectangle,2,0)	= 2;
  wzCellTypeSPoint1(wzCellType3Rectangle,2,1)	= 3;
  wzCellTypeSPoint1(wzCellType3Rectangle,3,0)	= 3;
  wzCellTypeSPoint1(wzCellType3Rectangle,3,1)	= 0;
  // <<< wzCellType3Tetrahedron >>> 
  wzCellTypeName(wzCellType3Tetrahedron)	= "Tetrahedron 3D";
  wzCellTypePoints(wzCellType3Tetrahedron)	= 4;
  wzCellTypeLines(wzCellType3Tetrahedron)	= 6;
  wzCellTypeSides(wzCellType3Tetrahedron)	= 4;
  wzCellTypeCODIM(wzCellType3Tetrahedron)	= wzIsRegion;
  wzCellTypeGDIM(wzCellType3Tetrahedron)	= 3;
  wzCellTypeCDIM(wzCellType3Tetrahedron)	= 3;
  wzCellTypeEPoint1(wzCellType3Tetrahedron,0)	= 2;	// Line 0 
  wzCellTypeEPoint2(wzCellType3Tetrahedron,0)	= 3;
  wzCellTypeESideL(wzCellType3Tetrahedron,0)	= 0;
  wzCellTypeESideR(wzCellType3Tetrahedron,0)	= 1;
  wzCellTypeEPoint1(wzCellType3Tetrahedron,1)	= 3;	// Line 1 
  wzCellTypeEPoint2(wzCellType3Tetrahedron,1)	= 1;
  wzCellTypeESideL(wzCellType3Tetrahedron,1)	= 0;
  wzCellTypeESideR(wzCellType3Tetrahedron,1)	= 2;
  wzCellTypeEPoint1(wzCellType3Tetrahedron,2)	= 1;	// Line 2 
  wzCellTypeEPoint2(wzCellType3Tetrahedron,2)   	= 2;
  wzCellTypeESideL(wzCellType3Tetrahedron,2)	= 0;
  wzCellTypeESideR(wzCellType3Tetrahedron,2)	= 3;
  wzCellTypeEPoint1(wzCellType3Tetrahedron,3)	= 0;	// Line 3 
  wzCellTypeEPoint2(wzCellType3Tetrahedron,3)   	= 1;
  wzCellTypeESideL(wzCellType3Tetrahedron,3)	= 2;
  wzCellTypeESideR(wzCellType3Tetrahedron,3)	= 3;
  wzCellTypeEPoint1(wzCellType3Tetrahedron,4)	= 0;	// Line 4 
  wzCellTypeEPoint2(wzCellType3Tetrahedron,4)   	= 2;
  wzCellTypeESideL(wzCellType3Tetrahedron,4)	= 3;
  wzCellTypeESideR(wzCellType3Tetrahedron,4)	= 1;
  wzCellTypeEPoint1(wzCellType3Tetrahedron,5)	= 0;	// Line 5 
  wzCellTypeEPoint2(wzCellType3Tetrahedron,5)   	= 3;
  wzCellTypeESideL(wzCellType3Tetrahedron,5)	= 1;
  wzCellTypeESideR(wzCellType3Tetrahedron,5)	= 2;
  wzCellTypeSSize(wzCellType3Tetrahedron,0)	= 3;	// Side 0 
  wzCellTypeSSide(wzCellType3Tetrahedron,0,0)	= 3;
  wzCellTypeSSide(wzCellType3Tetrahedron,0,1)	= 1;
  wzCellTypeSSide(wzCellType3Tetrahedron,0,2)	= 2;
  wzCellTypeSLine(wzCellType3Tetrahedron,0,0)	= 2;
  wzCellTypeSLine(wzCellType3Tetrahedron,0,1)	= 0;
  wzCellTypeSLine(wzCellType3Tetrahedron,0,2)	= 1;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,0,0)	= 1;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,0,1)	= 2;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,0,2)	= 3;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,0,3)	= 1;
  wzCellTypeSSize(wzCellType3Tetrahedron,1)	= 3;	// Side 1 
  wzCellTypeSSide(wzCellType3Tetrahedron,1,0)	= 0;
  wzCellTypeSSide(wzCellType3Tetrahedron,1,1)	= 3;
  wzCellTypeSSide(wzCellType3Tetrahedron,1,2)	= 2;
  wzCellTypeSLine(wzCellType3Tetrahedron,1,0)	= 0;
  wzCellTypeSLine(wzCellType3Tetrahedron,1,1)	= 4;
  wzCellTypeSLine(wzCellType3Tetrahedron,1,2)	= 5;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,1,0)	= 3;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,1,1)	= 2;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,1,2)	= 0;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,1,3)	= 3;
  wzCellTypeSSize(wzCellType3Tetrahedron,2)	= 3;	// Side 2 
  wzCellTypeSSide(wzCellType3Tetrahedron,2,0)	= 3;
  wzCellTypeSSide(wzCellType3Tetrahedron,2,1)	= 0;
  wzCellTypeSSide(wzCellType3Tetrahedron,2,2)	= 1;
  wzCellTypeSLine(wzCellType3Tetrahedron,2,0)	= 3;
  wzCellTypeSLine(wzCellType3Tetrahedron,2,1)	= 1;
  wzCellTypeSLine(wzCellType3Tetrahedron,2,2)	= 5;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,2,0)	= 0;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,2,1)	= 1;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,2,2)	= 3;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,2,3)	= 0;
  wzCellTypeSSize(wzCellType3Tetrahedron,3)	= 3;	// Side 3 
  wzCellTypeSSide(wzCellType3Tetrahedron,3,0)	= 2;
  wzCellTypeSSide(wzCellType3Tetrahedron,3,1)	= 1;
  wzCellTypeSSide(wzCellType3Tetrahedron,3,2)	= 0;
  wzCellTypeSLine(wzCellType3Tetrahedron,3,0)	= 3;
  wzCellTypeSLine(wzCellType3Tetrahedron,3,1)	= 4;
  wzCellTypeSLine(wzCellType3Tetrahedron,3,2)	= 2;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,3,0)	= 1;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,3,1)	= 0;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,3,2)	= 2;
  wzCellTypeSPoint1(wzCellType3Tetrahedron,3,3)	= 1;

  for(t=0;t<=(int)wzCellTypeTypes;t++){
    switch(wzCellTypeCODIM(t)){
    case 0:
      wzCellTypesize[t] = wzCellTypepoints[t] + wzCellTypesides[t];
      break;
    case 1:
      wzCellTypesize[t] = wzCellTypepoints[t] + wzCellTypesides[t] + 2;
      break;
    default:
      wzCellTypesize[t] = wzCellTypepoints[t] + wzCellTypesides[t] + 1;
      break;
    }
    switch(wzCellTypeCDIM(t)){
    case 3:
      for(s=0;s<wzCellTypeSides(t);s++){
	for(e=0;e<wzCellTypeSSize(t,s);e++){
	  ee=wzCellTypeSLine(t,s,e);
	  if(wzCellTypeESideL(t,ee)==s){
	    wzAssert(wzCellTypeEPoint1(t,ee)==wzCellTypeSPoint1(t,s,e));
	    wzAssert(wzCellTypeEPoint2(t,ee)==wzCellTypeSPoint2(t,s,e));
	    wzAssert(wzCellTypeESideR(t,ee)==wzCellTypeSSide(t,s,e));
	  }else if(wzCellTypeESideR(t,ee)==s){
	    wzAssert(wzCellTypeEPoint2(t,ee)==wzCellTypeSPoint1(t,s,e));
	    wzAssert(wzCellTypeEPoint1(t,ee)==wzCellTypeSPoint2(t,s,e));
	    wzAssert(wzCellTypeESideL(t,ee)==wzCellTypeSSide(t,s,e));
	  }else	wzAssert(0);
	}}
      wzAssert(wzCellTypeSides(t)-wzCellTypeLines(t)+wzCellTypePoints(t)==2);
      break;
    case 2:
      wzAssert(wzCellTypeSides(t)==wzCellTypeLines(t));
      wzAssert(wzCellTypeSides(t)==wzCellTypePoints(t));
      break;
    case 1:
      break;
	case 0:
	  break;
    }
    wzAssert(wzCellTypeGDIM(t)-wzCellTypeCDIM(t)==(int)wzCellTypeCODIM(t));
  }
}