complementation.cc 7.82 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
#include <iomanip>
#include <iostream>
#include "tgbaalgos/dotty.hh"
#include "tgbaparse/public.hh"
#include "tgba/tgbaproduct.hh"
#include "tgbaalgos/gtec/gtec.hh"
#include "tgbaalgos/ltl2tgba_fm.hh"
#include "ltlparse/public.hh"
#include "tgbaalgos/stats.hh"
#include "tgbaalgos/emptiness.hh"
#include "ltlast/unop.hh"
#include "tgbaalgos/stats.hh"
#include "tgbaalgos/emptiness_stats.hh"
#include "ltlvisit/destroy.hh"
#include "ltlvisit/clone.hh"
#include "tgba/tgbatba.hh"

18
#include "tgba/tgbasafracomplement.hh"
19
#include "tgba/tgbacomplement.hh"
20
21
22
23
24

void usage(const char* prog)
{
  std::cout << "usage: " << prog << " [options]" << std::endl;
  std::cout << "with options" << std::endl
25
26
            << "-S                      Use safra complementation"
            << std::endl
27
28
29
30
31
32
            << "-s     buchi_automaton  display the safra automaton"
            << std::endl
            << "-a     buchi_automaton  display the complemented automaton"
            << std::endl
            << "-astat buchi_automaton  statistics for !a" << std::endl
            << "-fstat formula          statistics for !A_f" << std::endl
33
34
            << "-f     formula          test !A_f and !A_!f" << std::endl
            << "-p     formula          print the automaton for f" << std::endl;
35
36
37
38
39
40
41
42
43
44
45
46
}

