Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Use pointer arithmetic for accessing to dynar elements.
[simgrid.git] / src / xbt / xbt_str.c
1 /* xbt_str.c - various helping functions to deal with strings               */
2
3 /* Copyright (c) 2007, 2008, 2009, 2010. The SimGrid Team.
4  * All rights reserved.                                                     */
5
6 /* This program is free software; you can redistribute it and/or modify it
7  * under the terms of the license (GNU LGPL) which comes with this package. */
8
9 #include "portable.h"
10 #include "xbt/misc.h"
11 #include "xbt/sysdep.h"
12 #include "xbt/str.h"            /* headers of these functions */
13 #include "xbt/strbuff.h"
14 #include "xbt/matrix.h"         /* for the diff */
15
16 /**  @brief Strip whitespace (or other characters) from the end of a string.
17  *
18  * Strips the whitespaces from the end of s.
19  * By default (when char_list=NULL), these characters get stripped:
20  *
21  *      - " "           (ASCII 32       (0x20)) space.
22  *      - "\t"          (ASCII 9        (0x09)) tab.
23  *      - "\n"          (ASCII 10       (0x0A)) line feed.
24  *      - "\r"          (ASCII 13       (0x0D)) carriage return.
25  *      - "\0"          (ASCII 0        (0x00)) NULL.
26  *      - "\x0B"        (ASCII 11       (0x0B)) vertical tab.
27  *
28  * @param s The string to strip. Modified in place.
29  * @param char_list A string which contains the characters you want to strip.
30  *
31  */
32 void xbt_str_rtrim(char *s, const char *char_list)
33 {
34   char *cur = s;
35   const char *__char_list = " \t\n\r\x0B";
36   char white_char[256] = { 1, 0 };
37
38   if (!s)
39     return;
40
41   if (!char_list) {
42     while (*__char_list) {
43       white_char[(unsigned char) *__char_list++] = 1;
44     }
45   } else {
46     while (*char_list) {
47       white_char[(unsigned char) *char_list++] = 1;
48     }
49   }
50
51   while (*cur)
52     ++cur;
53
54   while ((cur >= s) && white_char[(unsigned char) *cur])
55     --cur;
56
57   *++cur = '\0';
58 }
59
60 /**  @brief Strip whitespace (or other characters) from the beginning of a string.
61  *
62  * Strips the whitespaces from the begining of s.
63  * By default (when char_list=NULL), these characters get stripped:
64  *
65  *      - " "           (ASCII 32       (0x20)) space.
66  *      - "\t"          (ASCII 9        (0x09)) tab.
67  *      - "\n"          (ASCII 10       (0x0A)) line feed.
68  *      - "\r"          (ASCII 13       (0x0D)) carriage return.
69  *      - "\0"          (ASCII 0        (0x00)) NULL.
70  *      - "\x0B"        (ASCII 11       (0x0B)) vertical tab.
71  *
72  * @param s The string to strip. Modified in place.
73  * @param char_list A string which contains the characters you want to strip.
74  *
75  */
76 void xbt_str_ltrim(char *s, const char *char_list)
77 {
78   char *cur = s;
79   const char *__char_list = " \t\n\r\x0B";
80   char white_char[256] = { 1, 0 };
81
82   if (!s)
83     return;
84
85   if (!char_list) {
86     while (*__char_list) {
87       white_char[(unsigned char) *__char_list++] = 1;
88     }
89   } else {
90     while (*char_list) {
91       white_char[(unsigned char) *char_list++] = 1;
92     }
93   }
94
95   while (*cur && white_char[(unsigned char) *cur])
96     ++cur;
97
98   memmove(s, cur, strlen(cur) + 1);
99 }
100
101 /**  @brief Strip whitespace (or other characters) from the end and the begining of a string.
102  *
103  * Strips the whitespaces from both the beginning and the end of s.
104  * By default (when char_list=NULL), these characters get stripped:
105  *
106  *      - " "           (ASCII 32       (0x20)) space.
107  *      - "\t"          (ASCII 9        (0x09)) tab.
108  *      - "\n"          (ASCII 10       (0x0A)) line feed.
109  *      - "\r"          (ASCII 13       (0x0D)) carriage return.
110  *      - "\0"          (ASCII 0        (0x00)) NULL.
111  *      - "\x0B"        (ASCII 11       (0x0B)) vertical tab.
112  *
113  * @param s The string to strip.
114  * @param char_list A string which contains the characters you want to strip.
115  *
116  */
117 void xbt_str_trim(char *s, const char *char_list)
118 {
119
120   if (!s)
121     return;
122
123   xbt_str_rtrim(s, char_list);
124   xbt_str_ltrim(s, char_list);
125 }
126
127 /**  @brief Replace double whitespaces (but no other characters) from the string.
128  *
129  * The function modifies the string so that each time that several spaces appear,
130  * they are replaced by a single space. It will only do so for spaces (ASCII 32, 0x20).
131  *
132  * @param s The string to strip. Modified in place.
133  *
134  */
135 void xbt_str_strip_spaces(char *s)
136 {
137   char *p = s;
138   int e = 0;
139
140   if (!s)
141     return;
142
143   while (1) {
144     if (!*p)
145       goto end;
146
147     if (*p != ' ')
148       break;
149
150     p++;
151   }
152
153   e = 1;
154
155   do {
156     if (e)
157       *s++ = *p;
158
159     if (!*++p)
160       goto end;
161
162     if (e ^ (*p != ' '))
163       if ((e = !e))
164         *s++ = ' ';
165   } while (1);
166
167 end:
168   *s = '\0';
169 }
170
171 /** @brief Substitutes a char for another in a string
172  *
173  * @param str the string to modify
174  * @param from char to search
175  * @param to char to put instead
176  * @param occurence number of changes to do (=0 means all)
177  */
178 void xbt_str_subst(char *str, char from, char to, int occurence)
179 {
180   char *p = str;
181   while (*p != '\0') {
182     if (*p == from) {
183       *p = to;
184       if (occurence == 1)
185         return;
186       occurence--;
187     }
188     p++;
189   }
190 }
191
192 /** @brief Replaces a set of variables by their values
193  *
194  * @param str The input of the replacement process
195  * @param patterns The changes to apply
196  * @return The string modified
197  *
198  * Both '$toto' and '${toto}' are valid (and the two writing are equivalent).
199  *
200  * If the variable name contains spaces, use the brace version (ie, ${toto tutu})
201  *
202  * You can provide a default value to use if the variable is not set in the dict by using
203  * '${var:=default}' or '${var:-default}'. These two forms are equivalent, even if they
204  * shouldn't to respect the shell standard (:= form should set the value in the dict,
205  * but does not) (BUG).
206  */
207
208 char *xbt_str_varsubst(const char *str, xbt_dict_t patterns)
209 {
210   xbt_strbuff_t buff = xbt_strbuff_new_from(str);
211   char *res;
212   xbt_strbuff_varsubst(buff, patterns);
213   res = buff->data;
214   xbt_strbuff_free_container(buff);
215   return res;
216 }
217
218
219 /** @brief Splits a string into a dynar of strings
220  *
221  * @param s: the string to split
222  * @param sep: a string of all chars to consider as separator.
223  *
224  * By default (with sep=NULL), these characters are used as separator:
225  *
226  *      - " "           (ASCII 32       (0x20)) space.
227  *      - "\t"          (ASCII 9        (0x09)) tab.
228  *      - "\n"          (ASCII 10       (0x0A)) line feed.
229  *      - "\r"          (ASCII 13       (0x0D)) carriage return.
230  *      - "\0"          (ASCII 0        (0x00)) NULL.
231  *      - "\x0B"        (ASCII 11       (0x0B)) vertical tab.
232  */
233
234 xbt_dynar_t xbt_str_split(const char *s, const char *sep)
235 {
236   xbt_dynar_t res = xbt_dynar_new(sizeof(char *), &xbt_free_ref);
237   const char *p, *q;
238   int done;
239   const char *sep_dflt = " \t\n\r\x0B";
240   char is_sep[256] = { 1, 0 };
241
242   /* check what are the separators */
243   memset(is_sep, 0, sizeof(is_sep));
244   if (!sep) {
245     while (*sep_dflt)
246       is_sep[(unsigned char) *sep_dflt++] = 1;
247   } else {
248     while (*sep)
249       is_sep[(unsigned char) *sep++] = 1;
250   }
251   is_sep[0] = 1;                /* End of string is also separator */
252
253   /* Do the job */
254   p = q = s;
255   done = 0;
256
257   if (s[0] == '\0')
258     return res;
259
260   while (!done) {
261     char *topush;
262     while (!is_sep[(unsigned char) *q]) {
263       q++;
264     }
265     if (*q == '\0')
266       done = 1;
267
268     topush = xbt_malloc(q - p + 1);
269     memcpy(topush, p, q - p);
270     topush[q - p] = '\0';
271     xbt_dynar_push(res, &topush);
272     p = ++q;
273   }
274
275   return res;
276 }
277
278 /**
279  * \brief This functions splits a string after using another string as separator
280  * For example A!!B!!C splitted after !! will return the dynar {A,B,C}
281  * \return An array of dynars containing the string tokens
282  */
283 xbt_dynar_t xbt_str_split_str(const char *s, const char *sep)
284 {
285   xbt_dynar_t res = xbt_dynar_new(sizeof(char *), &xbt_free_ref);
286   int done;
287   const char *p, *q;
288
289   p = q = s;
290   done = 0;
291
292   if (s[0] == '\0')
293     return res;
294   if (sep[0] == '\0') {
295     s = xbt_strdup(s);
296     xbt_dynar_push(res, &s);
297     return res;
298   }
299
300   while (!done) {
301     char *to_push;
302     int v = 0;
303     //get the start of the first occurence of the substring
304     q = strstr(p, sep);
305     //if substring was not found add the entire string
306     if (NULL == q) {
307       v = strlen(p);
308       to_push = malloc(v + 1);
309       memcpy(to_push, p, v);
310       to_push[v] = '\0';
311       xbt_dynar_push(res, &to_push);
312       done = 1;
313     } else {
314       //get the appearance
315       to_push = malloc(q - p + 1);
316       memcpy(to_push, p, q - p);
317       //add string terminator
318       to_push[q - p] = '\0';
319       xbt_dynar_push(res, &to_push);
320       p = q + strlen(sep);
321     }
322   }
323   return res;
324 }
325
326 /** @brief Just like @ref xbt_str_split_quoted (Splits a string into a dynar of strings), but without memory allocation
327  *
328  * The string passed as argument must be writable (not const)
329  * The elements of the dynar are just parts of the string passed as argument.
330  *
331  * To free the structure constructed by this function, free the first element and free the dynar:
332  *
333  * free(xbt_dynar_get_ptr(dynar,0));
334  * xbt_dynar_free(&dynar);
335  */
336 xbt_dynar_t xbt_str_split_quoted_in_place(char *s) {
337   xbt_dynar_t res = xbt_dynar_new(sizeof(char *), NULL);
338   char *beg, *end;              /* pointers around the parsed chunk */
339   int in_simple_quote = 0, in_double_quote = 0;
340   int done = 0;
341   int ctn = 0;                  /* Got something in this block */
342
343   if (s[0] == '\0')
344     return res;
345
346   beg = s;
347
348   /* do not trim leading spaces: caller responsability to clean his cruft */
349   end = beg;
350
351   while (!done) {
352
353
354     switch (*end) {
355     case '\\':
356       ctn = 1;
357       /* Protected char; move it closer */
358       memmove(end, end + 1, strlen(end));
359       if (*end == '\0')
360         THROWF(arg_error, 0, "String ends with \\");
361       end++;                    /* Pass the protected char */
362       break;
363
364     case '\'':
365       ctn = 1;
366       if (!in_double_quote) {
367         in_simple_quote = !in_simple_quote;
368         memmove(end, end + 1, strlen(end));
369       } else {
370         /* simple quote protected by double ones */
371         end++;
372       }
373       break;
374     case '"':
375       ctn = 1;
376       if (!in_simple_quote) {
377         in_double_quote = !in_double_quote;
378         memmove(end, end + 1, strlen(end));
379       } else {
380         /* double quote protected by simple ones */
381         end++;
382       }
383       break;
384
385     case ' ':
386     case '\t':
387     case '\n':
388     case '\0':
389       if (*end == '\0' && (in_simple_quote || in_double_quote)) {
390         THROWF(arg_error, 0,
391                "End of string found while searching for %c in %s",
392                (in_simple_quote ? '\'' : '"'), s);
393       }
394       if (in_simple_quote || in_double_quote) {
395         end++;
396       } else {
397         if (*end == '\0')
398           done = 1;
399
400         *end = '\0';
401         if (ctn) {
402           /* Found a separator. Push the string if contains something */
403           xbt_dynar_push(res, &beg);
404         }
405         ctn = 0;
406
407         if (done)
408           break;
409
410         beg = ++end;
411         /* trim within the string, manually to speed things up */
412         while (*beg == ' ')
413           beg++;
414         end = beg;
415       }
416       break;
417
418     default:
419       ctn = 1;
420       end++;
421     }
422   }
423   return res;
424 }
425
426 /** @brief Splits a string into a dynar of strings, taking quotes into account
427  *
428  * It basically does the same argument separation than the shell, where white
429  * spaces can be escaped and where arguments are never split within a
430  * quote group.
431  * Several subsequent spaces are ignored (unless within quotes, of course).
432  * You may want to trim the input string, if you want to avoid empty entries
433  *
434  */
435
436 xbt_dynar_t xbt_str_split_quoted(const char *s)
437 {
438   xbt_dynar_t res = xbt_dynar_new(sizeof(char *), &xbt_free_ref);
439   xbt_dynar_t parsed;
440   char *str_to_free;            /* we have to copy the string before, to handle backslashes */
441   unsigned int cursor;
442   char *p;
443
444   if (s[0] == '\0')
445     return res;
446   str_to_free = xbt_strdup(s);
447
448   parsed = xbt_str_split_quoted_in_place(str_to_free);
449   xbt_dynar_foreach(parsed,cursor,p) {
450     char *q=xbt_strdup(p);
451     xbt_dynar_push(res,&q);
452   }
453   free(str_to_free);
454   xbt_dynar_shrink(res, 0);
455   xbt_dynar_free(&parsed);
456   return res;
457 }
458
459 /** @brief Join a set of strings as a single string */
460 char *xbt_str_join(xbt_dynar_t dyn, const char *sep)
461 {
462   int len = 1, dyn_len = xbt_dynar_length(dyn);
463   unsigned int cpt;
464   char *cursor;
465   char *res, *p;
466
467   if (!dyn_len)
468     return xbt_strdup("");
469
470   /* compute the length */
471   xbt_dynar_foreach(dyn, cpt, cursor) {
472     len += strlen(cursor);
473   }
474   len += strlen(sep) * dyn_len;
475   /* Do the job */
476   res = xbt_malloc(len);
477   p = res;
478   xbt_dynar_foreach(dyn, cpt, cursor) {
479     if ((int) cpt < dyn_len - 1)
480       p += sprintf(p, "%s%s", cursor, sep);
481     else
482       p += sprintf(p, "%s", cursor);
483   }
484   return res;
485 }
486 /** @brief Join a set of strings as a single string
487  *
488  * The parameter must be a NULL-terminated array of chars,
489  * just like xbt_dynar_to_array() produces
490  */
491 char *xbt_str_join_array(const char *const *strs, const char *sep)
492 {
493   char *res,*q;
494   int amount_strings=0;
495   int len=0;
496   int i;
497
498   if ((!strs) || (!strs[0]))
499     return xbt_strdup("");
500
501   /* compute the length before malloc */
502   for (i=0;strs[i];i++) {
503     len += strlen(strs[i]);
504     amount_strings++;
505   }
506   len += strlen(sep) * amount_strings;
507
508   /* Do the job */
509   q = res = xbt_malloc(len);
510   for (i=0;strs[i];i++) {
511     if (i!=0) { // not first loop
512       q += sprintf(q, "%s%s", sep, strs[i]);
513     } else {
514       q += sprintf(q,"%s",strs[i]);
515     }
516   }
517   return res;
518 }
519
520 #if defined(SIMGRID_NEED_GETLINE) || defined(DOXYGEN)
521 /** @brief Get a single line from the stream (reimplementation of the GNU getline)
522  *
523  * This is a redefinition of the GNU getline function, used on platforms where it does not exists.
524  *
525  * getline() reads an entire line from stream, storing the address of the buffer
526  * containing the text into *buf.  The buffer is null-terminated and includes
527  * the newline character, if one was found.
528  *
529  * If *buf is NULL, then getline() will allocate a buffer for storing the line,
530  * which should be freed by the user program.  Alternatively, before calling getline(),
531  * *buf can contain a pointer to a malloc()-allocated buffer *n bytes in size.  If the buffer
532  * is not large enough to hold the line, getline() resizes it with realloc(), updating *buf and *n
533  * as necessary.  In either case, on a successful call, *buf and *n will be updated to
534  * reflect the buffer address and allocated size respectively.
535  */
536 long getline(char **buf, size_t * n, FILE * stream)
537 {
538
539   size_t i;
540   int ch;
541
542   if (!*buf) {
543     *buf = xbt_malloc(512);
544     *n = 512;
545   }
546
547   if (feof(stream))
548     return (ssize_t) - 1;
549
550   for (i = 0; (ch = fgetc(stream)) != EOF; i++) {
551
552     if (i >= (*n) + 1)
553       *buf = xbt_realloc(*buf, *n += 512);
554
555     (*buf)[i] = ch;
556
557     if ((*buf)[i] == '\n') {
558       i++;
559       (*buf)[i] = '\0';
560       break;
561     }
562   }
563
564   if (i == *n)
565     *buf = xbt_realloc(*buf, *n += 1);
566
567   (*buf)[i] = '\0';
568
569   return (ssize_t) i;
570 }
571
572 #endif                          /* HAVE_GETLINE */
573
574 /*
575  * Diff related functions
576  *
577  * Implementation of the algorithm described in "An O(NP) Sequence
578  * Comparison Algorithm", by Sun Wu, Udi Manber, Gene Myers, and Webb
579  * Miller.
580  */
581
582 struct subsequence {
583     int x, y;                   /* starting coordinates */
584     int len;                    /* length */
585 };
586
587 static XBT_INLINE
588 void diff_snake(const char *vec_a[], int a0, int len_a,
589                 const char *vec_b[], int b0, int len_b,
590                 struct subsequence *seqs, int *fp, int k, int limit)
591 {
592   int record_seq;
593   int x, y;
594   int fp_left = fp[k - 1] + 1;
595   int fp_right = fp[k + 1];
596   if (fp_left > fp_right) {
597     x = fp_left;
598     record_seq = k - 1;
599   } else {
600     x = fp_right;
601     record_seq = k + 1;
602   }
603   y = x - k;
604   if (x + y <= limit) {
605     seqs[k].x = x;
606     seqs[k].y = y;
607     record_seq = k;
608   } else {
609     seqs[k] = seqs[record_seq];
610   }
611   while (x < len_a && y < len_b && !strcmp(vec_a[a0 + x], vec_b[b0 + y]))
612     ++x, ++y;
613   fp[k] = x;
614   if (record_seq == k)
615     seqs[k].len = x - seqs[k].x;
616 }
617
618 /* Returns the length of a shortest edit script, and a common
619  * subsequence from the middle.
620  */
621 static int diff_middle_subsequence(const char *vec_a[], int a0,  int len_a,
622                                    const char *vec_b[], int b0,  int len_b,
623                                    struct subsequence *subseq,
624                                    struct subsequence *seqs, int *fp)
625 {
626   const int delta = len_a - len_b;
627   const int limit = (len_a + len_b) / 2;
628   int kmin;
629   int kmax;
630   int k;
631   int p = -1;
632
633   if (delta >= 0) {
634     kmin = 0;
635     kmax = delta;
636   } else {
637     kmin = delta;
638     kmax = 0;
639   }
640   for (k = kmin; k <= kmax; k++)
641     fp[k] = -1;
642   do {
643     p++;
644     fp[kmin - p - 1] = fp[kmax + p + 1] = -1;
645     for (k = kmax + p; k > delta; k--)
646       diff_snake(vec_a, a0, len_a, vec_b, b0, len_b, seqs, fp, k, limit);
647     for (k = kmin - p; k <= delta; k++)
648       diff_snake(vec_a, a0, len_a, vec_b, b0, len_b, seqs, fp, k, limit);
649   } while (fp[delta] != len_a);
650
651   subseq->x = a0 + seqs[delta].x;
652   subseq->y = b0 + seqs[delta].y;
653   subseq->len = seqs[delta].len;
654   return abs(delta) + 2 * p;;
655 }
656
657 /* Finds a longest common subsequence.
658  */
659 static void diff_compute_lcs(const char *vec_a[], int a0, int len_a,
660                              const char *vec_b[], int b0, int len_b,
661                              xbt_dynar_t common_sequence,
662                              struct subsequence *seqs, int *fp)
663 {
664   if (len_a > 0 && len_b > 0) {
665     struct subsequence subseq;
666     int ses_len = diff_middle_subsequence(vec_a, a0, len_a, vec_b, b0, len_b,
667                                           &subseq, seqs, fp);
668     if (ses_len > 1) {
669       int u = subseq.x + subseq.len;
670       int v = subseq.y + subseq.len;
671       diff_compute_lcs(vec_a, a0, subseq.x - a0, vec_b, b0, subseq.y - b0,
672                        common_sequence, seqs, fp);
673       if (subseq.len > 0)
674         xbt_dynar_push(common_sequence, &subseq);
675       diff_compute_lcs(vec_a, u, a0 + len_a - u, vec_b, v, b0 + len_b - v,
676                        common_sequence, seqs, fp);
677     } else {
678       int len = MIN(len_a, len_b) - subseq.len;
679       if (subseq.x == a0 && subseq.y == b0) {
680         if (subseq.len > 0)
681           xbt_dynar_push(common_sequence, &subseq);
682         if (len > 0) {
683           struct subsequence subseq0 = {a0 + len_a - len,
684                                         b0 + len_b - len, len};
685           xbt_dynar_push(common_sequence, &subseq0);
686         }
687       } else {
688         if (len > 0) {
689           struct subsequence subseq0 = {a0, b0, len};
690           xbt_dynar_push(common_sequence, &subseq0);
691         }
692         if (subseq.len > 0)
693           xbt_dynar_push(common_sequence, &subseq);
694       }
695     }
696   }
697 }
698
699 static int diff_member(const char *s, const char *vec[], int from, int to)
700 {
701   for ( ; from < to ; from++)
702     if (!strcmp(s, vec[from]))
703       return 1;
704   return 0;
705 }
706
707 /* Extract common prefix.
708  */
709 static void diff_easy_prefix(const char *vec_a[], int *a0_p, int *len_a_p,
710                              const char *vec_b[], int *b0_p, int *len_b_p,
711                              xbt_dynar_t common_sequence)
712 {
713   int a0 = *a0_p;
714   int b0 = *b0_p;
715   int len_a = *len_a_p;
716   int len_b = *len_b_p;
717
718   while (len_a > 0 && len_b > 0) {
719     struct subsequence subseq = {a0, b0, 0};
720     while (len_a > 0 && len_b > 0 && !strcmp(vec_a[a0], vec_b[b0])) {
721       a0++, len_a--;
722       b0++, len_b--;
723       subseq.len++;
724     }
725     if (subseq.len > 0)
726       xbt_dynar_push(common_sequence, &subseq);
727     if (len_a > 0 && len_b > 0 &&
728         !diff_member(vec_a[a0], vec_b, b0 + 1, b0 + len_b)) {
729       a0++, len_a--;
730     } else {
731       break;
732     }
733   }
734
735   *a0_p = a0;
736   *b0_p = b0;
737   *len_a_p = len_a;
738   *len_b_p = len_b;
739 }
740
741 /* Extract common suffix.
742  */
743 static void diff_easy_suffix(const char *vec_a[], int *a0_p, int *len_a_p,
744                              const char *vec_b[], int *b0_p, int *len_b_p,
745                              xbt_dynar_t common_suffix)
746 {
747   int a0 = *a0_p;
748   int b0 = *b0_p;
749   int len_a = *len_a_p;
750   int len_b = *len_b_p;
751
752   while (len_a > 0 && len_b > 0){
753     struct subsequence subseq;
754     subseq.len = 0;
755     while (len_a > 0 && len_b > 0 &&
756            !strcmp(vec_a[a0 + len_a - 1], vec_b[b0 + len_b - 1])) {
757       len_a--;
758       len_b--;
759       subseq.len++;
760     }
761     if (subseq.len > 0) {
762       subseq.x = a0 + len_a;
763       subseq.y = b0 + len_b;
764       xbt_dynar_push(common_suffix, &subseq);
765     }
766     if (len_a > 0 && len_b > 0 &&
767         !diff_member(vec_b[b0 + len_b - 1], vec_a, a0, a0 + len_a - 1)) {
768       len_b--;
769     } else {
770       break;
771     }
772   }
773
774   *a0_p = a0;
775   *b0_p = b0;
776   *len_a_p = len_a;
777   *len_b_p = len_b;
778 }
779
780 /** @brief Compute the unified diff of two strings */
781 char *xbt_str_diff(const char *a, const char *b)
782 {
783   xbt_dynar_t da = xbt_str_split(a, "\n");
784   xbt_dynar_t db = xbt_str_split(b, "\n");
785   xbt_dynar_t common_sequence, common_suffix;
786   size_t len;
787   const char **vec_a, **vec_b;
788   int a0, b0;
789   int len_a, len_b;
790   int max;
791   int *fp_base, *fp;
792   struct subsequence *seqs_base, *seqs;
793   struct subsequence subseq;
794   xbt_dynar_t diff;
795   char *res;
796   int x, y;
797   unsigned s;
798
799   /* Clean empty lines at the end of da and db */
800   len = strlen(a);
801   if (len > 0 && a[len - 1] == '\n')
802     xbt_dynar_pop(da, NULL);
803   len = strlen(b);
804   if (len > 0 && b[len - 1] == '\n')
805     xbt_dynar_pop(db, NULL);
806
807   /* Various initializations */
808   /* Assume that dynar's content is contiguous */
809   a0 = 0;
810   len_a = xbt_dynar_length(da);
811   vec_a = len_a ? xbt_dynar_get_ptr(da, 0) : NULL;
812   b0 = 0;
813   len_b = xbt_dynar_length(db);
814   vec_b = len_b ? xbt_dynar_get_ptr(db, 0) : NULL;
815   max = MAX(len_a, len_b) + 1;
816   fp_base = xbt_new(int, 2 * max + 1);
817   fp = fp_base + max;           /* indexes in [-max..max] */
818   seqs_base = xbt_new(struct subsequence, 2 * max + 1);
819   seqs = seqs_base + max;       /* indexes in [-max..max] */
820   common_sequence = xbt_dynar_new(sizeof(struct subsequence), NULL);
821   common_suffix = xbt_dynar_new(sizeof(struct subsequence), NULL);
822
823   /* Add a sentinel a the end of the sequence */
824   subseq.x = len_a;
825   subseq.y = len_b;
826   subseq.len = 0;
827   xbt_dynar_push(common_suffix, &subseq);
828
829   /* Compute the Longest Common Subsequence */
830   diff_easy_prefix(vec_a, &a0, &len_a, vec_b, &b0, &len_b, common_sequence);
831   diff_easy_suffix(vec_a, &a0, &len_a, vec_b, &b0, &len_b, common_suffix);
832   diff_compute_lcs(vec_a, a0, len_a, vec_b, b0, len_b, common_sequence, seqs, fp);
833   while (!xbt_dynar_is_empty(common_suffix)) {
834     xbt_dynar_pop(common_suffix, &subseq);
835     xbt_dynar_push(common_sequence, &subseq);
836   }
837
838   /* Build a Shortest Edit Script, and the final result */
839   diff = xbt_dynar_new(sizeof(char *), &xbt_free_ref);
840   x = 0;
841   y = 0;
842   xbt_dynar_foreach(common_sequence, s, subseq) {
843     while (x < subseq.x) {
844       char *topush = bprintf("- %s", vec_a[x++]);
845       xbt_dynar_push_as(diff, char*, topush);
846     }
847     while (y < subseq.y) {
848       char *topush = bprintf("+ %s", vec_b[y++]);
849       xbt_dynar_push_as(diff, char*, topush);
850     }
851     while (x < subseq.x + subseq.len) {
852       char *topush = bprintf("  %s", vec_a[x++]);
853       xbt_dynar_push_as(diff, char*, topush);
854       y++;
855     }
856   }
857   res = xbt_str_join(diff, "\n");
858
859   xbt_free(fp_base);
860   xbt_free(seqs_base);
861   xbt_dynar_free(&db);
862   xbt_dynar_free(&da);
863   xbt_dynar_free(&common_sequence);
864   xbt_dynar_free(&common_suffix);
865   xbt_dynar_free(&diff);
866
867   return res;
868 }
869
870
871 /** @brief creates a new string containing what can be read on a fd
872  *
873  */
874 char *xbt_str_from_file(FILE * file)
875 {
876   xbt_strbuff_t buff = xbt_strbuff_new();
877   char *res;
878   char bread[1024];
879   memset(bread, 0, 1024);
880
881   while (!feof(file)) {
882     int got = fread(bread, 1, 1023, file);
883     bread[got] = '\0';
884     xbt_strbuff_append(buff, bread);
885   }
886
887   res = buff->data;
888   xbt_strbuff_free_container(buff);
889   return res;
890 }
891
892 #ifdef SIMGRID_TEST
893 #include "xbt/str.h"
894
895 #define mytest(name, input, expected) \
896   xbt_test_add(name); \
897   d=xbt_str_split_quoted(input); \
898   s=xbt_str_join(d,"XXX"); \
899   xbt_test_assert(!strcmp(s,expected),\
900                    "Input (%s) leads to (%s) instead of (%s)", \
901                    input,s,expected);\
902                    free(s); \
903                    xbt_dynar_free(&d);
904
905 XBT_TEST_SUITE("xbt_str", "String Handling");
906 XBT_TEST_UNIT("xbt_str_split_quoted", test_split_quoted, "test the function xbt_str_split_quoted")
907 {
908   xbt_dynar_t d;
909   char *s;
910
911   mytest("Empty", "", "");
912   mytest("Basic test", "toto tutu", "totoXXXtutu");
913   mytest("Useless backslashes", "\\t\\o\\t\\o \\t\\u\\t\\u",
914          "totoXXXtutu");
915   mytest("Protected space", "toto\\ tutu", "toto tutu");
916   mytest("Several spaces", "toto   tutu", "totoXXXtutu");
917   mytest("LTriming", "  toto tatu", "totoXXXtatu");
918   mytest("Triming", "  toto   tutu  ", "totoXXXtutu");
919   mytest("Single quotes", "'toto tutu' tata", "toto tutuXXXtata");
920   mytest("Double quotes", "\"toto tutu\" tata", "toto tutuXXXtata");
921   mytest("Mixed quotes", "\"toto' 'tutu\" tata", "toto' 'tutuXXXtata");
922   mytest("Backslashed quotes", "\\'toto tutu\\' tata",
923          "'totoXXXtutu'XXXtata");
924   mytest("Backslashed quotes + quotes", "'toto \\'tutu' tata",
925          "toto 'tutuXXXtata");
926
927 }
928
929 #define mytest_str(name, input, separator, expected) \
930   xbt_test_add(name); \
931   d=xbt_str_split_str(input, separator); \
932   s=xbt_str_join(d,"XXX"); \
933   xbt_test_assert(!strcmp(s,expected),\
934                    "Input (%s) leads to (%s) instead of (%s)", \
935                    input,s,expected);\
936                    free(s); \
937                    xbt_dynar_free(&d);
938
939 XBT_TEST_UNIT("xbt_str_split_str", test_split_str, "test the function xbt_str_split_str")
940 {
941   xbt_dynar_t d;
942   char *s;
943
944   mytest_str("Empty string and separator", "", "", "");
945   mytest_str("Empty string", "", "##", "");
946   mytest_str("Empty separator", "toto", "", "toto");
947   mytest_str("String with no separator in it", "toto", "##", "toto");
948   mytest_str("Basic test", "toto##tutu", "##", "totoXXXtutu");
949 }
950
951 /* Last args are format string and parameters for xbt_test_add */
952 #define mytest_diff(s1, s2, diff, ...)                                  \
953   do {                                                                  \
954     char *mytest_diff_res;                                              \
955     xbt_test_add(__VA_ARGS__);                                          \
956     mytest_diff_res = xbt_str_diff(s1, s2);                             \
957     xbt_test_assert(!strcmp(mytest_diff_res, diff),                     \
958                     "Wrong output:\n--- got:\n%s\n--- expected:\n%s\n---", \
959                     mytest_diff_res, diff);                             \
960     free(mytest_diff_res);                                              \
961   } while (0)
962
963 XBT_TEST_UNIT("xbt_str_diff", test_diff, "test the function xbt_str_diff")
964 {
965   unsigned i;
966
967   /* Trivial cases */
968   mytest_diff("a", "a", "  a", "1 word, no difference");
969   mytest_diff("a", "A", "- a\n+ A", "1 word, different");
970   mytest_diff("a\n", "a\n", "  a", "1 line, no difference");
971   mytest_diff("a\n", "A\n", "- a\n+ A", "1 line, different");
972
973   /* Empty strings */
974   mytest_diff("", "", "", "empty strings");
975   mytest_diff("", "a", "+ a", "1 word, added");
976   mytest_diff("a", "", "- a", "1 word, removed");
977   mytest_diff("", "a\n", "+ a", "1 line, added");
978   mytest_diff("a\n", "", "- a", "1 line, removed");
979   mytest_diff("", "a\nb\nc\n", "+ a\n+ b\n+ c", "4 lines, all added");
980   mytest_diff("a\nb\nc\n", "", "- a\n- b\n- c", "4 lines, all removed");
981
982   /* Empty lines */
983   mytest_diff("\n", "\n", "  ", "empty lines");
984   mytest_diff("", "\n", "+ ", "empty line, added");
985   mytest_diff("\n", "", "- ", "empty line, removed");
986
987   mytest_diff("a", "\na", "+ \n  a", "empty line added before word");
988   mytest_diff("a", "a\n\n", "  a\n+ ", "empty line added after word");
989   mytest_diff("\na", "a", "- \n  a", "empty line removed before word");
990   mytest_diff("a\n\n", "a", "  a\n- ", "empty line removed after word");
991
992   mytest_diff("a\n", "\na\n", "+ \n  a", "empty line added before line");
993   mytest_diff("a\n", "a\n\n", "  a\n+ ", "empty line added after line");
994   mytest_diff("\na\n", "a\n", "- \n  a", "empty line removed before line");
995   mytest_diff("a\n\n", "a\n", "  a\n- ", "empty line removed after line");
996
997   mytest_diff("a\nb\nc\nd\n", "\na\nb\nc\nd\n", "+ \n  a\n  b\n  c\n  d",
998               "empty line added before 4 lines");
999   mytest_diff("a\nb\nc\nd\n", "a\nb\nc\nd\n\n", "  a\n  b\n  c\n  d\n+ ",
1000               "empty line added after 4 lines");
1001   mytest_diff("\na\nb\nc\nd\n", "a\nb\nc\nd\n", "- \n  a\n  b\n  c\n  d",
1002               "empty line removed before 4 lines");
1003   mytest_diff("a\nb\nc\nd\n\n", "a\nb\nc\nd\n", "  a\n  b\n  c\n  d\n- ",
1004               "empty line removed after 4 lines");
1005
1006   /* Missing newline at the end of one of the strings */
1007   mytest_diff("a\n", "a", "  a", "1 line, 1 word, no difference");
1008   mytest_diff("a", "a\n", "  a", "1 word, 1 line, no difference");
1009   mytest_diff("a\n", "A", "- a\n+ A", "1 line, 1 word, different");
1010   mytest_diff("a", "A\n", "- a\n+ A", "1 word, 1 line, different");
1011
1012   mytest_diff("a\nb\nc\nd", "a\nb\nc\nd\n", "  a\n  b\n  c\n  d",
1013               "4 lines, no newline on first");
1014   mytest_diff("a\nb\nc\nd\n", "a\nb\nc\nd", "  a\n  b\n  c\n  d",
1015               "4 lines, no newline on second");
1016
1017   /* Four lines, all combinations of differences */
1018   for (i = 0 ; i < (1U << 4) ; i++) {
1019     char d2[4 + 1];
1020     char s2[4 * 2 + 1];
1021     char res[4 * 8 + 1];
1022     char *pd = d2;
1023     char *ps = s2;
1024     char *pr = res;
1025     unsigned j = 0;
1026     while (j < 4) {
1027       unsigned k;
1028       for (/* j */ ; j < 4 && !(i & (1U << j)) ; j++) {
1029         *pd++ = "abcd"[j];
1030         ps += sprintf(ps, "%c\n", "abcd"[j]);
1031         pr += sprintf(pr, "  %c\n", "abcd"[j]);
1032       }
1033       for (k = j ; k < 4 && (i & (1U << k)) ; k++) {
1034         *pd++ = "ABCD"[k];
1035         ps += sprintf(ps, "%c\n", "ABCD"[k]);
1036         pr += sprintf(pr, "- %c\n", "abcd"[k]);
1037       }
1038       for (/* j */ ; j < k ; j++) {
1039         pr += sprintf(pr, "+ %c\n", "ABCD"[j]);
1040       }
1041     }
1042     *pd = '\0';
1043     *--pr = '\0';               /* strip last '\n' from expected result */
1044     mytest_diff("a\nb\nc\nd\n", s2, res,
1045                 "compare (abcd) with changed (%s)", d2);
1046   }
1047
1048   /* Subsets of four lines, do not test for empty subset */
1049   for (i = 1 ; i < (1U << 4) ; i++) {
1050     char d2[4 + 1];
1051     char s2[4 * 2 + 1];
1052     char res[4 * 8 + 1];
1053     char *pd = d2;
1054     char *ps = s2;
1055     char *pr = res;
1056     unsigned j = 0;
1057     while (j < 4) {
1058       for (/* j */ ; j < 4 && (i & (1U << j)) ; j++) {
1059         *pd++ = "abcd"[j];
1060         ps += sprintf(ps, "%c\n", "abcd"[j]);
1061         pr += sprintf(pr, "  %c\n", "abcd"[j]);
1062       }
1063       for (/* j */; j < 4 && !(i & (1U << j)) ; j++) {
1064         pr += sprintf(pr, "- %c\n", "abcd"[j]);
1065       }
1066     }
1067     *pd = '\0';
1068     *--pr = '\0';               /* strip last '\n' from expected result */
1069     mytest_diff("a\nb\nc\nd\n", s2, res,
1070                 "compare (abcd) with subset (%s)", d2);
1071
1072     for (pr = res ; *pr != '\0' ; pr++)
1073       if (*pr == '-')
1074         *pr = '+';
1075     mytest_diff(s2, "a\nb\nc\nd\n", res,
1076                 "compare subset (%s) with (abcd)", d2);
1077   }
1078 }
1079
1080 #endif                          /* SIMGRID_TEST */