wzmaps.hxx

Delaunay三角形的网格剖分程序

#ifndef wzmaps_hxx
#define wzmaps_hxx

#include "wzmap.hxx"

// 1D coordinates:

class wz1Dexp;
class wz1Dlog;
class wz1Dsinh;
class wz1Datan;

// 2D coordinates:

class wz2Dpolar;
class wz2Dexp;
class wz2Dlog;
class wz2Dcosh;
class wz2Dacosh;
class wz2Dshukovski;
class wz2DroundedBox;

// 3D coordinates:

class wz3Dspherical;
class wz3Dcylinder;

// full declarations:

class wz1Dexp
  : public wz1Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz1Dexp(wz1dmap m=0);
};

class wz1Dlog
  : public wz1Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz1Dlog(wz1dmap m=0);
};

class wz1Dsinh
  : public wz1Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz1Dsinh(wz1dmap m=0);
};

class wz1Datan
  : public wz1Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz1Datan(wz1dmap m=0);
};

class wz2Dpolar
  : public wz2Dcoordinates
{
  wzFloat Epsilon;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Dpolar(wz2dmap m=0);
  void setMinimalRadius(wzFloat eps);
};

class wz2Delliptic
  : public wz2Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Delliptic(wz2dmap m=0);
};

class wz2Dexp
  : public wz2Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Dexp(wz2dmap m=0);
};

class wz2Dlog
  : public wz2Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Dlog(wz2dmap m=0);
};

class wz2Dcosh
  : public wz2Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Dcosh(wz2dmap m=0);
};
class wz2Dacosh
  : public wz2Dcoordinates
{
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Dacosh(wz2dmap m=0);
};

class wz2Dshukovski
  : public wz2Dcoordinates
{
  wzFloat yc;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  void setYofCircleCenter(wzFloat y=0);
  wz2Dshukovski(wz2dmap m=0,wzFloat y=wzInfty);
};

class wz2DshiftShukovski
  : public wz2Dcoordinates
{
  wzFloat lambda,rr,m0,m1,d0,d1,e0,e1,c0,c1;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  void setYofCircleCenter(wzFloat y=0);
  void setShiftPoints(wzFloat x0,wzFloat y0,wzFloat x1,wzFloat y1);
  wz2DshiftShukovski(wz2dmap m=0,wzFloat y=0);
};

class wz2DroundedBox
  : public wz2Dcoordinates
{
  wzFloat lx,ly,st1,st2,st3,st4,st5,up,phi;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2DroundedBox( wzFloat xl, wzFloat yl, wzFloat rr);
};

class wz3Dspherical
  : public wz3Dcoordinates
{
  wzFloat Epsilon;
  wzFloat Epsilon1;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz3Dspherical();
  void setMinimalRadius(wzFloat eps, wzFloat eps1);
};

class wz3Dcylinder
  : public wz3Dcoordinates
{
  wzFloat Epsilon;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz3Dcylinder();
  void setMinimalRadius(wzFloat eps);
};

class wz3Dcylindrical
  : public wz3Dcoordinates
{
  wz2dmap   xy;
  wz1dmap   z;
  wzMap*    XY;
  wzMap*    Z;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz3Dcylindrical(wz2dmap xy=0, wz1dmap z=0);
  void composeXY(wz2dmap z=0);
  void composeZ(wz1dmap z=0);
private:
  void compose(wz3dmap z=0);
};

class wz2Drectangular
  : public wz2Dcoordinates
{
  wz1dmap   x;
  wz1dmap   y;
  wzMap*    xx;
  wzMap*    yy;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz2Drectangular(wz1dmap x=0, wz1dmap y=0);
  void composeX(wz1dmap x=0);
  void composeY(wz1dmap y=0);
private:
  void compose(wz2dmap x=0);
};

class wz3Drectangular
  : public wz3Dcoordinates
{
  wz1dmap   x;
  wz1dmap   y;
  wz1dmap   z;
  wzMap*    xx;
  wzMap*    yy;
  wzMap*    zz;
public:
  void operator()(wzFloat *x, const wzFloat *u) const;
  void inverse   (wzFloat *u, const wzFloat *x) const;
  wz3Drectangular(wz1dmap x=0, wz1dmap y=0, wz1dmap z=0);
  void composeX(wz1dmap x=0);
  void composeY(wz1dmap y=0);
  void composeZ(wz1dmap z=0);
private:
  void compose(wz3dmap x=0);
};

#ifdef  wzUseObsolete
#define wzcoordinates wzmap
#define wzScaleExp wz1Dexp
#define wzScaleLog wz1Dlog
#define wzPolarCoordinates wz2Dpolar
#define wzCylinderCoordinates wz3Dcylinder
#define wzSphericalCoordinates wz3Dspherical
#endif

#endif