ibgotype.hxx

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

#ifndef ibgotype_hxx
#define ibgotype_hxx

#include "wzgrid.hxx"
#include "wzchart.hxx"
#include "wztensorgrid.hxx"

typedef ibgIndex ibgdLine;
typedef ibgIndex ibgdLineUp;
typedef ibgIndex ibgdLineSide;
typedef ibgIndex ibgdCorner;

const  ibgIndex ibgdMaxXDim    = 3;
const  ibgIndex ibgdMaxLines   = 6;
const  ibgIndex ibgdMaxLinesUp = 3;
const  ibgIndex ibgdMaxLineSides = 4;
const  ibgIndex ibgdMaxPlanes = 12;
const  ibgIndex ibgdMaxPlaneTypes = 3;
const  ibgIndex ibgdMaxPlaneCorners = 4;
const  ibgIndex ibgdMaxCorners = 8;
const  ibgIndex ibgdMaxCornerSides = 3;

class ibgOctreeType{
public:
  ibgOctreeType(ibgIndex dim);
  int	sOther[ibgdMaxLineSides+1];
  int	ptReverse[ibgdMaxPlaneCorners][ibgdMaxPlaneTypes];
  int	pReverse[ibgdMaxPlanes];
  int	pSide[3][ibgdMaxPlanes];
  int	pSide0[ibgdMaxPlanes];
  int	pSide1[ibgdMaxPlanes];
  int	pOrtho[ibgdMaxPlanes];
  int	csPlane[ibgdMaxCorners][ibgdMaxCornerSides];
  int	csPlaneSide[ibgdMaxCorners][ibgdMaxCornerSides];
  int	cReverse[ibgdMaxCorners];
  int	lsLine[ibgdMaxLines][ibgdMaxLineSides+1];
  int	llPlane[ibgdMaxLines][ibgdMaxLines];
  int	llOrtho[ibgdMaxLines][ibgdMaxLines];
  int	llSide[ibgdMaxLines][ibgdMaxLines];  // reverse to lsLine
  int	llCell[ibgdMaxLines][ibgdMaxLines];
  int	lsPlane[ibgdMaxLines][ibgdMaxLineSides+1];
  int	lsCell[ibgdMaxLines][ibgdMaxLineSides+1];
  ibgIndex	GridDim;
  static const ibgIndex	XDim;
  ibgIndex	Lines;
  ibgIndex	LinesUp;
  ibgIndex	LineSides;
  ibgIndex	PlaneSides;
  ibgIndex	PlaneCorners;
  ibgIndex	Corners;
  ibgIndex	Planes;
  ibgBoolean	isUp(ibgdLine e) 	const {return e<LinesUp;}
  ibgBoolean	isDown(ibgdLine e)	const {return e>=LinesUp;}
  ibgdLine	up(ibgdLine e)   	const {return e-LinesUp;}
  ibgdLine	down(ibgdLine e) 	const {return e+LinesUp;}
  ibgdLine 	Other(ibgdLine e) 	const {return isUp(e)?down(e):up(e);}
  ibgdLine 	Up(ibgdLine e) 		const {return isUp(e)?e:up(e);}
  ibgdLine 	Down(ibgdLine e) 	const {return isUp(e)?down(e):e;}
  ibgBoolean	llParallel(ibgdLine l1,ibgdLine l2)const 
  {return llOrtho[l1][l2]==Lines;}
  ibgBoolean	llOrthogonal(ibgdLine l1,ibgdLine l2)const 
  {return llOrtho[l1][l2]!=Lines;}
};

typedef const ibgOctreeType ibgConstOctreeType;

#endif