cog1d.hxx

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

#ifndef cog1d_hxx
#define cog1d_hxx

#ifndef cog_hxx
#include "cog.hxx"
#endif

class Cogeometry0D: public Cogeometry{
  cogRegion r;
public:
  Cogeometry0D(cogRegion r0=cogRegion(1))
    :r(r0){;}
  cogIndex   Point(cogPoint& p) const
  {p.segment()=r; return cogRCRegionFound;}
};

class Cogeometry1D: public Cogeometry{
  cogIndex num;
  cogIndex last;
  wzArray<cogFloat>  x;  // x[0],...,x[num-1]:	boundary positions;
  wzArray<cogRegion> r;  // r[0],...,r[num]:	regions before;
  wzArray<cogFace>   f;  // f[0],...,f[num-1]:	boundary faces;
  void setLast(cogIndex last) const;
public:
  Cogeometry1D(cogIndex num, cogFloat *x, cogRegion* r);
  Cogeometry1D(cogFloat xmin, cogFloat xmax, 
	       cogIndex num_of_regions_between = 1,
	       cogRegion first_inside  = cogRegion(2),
	       cogRegion outside       = cogRegion(1));
  Cogeometry1D(cogRegion r0=cogRegion(1),
	       cogFloat x1=0.0       ,cogRegion r1=cogRegion(2),
	       cogFloat x2=cogNoFloat,cogRegion r2=cogRegion(),
	       cogFloat x3=cogNoFloat,cogRegion r3=cogRegion(),
	       cogFloat x4=cogNoFloat,cogRegion r4=cogRegion(),
	       cogFloat x5=cogNoFloat,cogRegion r5=cogRegion(),
	       cogFloat x6=cogNoFloat,cogRegion r6=cogRegion(),
	       cogFloat x7=cogNoFloat,cogRegion r7=cogRegion(),
	       cogFloat x8=cogNoFloat,cogRegion r8=cogRegion(),
	       cogFloat x9=cogNoFloat,cogRegion r9=cogRegion());
  Cogeometry1D(cogIndex num, cogFloat *x, cogRegion* r, cogFace* f);
  Cogeometry1D(cogRegion s0,
	       cogFloat x1           ,cogRegion r1            ,cogFace f1          ,
	       cogFloat x2=cogNoFloat,cogRegion r2=cogRegion(),cogFace f2=cogFace(),
	       cogFloat x3=cogNoFloat,cogRegion r3=cogRegion(),cogFace f3=cogFace(),
	       cogFloat x4=cogNoFloat,cogRegion r4=cogRegion(),cogFace f4=cogFace(),
	       cogFloat x5=cogNoFloat,cogRegion r5=cogRegion(),cogFace f5=cogFace(),
	       cogFloat x6=cogNoFloat,cogRegion r6=cogRegion(),cogFace f6=cogFace(),
	       cogFloat x7=cogNoFloat,cogRegion r7=cogRegion(),cogFace f7=cogFace(),
	       cogFloat x8=cogNoFloat,cogRegion r8=cogRegion(),cogFace f8=cogFace(),
	       cogFloat x9=cogNoFloat,cogRegion r9=cogRegion(),cogFace f9=cogFace());
  // The region "rBefore" ends at x. After x we have the old region of x.
  void before(cogFloat x, cogRegion rBefore, cogFace face=cogFace());
  // The region "rAfter" begins at x. Before x we have the old region of x.
  void after( cogFloat x, cogRegion rAfter,  cogFace face=cogFace());
  // The replaced virtual functions of Cogeometry:
  cogIndex   Point(cogPoint& s) const;
  cogIndex   Line (cogFlag1& f, const cogLine& s) const;
};
#endif