b149507f0bdd54376987b7b3beeb975abaa61d25
[scilab.git] / scilab / modules / scicos_blocks / src / c / cscope.c
1 /*
2  *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3  *  Copyright (C) 2011-2012 - Scilab Enterprises - Clement DAVID
4  *
5  * Copyright (C) 2012 - 2016 - Scilab Enterprises
6  *
7  * This file is hereby licensed under the terms of the GNU GPL v2.0,
8  * pursuant to article 5.3.4 of the CeCILL v.2.1.
9  * This file was originally licensed under the terms of the CeCILL v2.1,
10  * and continues to be available under such terms.
11  * For more information, see the COPYING file which you should have received
12  * along with this program.
13  *
14  */
15
16 #include <string.h>
17
18 #include "dynlib_scicos_blocks.h"
19 #include "scoUtils.h"
20
21 #include "sci_malloc.h"
22 #include "elementary_functions.h"
23
24 #include "getGraphicObjectProperty.h"
25 #include "setGraphicObjectProperty.h"
26 #include "graphicObjectProperties.h"
27 #include "createGraphicObject.h"
28 #include "deleteGraphicObject.h"
29
30 #include "CurrentFigure.h"
31 #include "CurrentObject.h"
32
33 #include "scicos_block4.h"
34 #include "scicos.h"
35
36 #include "localization.h"
37 #include "os_string.h"
38
39 #include "FigureList.h"
40 #include "BuildObjects.h"
41 #include "AxesModel.h"
42
43 #define HISTORY_POINTS_THRESHOLD 4096
44
45 /*****************************************************************************
46  * Internal container structure
47  ****************************************************************************/
48
49 /**
50  * Container structure
51  */
52 typedef struct
53 {
54     struct
55     {
56         int numberOfPoints;
57         double ***bufferCoordinates;
58         int maxNumberOfPoints;
59         double ***historyCoordinates;
60     } internal;
61
62     struct
63     {
64         int periodCounter;
65
66         BOOL disableBufferUpdate;
67         int historyUpdateCounter;
68
69         int cachedFigureUID;
70         int cachedAxeUID;
71         int* cachedBufferPolylinesUIDs;
72         int* cachedHistoryPolylinesUIDs;
73     } scope;
74 } sco_data;
75
76 /**
77  * Get (and allocate on demand) the internal data used on this scope
78  * \param block the block
79  * \return the scope data
80  */
81 static sco_data *getScoData(scicos_block * block);
82
83 /**
84  * Release any internal data
85  *
86  * \param block the block
87  */
88 static void freeScoData(scicos_block * block);
89
90 /**
91  * Realloc the history buffer data
92  *
93  * \param block the block
94  * \param numberOfPoints realloc to store this number of points
95  */
96 static sco_data *reallocHistoryBuffer(scicos_block * block, int numberOfPoints);
97
98 /**
99  * Set values into the coordinates buffer.
100  *
101  * \param block the block
102  * \param coordinates the buffer
103  * \param numberOfPoints the number of points to set (actual)
104  * \param bufferPoints the buffer size (max)
105  * \param t the time to set
106  * \param value the value to set
107  */
108 static void setBuffersCoordinates(scicos_block * block, double* coordinates, const int numberOfPoints,
109                                   const int bufferPoints, const double t, const double value);
110
111 /**
112  * Append the data to the current data
113  *
114  * \param block the block
115  * \param input the input (0-indexed)
116  * \param t the current time
117  * \param data the data to append
118  */
119 static void appendData(scicos_block * block, int input, double t, double *data);
120
121 /**
122  * Push the block data to the polyline
123  *
124  * \param block the block
125  * \param input the selected input
126  * \param row the selected row
127  * \param iPolylineUID the polyline uid
128  *
129  */
130 static BOOL pushData(scicos_block * block, int input, int row);
131
132 /*****************************************************************************
133  * Graphics utils
134  ****************************************************************************/
135
136 /**
137  * Get (and allocate on demand) the figure associated with the block
138  * \param block the block
139  * \return a valid figure UID or NULL on error
140  */
141 static int getFigure(scicos_block * block);
142
143 /**
144  * Get (and allocate on demand) the axe associated with the input
145  *
146  * \param iFigureUID the parent figure UID
147  * \param block the block
148  * \param input the current input index (0-indexed)
149  * \return a valid axe UID or NULL on error
150  */
151 static int getAxe(int iFigureUID, scicos_block * block, int input);
152
153 /**
154  * Get (and allocate on demand) the polyline associated with the row
155  *
156  * \param iAxeUID the parent axe UID
157  * \param block the block
158  * \param row the current row index (0-indexed)
159  * \param history get the history polyline
160  * \return a valid polyline UID or NULL on error
161  */
162 static int getPolyline(int iAxeUID, scicos_block * block, int row, BOOL history);
163
164 /**
165  * Delete all the buffer polylines.
166  *
167  * \param block the block
168  */
169 static void deleteBufferPolylines(scicos_block * block);
170
171 /**
172  * Set the polylines history size and push the history buffer
173  *
174  * \param block the block
175  * \param input the input port index
176  * \param maxNumberOfPoints the size of the buffer
177  */
178 static BOOL pushHistory(scicos_block * block, int input, int maxNumberOfPoints);
179
180 /**
181  * Set the polylines bounds
182  *
183  * \param block the current block
184  * \param iAxeUID the axe uid
185  * \param periodCounter number of past periods since startup
186  */
187 static BOOL setPolylinesBounds(scicos_block * block, int iAxeUID, int periodCounter);
188
189 /*****************************************************************************
190  * Simulation function
191  ****************************************************************************/
192
193 /** \fn void cscope(scicos_block * block,int flag)
194     \brief the computational function
195     \param block A pointer to a scicos_block
196     \param flag An int which indicates the state of the block (init, update, ending)
197 */
198 SCICOS_BLOCKS_IMPEXP void cscope(scicos_block * block, scicos_flag flag)
199 {
200     int iFigureUID;
201
202     double t;
203     double *u;
204     sco_data *sco;
205
206     int i;
207     BOOL result;
208
209     switch (flag)
210     {
211
212         case Initialization:
213             sco = getScoData(block);
214             if (sco == NULL)
215             {
216                 set_block_error(-5);
217                 break;
218             }
219             iFigureUID = getFigure(block);
220             if (iFigureUID == 0)
221             {
222                 // allocation error
223                 set_block_error(-5);
224                 break;
225             }
226             break;
227
228         case StateUpdate:
229             iFigureUID = getFigure(block);
230             if (iFigureUID == 0)
231             {
232                 // allocation error
233                 set_block_error(-5);
234                 break;
235             }
236
237             t = get_scicos_time();
238             u = GetRealInPortPtrs(block, 1);
239
240             appendData(block, 0, t, u);
241
242             for (i = 0; i < block->insz[0]; i++)
243             {
244                 result = pushData(block, 0, i);
245                 if (result == FALSE)
246                 {
247                     Coserror("%s: unable to push some data.", "cscope");
248                     break;
249                 }
250             }
251             break;
252
253         case Ending:
254             sco = getScoData(block);
255             sco = reallocHistoryBuffer(block, sco->internal.maxNumberOfPoints + sco->internal.numberOfPoints);
256             sco->scope.disableBufferUpdate = FALSE;
257             sco->scope.historyUpdateCounter = 0;
258 #ifdef WITH_GUI
259             pushHistory(block, 0, sco->internal.maxNumberOfPoints);
260             deleteBufferPolylines(block);
261 #endif
262             freeScoData(block);
263             break;
264
265         default:
266             break;
267     }
268 }
269
270 /*-------------------------------------------------------------------------*/
271
272 /*****************************************************************************
273  *
274  * Container management
275  *
276  ****************************************************************************/
277
278 static sco_data *getScoData(scicos_block * block)
279 {
280     sco_data *sco = (sco_data *) * (block->work);
281     int i, j;
282
283     if (sco == NULL)
284     {
285         /*
286          * Data allocation
287          */
288
289         sco = (sco_data *) MALLOC(sizeof(sco_data));
290         if (sco == NULL)
291         {
292             goto error_handler_sco;
293         }
294
295         // 0 points out of a block->ipar[2] points buffer
296         sco->internal.numberOfPoints = 0;
297
298         sco->internal.bufferCoordinates = (double ***)CALLOC(block->nin, sizeof(double **));
299         if (sco->internal.bufferCoordinates == NULL)
300         {
301             goto error_handler_bufferCoordinates;
302         }
303
304         for (i = 0; i < block->nin; i++)
305         {
306             sco->internal.bufferCoordinates[i] = (double **)CALLOC(block->insz[i], sizeof(double *));
307             if (sco->internal.bufferCoordinates[i] == NULL)
308             {
309                 goto error_handler_bufferCoordinates_i;
310             }
311         }
312         for (i = 0; i < block->nin; i++)
313         {
314             for (j = 0; j < block->insz[i]; j++)
315             {
316                 sco->internal.bufferCoordinates[i][j] = (double *)CALLOC(3 * block->ipar[2], sizeof(double));
317
318                 if (sco->internal.bufferCoordinates[i][j] == NULL)
319                 {
320                     goto error_handler_bufferCoordinates_ij;
321                 }
322             }
323         }
324
325         // 0 points out of a 0 points history
326         sco->internal.maxNumberOfPoints = 0;
327
328         sco->internal.historyCoordinates = (double ***)CALLOC(block->nin, sizeof(double **));
329         if (sco->internal.historyCoordinates == NULL)
330         {
331             goto error_handler_historyCoordinates;
332         }
333
334         for (i = 0; i < block->nin; i++)
335         {
336             sco->internal.historyCoordinates[i] = (double **)CALLOC(block->insz[i], sizeof(double *));
337             if (sco->internal.historyCoordinates[i] == NULL)
338             {
339                 goto error_handler_historyCoordinates_i;
340             }
341         }
342
343         sco->scope.periodCounter = 0;
344
345         // flag to avoid pushing the buffer each time
346         sco->scope.disableBufferUpdate = FALSE;
347         // counter use to delay the history push
348         sco->scope.historyUpdateCounter = 0;
349
350         sco->scope.cachedFigureUID = 0;
351         sco->scope.cachedAxeUID = 0;
352         sco->scope.cachedBufferPolylinesUIDs = (int*)CALLOC(block->insz[0], sizeof(int));
353         sco->scope.cachedHistoryPolylinesUIDs = (int*)CALLOC(block->insz[0], sizeof(int));
354
355         *(block->work) = sco;
356     }
357
358     return sco;
359
360     /*
361      * Error management (out of normal flow)
362      */
363
364 error_handler_historyCoordinates_i:
365     for (j = 0; j < i; j++)
366     {
367         FREE(sco->internal.historyCoordinates[j]);
368     }
369     FREE(sco->internal.historyCoordinates);
370 error_handler_historyCoordinates:
371 error_handler_bufferCoordinates_ij:
372     for (i = 0; i < block->nin - 1; i++)
373     {
374         for (j = 0; j < block->insz[i] - 1; j++)
375         {
376             double* ptr = sco->internal.bufferCoordinates[i][j];
377             if (ptr != NULL)
378             {
379                 FREE(ptr);
380             }
381         }
382     }
383     i = block->nin - 1;
384 error_handler_bufferCoordinates_i:
385     for (j = 0; j < i; j++)
386     {
387         FREE(sco->internal.bufferCoordinates[j]);
388     }
389     FREE(sco->internal.bufferCoordinates);
390 error_handler_bufferCoordinates:
391     FREE(sco);
392 error_handler_sco:
393     // allocation error
394     set_block_error(-5);
395     return NULL;
396 }
397
398 static void freeScoData(scicos_block * block)
399 {
400     sco_data *sco = (sco_data *) * (block->work);
401     int i, j;
402
403     if (sco != NULL)
404     {
405         for (i = 0; i < block->nin; i++)
406         {
407             for (j = 0; j < block->insz[i]; j++)
408             {
409                 if (sco->internal.historyCoordinates[i][j] != NULL)
410                 {
411                     FREE(sco->internal.historyCoordinates[i][j]);
412                 }
413                 FREE(sco->internal.bufferCoordinates[i][j]);
414             }
415             FREE(sco->internal.historyCoordinates[i]);
416             FREE(sco->internal.bufferCoordinates[i]);
417         }
418         FREE(sco->internal.historyCoordinates);
419         FREE(sco->internal.bufferCoordinates);
420
421         FREE(sco->scope.cachedHistoryPolylinesUIDs);
422         FREE(sco->scope.cachedBufferPolylinesUIDs);
423
424         FREE(sco);
425         *(block->work) = NULL;
426     }
427 }
428
429 static sco_data *reallocHistoryBuffer(scicos_block * block, int numberOfPoints)
430 {
431     sco_data *sco = (sco_data *) * (block->work);
432     int i;
433
434     double *ptr;
435     int allocatedNumberOfPoints;
436
437     int previousNumberOfPoints = sco->internal.maxNumberOfPoints;
438     int numberOfCopiedPoints = numberOfPoints - sco->internal.maxNumberOfPoints;
439
440     double *buffer;
441     int bufferNumberOfPoints = block->ipar[2];
442     int bufferNewPointInc;
443
444     if (previousNumberOfPoints == 0)
445     {
446         allocatedNumberOfPoints = numberOfPoints;
447         bufferNewPointInc = 0;
448     }
449     else
450     {
451         allocatedNumberOfPoints = numberOfPoints - 1;
452         bufferNewPointInc = 1;
453     }
454
455     if (sco->scope.historyUpdateCounter <= 0)
456     {
457         if (numberOfPoints > HISTORY_POINTS_THRESHOLD)
458         {
459             sco->scope.disableBufferUpdate = TRUE;
460             sco->scope.historyUpdateCounter = numberOfPoints / HISTORY_POINTS_THRESHOLD;
461         }
462         else
463         {
464             sco->scope.disableBufferUpdate = FALSE;
465             sco->scope.historyUpdateCounter = 0;
466         }
467     }
468
469     for (i = 0; i < block->insz[0]; i++)
470     {
471         ptr = (double *)MALLOC(3 * allocatedNumberOfPoints * sizeof(double));
472         if (ptr == NULL)
473         {
474             goto error_handler;
475         }
476
477         // memcpy existing X-axis values from the history
478         memcpy(ptr, sco->internal.historyCoordinates[0][i], previousNumberOfPoints * sizeof(double));
479         // memcpy existing Y-axis values from the history
480         memcpy(ptr + allocatedNumberOfPoints, sco->internal.historyCoordinates[0][i] + previousNumberOfPoints, previousNumberOfPoints * sizeof(double));
481         // clear the last points, the Z-axis values
482         memset(ptr + 2 * allocatedNumberOfPoints, 0, allocatedNumberOfPoints * sizeof(double));
483
484         // then set the last points to the last values for X-axis and Y-axis values from the buffer points
485         buffer = sco->internal.bufferCoordinates[0][i];
486         memcpy(ptr + previousNumberOfPoints, buffer + bufferNewPointInc, (numberOfCopiedPoints - bufferNewPointInc) * sizeof(double));
487         memcpy(ptr + allocatedNumberOfPoints + previousNumberOfPoints, buffer + bufferNumberOfPoints + bufferNewPointInc, (numberOfCopiedPoints - bufferNewPointInc) * sizeof(double));
488
489         FREE(sco->internal.historyCoordinates[0][i]);
490         sco->internal.historyCoordinates[0][i] = ptr;
491     }
492
493     sco->internal.maxNumberOfPoints = allocatedNumberOfPoints;
494     return sco;
495
496 error_handler:
497     freeScoData(block);
498     // allocation error
499     set_block_error(-5);
500     return NULL;
501 }
502
503 static void setBuffersCoordinates(scicos_block* block, double* coordinates, const int numberOfPoints,
504                                   const int bufferPoints, const double t, const double value)
505 {
506     int setLen;
507     sco_data *sco = (sco_data *) * (block->work);
508
509     if (sco->scope.disableBufferUpdate == TRUE)
510     {
511         coordinates[numberOfPoints] = t;
512         coordinates[bufferPoints + numberOfPoints] = value;
513         return;
514     }
515
516     // X-axis values first
517     for (setLen = numberOfPoints; setLen < bufferPoints; setLen++)
518     {
519         coordinates[setLen] = t;
520     }
521     // then Y-axis values
522     for (setLen = numberOfPoints; setLen < bufferPoints; setLen++)
523     {
524         coordinates[bufferPoints + setLen] = value;
525     }
526     // then Z-axis values (always clear'ed)
527 }
528
529 static void appendData(scicos_block * block, int input, double t, double *data)
530 {
531     int i;
532
533     sco_data *sco = (sco_data *) * (block->work);
534
535     /*
536      * Handle the case where the t is greater than the data_bounds
537      */
538     if (t > ((sco->scope.periodCounter + 1) * block->rpar[3]))
539     {
540         sco->scope.periodCounter++;
541
542         // set the buffer coordinates to the last point
543         for (i = 0; i < block->insz[input]; i++)
544         {
545             sco->internal.bufferCoordinates[input][i][0] = sco->internal.bufferCoordinates[input][i][sco->internal.numberOfPoints - 1];
546             sco->internal.bufferCoordinates[input][i][block->ipar[2]] = sco->internal.bufferCoordinates[input][i][block->ipar[2] + sco->internal.numberOfPoints - 1];
547         }
548         sco->internal.numberOfPoints = 1;
549
550         // clear the history coordinates
551         sco->internal.maxNumberOfPoints = 0;
552         for (i = 0; i < block->insz[input]; i++)
553         {
554             if (sco->internal.historyCoordinates[input][i] != NULL)
555             {
556                 FREE(sco->internal.historyCoordinates[input][i]);
557                 sco->internal.historyCoordinates[input][i] = NULL;
558             }
559         }
560
561         // configure scope setting
562         if (setPolylinesBounds(block, getAxe(getFigure(block), block, input), sco->scope.periodCounter) == FALSE)
563         {
564             set_block_error(-5);
565             freeScoData(block);
566             sco = NULL;
567         }
568     }
569
570     /*
571      * Handle the case where the scope has more points than maxNumberOfPoints
572      */
573     if (sco != NULL && sco->internal.numberOfPoints >= block->ipar[2])
574     {
575         int maxNumberOfPoints = sco->internal.maxNumberOfPoints;
576
577         // on a full scope, re-alloc history coordinates
578         maxNumberOfPoints = maxNumberOfPoints + block->ipar[2];
579         sco = reallocHistoryBuffer(block, maxNumberOfPoints);
580
581         // set the buffer coordinates to the last point
582         for (i = 0; i < block->insz[input]; i++)
583         {
584             sco->internal.bufferCoordinates[input][i][0] = sco->internal.bufferCoordinates[input][i][block->ipar[2] - 1];
585             sco->internal.bufferCoordinates[input][i][block->ipar[2]] = sco->internal.bufferCoordinates[input][i][2 * block->ipar[2] - 1];
586         }
587         sco->internal.numberOfPoints = 1;
588
589         // reconfigure related graphic objects
590         if (pushHistory(block, input, sco->internal.maxNumberOfPoints) == FALSE)
591         {
592             set_block_error(-5);
593             freeScoData(block);
594             sco = NULL;
595         }
596     }
597
598     /*
599      * Update data
600      */
601     if (sco != NULL)
602     {
603         for (i = 0; i < block->insz[input]; i++)
604         {
605             const double value = data[i];
606             setBuffersCoordinates(block, sco->internal.bufferCoordinates[input][i], sco->internal.numberOfPoints, block->ipar[2], t, value);
607         }
608
609         sco->internal.numberOfPoints++;
610     }
611 }
612
613 static BOOL pushData(scicos_block * block, int input, int row)
614 {
615     int iFigureUID;
616     int iAxeUID;
617     int iPolylineUID;
618
619     double *data;
620     sco_data *sco;
621
622     iFigureUID = getFigure(block);
623     iAxeUID = getAxe(iFigureUID, block, input);
624     iPolylineUID = getPolyline(iAxeUID, block, row, FALSE);
625
626     sco = getScoData(block);
627     if (sco == NULL)
628     {
629         return FALSE;
630     }
631
632     // do not push any data if disabled
633     if (sco->scope.disableBufferUpdate == TRUE)
634     {
635         return TRUE;
636     }
637
638     // select the right input and row
639     data = sco->internal.bufferCoordinates[input][row];
640
641     return setGraphicObjectProperty(iPolylineUID, __GO_DATA_MODEL_COORDINATES__, data, jni_double_vector, block->ipar[2]);
642 }
643
644 /*****************************************************************************
645  *
646  * Graphic utils
647  *
648  ****************************************************************************/
649
650 /**
651  * Set properties on the figure.
652  *
653  * \param iFigureUID the figure uid
654  * \param block the current block
655  */
656 static void setFigureSettings(int iFigureUID, scicos_block * block)
657 {
658     char *label = NULL;
659
660     int nipar = GetNipar(block);
661     int *ipar = GetIparPtrs(block);
662
663     int win_pos[2];
664     int win_dim[2];
665
666     win_pos[0] = ipar[(nipar - 1) - 3];
667     win_pos[1] = ipar[(nipar - 1) - 2];
668     win_dim[0] = ipar[(nipar - 1) - 1];
669     win_dim[1] = ipar[nipar - 1];
670
671     if (win_pos[0] > 0 && win_pos[1] > 0)
672     {
673         setGraphicObjectProperty(iFigureUID, __GO_POSITION__, &win_pos, jni_int_vector, 2);
674     }
675
676     if (win_dim[0] > 0 && win_dim[1] > 0)
677     {
678         setGraphicObjectProperty(iFigureUID, __GO_SIZE__, &win_dim, jni_int_vector, 2);
679     }
680
681     label = GetLabelPtrs(block);
682     if (label != NULL)
683     {
684         if (strlen(label) > 0)
685         {
686             setGraphicObjectProperty(iFigureUID, __GO_NAME__, label, jni_string, 1);
687         }
688     }
689 };
690
691 /*****************************************************************************
692  *
693  * Graphic
694  *
695  ****************************************************************************/
696
697 static int getFigure(scicos_block * block)
698 {
699 #ifdef WITH_GUI
700     signed int figNum;
701     int iFigureUID = 0;
702     int iAxe = 0;
703     int i__1 = 1;
704     sco_data *sco = (sco_data *) * (block->work);
705
706     int i;
707
708     // assert the sco is not NULL
709     if (sco == NULL)
710     {
711         return 0;
712     }
713
714     // fast path for an existing object
715     if (sco->scope.cachedFigureUID != 0)
716     {
717         return sco->scope.cachedFigureUID;
718     }
719
720     figNum = block->ipar[0];
721
722     // with a negative id, use the block number indexed from a constant.
723     if (figNum < 0)
724     {
725         figNum = 20000 + get_block_number();
726     }
727
728     iFigureUID = getFigureFromIndex(figNum);
729     // create on demand
730     if (iFigureUID == 0)
731     {
732         iFigureUID = createNewFigureWithAxes();
733         setGraphicObjectProperty(iFigureUID, __GO_ID__, &figNum, jni_int, 1);
734
735         // the stored uid is a reference to the figure map, not to the current figure
736         iFigureUID = getFigureFromIndex(figNum);
737         sco->scope.cachedFigureUID = iFigureUID;
738
739         // set configured parameters
740         setFigureSettings(iFigureUID, block);
741
742         // allocate the axes through the getter
743         for (i = 0; i < GetNin(block); i++)
744         {
745             iAxe = getAxe(iFigureUID, block, i);
746
747             /*
748              * Setup according to block settings
749              */
750             setLabel(iAxe, __GO_X_AXIS_LABEL__, "t");
751             setLabel(iAxe, __GO_Y_AXIS_LABEL__, "y");
752
753             setGraphicObjectProperty(iAxe, __GO_X_AXIS_VISIBLE__, &i__1, jni_bool, 1);
754             setGraphicObjectProperty(iAxe, __GO_Y_AXIS_VISIBLE__, &i__1, jni_bool, 1);
755
756             setPolylinesBounds(block, iAxe, 0);
757         }
758     }
759     else
760     {
761         // set configured parameters
762         setFigureSettings(iFigureUID, block);
763     }
764
765     if (sco->scope.cachedFigureUID == 0)
766     {
767         sco->scope.cachedFigureUID = iFigureUID;
768     }
769     return iFigureUID;
770 #else
771     Coserror("%s: Scilab is compiled without GUI, can not use Scope.", "cscope");
772     return 0;
773 #endif
774 }
775
776 static int getAxe(int iFigureUID, scicos_block * block, int input)
777 {
778     int iAxe;
779     int i;
780     sco_data *sco = (sco_data *) * (block->work);
781
782     // assert the sco is not NULL
783     if (sco == NULL)
784     {
785         return 0;
786     }
787
788     // fast path for an existing object
789     if (sco->scope.cachedAxeUID != 0)
790     {
791         return sco->scope.cachedAxeUID;
792     }
793
794     iAxe = findChildWithKindAt(iFigureUID, __GO_AXES__, input);
795
796     /*
797      * Allocate if necessary
798      */
799     if (iAxe == 0)
800     {
801         cloneAxesModel(iFigureUID);
802         iAxe = findChildWithKindAt(iFigureUID, __GO_AXES__, input);
803     }
804
805     /*
806      * Setup on first access
807      */
808     if (iAxe != 0)
809     {
810         // allocate the polylines through the getter
811         for (i = 0; i < block->insz[input]; i++)
812         {
813             getPolyline(iAxe, block, i, TRUE);
814         }
815         for (i = 0; i < block->insz[input]; i++)
816         {
817             getPolyline(iAxe, block, i, FALSE);
818         }
819     }
820     else
821     {
822         return 0;
823     }
824
825     /*
826      * then cache with local storage
827      */
828     sco->scope.cachedAxeUID = iAxe;
829     return sco->scope.cachedAxeUID;
830 }
831
832 static int getPolyline(int iAxeUID, scicos_block * block, int row, BOOL history)
833 {
834     int iPolyline;
835     BOOL b__true = TRUE;
836
837     int color;
838
839     int* polylinesUIDs;
840     int polylineIndex;
841
842     sco_data *sco = (sco_data *) * (block->work);
843
844     // assert the sco is not NULL
845     if (sco == NULL)
846     {
847         return 0;
848     }
849
850     if (!history)
851     {
852         polylinesUIDs = sco->scope.cachedBufferPolylinesUIDs;
853         polylineIndex = block->insz[0] + row;
854     }
855     else
856     {
857         polylinesUIDs = sco->scope.cachedHistoryPolylinesUIDs;
858         polylineIndex = row;
859     }
860
861     // assert that the structure is in a good shape
862     if (polylinesUIDs == NULL)
863     {
864         return 0;
865     }
866
867     // fast path for an existing object
868     if (polylinesUIDs[row] != 0)
869     {
870         return polylinesUIDs[row];
871     }
872
873     iPolyline = findChildWithKindAt(iAxeUID, __GO_POLYLINE__, polylineIndex);
874
875     /*
876      * Allocate if necessary
877      */
878     if (iPolyline == 0)
879     {
880         iPolyline = createGraphicObject(__GO_POLYLINE__);
881
882         if (iPolyline != 0)
883         {
884             createDataObject(iPolyline, __GO_POLYLINE__);
885             setGraphicObjectRelationship(iAxeUID, iPolyline);
886         }
887         else
888         {
889             return 0;
890         }
891     }
892
893     /*
894      * Setup on first access
895      */
896
897     /*
898      * Default setup of the nGons property
899      */
900     {
901         int nGons = 1;
902         setGraphicObjectProperty(iPolyline, __GO_DATA_MODEL_NUM_GONS__, &nGons, jni_int, 1);
903     }
904
905     color = block->ipar[3 + row];
906     if (color > 0)
907     {
908         setGraphicObjectProperty(iPolyline, __GO_LINE_MODE__, &b__true, jni_bool, 1);
909         setGraphicObjectProperty(iPolyline, __GO_LINE_COLOR__, &color, jni_int, 1);
910     }
911     else
912     {
913         int iMarkSize = 4;
914         color = -color;
915         setGraphicObjectProperty(iPolyline, __GO_MARK_MODE__, &b__true, jni_bool, 1);
916         setGraphicObjectProperty(iPolyline, __GO_MARK_STYLE__, &color, jni_int, 1);
917         setGraphicObjectProperty(iPolyline, __GO_MARK_SIZE__, &iMarkSize, jni_int, 1);
918     }
919
920     {
921         int iClipState = 1; //on
922         setGraphicObjectProperty(iPolyline, __GO_CLIP_STATE__, &iClipState, jni_int, 1);
923     }
924
925     /*
926      * then cache with local storage
927      */
928     polylinesUIDs[row] = iPolyline;
929     return polylinesUIDs[row];
930 }
931
932 static void deleteBufferPolylines(scicos_block * block)
933 {
934     int i, j;
935
936     int iPolylineUID;
937
938     sco_data *sco;
939
940     sco = getScoData(block);
941     for (i = 0; i < block->nin; i++)
942     {
943         for (j = 0; j < block->insz[i]; j++)
944         {
945             iPolylineUID = sco->scope.cachedBufferPolylinesUIDs[j];
946             deleteGraphicObject(iPolylineUID);
947         }
948     }
949 }
950
951 static BOOL pushHistory(scicos_block * block, int input, int maxNumberOfPoints)
952 {
953     int i;
954
955     int iFigureUID;
956     int iAxeUID;
957     int iPolylineUID;
958
959     double *data;
960     sco_data *sco;
961
962     BOOL result = TRUE;
963
964     sco = getScoData(block);
965     iFigureUID = getFigure(block);
966     iAxeUID = getAxe(iFigureUID, block, input);
967
968     // push the data only if the counter == 0, decrement the counter if positive
969     if (sco->scope.historyUpdateCounter > 0)
970     {
971         sco->scope.historyUpdateCounter--;
972     }
973     if (sco->scope.historyUpdateCounter > 0)
974     {
975         return result;
976     }
977
978     for (i = 0; i < block->insz[input]; i++)
979     {
980         iPolylineUID = getPolyline(iAxeUID, block, i, TRUE);
981
982         data = sco->internal.historyCoordinates[input][i];
983         result = setGraphicObjectProperty(iPolylineUID, __GO_DATA_MODEL_COORDINATES__, data, jni_double_vector, maxNumberOfPoints);
984         if (result == FALSE)
985         {
986             return result;
987         }
988     }
989
990     return result;
991 }
992
993 static BOOL setPolylinesBounds(scicos_block * block, int iAxeUID, int periodCounter)
994 {
995     double dataBounds[6];
996     double period = block->rpar[3];
997
998     dataBounds[0] = periodCounter * period; // xMin
999     dataBounds[1] = (periodCounter + 1) * period;   // xMax
1000     dataBounds[2] = block->rpar[1]; // yMin
1001     dataBounds[3] = block->rpar[2]; // yMax
1002     dataBounds[4] = -1.0;       // zMin
1003     dataBounds[5] = 1.0;        // zMax
1004
1005     return setGraphicObjectProperty(iAxeUID, __GO_DATA_BOUNDS__, dataBounds, jni_double_vector, 6);
1006 }