Actual source code: lg.c

petsc-3.6.1 2015-07-22
Report Typos and Errors
  2: #include <../src/sys/classes/draw/utils/lgimpl.h>  /*I   "petscdraw.h"  I*/

  6: /*@
  7:    PetscDrawLGAddCommonPoint - Adds another point to each of the line graphs. All the points share
  8:       the same new X coordinate.  The new point must have an X coordinate larger than the old points.

 10:    Not Collective, but ignored by all processors except processor 0 in PetscDrawLG

 12:    Input Parameters:
 13: +  lg - the LineGraph data structure
 14: .   x - the common x coordiante point
 15: -   y - the new y coordinate point for each curve.

 17:    Level: intermediate

 19:    Concepts: line graph^adding points

 21: .seealso: PetscDrawLGAddPoints(), PetscDrawLGAddPoint()
 22: @*/
 23: PetscErrorCode  PetscDrawLGAddCommonPoint(PetscDrawLG lg,const PetscReal x,const PetscReal *y)
 24: {
 26:   PetscInt       i;

 29:   if (lg && ((PetscObject)lg)->classid == PETSC_DRAW_CLASSID) return(0);

 32:   if (lg->loc+lg->dim >= lg->len) { /* allocate more space */
 33:     PetscReal *tmpx,*tmpy;
 34:     PetscMalloc2(lg->len+lg->dim*CHUNCKSIZE,&tmpx,lg->len+lg->dim*CHUNCKSIZE,&tmpy);
 35:     PetscLogObjectMemory((PetscObject)lg,2*lg->dim*CHUNCKSIZE*sizeof(PetscReal));
 36:     PetscMemcpy(tmpx,lg->x,lg->len*sizeof(PetscReal));
 37:     PetscMemcpy(tmpy,lg->y,lg->len*sizeof(PetscReal));
 38:     PetscFree2(lg->x,lg->y);
 39:     lg->x    = tmpx;
 40:     lg->y    = tmpy;
 41:     lg->len += lg->dim*CHUNCKSIZE;
 42:   }
 43:   for (i=0; i<lg->dim; i++) {
 44:     if (x > lg->xmax) lg->xmax = x;
 45:     if (x < lg->xmin) lg->xmin = x;
 46:     if (y[i] > lg->ymax) lg->ymax = y[i];
 47:     if (y[i] < lg->ymin) lg->ymin = y[i];

 49:     lg->x[lg->loc]   = x;
 50:     lg->y[lg->loc++] = y[i];
 51:   }
 52:   lg->nopts++;
 53:   return(0);
 54: }

 58: /*@
 59:    PetscDrawLGAddPoint - Adds another point to each of the line graphs.
 60:    The new point must have an X coordinate larger than the old points.

 62:    Not Collective, but ignored by all processors except processor 0 in PetscDrawLG

 64:    Input Parameters:
 65: +  lg - the LineGraph data structure
 66: -  x, y - the points to two arrays containing the new x and y
 67:           point for each curve.

 69:    Level: intermediate

 71:    Concepts: line graph^adding points

 73: .seealso: PetscDrawLGAddPoints(), PetscDrawLGAddCommonPoint()
 74: @*/
 75: PetscErrorCode  PetscDrawLGAddPoint(PetscDrawLG lg,const PetscReal *x,const PetscReal *y)
 76: {
 78:   PetscInt       i;
 79:   PetscReal      xx;

 82:   if (lg && ((PetscObject)lg)->classid == PETSC_DRAW_CLASSID) return(0);

 85:   if (lg->loc+lg->dim >= lg->len) { /* allocate more space */
 86:     PetscReal *tmpx,*tmpy;
 87:     PetscMalloc2(lg->len+lg->dim*CHUNCKSIZE,&tmpx,lg->len+lg->dim*CHUNCKSIZE,&tmpy);
 88:     PetscLogObjectMemory((PetscObject)lg,2*lg->dim*CHUNCKSIZE*sizeof(PetscReal));
 89:     PetscMemcpy(tmpx,lg->x,lg->len*sizeof(PetscReal));
 90:     PetscMemcpy(tmpy,lg->y,lg->len*sizeof(PetscReal));
 91:     PetscFree2(lg->x,lg->y);
 92:     lg->x    = tmpx;
 93:     lg->y    = tmpy;
 94:     lg->len += lg->dim*CHUNCKSIZE;
 95:   }
 96:   for (i=0; i<lg->dim; i++) {
 97:     if (!x) {
 98:       xx = lg->nopts;
 99:     } else {
100:       xx = x[i];
101:     }
102:     if (xx > lg->xmax) lg->xmax = xx;
103:     if (xx < lg->xmin) lg->xmin = xx;
104:     if (y[i] > lg->ymax) lg->ymax = y[i];
105:     if (y[i] < lg->ymin) lg->ymin = y[i];

107:     lg->x[lg->loc]   = xx;
108:     lg->y[lg->loc++] = y[i];
109:   }
110:   lg->nopts++;
111:   return(0);
112: }

