xginitc.c

xgrafix 是PTSG模拟程序中的图形截面库 改版本是最新版本

  else if (!strcmp(PlotType, "loglinlin")) plottype = LOG_LIN_LIN;
  else if (!strcmp(PlotType, "loglinlog")) plottype = LOG_LIN_LOG;
  else if (!strcmp(PlotType, "logloglin")) plottype = LOG_LOG_LIN;
  else if (!strcmp(PlotType, "logloglog")) plottype = LOG_LOG_LOG;
  else {
    printf("Unrecognized plot string '%s' for Window '%s'.\n", PlotType, Z_Label);
    exit(-1);
  }

  XGSetupWindow(Z_Label, State, ulx, uly, THREED, data_ptr, label_ptr, 
		plottype);
}

void XGSet3DFlag(char *PlotType, char *X_Label, char *Y_Label, char *Z_Label,
		 SCALAR Theta, SCALAR Phi, char *State, int ulx, int uly,
		 SCALAR X_Scale, SCALAR Y_Scale, SCALAR Z_Scale, 
		 int X_Auto_Rescale, int Y_Auto_Rescale, int Z_Auto_Rescale,
		 SCALAR X_Min, SCALAR X_Max, SCALAR Y_Min, SCALAR Y_Max, 
		 SCALAR Z_Min, SCALAR Z_Max, int *openFlag)
{
  XGSet3D(PlotType, X_Label, Y_Label, Z_Label, Theta, Phi, State, ulx, uly,
	  X_Scale, Y_Scale, Z_Scale, X_Auto_Rescale, Y_Auto_Rescale,
	  Z_Auto_Rescale, X_Min, X_Max, Y_Min, Y_Max, Z_Min, Z_Max);
  
  theWindowArray[numberOfWindows-1]->openFlag = openFlag;
}

void XGSet2D(char *PlotType, char *X_Label, char *Y_Label, char *State,
	     int ulx, int uly, SCALAR X_Scale, SCALAR Y_Scale, 
	     int X_Auto_Rescale, int Y_Auto_Rescale, SCALAR X_Min, 
	     SCALAR X_Max, SCALAR Y_Min, SCALAR Y_Max)
{
  int             plottype;
  
  DataType        data_ptr = NULL;
  LabelType       label_ptr;
  
  label_ptr = SetupLabelStruct();
  label_ptr->Y_Label = strdup(Y_Label);
  label_ptr->Y_Min = Y_Min;
  label_ptr->Y_Max = Y_Max;
  label_ptr->Y_Scale = Y_Scale;
  label_ptr->Y_Auto_Rescale = Y_Auto_Rescale;

  label_ptr->X_Label = strdup(X_Label);
  label_ptr->X_Min = X_Min;
  label_ptr->X_Max = X_Max;
  label_ptr->X_Scale = X_Scale;
  label_ptr->X_Auto_Rescale = X_Auto_Rescale;
  
  if (!strcmp(PlotType, "linlin"))      plottype = LIN_LIN;
  else if (!strcmp(PlotType, "linlog")) plottype = LIN_LOG;
  else if (!strcmp(PlotType, "loglog")) plottype = LOG_LOG;
  else if (!strcmp(PlotType, "loglin")) plottype = LOG_LIN;
  else {
    printf("Unrecognized plot string '%s' for Window '%s'.\n", PlotType, Y_Label);
    exit(-1);
  }
  
  XGSetupWindow(Y_Label, State, ulx, uly, TWOD, data_ptr, label_ptr, plottype);
}

void XGSet2D_OS(char *PlotType, char *X_Label, char *Y_Label, char *State,
	     int ulx, int uly, SCALAR X_Offset, SCALAR Y_Offset,
             SCALAR X_Scale, SCALAR Y_Scale, 
	     int X_Auto_Rescale, int Y_Auto_Rescale, SCALAR X_Min, 
	     SCALAR X_Max, SCALAR Y_Min, SCALAR Y_Max)
{
  int             plottype;
  
  DataType        data_ptr = NULL;
  LabelType       label_ptr;
  
  label_ptr = SetupLabelStruct();
  label_ptr->Y_Label = strdup(Y_Label);
  label_ptr->Y_Min = Y_Min;
  label_ptr->Y_Max = Y_Max;
  label_ptr->Y_Scale = Y_Scale;
  label_ptr->Y_Offset = Y_Offset;
  label_ptr->Y_Auto_Rescale = Y_Auto_Rescale;

  label_ptr->X_Label = strdup(X_Label);
  label_ptr->X_Min = X_Min;
  label_ptr->X_Max = X_Max;
  label_ptr->X_Scale = X_Scale;
  label_ptr->X_Offset = X_Offset;
  label_ptr->X_Auto_Rescale = X_Auto_Rescale;
  
  if (!strcmp(PlotType, "linlin"))      plottype = LIN_LIN;
  else if (!strcmp(PlotType, "linlog")) plottype = LIN_LOG;
  else if (!strcmp(PlotType, "loglog")) plottype = LOG_LOG;
  else if (!strcmp(PlotType, "loglin")) plottype = LOG_LIN;
  else {
    fprintf(stderr, "Unrecognized plot string '%s' for Window '%s'.\n", PlotType, Y_Label);
    exit(-1);
  }
  
  XGSetupWindow(Y_Label, State, ulx, uly, TWOD, data_ptr, label_ptr, plottype);
}

