Changeset 3882
- Timestamp:
- Aug 21, 2013 11:47:45 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/gui/src/RpDxToVtk.c
r3871 r3882 24 24 25 25 #define DO_WEDGES 26 //#define CHECK_WINDINGS 26 27 27 28 static INLINE char * … … 57 58 const char *p; 58 59 char mesg[2000]; 59 double *array , scale, vmin, vmax;60 double *array; 60 61 int iX, iY, iZ; 61 62 … … 65 66 return TCL_ERROR; 66 67 } 67 vmin = DBL_MAX, vmax = -DBL_MAX;68 68 iX = iY = iZ = 0; 69 69 for (i = 0; i < nPoints; i++) { … … 93 93 } 94 94 array[loc] = value; 95 if (value < vmin) { 96 vmin = value; 97 } 98 if (value > vmax) { 99 vmax = value; 100 } 101 } 102 scale = 1.0 / (vmax - vmin); 95 } 103 96 for (i = 0; i < nPoints; i++) { 104 #ifdef notdef105 sprintf(mesg, "%g\n", (array[i] - vmin) * scale);106 #endif107 97 sprintf(mesg, "%g\n", array[i]); 108 98 Tcl_AppendToObj(objPtr, mesg, -1); … … 120 110 const char *p; 121 111 char mesg[2000]; 122 double *array , scale, vmin, vmax;112 double *array; 123 113 124 114 p = *stringPtr; … … 127 117 return TCL_ERROR; 128 118 } 129 vmin = DBL_MAX, vmax = -DBL_MAX;130 119 for (i = 0; i < nPoints; i++) { 131 120 double value; … … 145 134 p = nextPtr; 146 135 array[i] = value; 147 if (value < vmin) { 148 vmin = value; 149 } 150 if (value > vmax) { 151 vmax = value; 152 } 153 } 154 scale = 1.0 / (vmax - vmin); 136 } 155 137 for (i = 0; i < nPoints; i++) { 156 #ifdef notdef157 sprintf(mesg, "%g\n", (array[i] - vmin) * scale);158 #endif159 138 sprintf(mesg, "%g\n", array[i]); 160 139 Tcl_AppendToObj(objPtr, mesg, -1); … … 172 151 const char *p; 173 152 char mesg[2000]; 174 double *array , scale, vmin, vmax;153 double *array; 175 154 int iXY, iZ; 176 155 int nPoints; … … 183 162 return TCL_ERROR; 184 163 } 185 vmin = DBL_MAX, vmax = -DBL_MAX;186 164 iXY = iZ = 0; 187 165 for (i = 0; i < nPoints; i++) { … … 208 186 } 209 187 array[loc] = value; 210 if (value < vmin) { 211 vmin = value; 212 } 213 if (value > vmax) { 214 vmax = value; 215 } 216 } 217 scale = 1.0 / (vmax - vmin); 188 } 218 189 for (i = 0; i < nPoints; i++) { 219 #ifdef notdef220 sprintf(mesg, "%g\n", (array[i] - vmin) * scale);221 #endif222 190 sprintf(mesg, "%g\n", array[i]); 223 191 Tcl_AppendToObj(objPtr, mesg, -1); … … 278 246 return TCL_OK; 279 247 } 280 248 #if defined(DO_WEDGES) && defined(CHECK_WINDINGS) 281 249 static void 282 250 Normalize(double v[3]) … … 301 269 vout[2] = v1[0]*v2[1] - v1[1]*v2[0]; 302 270 } 303 271 #endif 304 272 /* 305 273 * DxToVtk string … … 565 533 for (iz = 0; iz < count[2]; iz++) { 566 534 for (i = 0; i < nXYPoints; i++) { 567 sprintf(mesg, "%g %g %g\n", points[i*2], points[i*2+1], origin[2] + dz * iz);535 sprintf(mesg, "%g %g %g\n", points[i*2], points[i*2+1], (origin[2] + dz * iz)); 568 536 Tcl_AppendToObj(pointsObjPtr, mesg, -1); 569 537 } … … 616 584 } 617 585 fclose(ftmp); 618 #ifdef notdef586 #ifdef CHECK_WINDINGS 619 587 double normal[3]; 620 588 normal[0] = normal[1] = 0.0; … … 637 605 // skip boundary points we added 638 606 if (c0 >= 0 && c1 >= 0 && c2 >= 0) { 639 #ifdef notdef607 #ifdef CHECK_WINDINGS 640 608 /* Winding of base triangle should point to 641 609 outside of cell using right hand rule */ … … 675 643 unlink(fcells); 676 644 } 677 free(points);678 645 sprintf(mesg, "CELLS %d %d\n", nCells, 7*nCells); 679 646 Tcl_AppendToObj(objPtr, mesg, -1); … … 686 653 } 687 654 #endif 655 if (points != NULL) { 656 free(points); 657 } 688 658 sprintf(mesg, "POINT_DATA %d\n", nPoints); 689 659 Tcl_AppendToObj(objPtr, mesg, -1);
Note: See TracChangeset
for help on using the changeset viewer.