Plo2dn.c File Reference

#include <string.h>
#include <stdio.h>
#include <math.h>
#include "math_graphics.h"
#include "PloEch.h"
#include "../../sparse/includes/spConfig.h"
#include "SetProperty.h"
#include "GetProperty.h"
#include "InitObjects.h"
#include "DrawObjects.h"
#include "BuildObjects.h"
#include "Axes.h"
#include "BasicAlgos.h"
#include "clipping.h"
#include "sciprint.h"
#include "CurrentObjectsManagement.h"
#include "MALLOC.h"

Include dependency graph for Plo2dn.c:

Go to the source code of this file.

Defines
#define	spINSIDE_SPARSE
Functions
void	initsubwin ()
void	compute_data_bounds (int cflag, char dataflag, double *x, double *y, integer n1, integer n2, double *drect)
void	compute_data_bounds2 (int cflag, char dataflag, char *logflags, double *x, double *y, integer n1, integer n2, double *drect)
BOOL	update_specification_bounds (sciPointObj *psubwin, double *rect, int flag)
int	re_index_brect (double *brect, double *drect)
BOOL	strflag2axes_properties (sciPointObj *psubwin, char *strflag)
char **	FreeUserLabels (char **u_xlabels, int *u_nxgrads)
double *	FreeUserGrads (double *u_xgrads)
double *	AllocUserGrads (double *u_xgrads, int nb)
int	CopyUserGrads (double *u_xgrad_SRC, double *u_xgrad_DEST, int dim)
char **	AllocAndSetUserLabels (char **u_xlabels, double *u_xgrads, int u_nxgrads, char logflag)
char **	AllocAndSetUserLabelsFromMdl (char **u_xlabels, char **u_xlabels_MDL, int u_nxgrads)
int	CreatePrettyGradsFromNax (sciPointObj *psubwin, int *Nax)
int	GraduateWithNax (sciSubWindow *ppsubwin, double *min, double *max, int nbtics, double *grads)
int	ChooseGoodFormat (char *c_format, char logflag, double *_grads, int n_grads)
int	plot2dn (integer ptype, char *logflags, double *x, double *y, integer *n1, integer *n2, integer *style, char *strflag, char *legend, double *brect, integer *aaint, BOOL flagNax, integer lstr1, integer lstr2)
BOOL	update_specification_bounds (sciPointObj *psubwin, rect, int flag)

---

Define Documentation

#define spINSIDE_SPARSE

Definition at line 12 of file Plo2dn.c.

---

Function Documentation

char ** AllocAndSetUserLabels	(	char **	u_xlabels,
		double *	u_xgrads,
		int	u_nxgrads,
		char	logflag
	)

Definition at line 732 of file Plo2dn.c.

References ChooseGoodFormat(), foo(), i, MALLOC, NULL, and sciprint().

Referenced by CreatePrettyGradsFromNax().

{
  int i;
  char c_format[5]; 
  int nbtics = u_nxgrads;

  if(u_xgrads == NULL)
    return (char **) NULL;
  
  if(u_xlabels != NULL)
    {
      sciprint("Impossible: u_xlabels must be freed before re-allocating");
      return (char **) NULL;
    }
  

  if((u_xlabels=(char **)MALLOC(u_nxgrads*sizeof(char *)))==NULL){
    sciprint("No more place for allocating user labels using Nax");
    return (char **) NULL;
  }

  ChooseGoodFormat(c_format,logflag,u_xgrads,u_nxgrads);
  
  for(i=0;i<nbtics;i++)
    {  
      char foo[100];
      
      sprintf(foo,c_format, u_xgrads[i]);
      
      if((u_xlabels[i]=(char *)MALLOC((strlen(foo)+1)*sizeof(char )))==NULL){
        sciprint("No more place for allocating u_xlabels");
        return (char **) NULL;
      }
      
      strcpy(u_xlabels[i],foo);
    }
  
  return u_xlabels;
}

Here is the call graph for this function:

Here is the caller graph for this function:

char** AllocAndSetUserLabelsFromMdl	(	char **	u_xlabels,
		char **	u_xlabels_MDL,
		int	u_nxgrads
	)

double* AllocUserGrads	(	double *	u_xgrads,
		int	nb
	)

int ChooseGoodFormat	(	char *	c_format,
		char	logflag,
		double *	_grads,
		int	n_grads
	)

Definition at line 1199 of file Format.c.

References ChoixFormatE().