int main(int argc, char* argv[])
{
  char *file = 0;
  bool print_safra = false;
  bool print_automaton = false;
  bool check = false;
  int return_value = 0;
  bool stats = false;
  bool formula = false;
  bool automaton = false;
47
48
  bool safra = false;
  bool print_formula = false;
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75

  if (argc < 3)
  {
    usage(argv[0]);
    return 1;
  }

  for (int i = 1; i < argc; ++i)
  {
    if (argv[i][0] == '-')
    {
      if (strcmp(argv[i] + 1, "astat") == 0)
      {
        stats = true;
        automaton = true;
        continue;
      }

      if (strcmp(argv[i] + 1, "fstat") == 0)
      {
        stats = true;
        formula = true;
        continue;
      }

      switch (argv[i][1])
      {
76
77
        case 'S':
          safra = true; break;
78
        case 's':
79
          safra = true; print_safra = true; break;
80
81
82
83
        case 'a':
          print_automaton = true; break;
        case 'f':
          check = true; break;
84
85
        case 'p':
          print_formula = true; break;
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
        default:
          std::cerr << "unrecognized option `-" << argv[i][1]
                    << "'" << std::endl;
          return 2;
      }
    }
    else
      file = argv[i];
  }

  if (file == 0)
  {
    usage(argv[0]);
    return 1;
  }

  spot::bdd_dict* dict = new spot::bdd_dict();
  if (print_automaton || print_safra)
  {
    spot::ltl::environment& env(spot::ltl::default_environment::instance());
    spot::tgba_parse_error_list pel;
    spot::tgba_explicit* a = spot::tgba_parse(file, pel, dict, env);
    if (spot::format_tgba_parse_errors(std::cerr, file, pel))
      return 2;

111
112
113
114
115
116
    spot::tgba* complement = 0;

    if (safra)
      complement = new spot::tgba_safra_complement(a);
    else
      complement = new spot::tgba_complement(a);
117
118
119
120
121

    if (print_automaton)
      spot::dotty_reachable(std::cout, complement);

    if (print_safra)
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
    {
      spot::tgba_safra_complement* safra_complement =
        dynamic_cast<spot::tgba_safra_complement*>(complement);
      spot::display_safra(safra_complement);
    }
    delete complement;
    delete a;
  }
  else if (print_formula)
  {
    spot::tgba* a;
    spot::ltl::formula* f1 = 0;

    spot::ltl::parse_error_list p1;
    f1 = spot::ltl::parse(file, p1);

    if (spot::ltl::format_parse_errors(std::cerr, file, p1))
      return 2;

    a = spot::ltl_to_tgba_fm(f1, dict);

    spot::tgba* complement = 0;
    if (safra)
      complement = new spot::tgba_safra_complement(a);
    else
      complement = new spot::tgba_complement(a);

    spot::dotty_reachable(std::cout, complement);
    spot::ltl::destroy(f1);
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
    delete complement;
    delete a;
  }
  else if (stats)
  {
    spot::tgba* a;
    spot::ltl::formula* f1 = 0;

    if (formula)
    {
      spot::ltl::parse_error_list p1;
      f1 = spot::ltl::parse(file, p1);

      if (spot::ltl::format_parse_errors(std::cerr, file, p1))
        return 2;

      a = spot::ltl_to_tgba_fm(f1, dict);
    }
    else
    {
      spot::tgba_parse_error_list pel;
      spot::ltl::environment& env(spot::ltl::default_environment::instance());
      a = spot::tgba_parse(file, pel, dict, env);
      if (spot::format_tgba_parse_errors(std::cerr, file, pel))
        return 2;
    }

178
    spot::tgba_safra_complement* safra_complement =
179
      new spot::tgba_safra_complement(a);
180
181
182
183

    spot::tgba_statistics a_size =  spot::stats_reachable(a);
    std::cout << "Original: "
              << a_size.states << ", "
184
185
186
              << a_size.transitions << ", "
              << a->number_of_acceptance_conditions()
              << std::endl;
187
188
189
190
191

    spot::tgba *buchi = new spot::tgba_sba_proxy(a);
    a_size =  spot::stats_reachable(buchi);
    std::cout << "Buchi: "
              << a_size.states << ", "
192
193
194
              << a_size.transitions << ", "
              << buchi->number_of_acceptance_conditions()
              << std::endl;
195
196
    delete buchi;

197
198
199
200
201
202
203
204
205
206
207
208
    spot::tgba_statistics b_size =  spot::stats_reachable(safra_complement);
    std::cout << "Safra Complement: "
              << b_size.states << ", "
              << b_size.transitions << ", "
              << safra_complement->number_of_acceptance_conditions()
              << std::endl;

    spot::tgba_complement* complement =
      new spot::tgba_complement(a);

    b_size =  spot::stats_reachable(complement);
    std::cout << "GBA Complement: "
209
              << b_size.states << ", "
210
211
212
              << b_size.transitions << ", "
              << complement->number_of_acceptance_conditions()
              << std::endl;
213
214
215
216
217
218
219
220
221
222
223
224

    delete complement;
    delete a;
    if (formula)
    {
      spot::ltl::formula* nf1 =
        spot::ltl::unop::instance(spot::ltl::unop::Not,
                                  spot::ltl::clone(f1));
      spot::tgba* a2 = spot::ltl_to_tgba_fm(nf1, dict);
      spot::tgba_statistics a_size =  spot::stats_reachable(a2);
      std::cout << "Not Formula: "
                << a_size.states << ", "
225
226
227
                << a_size.transitions << ", "
                << a2->number_of_acceptance_conditions()
                << std::endl;
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245

      delete a2;
      spot::ltl::destroy(f1);
      spot::ltl::destroy(nf1);
    }
  }
  else
  {
    spot::ltl::parse_error_list p1;
    spot::ltl::formula* f1 = spot::ltl::parse(file, p1);

    if (spot::ltl::format_parse_errors(std::cerr, file, p1))
      return 2;

    spot::tgba* Af = spot::ltl_to_tgba_fm(f1, dict);
    spot::ltl::formula* nf1 = spot::ltl::unop::instance(spot::ltl::unop::Not,
                                                        spot::ltl::clone(f1));
    spot::tgba* Anf = spot::ltl_to_tgba_fm(nf1, dict);
246
247
248
249
250
251
252
253
254
255
256
257
258

    spot::tgba* nAf;
    spot::tgba* nAnf;
    if (safra)
    {
      nAf = new spot::tgba_safra_complement(Af);
      nAnf = new spot::tgba_safra_complement(Anf);
    }
    else
    {
      nAf = new spot::tgba_complement(Af);
      nAnf = new spot::tgba_complement(Anf);
    }
259
260
261
262
    spot::tgba* prod = new spot::tgba_product(nAf, nAnf);
    spot::emptiness_check* ec = spot::couvreur99(prod);
    spot::emptiness_check_result* res = ec->check();
    spot::tgba_statistics a_size =  spot::stats_reachable(ec->automaton());
263
264
265
266
267
268
269
    std::cout << "States: "
              << a_size.states << std::endl
              << "Transitions: "
              << a_size.transitions << std::endl
              << "Acc Cond: "
              << ec->automaton()->number_of_acceptance_conditions()
              << std::endl;
270
271
    if (res)
    {
272
      std::cout << "FAIL";
273
274
275
      return_value = 1;
    }
    else
276
      std::cout << "OK";
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
    std::cout << std::endl;

    delete res;
    delete ec;
    delete prod;
    delete nAf;
    delete Af;
    delete nAnf;
    delete Anf;

    spot::ltl::destroy(nf1);
    spot::ltl::destroy(f1);

  }

  delete dict;

  return return_value;
}