wzview.cxx

Delaunay三角形的网格剖分程序

//#include "ibgoutput.hxx"
#include "wzview.hxx"
#include "wzdefault.hxx"

// color settings;

static wzDefaultSegmentDescriptionClass black_background;
static wzDefaultSegmentDescriptionClass white_background;

static void defineColors()
{
  black_background.use();
  white_background.use();
  wzRegion(  ).color()="#aaa000000";
  wzRegion( 1).color()="#099011811";
  wzRegion( 2).color()="#302000";
  wzRegion( 3).color()="#000300a11";
  wzRegion( 4).color()="#700900900";
  wzRegion( 5).color()="#d11a00155";
  wzRegion( 6).color()="#311200211";
  wzRegion( 7).color()="#000d11100";
  wzRegion( 8).color()="#100300111";
  wzRegion( 9).color()="#111111111";
  wzRegion(10).color()="#100011211";
  wzRegion(11).color()="#400711911";
  wzRegion(12).color()="#400511911";
  wzRegion(13).color()="#400211911";
  wzRegion(14).color()="#600711911";
}

static int iso = 0;

static void defineTransformation(wzWindow win)
{
  wzFloat xmin[3],xmax[3],dd,dx;
  wzIndex i;
  if(wzCurrentData::grid){
    for(i=0;i<3;i++){
      xmin[i] = wzCurrentData::grid->cgmin(i);
      xmax[i] = wzCurrentData::grid->cgmax(i);
    }
  }else if (wzCurrentData::octree){
    for(i=0;i<3;i++){
      xmin[i] = wzCurrentData::octree->min[i];
      xmax[i] = wzCurrentData::octree->max[i];
    }
  }else{
    for(i=0;i<3;i++){
      xmin[i] = 0;
      xmax[i] = 1;
    }
  }    
  if(iso){
    for(i=0;i<3;i++){
      wzCurrentData::shift[i]  = xmin[i];
      wzCurrentData::factor[i] = dx = xmax[i] - xmin[i];
      xmax[i] = 1;
	if(dx<=0) {wzCurrentData::factor[i] = 1; xmax[i] = 0;}
    }
    wzwSetBaseVolume(win,0,xmax[0],0,xmax[1],0,xmax[2]);
  }else{
    dd = 0;
    for(i=0;i<3;i++){
      wzCurrentData::shift[i] = xmin[i];
      dx = xmax[i] - xmin[i];
      if(dd<dx) dd=dx;
    }
    for(i=0;i<3;i++){
      wzCurrentData::factor[i] = dd;
    }
    wzwSetBaseVolume(win,
		     0,(xmax[0]-xmin[0])/dd,
		     0,(xmax[1]-xmin[1])/dd,
		     0,(xmax[2]-xmin[2])/dd);
  }
}

static  wzWindow window = 0;

void v(){wzView();}

void wzView()
{
  char *title = "debug";
  int i=1;
  if(!window) window = wzwCreate(title,"help: F1",
				 "wzgridview.html#intro",0,0,600,600);
  defineColors();
  defineTransformation(window); 
  if(wzCurrentData::grid)
    wzGridView(window,i++,wzCurrentData::grid,"current grid","kikifax");
  if(wzCurrentData::geometry)
    cogeometryView(window,i++,wzCurrentData::geometry,"geometry","kikifax");
  if(wzCurrentData::octree)
    ibgOctreeView(window,i++,wzCurrentData::octree,"octree","kikifax");

  wzwSetDataValid(window,0);
  wzwView(window);
  //  wzwDestroy(window);
  //  window = 0;
}

void wzView(wzgrid g)
{
  char *title = "grid";
  if(!window) window = wzwCreate(title,"help: F1",
				 "wzgridview.html#intro",0,0,600,600);
  defineColors();
  defineTransformation(window);
  wzGridView(window,1,g,"grid","kikifax");
  wzwSetDataValid(window,0);
  wzwView(window);
  wzwDestroy(window);
  window = 0;
}

void wzView(cogeometry g)
{
  char *title = "geometry";
  if(!window) window = wzwCreate(title,"help: F1",
				 "cogview.html#intro",0,0,600,600);
  defineColors();
  defineTransformation(window);
  cogeometryView(window,1,g,"geometry","kikifax");
  wzwSetDataValid(window,0);
  wzwView(window);
  wzwDestroy(window);
  window = 0;
}