147 {
149 long iArg;
150 SCANNED_ARG *scArg;
151 char *input = NULL, *output = NULL, *columnName = NULL;
152 long mode = -1, matchCode, tmpfile_used;
153 int64_t i, j, rows, rowLimit, pagesPerPage=0;
154 long breakNext;
155 double gapAmount = 0, gapFactor = 0, changeAmount, changeBase, *columnData = NULL;
156 char *matchPattern = NULL;
157 long matchPatternAfter = 0;
158 double increaseOfAmount = -1, decreaseOfAmount = -1;
159 long retval;
160 int32_t dataType;
161 long overlap = 0;
162 unsigned long flags = 0, pipeFlags = 0, changeFlags = 0, decreaseOfFlags = 0, increaseOfFlags = 0;
163 unsigned long majorOrderFlag;
164 char **stringData = NULL;
165 short columnMajorOrder = -1;
166
168 argc =
scanargs(&scArg, argc, argv);
169 if (argc < 2) {
170 fprintf(stderr, "%s", USAGE);
171 return EXIT_FAILURE;
172 }
173
174 for (iArg = 1; iArg < argc; iArg++) {
175 if (scArg[iArg].arg_type == OPTION) {
176 switch (matchCode =
match_string(scArg[iArg].list[0], option, N_OPTIONS, 0)) {
177 case SET_MAJOR_ORDER:
178 majorOrderFlag = 0;
179 scArg[iArg].n_items -= 1;
180 if (scArg[iArg].n_items > 0 &&
181 !
scanItemList(&majorOrderFlag, scArg[iArg].list + 1, &scArg[iArg].n_items, 0,
182 "row", -1, NULL, 0, SDDS_ROW_MAJOR_ORDER,
183 "column", -1, NULL, 0, SDDS_COLUMN_MAJOR_ORDER, NULL)) {
184 SDDS_Bomb(
"invalid -majorOrder syntax/values");
185 }
186 if (majorOrderFlag & SDDS_COLUMN_MAJOR_ORDER)
187 columnMajorOrder = 1;
188 else if (majorOrderFlag & SDDS_ROW_MAJOR_ORDER)
189 columnMajorOrder = 0;
190 break;
191
192 case SET_GAPIN:
193 if ((scArg[iArg].n_items -= 2) < 0 ||
194 !
scanItemList(&flags, scArg[iArg].list + 2, &scArg[iArg].n_items, 0,
195 "amount",
SDDS_DOUBLE, &gapAmount, 1, GAPIN_AMOUNT,
196 "factor",
SDDS_DOUBLE, &gapFactor, 1, GAPIN_FACTOR, NULL) ||
197 (flags & GAPIN_AMOUNT && gapAmount <= 0) ||
198 (flags & GAPIN_FACTOR && gapFactor <= 0)) {
199 fprintf(stderr, "Error: invalid -gapin syntax/values\n");
200 return EXIT_FAILURE;
201 }
202 columnName = scArg[iArg].list[1];
203 mode = matchCode;
204 break;
205
206 case SET_INCREASEOF:
207 increaseOfFlags = 0;
208 if (scArg[iArg].n_items < 2) {
209 fprintf(stderr, "Error: invalid -increaseOf syntax\n");
210 return EXIT_FAILURE;
211 }
212 scArg[iArg].n_items -= 2;
213 if (!
scanItemList(&increaseOfFlags, scArg[iArg].list + 2, &scArg[iArg].n_items, 0,
214 "amount",
SDDS_DOUBLE, &increaseOfAmount, 1, INCREASEOF_AMOUNT,
215 "cumulative", -1, NULL, 0, INCREASEOF_CUMULATIVE,
216 "reset", -1, NULL, 0, INCREASEOF_RESET, NULL) ||
217 ((flags & INCREASEOF_AMOUNT) && increaseOfAmount <= 0)) {
218 fprintf(stderr, "Error: invalid -increaseOf syntax\n");
219 return EXIT_FAILURE;
220 }
221 columnName = scArg[iArg].list[1];
222 mode = matchCode;
223 break;
224
225 case SET_DECREASEOF:
226 decreaseOfFlags = 0;
227 if (scArg[iArg].n_items < 2) {
228 fprintf(stderr, "Error: invalid -decreaseOf syntax\n");
229 return EXIT_FAILURE;
230 }
231 scArg[iArg].n_items -= 2;
232 if (!
scanItemList(&decreaseOfFlags, scArg[iArg].list + 2, &scArg[iArg].n_items, 0,
233 "amount",
SDDS_DOUBLE, &decreaseOfAmount, 1, DECREASEOF_AMOUNT,
234 "cumulative", -1, NULL, 0, DECREASEOF_CUMULATIVE,
235 "reset", -1, NULL, 0, DECREASEOF_RESET, NULL) ||
236 ((flags & DECREASEOF_AMOUNT) && decreaseOfAmount <= 0)) {
237 fprintf(stderr, "Error: invalid -decreaseOf syntax\n");
238 return EXIT_FAILURE;
239 }
240 columnName = scArg[iArg].list[1];
241 mode = matchCode;
242 break;
243
244 case SET_CHANGEOF:
245 if ((scArg[iArg].n_items -= 2) < 0 ||
246 !
scanItemList(&changeFlags, scArg[iArg].list + 2, &scArg[iArg].n_items, 0,
247 "amount",
SDDS_DOUBLE, &changeAmount, 1, CHANGEOF_AMOUNT,
248 "base",
SDDS_DOUBLE, &changeBase, 1, CHANGEOF_BASE, NULL) ||
249 (changeFlags & CHANGEOF_AMOUNT && changeAmount <= 0)) {
250 fprintf(stderr, "Error: invalid -changeof syntax/values\n");
251 return EXIT_FAILURE;
252 }
253 columnName = scArg[iArg].list[1];
254 mode = matchCode;
255 break;
256
257 case SET_ROWLIMIT:
258 if (scArg[iArg].n_items < 2) {
259 fprintf(stderr, "Error: invalid -rowlimit syntax\n");
260 return EXIT_FAILURE;
261 }
262 if (sscanf(scArg[iArg].list[1], "%" SCNd64, &rowLimit) != 1 || rowLimit <= 0) {
263 fprintf(stderr, "Error: invalid -rowlimit syntax\n");
264 return EXIT_FAILURE;
265 }
266 if (scArg[iArg].n_items > 2) {
267 scArg[iArg].n_items -= 2;
268 if (!
scanItemList(&flags, scArg[iArg].list + 2, &scArg[iArg].n_items, 0,
269 "overlap",
SDDS_LONG, &overlap, 1, ROWLIMIT_OVERLAP, NULL) ||
270 overlap < 0) {
271 fprintf(stderr, "Error: invalid overlap given in -rowlimit syntax\n");
272 return EXIT_FAILURE;
273 }
274 }
275 mode = matchCode;
276 break;
277
278 case SET_PIPE:
279 if (!
processPipeOption(scArg[iArg].list + 1, scArg[iArg].n_items - 1, &pipeFlags)) {
280 fprintf(stderr, "Error: invalid -pipe syntax\n");
281 return EXIT_FAILURE;
282 }
283 break;
284
285 case SET_MATCHTO:
286 if ((scArg[iArg].n_items != 3 && scArg[iArg].n_items != 4) ||
287 strlen(columnName = scArg[iArg].list[1]) == 0 ||
288 strlen(matchPattern = scArg[iArg].list[2]) == 0) {
289 fprintf(stderr, "Error: invalid -matchTo syntax\n");
290 return EXIT_FAILURE;
291 }
292 if (scArg[iArg].n_items == 4) {
293 if (strncmp(scArg[iArg].list[3], "after", strlen(scArg[iArg].list[3])) == 0)
294 matchPatternAfter = 1;
295 else {
296 fprintf(stderr, "Error: invalid -matchTo syntax\n");
297 return EXIT_FAILURE;
298 }
299 }
300 mode = matchCode;
301 break;
302
303 case SET_PAGES_PER_PAGE:
304 if (scArg[iArg].n_items != 2) {
305 fprintf(stderr, "Error: invalid -pagesPerPage syntax\n");
306 return EXIT_FAILURE;
307 }
308 if (sscanf(scArg[iArg].list[1], "%" SCNd64, &pagesPerPage) != 1 || pagesPerPage <= 0) {
309 fprintf(stderr, "Error: invalid -pagesPerPage syntax\n");
310 return EXIT_FAILURE;
311 }
312 mode = matchCode;
313 break;
314
315 default:
316 fprintf(stderr, "Error: unknown switch: %s\n", scArg[iArg].list[0]);
317 fprintf(stderr, "%s", USAGE);
318 return EXIT_FAILURE;
319 }
320 } else {
321 if (input == NULL)
322 input = scArg[iArg].list[0];
323 else if (output == NULL)
324 output = scArg[iArg].list[0];
325 else {
326 fprintf(stderr, "Error: too many filenames given\n");
327 return EXIT_FAILURE;
328 }
329 }
330 }
331
332 processFilenames(
"sddsbreak", &input, &output, pipeFlags, 0, &tmpfile_used);
333
334 if (mode == -1) {
335 fprintf(stderr, "Error: no break mode specified\n");
336 return EXIT_FAILURE;
337 }
338
342 return EXIT_FAILURE;
343 }
344
345 SDDSnew.layout.data_mode.no_row_counts = 0;
346 if (columnMajorOrder != -1)
347 SDDSnew.layout.data_mode.column_major = columnMajorOrder;
348 else
349 SDDSnew.layout.data_mode.column_major = SDDSold.layout.data_mode.column_major;
350
353 return EXIT_FAILURE;
354 }
355
356 if (mode != SET_ROWLIMIT && mode!=SET_PAGES_PER_PAGE) {
358 SDDS_SetError(
"Problem getting type information on given column");
360 return EXIT_FAILURE;
361 }
362 if (mode == SET_MATCHTO) {
364 fprintf(stderr, "Error: given column does not contain string data\n");
365 return EXIT_FAILURE;
366 }
368 if (!(mode == SET_CHANGEOF && !(changeFlags & CHANGEOF_AMOUNT) && !(changeFlags & CHANGEOF_BASE))) {
369 fprintf(stderr, "Error: given column does not contain numeric data\n");
370 return EXIT_FAILURE;
371 }
372 }
373 }
374
377 SDDS_SetError(
"Problem getting number of rows of tabular data");
379 return EXIT_FAILURE;
380 }
385 return EXIT_FAILURE;
386 }
387 if (rows == 0) {
390 return EXIT_FAILURE;
391 }
392 continue;
393 }
394
395 switch (mode) {
396 case SET_GAPIN:
398 if (!columnData) {
401 return EXIT_FAILURE;
402 }
403 if (!gapAmount && rows > 1) {
404 double *gap =
tmalloc(
sizeof(*gap) * rows);
405 for (i = 1; i < rows; i++)
406 gap[i - 1] = fabs(columnData[i] - columnData[i - 1]);
408 fprintf(stderr, "Error: unable to determine default gap amount--couldn't find median gap\n");
409 free(gap);
410 return EXIT_FAILURE;
411 }
412 gapAmount *= (gapFactor ? gapFactor : 2);
413 free(gap);
414 }
415 {
416 int64_t newStart = 0;
417 for (i = 1; i <= rows; i++) {
418 if (i != rows && fabs(columnData[i] - columnData[i - 1]) < gapAmount)
419 continue;
423 free(columnData);
424 return EXIT_FAILURE;
425 }
426 newStart = i;
427 }
428 }
429 free(columnData);
430 break;
431
432 case SET_INCREASEOF:
434 if (!columnData) {
437 return EXIT_FAILURE;
438 }
439 {
440 int64_t newStart = 0;
441 if (increaseOfAmount <= 0) {
442 for (i = 1; i <= rows; i++) {
443 if (i != rows && columnData[i] <= columnData[i - 1])
444 continue;
448 free(columnData);
449 return EXIT_FAILURE;
450 }
451 newStart = i;
452 }
453 } else {
454 if (increaseOfFlags & INCREASEOF_CUMULATIVE) {
455 long iref = 0;
456 for (i = 1; i <= rows; i++) {
457 if ((increaseOfFlags & INCREASEOF_RESET) && columnData[i] < columnData[iref])
458 iref = i;
459 if (i != rows && (columnData[i] - columnData[iref]) < increaseOfAmount)
460 continue;
464 free(columnData);
465 return EXIT_FAILURE;
466 }
467 newStart = i;
468 if (increaseOfFlags & INCREASEOF_CUMULATIVE)
469 iref = i;
470 }
471 } else {
472 for (i = 1; i <= rows; i++) {
473 if (i != rows && (columnData[i] - columnData[i - 1]) < increaseOfAmount)
474 continue;
478 free(columnData);
479 return EXIT_FAILURE;
480 }
481 newStart = i;
482 }
483 }
484 }
485 }
486 free(columnData);
487 break;
488
489 case SET_DECREASEOF:
491 if (!columnData) {
494 return EXIT_FAILURE;
495 }
496 {
497 int64_t newStart = 0;
498 if (decreaseOfAmount <= 0) {
499 for (i = 1; i <= rows; i++) {
500 if (i != rows && columnData[i] >= columnData[i - 1])
501 continue;
505 free(columnData);
506 return EXIT_FAILURE;
507 }
508 newStart = i;
509 }
510 } else {
511 if (decreaseOfFlags & DECREASEOF_CUMULATIVE) {
512 long iref = 0;
513 for (i = 1; i <= rows; i++) {
514 if ((decreaseOfFlags & DECREASEOF_RESET) && columnData[i] > columnData[iref])
515 iref = i;
516 if (i != rows && (columnData[iref] - columnData[i]) < decreaseOfAmount)
517 continue;
521 free(columnData);
522 return EXIT_FAILURE;
523 }
524 newStart = i;
525 if (decreaseOfFlags & DECREASEOF_CUMULATIVE)
526 iref = i;
527 }
528 } else {
529 for (i = 1; i <= rows; i++) {
530 if (i != rows && (columnData[i - 1] - columnData[i]) < decreaseOfAmount)
531 continue;
535 free(columnData);
536 return EXIT_FAILURE;
537 }
538 newStart = i;
539 }
540 }
541 }
542 }
543 free(columnData);
544 break;
545
546 case SET_CHANGEOF:
549 if (!columnData) {
552 return EXIT_FAILURE;
553 }
554 } else {
556 if (!stringData) {
559 return EXIT_FAILURE;
560 }
561 }
562 {
563 int64_t newStart = 0;
565 for (i = 1; i <= rows; i++) {
566 if (i != rows &&
567 ((dataType ==
SDDS_STRING && strcmp(stringData[i], stringData[i - 1]) == 0) ||
568 (dataType !=
SDDS_STRING && columnData[i] == columnData[i - 1])))
569 continue;
574 free(columnData);
575 else
577 return EXIT_FAILURE;
578 }
579 newStart = i;
580 }
581 } else {
582 long region = 0, lastRegion = 0;
583 if (!(changeFlags & CHANGEOF_BASE) && rows >= 1)
584 changeBase = columnData[0];
585 if (rows > 1)
586 lastRegion = (columnData[0] - changeBase) / changeAmount;
587
588 newStart = 0;
589 for (i = 1; i <= rows; i++) {
590 if (i != rows)
591 region = (columnData[i] - changeBase) / changeAmount;
592 if (i != rows && region == lastRegion)
593 continue;
598 free(columnData);
599 else
601 return EXIT_FAILURE;
602 }
603 newStart = i;
604 lastRegion = region;
605 }
606 }
607 }
609 free(columnData);
610 else
612 break;
613
614 case SET_MATCHTO:
616 if (!stringData) {
619 return EXIT_FAILURE;
620 }
621 {
622 int64_t newStart = 0;
623 breakNext = 0;
624 for (i = 1; i <= rows; i++) {
625 if (i != rows && !breakNext) {
626 if (
wild_match(stringData[i], matchPattern)) {
627 if (matchPatternAfter) {
628 breakNext = 1;
629 continue;
630 }
631 } else {
632 continue;
633 }
634 }
639 return EXIT_FAILURE;
640 }
641 breakNext = 0;
642 newStart = i;
643 }
644 }
646 break;
647
648 case SET_ROWLIMIT:
649 for (i = 0; i < rows; i += rowLimit - overlap) {
650 if ((j = i + rowLimit - 1) >= rows)
651 j = rows - 1;
654 return EXIT_FAILURE;
655 }
658 return EXIT_FAILURE;
659 }
660 if (j == rows - 1)
661 break;
662 }
663 break;
664
665 case SET_PAGES_PER_PAGE:
666 rowLimit = rows/pagesPerPage;
667 for (j=0; j<pagesPerPage; j++) {
668 int64_t iStart, iEnd;
669 iStart = j*rowLimit;
670 iEnd = j==(pagesPerPage-1)? rows-1 : (j+1)*rowLimit-1;
673 return EXIT_FAILURE;
674 }
677 return EXIT_FAILURE;
678 }
679 }
680 break;
681
682 default:
683 fprintf(stderr, "Error: unknown break mode code seen---this can't happen\n");
684 return EXIT_FAILURE;
685 }
686 }
687
688 if (retval == 0) {
690 return EXIT_FAILURE;
691 }
692
695 return EXIT_FAILURE;
696 }
697
700 return EXIT_FAILURE;
701 }
702
704 return EXIT_FAILURE;
705 }
706
707 return EXIT_SUCCESS;
708}
int32_t SDDS_CopyParameters(SDDS_DATASET *SDDS_target, SDDS_DATASET *SDDS_source)
int32_t SDDS_InitializeCopy(SDDS_DATASET *SDDS_target, SDDS_DATASET *SDDS_source, char *filename, char *filemode)
int32_t SDDS_CopyArrays(SDDS_DATASET *SDDS_target, SDDS_DATASET *SDDS_source)
int32_t SDDS_CopyRows(SDDS_DATASET *SDDS_target, SDDS_DATASET *SDDS_source, int64_t firstRow, int64_t lastRow)
int32_t SDDS_StartPage(SDDS_DATASET *SDDS_dataset, int64_t expected_n_rows)
int32_t SDDS_GetColumnInformation(SDDS_DATASET *SDDS_dataset, char *field_name, void *memory, int32_t mode,...)
Retrieves information about a specified column in the SDDS dataset.
int32_t SDDS_WritePage(SDDS_DATASET *SDDS_dataset)
Writes the current data table to the output file.
int32_t SDDS_WriteLayout(SDDS_DATASET *SDDS_dataset)
Writes the SDDS layout header to the output file.
int32_t SDDS_FreeStringArray(char **string, int64_t strings)
Frees an array of strings by deallocating each individual string.
void SDDS_SetError(char *error_text)
Records an error message in the SDDS error stack.
void SDDS_PrintErrors(FILE *fp, int32_t mode)
Prints recorded error messages to a specified file stream.
void SDDS_RegisterProgramName(const char *name)
Registers the executable program name for use in error messages.
void SDDS_Bomb(char *message)
Terminates the program after printing an error message and recorded errors.
#define SDDS_STRING
Identifier for the string data type.
#define SDDS_LONG
Identifier for the signed 32-bit integer data type.
#define SDDS_DOUBLE
Identifier for the double data type.
#define SDDS_NUMERIC_TYPE(type)
Checks if the given type identifier corresponds to any numeric type.
void * tmalloc(uint64_t size_of_block)
Allocates a memory block of the specified size with zero initialization.
long match_string(char *string, char **option, long n_options, long mode)
Matches a given string against an array of option strings based on specified modes.
long replaceFileAndBackUp(char *file, char *replacement)
Replaces a file with a replacement file and creates a backup of the original.
int scanargs(SCANNED_ARG **scanned, int argc, char **argv)
long processPipeOption(char **item, long items, unsigned long *flags)
void processFilenames(char *programName, char **input, char **output, unsigned long pipeFlags, long noWarnings, long *tmpOutputUsed)
long scanItemList(unsigned long *flags, char **item, long *items, unsigned long mode,...)
Scans a list of items and assigns values based on provided keywords and types.
int wild_match(char *string, char *template)
Determine whether one string is a wildcard match for another.