wz1dmaps.cxx

Delaunay三角形的网格剖分程序

#include "wzoutput.hxx"
#include "wzmaps.hxx"
#include <math.h>

void wz1Dcoordinates::operator()(wzFloat *x, const wzFloat *u) const
{
  wzMapPreprocessing(x,u);
  X[0]=U[0];
  wzMapPostprocessing(x);
}

void wz1Dcoordinates::inverse(wzFloat *u, const wzFloat *x) const
{
  wzMapInversePreprocessing(u,x);
  U[0]=X[0];
  wzMapInversePostprocessing(u);
}

void wz1Dlog::inverse(wzFloat *u, const wzFloat *x) const
{
  wzMapInversePreprocessing(u,x);
  U[0] = exp(X[0]);
  wzMapInversePostprocessing(u);
}
void wz1Dlog::operator()(wzFloat *x, const wzFloat *u) const
{
  wzMapPreprocessing(x,u);
  X[0] = log(U[0]);
  wzMapPostprocessing(x);
}

void wz1Dexp::operator()(wzFloat *x, const wzFloat *u) const
{
  wzMapPreprocessing(x,u);
  X[0] = exp(U[0]);
  wzMapPostprocessing(x);
}
void wz1Dexp::inverse(wzFloat *u, const wzFloat *x) const
{
  wzMapInversePreprocessing(u,x);
  U[0] = log(X[0]);
  wzMapInversePostprocessing(u);
}
void wz1Dsinh::operator()(wzFloat *x, const wzFloat *u) const
{
  wzMapPreprocessing(x,u);
  X[0] = sinh(U[0]);
  wzMapPostprocessing(x);
}

void wz1Dsinh::inverse(wzFloat *u, const wzFloat *x) const
{
  wzMapInversePreprocessing(u,x);
  U[0] = asinh(X[0]);
  wzMapInversePostprocessing(u);
}

void wz1Datan::operator()(wzFloat *x, const wzFloat *u) const
{
  wzMapPreprocessing(x,u);
  X[0] = atan(U[0]);
  wzMapPostprocessing(x);
}

void wz1Datan::inverse(wzFloat *u, const wzFloat *x) const
{
  wzMapInversePreprocessing(u,x);
  U[0] = tan(X[0]);
  wzMapInversePostprocessing(u);
}

wz1Dexp::wz1Dexp(wz1dmap m)
  :wz1Dcoordinates(m)
{
  name = "exponential function";
  image.setMinimum(0);
}
wz1Dsinh::wz1Dsinh(wz1dmap m)
  :wz1Dcoordinates(m)
{
  name = "sinus hyperbolicus";
}

wz1Datan::wz1Datan(wz1dmap m)
  :wz1Dcoordinates(m)
{
  name = "arcus tangens";
  image.setBorder(-wzPi/2,wzPi/2);
}

wz1Dlog::wz1Dlog(wz1dmap m)
  :wz1Dcoordinates(m)
{
  name = "logarithm";
  full.setMinimum(0);
  chart.setMinimum(0);
}