Actual source code: lg.c

      2: #include <../src/sys/classes/draw/utils/lgimpl.h>

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

      8:    Logically Collective on PetscDrawLG

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

     15:    Level: intermediate

     17:    Note: You must call PetscDrawLGDraw() to display any added points
     18:          Call PetscDrawLGReset() to remove all points

     20: .seealso: PetscDrawLGCreate(), PetscDrawLGAddPoints(), PetscDrawLGAddPoint(), PetscDrawLGReset(), PetscDrawLGDraw()
     21: @*/
     22: PetscErrorCode  PetscDrawLGAddCommonPoint(PetscDrawLG lg,const PetscReal x,const PetscReal *y)
     23: {
     25:   PetscInt       i;

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

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

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

     58:    Logically Collective on PetscDrawLG

     60:    Input Parameters:
     61: +  lg - the LineGraph data structure
     62: -  x, y - the points to two arrays containing the new x and y
     63:           point for each curve.

     65:    Note: You must call PetscDrawLGDraw() to display any added points
     66:          Call PetscDrawLGReset() to remove all points

     68:    Level: intermediate

     70: .seealso: PetscDrawLGCreate(), PetscDrawLGAddPoints(), PetscDrawLGAddCommonPoint(), PetscDrawLGReset(), PetscDrawLGDraw()
     71: @*/
     72: PetscErrorCode  PetscDrawLGAddPoint(PetscDrawLG lg,const PetscReal *x,const PetscReal *y)
     73: {
     75:   PetscInt       i;
     76:   PetscReal      xx;

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

    103:     lg->x[lg->loc]   = xx;
    104:     lg->y[lg->loc++] = y[i];
    105:   }
    106:   lg->nopts++;
    107:   return(0);
    108: }

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

    114:    Logically Collective on PetscDrawLG

    116:    Input Parameters:
    117: +  lg - the LineGraph data structure
    118: .  xx,yy - points to two arrays of pointers that point to arrays
    119:            containing the new x and y points for each curve.
    120: -  n - number of points being added

    122:    Level: intermediate

    124:    Note: You must call PetscDrawLGDraw() to display any added points
    125:          Call PetscDrawLGReset() to remove all points

    127: .seealso: PetscDrawLGCreate(), PetscDrawLGAddPoint(), PetscDrawLGAddCommonPoint(), PetscDrawLGReset(), PetscDrawLGDraw()
    128: @*/
    129: PetscErrorCode  PetscDrawLGAddPoints(PetscDrawLG lg,PetscInt n,PetscReal **xx,PetscReal **yy)
    130: {
    132:   PetscInt       i,j,k;
    133:   PetscReal      *x,*y;

    138:   if (lg->loc+n*lg->dim >= lg->len) { /* allocate more space */
    139:     PetscReal *tmpx,*tmpy;
    140:     PetscInt  chunk = CHUNCKSIZE;

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

    161:       lg->x[k] = x[i];
    162:       lg->y[k] = y[i];
    163:       k       += lg->dim;
    164:     }
    165:   }
    166:   lg->loc   += n*lg->dim;
    167:   lg->nopts += n;
    168:   return(0);
    169: }