Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
[xbt] Add unit test got xbt_dict with arbitrary (not a NULL temrinated string) key
[simgrid.git] / src / xbt / xbt_str.c
1 /* xbt_str.c - various helping functions to deal with strings               */
2
3 /* Copyright (c) 2007-2014. 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 = xbt_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 = xbt_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  * So if you don't store that argument elsewhere, you should free it in addition
331  * to freeing the dynar. This can be done by simply freeing the first argument
332  * of the dynar:
333  *  free(xbt_dynar_get_ptr(dynar,0));
334  *
335  * Actually this function puts a bunch of \0 in the memory area you passed as
336  * argument to separate the elements, and pushes the address of each chunk
337  * in the resulting dynar. Yes, that's uneven. Yes, that's gory. But that's efficient.
338  */
339 xbt_dynar_t xbt_str_split_quoted_in_place(char *s) {
340   xbt_dynar_t res = xbt_dynar_new(sizeof(char *), NULL);
341   char *beg, *end;              /* pointers around the parsed chunk */
342   int in_simple_quote = 0, in_double_quote = 0;
343   int done = 0;
344   int ctn = 0;                  /* Got something in this block */
345
346   if (s[0] == '\0')
347     return res;
348
349   beg = s;
350
351   /* do not trim leading spaces: caller responsibility to clean his cruft */
352   end = beg;
353
354   while (!done) {
355
356
357     switch (*end) {
358     case '\\':
359       ctn = 1;
360       /* Protected char; move it closer */
361       memmove(end, end + 1, strlen(end));
362       if (*end == '\0')
363         THROWF(arg_error, 0, "String ends with \\");
364       end++;                    /* Pass the protected char */
365       break;
366
367     case '\'':
368       ctn = 1;
369       if (!in_double_quote) {
370         in_simple_quote = !in_simple_quote;
371         memmove(end, end + 1, strlen(end));
372       } else {
373         /* simple quote protected by double ones */
374         end++;
375       }
376       break;
377     case '"':
378       ctn = 1;
379       if (!in_simple_quote) {
380         in_double_quote = !in_double_quote;
381         memmove(end, end + 1, strlen(end));
382       } else {
383         /* double quote protected by simple ones */
384         end++;
385       }
386       break;
387
388     case ' ':
389     case '\t':
390     case '\n':
391     case '\0':
392       if (*end == '\0' && (in_simple_quote || in_double_quote)) {
393         THROWF(arg_error, 0,
394                "End of string found while searching for %c in %s",
395                (in_simple_quote ? '\'' : '"'), s);
396       }
397       if (in_simple_quote || in_double_quote) {
398         end++;
399       } else {
400         if (*end == '\0')
401           done = 1;
402
403         *end = '\0';
404         if (ctn) {
405           /* Found a separator. Push the string if contains something */
406           xbt_dynar_push(res, &beg);
407         }
408         ctn = 0;
409
410         if (done)
411           break;
412
413         beg = ++end;
414         /* trim within the string, manually to speed things up */
415         while (*beg == ' ')
416           beg++;
417         end = beg;
418       }
419       break;
420
421     default:
422       ctn = 1;
423       end++;
424     }
425   }
426   return res;
427 }
428
429 /** @brief Splits a string into a dynar of strings, taking quotes into account
430  *
431  * It basically does the same argument separation than the shell, where white
432  * spaces can be escaped and where arguments are never split within a
433  * quote group.
434  * Several subsequent spaces are ignored (unless within quotes, of course).
435  * You may want to trim the input string, if you want to avoid empty entries
436  *
437  */
438
439 xbt_dynar_t xbt_str_split_quoted(const char *s)
440 {
441   xbt_dynar_t res = xbt_dynar_new(sizeof(char *), &xbt_free_ref);
442   xbt_dynar_t parsed;
443   char *str_to_free;            /* we have to copy the string before, to handle backslashes */
444   unsigned int cursor;
445   char *p;
446
447   if (s[0] == '\0')
448     return res;
449   str_to_free = xbt_strdup(s);
450
451   parsed = xbt_str_split_quoted_in_place(str_to_free);
452   xbt_dynar_foreach(parsed,cursor,p) {
453     char *q=xbt_strdup(p);
454     xbt_dynar_push(res,&q);
455   }
456   free(str_to_free);
457   xbt_dynar_shrink(res, 0);
458   xbt_dynar_free(&parsed);
459   return res;
460 }
461
462 /** @brief Join a set of strings as a single string */
463 char *xbt_str_join(xbt_dynar_t dyn, const char *sep)
464 {
465   int len = 1, dyn_len = xbt_dynar_length(dyn);
466   unsigned int cpt;
467   char *cursor;
468   char *res, *p;
469
470   if (!dyn_len)
471     return xbt_strdup("");
472
473   /* compute the length */
474   xbt_dynar_foreach(dyn, cpt, cursor) {
475     len += strlen(cursor);
476   }
477   len += strlen(sep) * dyn_len;
478   /* Do the job */
479   res = xbt_malloc(len);
480   p = res;
481   xbt_dynar_foreach(dyn, cpt, cursor) {
482     if ((int) cpt < dyn_len - 1)
483       p += sprintf(p, "%s%s", cursor, sep);
484     else
485       p += sprintf(p, "%s", cursor);
486   }
487   return res;
488 }
489 /** @brief Join a set of strings as a single string
490  *
491  * The parameter must be a NULL-terminated array of chars,
492  * just like xbt_dynar_to_array() produces
493  */
494 char *xbt_str_join_array(const char *const *strs, const char *sep)
495 {
496   char *res,*q;
497   int amount_strings=0;
498   int len=0;
499   int i;
500
501   if ((!strs) || (!strs[0]))
502     return xbt_strdup("");
503
504   /* compute the length before malloc */
505   for (i=0;strs[i];i++) {
506     len += strlen(strs[i]);
507     amount_strings++;
508   }
509   len += strlen(sep) * amount_strings;
510
511   /* Do the job */
512   q = res = xbt_malloc(len);
513   for (i=0;strs[i];i++) {
514     if (i!=0) { // not first loop
515       q += sprintf(q, "%s%s", sep, strs[i]);
516     } else {
517       q += sprintf(q,"%s",strs[i]);
518     }
519   }
520   return res;
521 }
522
523 /** @brief Get a single line from the stream (reimplementation of the GNU getline)
524  *
525  * This is a reimplementation of the GNU getline function, so that our code don't depends on the GNU libc.
526  *
527  * xbt_getline() reads an entire line from stream, storing the address of the
528  * buffer containing the text into *buf.  The buffer is null-terminated and
529  * includes the newline character, if one was found.
530  *
531  * If *buf is NULL, then xbt_getline() will allocate a buffer for storing the
532  * line, which should be freed by the user program.
533  *
534  * Alternatively, before calling xbt_getline(), *buf can contain a pointer to a
535  * malloc()-allocated buffer *n bytes in size.  If the buffer is not large
536  * enough to hold the line, xbt_getline() resizes it with realloc(), updating
537  * *buf and *n as necessary.
538  *
539  * In either case, on a successful call, *buf and *n will be updated to reflect
540  * the buffer address and allocated size respectively.
541  */
542 ssize_t xbt_getline(char **buf, size_t *n, FILE *stream)
543 {
544   ssize_t i;
545   int ch;
546
547   ch = getc(stream);
548   if (ferror(stream) || feof(stream))
549     return -1;
550
551   if (!*buf) {
552     *n = 512;
553     *buf = xbt_malloc(*n);
554   }
555
556   i = 0;
557   do {
558     if (i == *n)
559       *buf = xbt_realloc(*buf, *n += 512);
560     (*buf)[i++] = ch;
561   } while (ch != '\n' && (ch = getc(stream)) != EOF);
562
563   if (i == *n)
564     *buf = xbt_realloc(*buf, *n += 1);
565   (*buf)[i] = '\0';
566
567   return i;
568 }
569
570 /*
571  * Diff related functions
572  *
573  * Implementation of the algorithm described in "An O(NP) Sequence Comparison
574  * Algorithm", by Sun Wu, Udi Manber, Gene Myers, and Webb Miller (Information
575  * Processing Letters 35(6):317-323, 1990), with the linear-space
576  * divide-and-conquer strategy described in "An O(ND) Difference Algorithm and
577  * Its Variations", by Eugene W. Myers (Algorithmica 1:251-266, 1986).
578  */
579
580 struct subsequence {
581     int x, y;                   /* starting coordinates */
582     int len;                    /* length */
583 };
584
585 static XBT_INLINE
586 void diff_snake(const char *vec_a[], int a0, int len_a,
587                 const char *vec_b[], int b0, int len_b,
588                 struct subsequence *seqs, int *fp, int k, int limit)
589 {
590   int record_seq;
591   int x, y;
592   int fp_left = fp[k - 1] + 1;
593   int fp_right = fp[k + 1];
594   if (fp_left > fp_right) {
595     x = fp_left;
596     record_seq = k - 1;
597   } else {
598     x = fp_right;
599     record_seq = k + 1;
600   }
601   y = x - k;
602   if (x + y <= limit) {
603     seqs[k].x = x;
604     seqs[k].y = y;
605     record_seq = k;
606   } else {
607     seqs[k] = seqs[record_seq];
608   }
609   while (x < len_a && y < len_b && !strcmp(vec_a[a0 + x], vec_b[b0 + y]))
610     ++x, ++y;
611   fp[k] = x;
612   if (record_seq == k)
613     seqs[k].len = x - seqs[k].x;
614 }
615
616 /* Returns the length of a shortest edit script, and a common
617  * subsequence from the middle.
618  */
619 static int diff_middle_subsequence(const char *vec_a[], int a0,  int len_a,
620                                    const char *vec_b[], int b0,  int len_b,
621                                    struct subsequence *subseq,
622                                    struct subsequence *seqs, int *fp)
623 {
624   const int delta = len_a - len_b;
625   const int limit = (len_a + len_b) / 2;
626   int kmin;
627   int kmax;
628   int k;
629   int p = -1;
630
631   if (delta >= 0) {
632     kmin = 0;
633     kmax = delta;
634   } else {
635     kmin = delta;
636     kmax = 0;
637   }
638   for (k = kmin; k <= kmax; k++)
639     fp[k] = -1;
640   do {
641     p++;
642     fp[kmin - p - 1] = fp[kmax + p + 1] = -1;
643     for (k = kmax + p; k > delta; k--)
644       diff_snake(vec_a, a0, len_a, vec_b, b0, len_b, seqs, fp, k, limit);
645     for (k = kmin - p; k <= delta; k++)
646       diff_snake(vec_a, a0, len_a, vec_b, b0, len_b, seqs, fp, k, limit);
647   } while (fp[delta] != len_a);
648
649   subseq->x = a0 + seqs[delta].x;
650   subseq->y = b0 + seqs[delta].y;
651   subseq->len = seqs[delta].len;
652   return abs(delta) + 2 * p;;
653 }
654
655 /* Finds a longest common subsequence.
656  * Returns its length.
657  */
658 static int diff_compute_lcs(const char *vec_a[], int a0, int len_a,
659                             const char *vec_b[], int b0, int len_b,
660                             xbt_dynar_t common_sequence,
661                             struct subsequence *seqs, int *fp)
662 {
663   if (len_a > 0 && len_b > 0) {
664     struct subsequence subseq;
665     int ses_len = diff_middle_subsequence(vec_a, a0, len_a, vec_b, b0, len_b,
666                                           &subseq, seqs, fp);
667     int lcs_len = (len_a + len_b - ses_len) / 2;
668     if (lcs_len == 0) {
669       return 0;
670     } else if (ses_len > 1) {
671       int lcs_len1 = subseq.len;
672       if (lcs_len1 < lcs_len)
673         lcs_len1 += diff_compute_lcs(vec_a, a0, subseq.x - a0,
674                                      vec_b, b0, subseq.y - b0,
675                                      common_sequence, seqs, fp);
676       if (subseq.len > 0)
677         xbt_dynar_push(common_sequence, &subseq);
678       if (lcs_len1 < lcs_len) {
679         int u = subseq.x + subseq.len;
680         int v = subseq.y + subseq.len;
681         diff_compute_lcs(vec_a, u, a0 + len_a - u, vec_b, v, b0 + len_b - v,
682                          common_sequence, seqs, fp);
683       }
684     } else {
685       int len = MIN(len_a, len_b) - subseq.len;
686       if (subseq.x == a0 && subseq.y == b0) {
687         if (subseq.len > 0)
688           xbt_dynar_push(common_sequence, &subseq);
689         if (len > 0) {
690           struct subsequence subseq0 = {a0 + len_a - len,
691                                         b0 + len_b - len, len};
692           xbt_dynar_push(common_sequence, &subseq0);
693         }
694       } else {
695         if (len > 0) {
696           struct subsequence subseq0 = {a0, b0, len};
697           xbt_dynar_push(common_sequence, &subseq0);
698         }
699         if (subseq.len > 0)
700           xbt_dynar_push(common_sequence, &subseq);
701       }
702     }
703     return lcs_len;
704   } else {
705     return 0;
706   }
707 }
708
709 static int diff_member(const char *s, const char *vec[], int from, int to)
710 {
711   for ( ; from < to ; from++)
712     if (!strcmp(s, vec[from]))
713       return 1;
714   return 0;
715 }
716
717 /* Extract common prefix.
718  */
719 static void diff_easy_prefix(const char *vec_a[], int *a0_p, int *len_a_p,
720                              const char *vec_b[], int *b0_p, int *len_b_p,
721                              xbt_dynar_t common_sequence)
722 {
723   int a0 = *a0_p;
724   int b0 = *b0_p;
725   int len_a = *len_a_p;
726   int len_b = *len_b_p;
727
728   while (len_a > 0 && len_b > 0) {
729     struct subsequence subseq = {a0, b0, 0};
730     while (len_a > 0 && len_b > 0 && !strcmp(vec_a[a0], vec_b[b0])) {
731       a0++, len_a--;
732       b0++, len_b--;
733       subseq.len++;
734     }
735     if (subseq.len > 0)
736       xbt_dynar_push(common_sequence, &subseq);
737     if (len_a > 0 && len_b > 0 &&
738         !diff_member(vec_a[a0], vec_b, b0 + 1, b0 + len_b)) {
739       a0++, len_a--;
740     } else {
741       break;
742     }
743   }
744
745   *a0_p = a0;
746   *b0_p = b0;
747   *len_a_p = len_a;
748   *len_b_p = len_b;
749 }
750
751 /* Extract common suffix.
752  */
753 static void diff_easy_suffix(const char *vec_a[], int *a0_p, int *len_a_p,
754                              const char *vec_b[], int *b0_p, int *len_b_p,
755                              xbt_dynar_t common_suffix)
756 {
757   int a0 = *a0_p;
758   int b0 = *b0_p;
759   int len_a = *len_a_p;
760   int len_b = *len_b_p;
761
762   while (len_a > 0 && len_b > 0){
763     struct subsequence subseq;
764     subseq.len = 0;
765     while (len_a > 0 && len_b > 0 &&
766            !strcmp(vec_a[a0 + len_a - 1], vec_b[b0 + len_b - 1])) {
767       len_a--;
768       len_b--;
769       subseq.len++;
770     }
771     if (subseq.len > 0) {
772       subseq.x = a0 + len_a;
773       subseq.y = b0 + len_b;
774       xbt_dynar_push(common_suffix, &subseq);
775     }
776     if (len_a > 0 && len_b > 0 &&
777         !diff_member(vec_b[b0 + len_b - 1], vec_a, a0, a0 + len_a - 1)) {
778       len_b--;
779     } else {
780       break;
781     }
782   }
783
784   *a0_p = a0;
785   *b0_p = b0;
786   *len_a_p = len_a;
787   *len_b_p = len_b;
788 }
789
790 /** @brief Compute the unified diff of two strings */
791 char *xbt_str_diff(const char *a, const char *b)
792 {
793   xbt_dynar_t da = xbt_str_split(a, "\n");
794   xbt_dynar_t db = xbt_str_split(b, "\n");
795   xbt_dynar_t common_sequence, common_suffix;
796   size_t len;
797   const char **vec_a, **vec_b;
798   int a0, b0;
799   int len_a, len_b;
800   int max;
801   int *fp_base, *fp;
802   struct subsequence *seqs_base, *seqs;
803   struct subsequence subseq;
804   xbt_dynar_t diff;
805   char *res;
806   int x, y;
807   unsigned s;
808
809   /* Clean empty lines at the end of da and db */
810   len = strlen(a);
811   if (len > 0 && a[len - 1] == '\n')
812     xbt_dynar_pop(da, NULL);
813   len = strlen(b);
814   if (len > 0 && b[len - 1] == '\n')
815     xbt_dynar_pop(db, NULL);
816
817   /* Various initializations */
818   /* Assume that dynar's content is contiguous */
819   a0 = 0;
820   len_a = xbt_dynar_length(da);
821   vec_a = len_a ? xbt_dynar_get_ptr(da, 0) : NULL;
822   b0 = 0;
823   len_b = xbt_dynar_length(db);
824   vec_b = len_b ? xbt_dynar_get_ptr(db, 0) : NULL;
825   max = MAX(len_a, len_b) + 1;
826   fp_base = xbt_new(int, 2 * max + 1);
827   fp = fp_base + max;           /* indexes in [-max..max] */
828   seqs_base = xbt_new(struct subsequence, 2 * max + 1);
829   seqs = seqs_base + max;       /* indexes in [-max..max] */
830   common_sequence = xbt_dynar_new(sizeof(struct subsequence), NULL);
831   common_suffix = xbt_dynar_new(sizeof(struct subsequence), NULL);
832
833   /* Add a sentinel a the end of the sequence */
834   subseq.x = len_a;
835   subseq.y = len_b;
836   subseq.len = 0;
837   xbt_dynar_push(common_suffix, &subseq);
838
839   /* Compute the Longest Common Subsequence */
840   diff_easy_prefix(vec_a, &a0, &len_a, vec_b, &b0, &len_b, common_sequence);
841   diff_easy_suffix(vec_a, &a0, &len_a, vec_b, &b0, &len_b, common_suffix);
842   diff_compute_lcs(vec_a, a0, len_a, vec_b, b0, len_b, common_sequence, seqs, fp);
843   while (!xbt_dynar_is_empty(common_suffix)) {
844     xbt_dynar_pop(common_suffix, &subseq);
845     xbt_dynar_push(common_sequence, &subseq);
846   }
847
848   /* Build a Shortest Edit Script, and the final result */
849   diff = xbt_dynar_new(sizeof(char *), &xbt_free_ref);
850   x = 0;
851   y = 0;
852   xbt_dynar_foreach(common_sequence, s, subseq) {
853     while (x < subseq.x) {
854       char *topush = bprintf("- %s", vec_a[x++]);
855       xbt_dynar_push_as(diff, char*, topush);
856     }
857     while (y < subseq.y) {
858       char *topush = bprintf("+ %s", vec_b[y++]);
859       xbt_dynar_push_as(diff, char*, topush);
860     }
861     while (x < subseq.x + subseq.len) {
862       char *topush = bprintf("  %s", vec_a[x++]);
863       xbt_dynar_push_as(diff, char*, topush);
864       y++;
865     }
866   }
867   res = xbt_str_join(diff, "\n");
868
869   xbt_free(fp_base);
870   xbt_free(seqs_base);
871   xbt_dynar_free(&db);
872   xbt_dynar_free(&da);
873   xbt_dynar_free(&common_sequence);
874   xbt_dynar_free(&common_suffix);
875   xbt_dynar_free(&diff);
876
877   return res;
878 }
879
880
881 /** @brief creates a new string containing what can be read on a fd
882  *
883  */
884 char *xbt_str_from_file(FILE * file)
885 {
886   xbt_strbuff_t buff = xbt_strbuff_new();
887   char *res;
888   char bread[1024];
889   memset(bread, 0, 1024);
890
891   while (!feof(file)) {
892     int got = fread(bread, 1, 1023, file);
893     bread[got] = '\0';
894     xbt_strbuff_append(buff, bread);
895   }
896
897   res = buff->data;
898   xbt_strbuff_free_container(buff);
899   return res;
900 }
901
902 /* @brief Retrun 1 if string 'str' starts with string 'start'
903  *
904  * \param str a string
905  * \param start the string to search in str
906  *
907  * \return 1 if 'str' starts with 'start'
908  */
909 int xbt_str_start_with(const char* str, const char* start)
910 {
911   int i;
912   size_t l_str = strlen(str);
913   size_t l_start = strlen(start);
914
915   if(l_start > l_str) return 0;
916
917   for(i = 0; i< l_start; i++){
918     if(str[i] != start[i]) return 0;
919   }
920
921   return 1;
922 }
923
924 #ifdef SIMGRID_TEST
925 #include "xbt/str.h"
926
927 #define mytest(name, input, expected) \
928   xbt_test_add(name); \
929   d=xbt_str_split_quoted(input); \
930   s=xbt_str_join(d,"XXX"); \
931   xbt_test_assert(!strcmp(s,expected),\
932                    "Input (%s) leads to (%s) instead of (%s)", \
933                    input,s,expected);\
934                    free(s); \
935                    xbt_dynar_free(&d);
936
937 XBT_TEST_SUITE("xbt_str", "String Handling");
938 XBT_TEST_UNIT("xbt_str_split_quoted", test_split_quoted, "test the function xbt_str_split_quoted")
939 {
940   xbt_dynar_t d;
941   char *s;
942
943   mytest("Empty", "", "");
944   mytest("Basic test", "toto tutu", "totoXXXtutu");
945   mytest("Useless backslashes", "\\t\\o\\t\\o \\t\\u\\t\\u",
946          "totoXXXtutu");
947   mytest("Protected space", "toto\\ tutu", "toto tutu");
948   mytest("Several spaces", "toto   tutu", "totoXXXtutu");
949   mytest("LTriming", "  toto tatu", "totoXXXtatu");
950   mytest("Triming", "  toto   tutu  ", "totoXXXtutu");
951   mytest("Single quotes", "'toto tutu' tata", "toto tutuXXXtata");
952   mytest("Double quotes", "\"toto tutu\" tata", "toto tutuXXXtata");
953   mytest("Mixed quotes", "\"toto' 'tutu\" tata", "toto' 'tutuXXXtata");
954   mytest("Backslashed quotes", "\\'toto tutu\\' tata",
955          "'totoXXXtutu'XXXtata");
956   mytest("Backslashed quotes + quotes", "'toto \\'tutu' tata",
957          "toto 'tutuXXXtata");
958
959 }
960
961 #define mytest_str(name, input, separator, expected) \
962   xbt_test_add(name); \
963   d=xbt_str_split_str(input, separator); \
964   s=xbt_str_join(d,"XXX"); \
965   xbt_test_assert(!strcmp(s,expected),\
966                    "Input (%s) leads to (%s) instead of (%s)", \
967                    input,s,expected);\
968                    free(s); \
969                    xbt_dynar_free(&d);
970
971 XBT_TEST_UNIT("xbt_str_split_str", test_split_str, "test the function xbt_str_split_str")
972 {
973   xbt_dynar_t d;
974   char *s;
975
976   mytest_str("Empty string and separator", "", "", "");
977   mytest_str("Empty string", "", "##", "");
978   mytest_str("Empty separator", "toto", "", "toto");
979   mytest_str("String with no separator in it", "toto", "##", "toto");
980   mytest_str("Basic test", "toto##tutu", "##", "totoXXXtutu");
981 }
982
983 /* Last args are format string and parameters for xbt_test_add */
984 #define mytest_diff(s1, s2, diff, ...)                                  \
985   do {                                                                  \
986     char *mytest_diff_res;                                              \
987     xbt_test_add(__VA_ARGS__);                                          \
988     mytest_diff_res = xbt_str_diff(s1, s2);                             \
989     xbt_test_assert(!strcmp(mytest_diff_res, diff),                     \
990                     "Wrong output:\n--- got:\n%s\n--- expected:\n%s\n---", \
991                     mytest_diff_res, diff);                             \
992     free(mytest_diff_res);                                              \
993   } while (0)
994
995 XBT_TEST_UNIT("xbt_str_diff", test_diff, "test the function xbt_str_diff")
996 {
997   unsigned i;
998
999   /* Trivial cases */
1000   mytest_diff("a", "a", "  a", "1 word, no difference");
1001   mytest_diff("a", "A", "- a\n+ A", "1 word, different");
1002   mytest_diff("a\n", "a\n", "  a", "1 line, no difference");
1003   mytest_diff("a\n", "A\n", "- a\n+ A", "1 line, different");
1004
1005   /* Empty strings */
1006   mytest_diff("", "", "", "empty strings");
1007   mytest_diff("", "a", "+ a", "1 word, added");
1008   mytest_diff("a", "", "- a", "1 word, removed");
1009   mytest_diff("", "a\n", "+ a", "1 line, added");
1010   mytest_diff("a\n", "", "- a", "1 line, removed");
1011   mytest_diff("", "a\nb\nc\n", "+ a\n+ b\n+ c", "4 lines, all added");
1012   mytest_diff("a\nb\nc\n", "", "- a\n- b\n- c", "4 lines, all removed");
1013
1014   /* Empty lines */
1015   mytest_diff("\n", "\n", "  ", "empty lines");
1016   mytest_diff("", "\n", "+ ", "empty line, added");
1017   mytest_diff("\n", "", "- ", "empty line, removed");
1018
1019   mytest_diff("a", "\na", "+ \n  a", "empty line added before word");
1020   mytest_diff("a", "a\n\n", "  a\n+ ", "empty line added after word");
1021   mytest_diff("\na", "a", "- \n  a", "empty line removed before word");
1022   mytest_diff("a\n\n", "a", "  a\n- ", "empty line removed after word");
1023
1024   mytest_diff("a\n", "\na\n", "+ \n  a", "empty line added before line");
1025   mytest_diff("a\n", "a\n\n", "  a\n+ ", "empty line added after line");
1026   mytest_diff("\na\n", "a\n", "- \n  a", "empty line removed before line");
1027   mytest_diff("a\n\n", "a\n", "  a\n- ", "empty line removed after line");
1028
1029   mytest_diff("a\nb\nc\nd\n", "\na\nb\nc\nd\n", "+ \n  a\n  b\n  c\n  d",
1030               "empty line added before 4 lines");
1031   mytest_diff("a\nb\nc\nd\n", "a\nb\nc\nd\n\n", "  a\n  b\n  c\n  d\n+ ",
1032               "empty line added after 4 lines");
1033   mytest_diff("\na\nb\nc\nd\n", "a\nb\nc\nd\n", "- \n  a\n  b\n  c\n  d",
1034               "empty line removed before 4 lines");
1035   mytest_diff("a\nb\nc\nd\n\n", "a\nb\nc\nd\n", "  a\n  b\n  c\n  d\n- ",
1036               "empty line removed after 4 lines");
1037
1038   /* Missing newline at the end of one of the strings */
1039   mytest_diff("a\n", "a", "  a", "1 line, 1 word, no difference");
1040   mytest_diff("a", "a\n", "  a", "1 word, 1 line, no difference");
1041   mytest_diff("a\n", "A", "- a\n+ A", "1 line, 1 word, different");
1042   mytest_diff("a", "A\n", "- a\n+ A", "1 word, 1 line, different");
1043
1044   mytest_diff("a\nb\nc\nd", "a\nb\nc\nd\n", "  a\n  b\n  c\n  d",
1045               "4 lines, no newline on first");
1046   mytest_diff("a\nb\nc\nd\n", "a\nb\nc\nd", "  a\n  b\n  c\n  d",
1047               "4 lines, no newline on second");
1048
1049   /* Four lines, all combinations of differences */
1050   for (i = 0 ; i < (1U << 4) ; i++) {
1051     char d2[4 + 1];
1052     char s2[4 * 2 + 1];
1053     char res[4 * 8 + 1];
1054     char *pd = d2;
1055     char *ps = s2;
1056     char *pr = res;
1057     unsigned j = 0;
1058     while (j < 4) {
1059       unsigned k;
1060       for (/* j */ ; j < 4 && !(i & (1U << j)) ; j++) {
1061         *pd++ = "abcd"[j];
1062         ps += sprintf(ps, "%c\n", "abcd"[j]);
1063         pr += sprintf(pr, "  %c\n", "abcd"[j]);
1064       }
1065       for (k = j ; k < 4 && (i & (1U << k)) ; k++) {
1066         *pd++ = "ABCD"[k];
1067         ps += sprintf(ps, "%c\n", "ABCD"[k]);
1068         pr += sprintf(pr, "- %c\n", "abcd"[k]);
1069       }
1070       for (/* j */ ; j < k ; j++) {
1071         pr += sprintf(pr, "+ %c\n", "ABCD"[j]);
1072       }
1073     }
1074     *pd = '\0';
1075     *--pr = '\0';               /* strip last '\n' from expected result */
1076     mytest_diff("a\nb\nc\nd\n", s2, res,
1077                 "compare (abcd) with changed (%s)", d2);
1078   }
1079
1080   /* Subsets of four lines, do not test for empty subset */
1081   for (i = 1 ; i < (1U << 4) ; i++) {
1082     char d2[4 + 1];
1083     char s2[4 * 2 + 1];
1084     char res[4 * 8 + 1];
1085     char *pd = d2;
1086     char *ps = s2;
1087     char *pr = res;
1088     unsigned j = 0;
1089     while (j < 4) {
1090       for (/* j */ ; j < 4 && (i & (1U << j)) ; j++) {
1091         *pd++ = "abcd"[j];
1092         ps += sprintf(ps, "%c\n", "abcd"[j]);
1093         pr += sprintf(pr, "  %c\n", "abcd"[j]);
1094       }
1095       for (/* j */; j < 4 && !(i & (1U << j)) ; j++) {
1096         pr += sprintf(pr, "- %c\n", "abcd"[j]);
1097       }
1098     }
1099     *pd = '\0';
1100     *--pr = '\0';               /* strip last '\n' from expected result */
1101     mytest_diff("a\nb\nc\nd\n", s2, res,
1102                 "compare (abcd) with subset (%s)", d2);
1103
1104     for (pr = res ; *pr != '\0' ; pr++)
1105       if (*pr == '-')
1106         *pr = '+';
1107     mytest_diff(s2, "a\nb\nc\nd\n", res,
1108                 "compare subset (%s) with (abcd)", d2);
1109   }
1110 }
1111
1112 #endif                          /* SIMGRID_TEST */