{
  int last_index = n_grads - 1;

  if(logflag == 'l')
  {
    ChoixFormatE(c_format,
      exp10(_grads[0]),
      exp10(_grads[last_index]),
      (( exp10(_grads[last_index]))-( exp10(_grads[0])))/(last_index));
  }
  else
  {
    ChoixFormatE(c_format,
      _grads[0],
      _grads[last_index],
      ((_grads[last_index])-(_grads[0]))/(last_index)); /* Adding F.Leray 06.05.04 */
  }

  return 0;

}

Here is the call graph for this function:

void compute_data_bounds	(	int	cflag,
		char	dataflag,
		double *	x,
		double *	y,
		integer	n1,
		integer	n2,
		double *	drect
	)

Definition at line 369 of file Plo2dn.c.

References Maxi(), Mini(), and xd.

{
  int size_x,size_y;
  double xd[2];
  double *x1;
  switch ( dataflag ) {
  case 'e' : 
    xd[0] = 1.0; xd[1] = (double)n2;
    x1 = xd;size_x = (n2 != 0) ? 2 : 0 ;
    break; 
  case 'o' : 
    x1 = x;size_x = n2;
    break;
  case 'g' : 
  default  : 
    x1 = x;size_x = (cflag == 1) ? n1 : (n1*n2) ;
    break; 
  }

  if (size_x != 0) {
    drect[0] =  Mini(x1, size_x); 
    drect[1] =  Maxi(x1,size_x); 
  }
  else {
    drect[0] = 0.0;
    drect[1] = 10.0;
  }

  size_y = (cflag == 1) ? n2 : (n1*n2) ;
  if (size_y != 0) {
    drect[2] =  Mini(y, size_y); 
    drect[3] =  Maxi(y,size_y); 
  }
  else {
    drect[2] = 0.0;
    drect[3] = 10.0;
  }
  /* back to default values for  x=[] and y = [] */
  if ( drect[2] == LARGEST_REAL ) { drect[2] = 0.0; drect[3] = 10.0 ;} 
  if ( drect[0] == LARGEST_REAL ) { drect[0] = 0.0; drect[1] = 10.0 ;} 

}

Here is the call graph for this function:

void compute_data_bounds2	(	int	cflag,
		char	dataflag,
		char *	logflags,
		double *	x,
		double *	y,
		integer	n1,
		integer	n2,
		double *	drect
	)

Definition at line 294 of file Plo2dn.c.

References Maxi(), Mini(), sciFindStPosMin(), and xd.

{
  int size_x,size_y;
  double xd[2];
  double *x1;
  switch ( dataflag ) {
  case 'e' : 
    xd[0] = 1.0; xd[1] = (double)n2;
    x1 = xd;size_x = (n2 != 0) ? 2 : 0 ;
    break; 
  case 'o' : 
    x1 = x;size_x = n2;
    break;
  case 'g' : 
  default  : 
    x1 = x;size_x = (cflag == 1) ? n1 : (n1*n2) ;
    break; 
  }

  if (size_x != 0) {
    if(logflags[0] != 'l'){
      drect[0] =  Mini(x1, size_x); 
      drect[1] =  Maxi(x1,size_x); 
    }
    else { /* log. case */
      drect[0] =  sciFindStPosMin(x1,size_x); 
      drect[1] =  Maxi(x1,size_x); 
    }
    
  }
  else {
    if(logflags[0] != 'l'){
      drect[0] = 0.0;
      drect[1] = 10.0;
    }
    else{/* log. case */
      drect[0] = 1.0;
      drect[1] = 10.0;
    }
  }

  size_y = (cflag == 1) ? n2 : (n1*n2) ;
  if (size_y != 0) {
    if(logflags[1] != 'l'){
      drect[2] =  Mini(y, size_y); 
      drect[3] =  Maxi(y,size_y); 
    }
    else{/* log. case */
      drect[2] =  sciFindStPosMin(y,size_y); 
      drect[3] =  Maxi(y,size_y); 
    }

  }
  else {
    if(logflags[1] != 'l'){
      drect[2] = 0.0;
      drect[3] = 10.0;
    }
    else{/* log. case */
      drect[2] = 1.0;
      drect[3] = 10.0;
    }
  }
  /* back to default values for  x=[] and y = [] */
  if ( drect[2] == LARGEST_REAL ) { drect[2] = 0.0; drect[3] = 10.0 ;} 
  if ( drect[0] == LARGEST_REAL ) { drect[0] = 0.0; drect[1] = 10.0 ;} 

}

Here is the call graph for this function:

int CopyUserGrads	(	double *	u_xgrad_SRC,
		double *	u_xgrad_DEST,
		int	dim
	)

