1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
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
76
77
78
79
80
81
|
/*
* Copyright 2017 Sven Verdoolaege
*
* Use of this software is governed by the MIT license
*
* Written by Sven Verdoolaege.
*/
#include <isl_multi_macro.h>
/* Compute the sum of "multi1" and "multi2" on the union of their domains,
* with the actual sum on the shared domain and
* the defined expression on the symmetric difference of the domains.
*
* We simply iterate over the elements in both arguments and
* call isl_union_pw_aff_union_add on each of them, if there is
* at least one element.
*
* Otherwise, the two expressions have an explicit domain and
* the union of these explicit domains is computed.
* This assumes that the explicit domains are either both in terms
* of specific domains elements or both in terms of parameters.
* However, if one of the expressions does not have any constraints
* on its explicit domain, then this is allowed as well and the result
* is the expression with no constraints on its explicit domain.
*/
__isl_give MULTI(BASE) *FN(MULTI(BASE),union_add)(
__isl_take MULTI(BASE) *multi1, __isl_take MULTI(BASE) *multi2)
{
isl_bool has_domain, is_params1, is_params2;
if (!multi1)
goto error;
if (multi1->n > 0)
return FN(MULTI(BASE),bin_op)(multi1, multi2,
&FN(EL,union_add));
FN(MULTI(BASE),align_params_bin)(&multi1, &multi2);
if (FN(MULTI(BASE),check_equal_space)(multi1, multi2) < 0)
goto error;
if (FN(MULTI(BASE),check_has_explicit_domain)(multi1) < 0 ||
FN(MULTI(BASE),check_has_explicit_domain)(multi2) < 0)
goto error;
has_domain = FN(MULTI(BASE),has_non_trivial_domain)(multi1);
if (has_domain < 0)
goto error;
if (!has_domain) {
FN(MULTI(BASE),free)(multi2);
return multi1;
}
has_domain = FN(MULTI(BASE),has_non_trivial_domain)(multi2);
if (has_domain < 0)
goto error;
if (!has_domain) {
FN(MULTI(BASE),free)(multi1);
return multi2;
}
is_params1 = FN(DOM,is_params)(multi1->u.dom);
is_params2 = FN(DOM,is_params)(multi2->u.dom);
if (is_params1 < 0 || is_params2 < 0)
goto error;
if (is_params1 != is_params2)
isl_die(FN(MULTI(BASE),get_ctx)(multi1),
isl_error_invalid,
"cannot compute union of concrete domain and "
"parameter constraints", goto error);
multi1 = FN(MULTI(BASE),cow)(multi1);
if (!multi1)
goto error;
multi1->u.dom = FN(DOM,union)(multi1->u.dom,
FN(DOM,copy)(multi2->u.dom));
if (!multi1->u.dom)
goto error;
FN(MULTI(BASE),free)(multi2);
return multi1;
error:
FN(MULTI(BASE),free)(multi1);
FN(MULTI(BASE),free)(multi2);
return NULL;
}
|