116: /*@C
117:    PetscDrawLGAddPoints - Adds several points to each of the line graphs.
118:    The new points must have an X coordinate larger than the old points.

120:    Not Collective, but ignored by all processors except processor 0 in PetscDrawLG

122:    Input Parameters:
123: +  lg - the LineGraph data structure
124: .  xx,yy - points to two arrays of pointers that point to arrays
125:            containing the new x and y points for each curve.
126: -  n - number of points being added

128:    Level: intermediate


131:    Concepts: line graph^adding points

133: .seealso: PetscDrawLGAddPoint()
134: @*/
135: PetscErrorCode  PetscDrawLGAddPoints(PetscDrawLG lg,PetscInt n,PetscReal **xx,PetscReal **yy)
136: {
138:   PetscInt       i,j,k;
139:   PetscReal      *x,*y;

142:   if (lg && ((PetscObject)lg)->classid == PETSC_DRAW_CLASSID) return(0);
144:   if (lg->loc+n*lg->dim >= lg->len) { /* allocate more space */
145:     PetscReal *tmpx,*tmpy;
146:     PetscInt  chunk = CHUNCKSIZE;

148:     if (n > chunk) chunk = n;
149:     PetscMalloc2(lg->len+lg->dim*chunk,&tmpx,lg->len+lg->dim*chunk,&tmpy);
150:     PetscLogObjectMemory((PetscObject)lg,2*lg->dim*chunk*sizeof(PetscReal));
151:     PetscMemcpy(tmpx,lg->x,lg->len*sizeof(PetscReal));
152:     PetscMemcpy(tmpy,lg->y,lg->len*sizeof(PetscReal));
153:     PetscFree2(lg->x,lg->y);
154:     lg->x    = tmpx;
155:     lg->y    = tmpy;
156:     lg->len += lg->dim*chunk;
157:   }
158:   for (j=0; j<lg->dim; j++) {
159:     x = xx[j]; y = yy[j];
160:     k = lg->loc + j;
161:     for (i=0; i<n; i++) {
162:       if (x[i] > lg->xmax) lg->xmax = x[i];
163:       if (x[i] < lg->xmin) lg->xmin = x[i];
164:       if (y[i] > lg->ymax) lg->ymax = y[i];
165:       if (y[i] < lg->ymin) lg->ymin = y[i];

167:       lg->x[k] = x[i];
168:       lg->y[k] = y[i];
169:       k       += lg->dim;
170:     }
171:   }
172:   lg->loc   += n*lg->dim;
173:   lg->nopts += n;
174:   return(0);
175: }

179: /*@
180:    PetscDrawLGSetLimits - Sets the axis limits for a line graph. If more
181:    points are added after this call, the limits will be adjusted to
182:    include those additional points.

184:    Not Collective, but ignored by all processors except processor 0 in PetscDrawLG

186:    Input Parameters:
187: +  xlg - the line graph context
188: -  x_min,x_max,y_min,y_max - the limits

190:    Level: intermediate

192:    Concepts: line graph^setting axis

194: @*/
195: PetscErrorCode  PetscDrawLGSetLimits(PetscDrawLG lg,PetscReal x_min,PetscReal x_max,PetscReal y_min,PetscReal y_max)
196: {
198:   if (lg && ((PetscObject)lg)->classid == PETSC_DRAW_CLASSID) return(0);
200:   (lg)->xmin = x_min;
201:   (lg)->xmax = x_max;
202:   (lg)->ymin = y_min;
203:   (lg)->ymax = y_max;
204:   return(0);
205: }