void XGSet2DFlag(char *PlotType, char *X_Label, char *Y_Label, char *State,
		 int ulx, int uly, SCALAR X_Scale, SCALAR Y_Scale, 
		 int X_Auto_Rescale, int Y_Auto_Rescale, SCALAR X_Min, 
		 SCALAR X_Max, SCALAR Y_Min, SCALAR Y_Max, int *openFlag)
{
  XGSet2D(PlotType, X_Label, Y_Label, State, ulx, uly, X_Scale, Y_Scale,
	  X_Auto_Rescale, Y_Auto_Rescale, X_Min, X_Max, Y_Min, Y_Max);

  theWindowArray[numberOfWindows-1]->openFlag = openFlag;
}

void XGSet2DFlag_OS(char *PlotType, char *X_Label, char *Y_Label, char *State,
		 int ulx, int uly, SCALAR X_Offset, SCALAR Y_Offset,
		 SCALAR X_Scale, SCALAR Y_Scale, 
		 int X_Auto_Rescale, int Y_Auto_Rescale, SCALAR X_Min, 
		 SCALAR X_Max, SCALAR Y_Min, SCALAR Y_Max, int *openFlag)
{
  XGSet2D_OS(PlotType, X_Label, Y_Label, State, ulx, uly, X_Offset, Y_Offset,
	  X_Scale, Y_Scale,
	  X_Auto_Rescale, Y_Auto_Rescale, X_Min, X_Max, Y_Min, Y_Max);

  theWindowArray[numberOfWindows-1]->openFlag = openFlag;
}

void XGSet2DC(char *PlotType, char *X_Label, char *Y_Label, char *Z_Label,
	      char *State, int ulx, int uly, SCALAR X_Scale, SCALAR Y_Scale,
	      SCALAR Z_Scale, int X_Auto_Rescale, int Y_Auto_Rescale,
	      int Z_Auto_Rescale, SCALAR X_Min, SCALAR X_Max, SCALAR Y_Min, 
	      SCALAR Y_Max, SCALAR Z_Min, SCALAR Z_Max)
{
  XGSet3D(PlotType, X_Label, Y_Label, Z_Label, 0.0, 0.0, State, ulx, uly,
	  X_Scale, Y_Scale, Z_Scale, X_Auto_Rescale, Y_Auto_Rescale,
	  Z_Auto_Rescale, X_Min, X_Max, Y_Min, Y_Max, Z_Min, Z_Max);
}

void XGSet2DCFlag(char *PlotType, char *X_Label, char *Y_Label, char *Z_Label,
		  char *State, int ulx, int uly, SCALAR X_Scale,SCALAR Y_Scale,
		  SCALAR Z_Scale, int X_Auto_Rescale, int Y_Auto_Rescale,
		  int Z_Auto_Rescale, SCALAR X_Min, SCALAR X_Max, SCALAR Y_Min,
		  SCALAR Y_Max, SCALAR Z_Min, SCALAR Z_Max, int *openFlag)
{
  XGSet3DFlag(PlotType, X_Label, Y_Label, Z_Label, 0.0, 0.0, State, ulx, uly,
	      X_Scale, Y_Scale, Z_Scale, X_Auto_Rescale, Y_Auto_Rescale,
	      Z_Auto_Rescale, X_Min, X_Max, Y_Min, Y_Max, Z_Min, Z_Max, 
	      openFlag);
}

void XGCurve(SCALAR *x_array, SCALAR *y_array, int *npoints, int color)
{
  DataType        data_prevptr, data_ptr;

  data_prevptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr = SetupDataStruct();
  data_ptr->y = y_array;
  data_ptr->x = x_array;
  data_ptr->npoints = npoints;
  data_ptr->color = color;
  data_ptr->next = data_prevptr;

  theWindowArray[numberOfWindows-1]->data = data_ptr;
}

void XGCurveVector(SCALAR *x_array, SCALAR *y_array, int *npoints, int color,
		    int xSize, int xOffset, int ySize, int yOffset)
{
  DataType data_ptr;

  XGCurve(x_array, y_array, npoints, color);

  data_ptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr->xArraySize = xSize;
  data_ptr->yArraySize = ySize;
  data_ptr->xArrayOffset = xOffset;
  data_ptr->yArrayOffset = yOffset;
}

void XGScat2D(SCALAR *x_array, SCALAR *y_array, int *npoints, int color)
{
  DataType        data_prevptr, data_ptr;
  
  data_prevptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr = SetupDataStruct();
  data_ptr->y = y_array;
  data_ptr->x = x_array;
  data_ptr->npoints = npoints;
  data_ptr->color = color;
  data_ptr->next = data_prevptr;

  Set_Scatter(data_ptr);

  theWindowArray[numberOfWindows-1]->data = data_ptr;
}