int CreatePrettyGradsFromNax	(	sciPointObj *	psubwin,
		int *	Nax
	)

double * FreeUserGrads

(

double *

u_xgrads

)

Definition at line 690 of file Plo2dn.c.

References FREE, and NULL.

Referenced by CreatePrettyGradsFromNax().

{
  FREE(u_xgrads); u_xgrads = NULL;
  return u_xgrads;
}

Here is the caller graph for this function:

char** FreeUserLabels	(	char **	u_xlabels,
		int *	u_nxgrads
	)

int GraduateWithNax	(	sciSubWindow *	ppsubwin,
		double *	min,
		double *	max,
		int	nbtics,
		double *	grads
	)

Definition at line 654 of file Plo2dn.c.

References i, and pas.

Referenced by CreatePrettyGradsFromNax().

{
  int i;
  double pas;
  
  if(nbtics == 1){
    pas = 0.;
    grads[0] = (*min);
  }
  else{
    pas = (*max - *min) / (nbtics -1);
    
    for(i=0;i<nbtics;i++) 
      grads[i] = (*min) + pas*i;
  }

  return 0;
}

Here is the caller graph for this function:

void initsubwin

(

)

Definition at line 431 of file InitObjects.c.

References sciSubWindow::axes, BT_OFF, clearUserData(), Cscale2default(), FALSE, initSubWinSize(), AXES::limits, pSUBWIN_FEATURE, AXES::rect, reinitSubWin(), sciGetCurrentSubWin(), AXES::subint, and AXES::ticscolor.

{
  sciPointObj  * psubwin  ;
  sciSubWindow * ppSubWin ;
  

  Cscale2default(); 
  psubwin = sciGetCurrentSubWin();
  ppSubWin = pSUBWIN_FEATURE (psubwin) ;


  initSubWinSize( psubwin  ) ;

  clearUserData( psubwin ) ;

  reinitSubWin( psubwin ) ;
  
  (ppSubWin->axes).axes_visible[0] = FALSE;
  (ppSubWin->axes).axes_visible[1] = FALSE;
  (ppSubWin->axes).axes_visible[2] = FALSE;

  (ppSubWin->axes).reverse[0] = FALSE;
  (ppSubWin->axes).reverse[1] = FALSE;
  (ppSubWin->axes).reverse[2] = FALSE;
  
  ppSubWin->axes.rect = BT_OFF;  
  ppSubWin->axes.ticscolor = -1;
  ppSubWin->axes.subint[0] =  1;
  ppSubWin->axes.subint[1] =  1;
  ppSubWin->axes.subint[2] =  1;
  ppSubWin->axes.limits[0]  = 0; 

}

Here is the call graph for this function:

int plot2dn	(	integer	ptype,
		char *	logflags,
		double *	x,
		double *	y,
		integer *	n1,
		integer *	n2,
		integer *	style,
		char *	strflag,
		char *	legend,
		double *	brect,
		integer *	aaint,
		BOOL	flagNax,
		integer	lstr1,
		integer	lstr2
	)

Definition at line 90 of file Plo2dn.c.

References sciSubWindow::alpha, sciSubWindow::alpha_kp, AXES::auto_ticks, sciSubWindow::axes, checkRedrawing(), compute_data_bounds2(), ConstructCompound(), ConstructLegend(), ConstructPolyline(), CreatePrettyGradsFromNax(), DrawAxesIfRequired(), FALSE, sciSubWindow::FirstPlot, sciSubWindow::flagNax, frame_clip_off(), frame_clip_on(), FREE, sciSubWindow::is3d, sciSubWindow::logflags, MALLOC, Max, Merge3d(), Min, NULL, PD0, sciSubWindow::project, pSUBWIN_FEATURE, re_index_brect(), sciDrawObj(), sciDrawObjIfRequired(), sciGetCurrentFigure(), sciGetCurrentObj(), sciGetCurrentScilabXgc(), sciGetCurrentSubWin(), sciGetGraphicMode(), sciGetHandle(), sciGetSurface(), sciprint(), sciSetCurrentObj(), sciSetIsClipping(), sciSubWindow::SRect, strflag2axes_properties(), sciSubWindow::theta, sciSubWindow::theta_kp, TRUE, update_specification_bounds(), UpdateSubwinScale(), and sciSubWindow::with_leg.

Referenced by Objplot2d().

