/*------------------------------------------------
* plot3d
*-----------------------------------------------*/
void Objplot3d (char * fname ,
int * isfac ,
int * izcol ,
double x[] ,
double y[] ,
double z[] ,
double * zcol ,
int * m ,
int * n ,
double * theta ,
double * alpha ,
char * legend,
int * iflag ,
double * ebox ,
int * m1 , /*Adding F.Leray 12.03.04 and 19.03.04*/
int * n1 ,
int * m2 ,
int * n2 ,
int * m3 ,
int * n3 ,
int * m3n ,
int * n3n)
/* F.Leray 25.04.05 : warning here legends means "X@Y@Z": it is labels writings!! */
/* legends has not the same meaning than inside plot2dn (there, it is really the legends of the plotted curves)*/
{
sciTypeOf3D typeof3d;
int flagcolor = 0;
int* pObj = NULL;
int i = 0;
int iSubwinUID = 0;
int firstPlot = 0;
int clipState = 0;
int iNewSurfaceUID = 0;
/* =================================================
* Force SubWindow properties according to arguments
* ================================================= */
iSubwinUID = getCurrentSubWin();
checkRedrawing();
initSubWinTo3d(iSubwinUID, legend, iflag, *alpha, *theta, ebox, x, *m1 **n1, y, *m2 **n2, z, *m3 **n3);
/* =================================================
* Analyze arguments to find entity type
* ================================================= */
if (*isfac == 1)
{
if (*izcol == 0)
{
if (strcmp(fname, "plot3d1") == 0)
{
typeof3d = SCI_FAC3D;
flagcolor = 1;
}
else
{
typeof3d = SCI_FAC3D;
flagcolor = 0;
}
}
else if (*izcol == 2)
{
typeof3d = SCI_FAC3D;
flagcolor = 3;
}
else
{
typeof3d = SCI_FAC3D;
flagcolor = 2;
}
}
else if (*isfac == 0)
{
if (strcmp(fname, "plot3d1") == 0)
{
typeof3d = SCI_PLOT3D;
flagcolor = 1;
}
else
{
typeof3d = SCI_PLOT3D;
flagcolor = 0;
}
}
else
{
typeof3d = SCI_PARAM3D1;
flagcolor = 1;
}
/* =================================================
* Construct the Entities
* ================================================= */
/*Distinction here between SCI_PARAM3D1 and others*/
if (typeof3d != SCI_PARAM3D1)
{
int dimvectx = -1;
int dimvecty = -1;
if (*isfac == 1)
{
/* x is considered as a matrix */
dimvectx = -1;
}
else if (*m1 == 1) /* x is a row vector */
{
dimvectx = *n1;
}
else if (*n1 == 1) /* x is a column vector */
{
dimvectx = *m1;
}
else /* x is a matrix */
{
dimvectx = -1;
}
if (dimvectx > 1)
{
int monotony = checkMonotony(x, dimvectx);
if (monotony == 0)
{
Scierror(999, _("%s: x vector is not monotonous.\n"), "Objplot3d");
return;
}
}
if (*isfac == 1)
{
/* y is considered as a matrix */
dimvecty = -1;
}
else if (*m2 == 1) /* y is a row vector */
{
dimvecty = *n2;
}
else if (*n2 == 1) /* y is a column vector */
{
dimvecty = *m2;
}
else /* y is a matrix */
{
dimvecty = -1;
}
if (dimvecty > 1)
{
/* test the monotony on y*/
int monotony = checkMonotony(y, dimvecty);
if (monotony == 0)
{
Scierror(999, _("%s: y vector is not monotonous.\n"), "Objplot3d");
return;
}
}
iNewSurfaceUID = ConstructSurface(iSubwinUID, typeof3d,
x, y, z, zcol, *izcol, *m, *n, iflag, ebox, flagcolor,
isfac, m1, n1, m2, n2, m3, n3, m3n, n3n);
if (iNewSurfaceUID == 0)
{
Scierror(999, _("%s: No more memory.\n"), "Objplot3d");
return;
}
setCurrentObject(iNewSurfaceUID);
/* Force clipping, 1: CLIPGRF */
clipState = 1;
setGraphicObjectProperty(iNewSurfaceUID, __GO_CLIP_STATE__, &clipState, jni_int, 1);
}
else
{
int iNewPolylineUID = 0;
int iCurrentSubwinUID = 0;
if ((pObj = (int*)MALLOC (*n * sizeof (int))) == NULL)
{
Scierror(999, "%s: No more memory.\n", fname);
return;
}
iCurrentSubwinUID = getCurrentSubWin();
for (i = 0; i < *n; ++i)
{
if ((*n > 0) && (zcol != (double *)NULL))
{
if ((int) zcol[i] > 0)
{
int intzcol = (int) zcol[i];
iNewPolylineUID = ConstructPolyline
(iCurrentSubwinUID,
&(x[*m * i]), &(y[*m * i]), &(z[*m * i]), 0, *m, 1,
&intzcol, NULL, NULL, NULL, NULL, TRUE, FALSE, FALSE, FALSE);
}
else
{
int intzcol = (int) - zcol[i];
iNewPolylineUID = ConstructPolyline
(iCurrentSubwinUID,
&(x[*m * i]), &(y[*m * i]), &(z[*m * i]), 0, *m, 1,
NULL, NULL, &intzcol, NULL, NULL, FALSE, FALSE, TRUE, FALSE);
}
}
else
{
/* default case, nothing is given */
int curcolor = 0;
int *piCurColor = &curcolor;
getGraphicObjectProperty(iCurrentSubwinUID, __GO_LINE_COLOR__, jni_int, (void**)&piCurColor);
iNewPolylineUID = ConstructPolyline(iCurrentSubwinUID,
&(x[*m * i]), &(y[*m * i]), &(z[*m * i]), 0, *m, 1,
&curcolor, NULL, NULL, NULL, NULL, TRUE, FALSE, FALSE, FALSE);
}
if (iNewPolylineUID == 0)
{
Scierror(999, _("%s: No more memory.\n"), fname);
FREE(pObj);
return;
}
setCurrentObject(iNewPolylineUID);
setGraphicObjectRelationship(iCurrentSubwinUID, iNewPolylineUID);
/* Force clipping, 1: CLIPGRF */
clipState = 1;
setGraphicObjectProperty(iNewPolylineUID, __GO_CLIP_STATE__, &clipState, jni_int, 1);
pObj[i] = iNewPolylineUID;
}
/** construct Compound and make it current object**/
if (*n > 1)
{
int o = createCompound (iCurrentSubwinUID, pObj, *n);
setCurrentObject(o);
}
FREE(pObj);
}
/* =================================================
* Redraw Figure
* ================================================= */
// subwin has been modified
firstPlot = 0;
setGraphicObjectProperty(iSubwinUID, __GO_FIRST_PLOT__, &firstPlot, jni_bool, 1);
}
/*------------------------------------------------
* plot3d
*-----------------------------------------------*/
void Objplot3d (char * fname ,
int * isfac ,
int * izcol ,
double x[] ,
double y[] ,
double z[] ,
double * zcol ,
int * m ,
int * n ,
double * theta ,
double * alpha ,
char * legend,
int * iflag ,
double * ebox ,
int * m1 , /*Adding F.Leray 12.03.04 and 19.03.04*/
int * n1 ,
int * m2 ,
int * n2 ,
int * m3 ,
int * n3 ,
int * m3n ,
int * n3n)
/* F.Leray 25.04.05 : warning here legends means "X@Y@Z": it is labels writings!! */
/* legends has not the same meaning than inside plot2dn (there, it is really the legends of the plotted curves)*/
{
sciTypeOf3D typeof3d;
int flagcolor = 0;
long *hdltab = NULL;
int i = 0;
char *psubwinUID = NULL;
char *pobjUID = NULL;
char *parentFigureUID = NULL;
double drect[6];
double rotationAngles[2];
double* dataBounds = NULL;
char * loc = NULL;
char * legx = NULL;
char * legy = NULL;
char * legz = NULL;
char* labelId = NULL;
/* char * buff = NULL; */
int flag_x = 1;
int flag_y = 1;
int dimvectx = -1;
int dimvecty = -1;
int view = 0;
int linLogFlag;
int firstPlot = 0;
int *piFirstPlot = &firstPlot;
int box = 0;
int axisVisible = 0;
int autoScale = 0;
int *piAutoScale = &autoScale;
int isoview = 0;
int clipState = 0;
char *pNewSurfaceUID = NULL;
/* Initialisation drect A.C pour debuggueur */
drect[0] = 0.0;
drect[1] = 0.0;
drect[2] = 0.0;
drect[3] = 0.0;
drect[4] = 0.0;
drect[5] = 0.0;
/* =================================================
* Force SubWindow properties according to arguments
* ================================================= */
parentFigureUID = (char*)getCurrentFigure();
psubwinUID = (char*)getCurrentSubWin();
checkRedrawing();
/* Force 3D view */
view = 1;
setGraphicObjectProperty(psubwinUID, __GO_VIEW__, &view, jni_int, 1);
if (legend != NULL)
{
int textDimensions[2] = {1, 1};
/* F.Leray 25.04.05 replace the default labels by the user labels if specified */
loc = (char *) MALLOC((strlen(legend) + 1) * sizeof(char));
if (loc == NULL)
{
Scierror(999, _("%s: No more memory.\n"), "Objplot3d");
}
strcpy(loc, legend);
/* legx=strtok_r(loc,"@",&buff); */
legx = strtok(loc, "@");
if (legx != NULL)
{
getGraphicObjectProperty(psubwinUID, __GO_X_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_TEXT_ARRAY_DIMENSIONS__, textDimensions, jni_int_vector, 2);
setGraphicObjectProperty(labelId, __GO_TEXT_STRINGS__, &legx, jni_string_vector, textDimensions[0]*textDimensions[1]);
}
/* legy=strtok_r((char *)0,"@",&buff); */
legy = strtok((char *)NULL, "@"); /* NULL to begin at the last read character */
if (legy != NULL)
{
getGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_TEXT_ARRAY_DIMENSIONS__, textDimensions, jni_int_vector, 2);
setGraphicObjectProperty(labelId, __GO_TEXT_STRINGS__, &legy, jni_string_vector, textDimensions[0]*textDimensions[1]);
}
/* legz=strtok_r((char *)0,"@",&buff); */
legz = strtok((char *)NULL, "@");
if (legz != NULL)
{
getGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_TEXT_ARRAY_DIMENSIONS__, textDimensions, jni_int_vector, 2);
setGraphicObjectProperty(labelId, __GO_TEXT_STRINGS__, &legz, jni_string_vector, textDimensions[0]*textDimensions[1]);
}
}
/* Force psubwin->logflags to linear */
linLogFlag = 0;
setGraphicObjectProperty(psubwinUID, __GO_X_AXIS_LOG_FLAG__, &linLogFlag, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_LOG_FLAG__, &linLogFlag, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_LOG_FLAG__, &linLogFlag, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_FIRST_PLOT__, jni_bool, (void **)&piFirstPlot);
if (firstPlot == 0 && (iflag[2] == 0 || iflag[2] == 1))
{
/* Nothing to do: we leave as before */
}
else
{
int labelVisible;
if (iflag[2] == 0 || iflag[2] == 1)
{
if (firstPlot)
{
/* 0: OFF */
box = 0;
axisVisible = 0;
setGraphicObjectProperty(psubwinUID, __GO_X_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_BOX_TYPE__, &box, jni_int, 1);
labelVisible = 0;
getGraphicObjectProperty(psubwinUID, __GO_X_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
}
/*else no changes : the axes visible properties are driven by the previous plot */
}
else if (iflag[2] == 2)
{
/* 2: HIDDEN_AXES */
box = 2;
/* for 2d use only (when switching to 2d mode) */
setGraphicObjectProperty(psubwinUID, __GO_BOX_TYPE__, &box, jni_int, 1);
axisVisible = 0;
setGraphicObjectProperty(psubwinUID, __GO_X_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
labelVisible = 0;
getGraphicObjectProperty(psubwinUID, __GO_X_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
}
else if (iflag[2] == 3)
{
/* 1: ON */
box = 1;
/* for 2d use only (when switching to 2d mode) */
setGraphicObjectProperty(psubwinUID, __GO_BOX_TYPE__, &box, jni_int, 1);
axisVisible = 0;
setGraphicObjectProperty(psubwinUID, __GO_X_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
labelVisible = 1;
getGraphicObjectProperty(psubwinUID, __GO_X_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
}
else if (iflag[2] == 4)
{
/* 1: ON */
box = 1;
setGraphicObjectProperty(psubwinUID, __GO_BOX_TYPE__, &box, jni_int, 1);
axisVisible = 1;
setGraphicObjectProperty(psubwinUID, __GO_X_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
setGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_VISIBLE__, &axisVisible, jni_bool, 1);
labelVisible = 1;
getGraphicObjectProperty(psubwinUID, __GO_X_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Y_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
getGraphicObjectProperty(psubwinUID, __GO_Z_AXIS_LABEL__, jni_string, (void **)&labelId);
setGraphicObjectProperty(labelId, __GO_VISIBLE__, &labelVisible, jni_bool, 1);
}
}
rotationAngles[0] = *alpha;
rotationAngles[1] = *theta;
setGraphicObjectProperty(psubwinUID, __GO_ROTATION_ANGLES__, rotationAngles, jni_double_vector, 2);
getGraphicObjectProperty(psubwinUID, __GO_DATA_BOUNDS__, jni_double_vector, (void **)&dataBounds);
getGraphicObjectProperty(psubwinUID, __GO_AUTO_SCALE__, jni_bool, (void **)&piAutoScale);
if (autoScale)
{
/* compute and merge new specified bounds with data bounds */
switch (iflag[1])
{
case 0: /* do not change data bounds */
break;
case 1 :
case 3 :
case 5 :
case 7 : /* Force data bounds=ebox */
drect[0] = ebox[0]; /*xmin*/
drect[2] = ebox[2]; /*ymin*/
drect[1] = ebox[1]; /*xmax*/
drect[3] = ebox[3]; /*ymax*/
drect[4] = ebox[4]; /*zmin*/
drect[5] = ebox[5]; /*zmax*/
break;
case 2 :
case 4 :
case 6 :
case 8:/* Force data bounds to the x and y bounds */
getDrect(x, (*m1) * (*n1), &drect[0], &drect[1], dataBounds[0], dataBounds[1]);
getDrect(y, (*m2) * (*n2), &drect[2], &drect[3], dataBounds[2], dataBounds[3]);
getDrect(z, (*m3) * (*n3), &drect[4], &drect[5], dataBounds[4], dataBounds[5]);
break;
}
/* merge data bounds and drect */
if (!firstPlot)
{
drect[0] = Min(dataBounds[0], drect[0]); /* xmin */
drect[1] = Max(dataBounds[1], drect[1]); /* xmax */
drect[2] = Min(dataBounds[2], drect[2]); /* ymin */
drect[3] = Max(dataBounds[3], drect[3]); /* ymax */
drect[4] = Min(dataBounds[4], drect[4]); /* zmin */
drect[5] = Max(dataBounds[5], drect[5]); /* zmax */
}
if (iflag[1] != 0)
{
setGraphicObjectProperty(psubwinUID, __GO_DATA_BOUNDS__, drect, jni_double_vector, 6);
}
}
if (iflag[1] != 0)
{
isoview = (iflag[1] == 3 || iflag[1] == 4 || iflag[1] == 5 || iflag[1] == 6);
setGraphicObjectProperty(psubwinUID, __GO_ISOVIEW__, &isoview, jni_bool, 1);
}
/* =================================================
* Analyze arguments to find entity type
* ================================================= */
if (*isfac == 1)
{
if (*izcol == 0)
{
if (strcmp(fname, "plot3d1") == 0)
{
typeof3d = SCI_FAC3D;
flagcolor = 1;
}
else
{
typeof3d = SCI_FAC3D;
flagcolor = 0;
}
}
else if (*izcol == 2)
{
typeof3d = SCI_FAC3D;
flagcolor = 3;
}
else
{
typeof3d = SCI_FAC3D;
flagcolor = 2;
}
}
else if (*isfac == 0)
{
if (strcmp(fname, "plot3d1") == 0)
{
typeof3d = SCI_PLOT3D;
flagcolor = 1;
}
else
{
typeof3d = SCI_PLOT3D;
flagcolor = 0;
}
}
else
{
typeof3d = SCI_PARAM3D1;
flagcolor = 1;
}
/* =================================================
* Construct the Entities
* ================================================= */
/*Distinction here between SCI_PARAM3D1 and others*/
if (typeof3d != SCI_PARAM3D1)
{
if (*isfac == 1)
{
/* x is considered as a matrix */
dimvectx = -1;
}
else if (*m1 == 1) /* x is a row vector */
{
dimvectx = *n1;
}
else if (*n1 == 1) /* x is a column vector */
{
dimvectx = *m1;
}
else /* x is a matrix */
{
dimvectx = -1;
}
if (dimvectx > 1)
{
int monotony = checkMonotony(x, dimvectx);
if (monotony == 0)
{
Scierror(999, _("%s: x vector is not monotonous.\n"), "Objplot3d");
return;
}
flag_x = monotony;
}
if (*isfac == 1)
{
/* y is considered as a matrix */
dimvecty = -1;
}
else if (*m2 == 1) /* y is a row vector */
{
dimvecty = *n2;
}
else if (*n2 == 1) /* y is a column vector */
{
dimvecty = *m2;
}
else /* y is a matrix */
{
dimvecty = -1;
}
if (dimvecty > 1)
{
/* test the monotony on y*/
int monotony = checkMonotony(y, dimvecty);
if (monotony == 0)
{
Scierror(999, _("%s: y vector is not monotonous.\n"), "Objplot3d");
return;
}
flag_y = monotony;
}
pNewSurfaceUID = ConstructSurface(psubwinUID, typeof3d,
x, y, z, zcol, *izcol, *m, *n, iflag, ebox, flagcolor,
isfac, m1, n1, m2, n2, m3, n3, m3n, n3n);
if (pNewSurfaceUID == NULL)
{
Scierror(999, _("%s: No more memory.\n"), "Objplot3d");
return;
}
setCurrentObject(pNewSurfaceUID);
/* Force clipping, 1: CLIPGRF */
clipState = 1;
setGraphicObjectProperty(pNewSurfaceUID, __GO_CLIP_STATE__, &clipState, jni_int, 1);
releaseGraphicObjectProperty(__GO_PARENT__, pNewSurfaceUID, jni_string, 1);
}
else
{
char* pNewPolylineUID = NULL;
char* currentSubwinUID = NULL;
if ((hdltab = MALLOC (*n * sizeof (long))) == NULL)
{
Scierror(999, "%s: No more memory.\n", fname);
return;
}
currentSubwinUID = (char*)getCurrentSubWin();
for (i = 0; i < *n; ++i)
{
/* F.Leray Pb here: In fact we do not create a Surface but one or several 3D Polylines
Pb comes when wanting to access the fields "surface_color" or "flag" for example
in function sciSet (cf. matdes.c).
Question 1: Are these properties accessible from a SCI_PARAM3D1 ?
Question 2: Is "flag" obsolete and replaced by "color_mode"?*/
if ((*n > 0) && (zcol != (double *)NULL))
{
if ((int) zcol[i] > 0)
{
int intzcol = (int) zcol[i];
pNewPolylineUID = ConstructPolyline
(currentSubwinUID,
&(x[*m * i]), &(y[*m * i]), &(z[*m * i]), 0, *m, 1,
&intzcol, NULL, NULL, NULL, NULL, TRUE, FALSE, FALSE, FALSE);
}
else
{
int intzcol = (int) - zcol[i];
pNewPolylineUID = ConstructPolyline
(currentSubwinUID,
&(x[*m * i]), &(y[*m * i]), &(z[*m * i]), 0, *m, 1,
NULL, NULL, &intzcol, NULL, NULL, FALSE, FALSE, TRUE, FALSE);
}
}
else
{
/* default case, nothing is given */
int curcolor = 0;
int *piCurColor = &curcolor;
getGraphicObjectProperty(currentSubwinUID, __GO_LINE_COLOR__, jni_int, (void**)&piCurColor);
pNewPolylineUID = ConstructPolyline(currentSubwinUID,
&(x[*m * i]), &(y[*m * i]), &(z[*m * i]), 0, *m, 1,
&curcolor, NULL, NULL, NULL, NULL, TRUE, FALSE, FALSE, FALSE);
}
if (pNewPolylineUID == NULL)
{
Scierror(999, _("%s: No more memory.\n"), fname);
FREE(hdltab);
return;
}
setCurrentObject(pNewPolylineUID);
setGraphicObjectRelationship(currentSubwinUID, pNewPolylineUID);
releaseGraphicObjectProperty(__GO_PARENT__, pNewPolylineUID, jni_string, 1);
pNewPolylineUID = NULL;
pobjUID = (char*)getCurrentObject();
/* Force clipping, 1: CLIPGRF */
clipState = 1;
setGraphicObjectProperty(pobjUID, __GO_CLIP_STATE__, &clipState, jni_int, 1);
hdltab[i] = getHandle(pobjUID);
}
/** construct Compound and make it current object**/
if (*n > 1)
{
char* o = ConstructCompound (hdltab, *n);
setCurrentObject(o);
releaseGraphicObjectProperty(__GO_PARENT__, o, jni_string, 1);
}
FREE(hdltab);
}
/* =================================================
* Redraw Figure
* ================================================= */
// subwin has been modified
firstPlot = 0;
setGraphicObjectProperty(psubwinUID, __GO_FIRST_PLOT__, &firstPlot, jni_bool, 1);
FREE(loc);
loc = NULL;
}
/* set3ddata(pobj,cstk(l2), &l3, &numrow3, &numcol3) */
int set3ddata( char* pobjUID, AssignedList * tlist )
{
char* type;
int m1, n1, m2, n2, m3, n3;
int m3n, n3n;
int isFac3d;
double * pvecx = NULL;
double * pvecy = NULL;
double * pvecz = NULL;
int dimvectx = 0;
int dimvecty = 0;
double* inputColors;
int nbInputColors;
// number of specified colors
int nc = 0;
int izcol;
/* no copy now we just perform tests on the matrices */
pvecx = getCurrentDoubleMatrixFromList( tlist, &m1, &n1 );
pvecy = getCurrentDoubleMatrixFromList( tlist, &m2, &n2 );
pvecz = getCurrentDoubleMatrixFromList( tlist, &m3, &n3 );
if ( m1 * n1 == m3 * n3 && m1 * n1 == m2 * n2 && m1 * n1 != 1 )
{
if ( !(m1 == m2 && m2 == m3 && n1 == n2 && n2 == n3) )
{
Scierror(999, _("%s: Wrong size for arguments #%d, #%d and #%d: Incompatible length.\n"),"Tlist",1,2,3);
return SET_PROPERTY_ERROR;
}
}
else
{
if ( m2 * n2 != n3 )
{
Scierror(999, _("%s: Wrong size for arguments #%d and #%d: Incompatible length.\n"),"Tlist",2,3);
return SET_PROPERTY_ERROR;
}
if ( m1 * n1 != m3 )
{
Scierror(999, _("%s: Wrong size for arguments #%d and #%d: Incompatible length.\n"),"Tlist",1,3);
return SET_PROPERTY_ERROR;
}
if ( m1 * n1 <= 1 || m2 * n2 <= 1 )
{
Scierror(999, _("%s: Wrong size for arguments #%d and #%d: Should be >= %d.\n"),"Tlist",1,2,2);
return SET_PROPERTY_ERROR;
}
}
if ( m1 * n1 == 0 || m2 * n2 == 0 || m3 * n3 == 0 )
{
return sciReturnEmptyMatrix();
}
/* get color size if exists */
if ( getAssignedListNbElement( tlist ) == 4 )
{
getCurrentDoubleMatrixFromList( tlist, &m3n, &n3n ) ;
if ( m3n * n3n == m3 * n3 )
{
/* the color is a matrix, with same size as Z */
izcol = 2;
}
else if (m3n * n3n == n3 && (m3n == 1 || n3n == 1))
{
/* a vector with as many colors as facets */
izcol = 1;
}
else
{
Scierror(999, _("Wrong size for %s element: A %d-by-%d matrix or a vector of size %d expected.\n"), "color", m3, n3, n3);
return SET_PROPERTY_ERROR;
}
}
else
{
m3n = 0;
n3n = 0;
izcol = 0;
}
getGraphicObjectProperty(pobjUID, __GO_TYPE__, jni_string, &type);
if (strcmp(type, __GO_FAC3D__) == 0)
{
isFac3d = 1;
}
else
{
isFac3d = 0;
}
if ( m1 * n1 == m3 * n3 && m1 * n1 == m2 * n2 && m1 * n1 != 1 ) /* NG beg */
{
/* case isfac=1;*/
if (isFac3d == 0)
{
Scierror(999, _("Can not change the %s of graphic object: its type is %s.\n"),"typeof3d","SCI_PLOT3D");
return SET_PROPERTY_ERROR;
}
}
else
{
/* case isfac=0;*/
if (isFac3d == 1)
{
Scierror(999, _("Can not change the %s of graphic object: its type is %s.\n"),"typeof3d","SCI_FAC3D");
return SET_PROPERTY_ERROR;
}
}
/* check the monotony on x and y */
if (isFac3d == 1)
{
/* x is considered as a matrix */
dimvectx = -1;
}
else if ( m1 == 1 ) /* x is a row vector */
{
dimvectx = n1;
}
else if ( n1 == 1 ) /* x is a column vector */
{
dimvectx = m1;
}
else /* x is a matrix */
{
dimvectx = -1;
}
if ( dimvectx > 1 )
{
int monotony = checkMonotony( pvecx, dimvectx );
if ( monotony == 0 )
{
Scierror(999, _("%s: Wrong value: Vector is not monotonous.\n"),"Objplot3d");
return SET_PROPERTY_ERROR;
}
/* To be implemented within the MVC */
#if 0
psurf->flag_x = monotony;
#endif
}
if (isFac3d == 1)
{
/* x is considered as a matrix */
dimvecty = -1;
}
else if(m2 == 1) /* y is a row vector */
{
dimvecty = n2;
}
else if(n2 == 1) /* y is a column vector */
{
dimvecty = m2;
}
else /* y is a matrix */
{
dimvecty = -1;
}
if( dimvecty > 1 )
{
int monotony = checkMonotony( pvecy, dimvecty );
if ( monotony == 0 )
{
Scierror(999, _("%s: Wrong value: Vector is not monotonous.\n"),"Objplot3d");
return SET_PROPERTY_ERROR;
}
/* To be implemented within the MVC */
#if 0
psurf->flag_y = monotony;
#endif
}
/* get the values now */
rewindAssignedList( tlist );
pvecx = getCurrentDoubleMatrixFromList( tlist, &m1, &n1 );
pvecy = getCurrentDoubleMatrixFromList( tlist, &m2, &n2 );
pvecz = getCurrentDoubleMatrixFromList( tlist, &m3, &n3 );
if (isFac3d == 1)
{
int numElementsArray[3];
int result;
numElementsArray[0] = n1;
numElementsArray[1] = m1;
numElementsArray[2] = m3n * n3n;
result = setGraphicObjectProperty(pobjUID, __GO_DATA_MODEL_NUM_ELEMENTS_ARRAY__, numElementsArray, jni_int_vector, 3);
if (result == 0)
{
Scierror(999, _("%s: No more memory.\n"), "set3ddata");
return SET_PROPERTY_ERROR;
}
}
else if (isFac3d == 0)
{
int gridSize[4];
int result;
gridSize[0] = m1;
gridSize[1] = n1;
gridSize[2] = m2;
gridSize[3] = n2;
result = setGraphicObjectProperty(pobjUID, __GO_DATA_MODEL_GRID_SIZE__, gridSize, jni_int_vector, 4);
if (result == 0)
{
Scierror(999, _("%s: No more memory.\n"), "set3ddata");
return SET_PROPERTY_ERROR;
}
}
setGraphicObjectProperty(pobjUID, __GO_DATA_MODEL_X__, pvecx, jni_double_vector, m1*n1);
setGraphicObjectProperty(pobjUID, __GO_DATA_MODEL_Y__, pvecy, jni_double_vector, m2*n2);
setGraphicObjectProperty(pobjUID, __GO_DATA_MODEL_Z__, pvecz, jni_double_vector, m3*n3);
if( getAssignedListNbElement( tlist ) == 4 ) /* F.Leray There is a color matrix */
{
inputColors = getCurrentDoubleMatrixFromList( tlist, &m3n, &n3n );
nbInputColors = m3n * n3n;
}
else
{
inputColors = NULL;
nbInputColors = 0;
}
/*
* Plot 3d case not treated for now
* To be implemented
*/
if (isFac3d == 1)
{
setGraphicObjectProperty(pobjUID, __GO_DATA_MODEL_COLORS__, inputColors, jni_double_vector, nbInputColors);
}
/* Color vector/matrix dimensions: to be checked for MVC implementation */
#if 0
psurf->m3n = m3n; /* If m3n and n3n are 0, then it means that no color matrix/vector was in input*/
psurf->n3n = n3n;
#endif
return SET_PROPERTY_SUCCEED;
}