void XGScat3D(SCALAR *x_array, SCALAR *y_array, SCALAR *z_array,
	      int *npoints, int color)
{
  DataType        data_prevptr, data_ptr;
  
  data_prevptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr = SetupDataStruct();
  data_ptr->zscat = z_array;
  data_ptr->y = y_array;
  data_ptr->x = x_array;
  data_ptr->npoints = npoints;
  data_ptr->mpoints = npoints;
  data_ptr->color = Color[color];
  data_ptr->next = data_prevptr;

  Set_Scatter(data_ptr);

  theWindowArray[numberOfWindows-1]->data = data_ptr;
}

void XGVector(SCALAR *x_array, SCALAR *y_array, SCALAR **z_array, 
	      SCALAR **w_array, int *mpoints, int *npoints, int color)
{
  DataType        data_prevptr, data_ptr;
  
  data_prevptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr = SetupDataStruct();
  data_ptr->w = z_array; /* VV: notice these are switched */
  data_ptr->z = w_array;
  data_ptr->y = y_array;
  data_ptr->x = x_array;
  data_ptr->npoints = npoints;
  data_ptr->mpoints = mpoints;
  data_ptr->color = color;
  data_ptr->next = data_prevptr;

  Set_Vector(data_ptr);

  theWindowArray[numberOfWindows-1]->data = data_ptr;
}

void XGVectorVector(SCALAR *x_array, SCALAR *y_array, SCALAR **z_array, 
		    SCALAR **w_array, int *mpoints, int *npoints, int color,
		    int xSize, int xOffset, int ySize, int yOffset, int zSize,
		    int zOffset, int wSize, int wOffset)
{
  DataType        data_ptr;
  
  XGVector(x_array, y_array, z_array, w_array, mpoints, npoints, color);

  data_ptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr->xArraySize = xSize;
  data_ptr->yArraySize = ySize;
  data_ptr->xArrayOffset = xOffset;
  data_ptr->yArrayOffset = yOffset;
  data_ptr->zArraySize = wSize; /*DKW: these are switched as well */
  data_ptr->wArraySize = zSize;
  data_ptr->zArrayOffset = wOffset;
  data_ptr->wArrayOffset = zOffset;
}

void XGSurf(SCALAR *x_array, SCALAR *y_array, SCALAR **z_array,
	    int *mpoints, int *npoints, int color)
{
  DataType        data_prevptr, data_ptr;
  
  data_prevptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr = SetupDataStruct();
  data_ptr->z = z_array;
  data_ptr->y = y_array;
  data_ptr->x = x_array;
  data_ptr->mpoints = mpoints;
  data_ptr->npoints = npoints;
  data_ptr->color = Color[color];
  data_ptr->next = data_prevptr;

  theWindowArray[numberOfWindows-1]->data = data_ptr;
}

void XGSurfVector(SCALAR *x_array, SCALAR *y_array, SCALAR **z_array,
		  int *mpoints, int *npoints, int color, int xSize, 
		  int xOffset, int ySize, int yOffset, int zSize, int zOffset)
{
  DataType        data_ptr;
  
  XGSurf(x_array, y_array, z_array, mpoints, npoints, color);

  data_ptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr->xArraySize = xSize;
  data_ptr->yArraySize = ySize;
  data_ptr->zArraySize = zSize;
  data_ptr->xArrayOffset = xOffset;
  data_ptr->yArrayOffset = yOffset;
  data_ptr->zArrayOffset = zOffset;
}

void XGCont(SCALAR *x_array, SCALAR *y_array, SCALAR **z_array,
	    int *mpoints, int *npoints, int color)
{ 
  XGSurf(x_array, y_array, z_array, mpoints, npoints, color);
}

void XGIRSurf(SCALAR **x_array, SCALAR **y_array, SCALAR **z_array,
	      int *mpoints, int *npoints, int color)
{
  DataType        data_prevptr, data_ptr;
  
  data_prevptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr = SetupDataStruct();
  data_ptr->z = z_array;
  data_ptr->y_irr = y_array;
  data_ptr->x_irr = x_array;
  data_ptr->mpoints = mpoints;
  data_ptr->npoints = npoints;
  data_ptr->color = Color[color];
  data_ptr->next = data_prevptr;

  Set_Irr_Surf(data_ptr);

  theWindowArray[numberOfWindows-1]->data = data_ptr;
}

void XGIRSurfVector(SCALAR **x_array, SCALAR **y_array, SCALAR **z_array,
		    int *mpoints, int *npoints, int color, int xSize, 
		    int xOffset, int ySize, int yOffset, int zSize, 
		    int zOffset)
{
  DataType        data_ptr;
  
  XGIRSurf(x_array, y_array, z_array, mpoints, npoints, color);

  data_ptr = theWindowArray[numberOfWindows-1]->data;
  data_ptr->xArraySize = xSize;
  data_ptr->yArraySize = ySize;
  data_ptr->zArraySize = zSize;
  data_ptr->xArrayOffset = xOffset;
  data_ptr->yArrayOffset = yOffset;
  data_ptr->zArrayOffset = zOffset;
}