{
  int closeflag = 0;
  int jj = 0;
  sciPointObj **pptabofpointobj;
  sciPointObj  *psubwin;
  long hdl;
  long *hdltab;
  int cmpt=0/*,i*/;
  BOOL containsSurface = FALSE ;
  int with_leg;
  double drect[6];
  char dataflag/*,frameflag*/;
  sciSubWindow * ppsubwin = NULL;
  BOOL bounds_changed = FALSE;
  BOOL axes_properties_changed = FALSE;

  psubwin = sciGetCurrentSubWin();
  ppsubwin = pSUBWIN_FEATURE(psubwin);

  /* check if the auto_clear property is on and then erase everything */
  checkRedrawing() ;
  
  if ( sciGetSurface(psubwin) !=  NULL ) 
  {
    containsSurface = TRUE ;
  }
  
  if (sciGetSurface(psubwin) == (sciPointObj *) NULL) /* F.Leray 18.05.04 */
    {
      ppsubwin->is3d = FALSE;
      ppsubwin->project[2]= 0;
    }
  else
    {
      ppsubwin->theta_kp=ppsubwin->theta;
      ppsubwin->alpha_kp=ppsubwin->alpha;
    }
  
  ppsubwin->alpha  = 0.0;
  ppsubwin->theta  = 270.0;
  
  if (sciGetSurface(psubwin) != (sciPointObj *) NULL){
    if(sciGetCurrentScilabXgc () != (struct BCG *) NULL)
      UpdateSubwinScale(psubwin);
    ppsubwin->is3d = FALSE;
  }


  /* Force psubwin->axes.aaint to those given by argument aaint*/
  /* F.Leray 07.10.04 REMOVE AAINT*/
 /*  for (i=0;i<4;i++) pSUBWIN_FEATURE(psubwin)->axes.aaint[i] = aaint[i];  */

  

  /* Force psubwin->logflags to those given by argument*/
  if ( ppsubwin->FirstPlot )
  {
    ppsubwin->logflags[0]=logflags[1];
    ppsubwin->logflags[1]=logflags[2];
  }

  /* Force "cligrf" clipping */
  sciSetIsClipping (psubwin,0); 

  /* Force  axes_visible property */
  /*pSUBWIN_FEATURE (psubwin)->isaxes  = TRUE;*/ /* WHY ??? */

  if (sciGetGraphicMode (psubwin)->autoscaling) {
    /* compute and merge new specified bounds with psubwin->Srect */
    switch (strflag[1])  {
    case '0': 
      /* do not change psubwin->Srect */
      break;
    case '1' : case '3' : case '5' : case '7':
      /* Force psubwin->Srect=brect */
      re_index_brect(brect, drect);
      break;
    case '2' : case '4' : case '6' : case '8': case '9':
      /* Force psubwin->Srect to the x and y bounds */
      if ( (int)strlen(logflags) < 1) dataflag='g' ; else dataflag=logflags[0];
      compute_data_bounds2(0,dataflag,ppsubwin->logflags,x,y,*n1,*n2,drect);
      break;
    }
    if (!pSUBWIN_FEATURE(psubwin)->FirstPlot && 
        (strflag[1] == '5' || strflag[1] == '7' || strflag[1] == '8' || strflag[1] == '9')) { /* merge psubwin->Srect and drect */
      
      drect[0] = Min(ppsubwin->SRect[0],drect[0]); /*xmin*/
      drect[2] = Min(ppsubwin->SRect[2],drect[2]); /*ymin*/
      drect[1] = Max(ppsubwin->SRect[1],drect[1]); /*xmax*/
      drect[3] = Max(ppsubwin->SRect[3],drect[3]); /*ymax*/
      
    }
    if (strflag[1] != '0')
      bounds_changed = update_specification_bounds(psubwin, drect,2);
  } 
  
  if(ppsubwin->FirstPlot == TRUE) bounds_changed = TRUE;
  
  axes_properties_changed = strflag2axes_properties(psubwin, strflag);
     
  ppsubwin->FirstPlot = FALSE; /* just after strflag2axes_properties */
   
  with_leg= (strflag[0] == '1');
  ppsubwin->with_leg = with_leg;

  /* F.Leray 07.10.04 : trigger algo to init. manual graduation u_xgrads and 
     u_ygrads if nax (in matdes.c which is == aaint HERE) was specified */
  
  ppsubwin->flagNax = flagNax; /* store new value for flagNax */

  if(ppsubwin->flagNax == TRUE){
    if(ppsubwin->logflags[0] == 'n' && ppsubwin->logflags[1] == 'n')
      {
        ppsubwin->axes.auto_ticks[0] = FALSE; /* x and y graduations are imposed by Nax */
        ppsubwin->axes.auto_ticks[1] = FALSE;
        
        CreatePrettyGradsFromNax(psubwin,aaint);
      }
    else{
      sciprint("Warning : Nax does not work with logarithmic scaling\n");}
  }
  
  if(bounds_changed == TRUE || axes_properties_changed == TRUE)
    sciDrawObj(sciGetCurrentFigure());
  /*     EraseAndOrRedraw(psubwin); /\* inhibit EraseAndOrRedraw for now F.Leray 20.12.04 *\/ */
  
  /*---- Drawing the curves and the legends ----*/
  if ( *n1 != 0 ) {
    frame_clip_on ();
    if ((hdltab = MALLOC ((*n1+2) * sizeof (long))) == NULL) {
      sciprint ("Running out of memory for plot2d\n");
      return 0;   
    }
    if (with_leg) {
      /* pptabofpointobj allocated for legends */
      if ((pptabofpointobj = MALLOC((*n1)*sizeof(sciPointObj*))) == NULL) {
        sciprint ("Running out of memory for plot2d\n");
        FREE(hdltab);
        return 0;
      }
    }
    for (jj = 0;jj < *n1; jj++) {/*A.Djalel 3D axes*/
      sciPointObj * pobj = NULL;
      if (style[jj] > 0) { 
        sciSetCurrentObj (ConstructPolyline
                          (sciGetCurrentSubWin(),&(x[jj*(*n2)]),
                           &(y[jj*(*n2)]),PD0,closeflag,*n2,*n1,ptype,
                           &style[jj],NULL,NULL,NULL,NULL,TRUE,FALSE,FALSE,FALSE));
      }
      else {
        int minusstyle = -style[jj];
        sciSetCurrentObj (ConstructPolyline
                          (sciGetCurrentSubWin(),&(x[jj*(*n2)]),
                           &(y[jj*(*n2)]),PD0,closeflag,*n2,*n1,ptype,
                           NULL,NULL,&minusstyle,NULL,NULL,FALSE,FALSE,TRUE,FALSE));
      }
      pobj = sciGetCurrentObj();
      
      if (with_leg) pptabofpointobj[jj] = pobj;
      sciDrawObjIfRequired(pobj);
      
      hdl=sciGetHandle(pobj);
      hdltab[cmpt]=hdl;
      cmpt++;
    }
    
    DrawAxesIfRequired(sciGetCurrentObj ()); /* force axes redrawing once is sufficient (F.Leray 10.01.05) */
    
    frame_clip_off ();

    /*---- Drawing the Legends ----*/
    if (with_leg) {
      sciSetCurrentObj (ConstructLegend
                        (sciGetCurrentSubWin(),
                         legend, strlen(legend), *n1, style, pptabofpointobj)); 
      hdl=sciGetHandle(sciGetCurrentObj ());   
      hdltab[cmpt]=hdl;
      cmpt++;
      FREE(pptabofpointobj);
    }

    /*---- construct Compound ----*/
    sciSetCurrentObj(ConstructCompound (hdltab, cmpt)); 
    FREE(hdltab);
    if ( containsSurface ) {
      Merge3d(psubwin);
      sciDrawObj(sciGetCurrentFigure ());
    }
    return(0);
  }
  

  /*  sciDrawObj(sciGetCurrentFigure ());*/
  return(0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

int re_index_brect	(	double *	brect,
		double *	drect
	)

BOOL strflag2axes_properties	(	sciPointObj *	psubwin,
		char *	strflag
	)

BOOL update_specification_bounds	(	sciPointObj *	psubwin,
		rect	,
		int	flag
	)

Definition at line 411 of file Plo2dn.c.

References FALSE, pSUBWIN_FEATURE, sci_update_frame_bounds_2d(), sci_update_frame_bounds_3d(), and sciSubWindow::SRect.

{
  sciSubWindow * ppsubwin = pSUBWIN_FEATURE (psubwin);
  BOOL haschanged = FALSE;
  
  ppsubwin->SRect[0] = rect[0];
  ppsubwin->SRect[1] = rect[1];
  ppsubwin->SRect[2] = rect[2];
  ppsubwin->SRect[3] = rect[3];
  if (flag==3) {
    ppsubwin->SRect[4] = rect[4];
    ppsubwin->SRect[5] = rect[5];
  }
  
  if(flag != 3)
    haschanged = sci_update_frame_bounds_2d(psubwin);
  else
    haschanged = sci_update_frame_bounds_3d(psubwin);
  
  
  return haschanged;
}

Here is the call graph for this function:

BOOL update_specification_bounds	(	sciPointObj *	psubwin,
		double *	rect,
		int	flag
	)

---