{'fraction_non_alphanumeric': 0.047503603721661646, 'fraction_numerical': 0.008583409775914035, 'mean_word_length': 4.32258064516129, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 4, 'https://': 1, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Invariant Checking of NRA Transition Systems via Incremental Reduction to LRA with EUF

Model checking invariant properties of designs, represented as transition systems, with non-linear real arithmetic (NRA), is an important though very hard problem. On the one hand NRA is a hard-to-solve theory; on the other hand most of the powerful model checking techniques lack support for NRA. In this paper, we present a counterexample-guided abstraction refinement (CEGAR) approach that leverages linearization techniques from differential calculus to enable the use of mature and efficient model checking algorithms for transition systems on linear real arithmetic (LRA) with uninterpreted functions (EUF). The results of an empirical evaluation confirm the validity and potential of this approach.


Introduction
Invariant checking for infinite-state transition systems is a fundamental research area. Based on the recent improvements of SMT technologies, effective approaches have been developed for the case of transition systems with dynamics over Linear Real Arithmetic [4,9,18,21]. However, many real-world industrial designs (e.g. aerospace, automotive) require modeling as transition systems over non-linear arithmetic (NRA). Although both problems are undecidable, proving properties of the NRA transition systems turns out to be much harder than the linear case, and has in fact received much less attention. Approaches based on BMC and k-induction [15,29] are possible, so that nonlinearity is handled at the SMT-level, by means of an SMT(NRA) solver (e.g. Z3 [13], nlSAT [20], Yices [14], SMT-RAT [1]). Their power is however limited. Consider the following simple transition system: initially, x ≥ 2 ∧ y ≥ 2 ∧ z = x * y; the transition relation is defined by x = x+1∧y = y+1∧z = x * y . The property "it is always the case that z ≥ x + y" is not k-inductive, not even for a very large value of k. Thus, the typical proving techniques that are based on k-induction using an SMT(NRA) solver will not be able to prove it. In principle, it is also possible to lift other approaches (e.g. interpolation, IC3 [9,26]) to handle non-linearities at the level of the solver. However, this requires the SMT(NRA) solver to carry out interpolation or quantifier elimination,

Related Work
There are not many tools that deal with NRA transition systems. The most relevant is the recently proposed ISAT3 [28], that uses an interpolation-based [23,24] approach to prove invariants. In addition to NRA, it also supports trascendental functions and some form of differential equations. ISAT3 is built on an SMT solver based on numeric techniques (interval arithmetic), and is able to provide results that are accurate up to the specified precision. In fact, in addition to "safe" and "unsafe" answers, ISAT3 may return "maybe unsafe" when it finds an envelope of given precision that may (but is not guaranteed to) contain a counterexample. Another relevant tool is DREACH [22], a bounded model checker implemented on top of the DREAL [16] SMT solver, that adopts numerical techniques similar to ISAT3. DREACH has an expressiveness similar to ISAT3, but being a bounded model checker it is unable to prove properties.
The work in [8] follows a reduction-based approach to check invariants of NRA transition systems. It over-approximates the non-linear terms with a coarse abstraction, encoding into LRA some weak properties of multiplication like identity and sign. Another reduction-based approach is presented in [25] in the context of program analysis. The idea is to find a (tighter) convex approximation of polynomials in form of polyhedra, thus obtaining a conservative linear transition system. The key differences of our approach with respect to [8,25] are that we iteratively refine the abstraction, and we adopt a reduction to LRA+EUF. Furthermore, to the best of our knowledge, there is no available implementation of the approach [25] in a program analysis tool -it has been only shown to work on SMT problems.
The idea of approximating a univariate function (in particular the natural logarithm ln) with tangent lines is used in [30]. Here we abstract a bivariate function (multiplication), and use tangent planes for the refinement. We also exploit other properties (e.g. monotonicity) to derive additional axioms. The idea of using tangent planes (spaces) has been explored in [27], limited to the case of SMT solving. Another key differences is that the tangent planes area used to under-approximate predicates, while we use them to refine the over-approximation of the multiplication function.

Background
Properties of the Multiplication Function. Geometrically, the surface generated by the multiplication function f (x, y) def = x * y is shown in Fig. 1a and 1b. This kind of surface is known in geometry as hyperbolic paraboloid. A hyperbolic paraboloid is a doubly-ruled surface, i.e. for every point on the surface, there are two distinct lines projected from the surface such that they pass through the point. In case of the multiplication surface, the projected lines basically lie on the surface. Tangent Plane. The tangent plane to a surface at a point of interest (a, b) is a plane that "just touches" the surface at the point. The tangent planes can be used to linearly approximate the surface at the point of interest. An important property of the tangent plane to a hyperbolic paraboliod is that the two projected lines from the surface are also in the tangent plane, and they define how the plane cuts the surface (see Fig. 1c and 1d). The tangent plane T mul a,b (x, y) to the multiplication function f (x, y) at point (a, b) is calculated as follows: where d dx f (x, y) | (a,b) and d dy f (x, y) | (a,b) are the first-order partial derivatives of f (x, y) w.r.t. x and y respectively, evaluated at (a, b). T mul a,b (x, y) simplifies to: Logic and Satisfiability. We assume the standard first-order quantifier-free logical setting and standard notions of theory, model, satisfiability, and logical consequence. If ϕ is a formula, we denote with vars(ϕ) the set of its variables, and with atoms(ϕ) the set of its atoms. We write ϕ(X) to denote that vars(ϕ) ⊆ X. If x and y are two variables, we denote with ϕ{x → y} the formula obtained by replacing all the occurrences of x in ϕ with y. We extend this notation to ordered sequences of variables in the natural way.
If µ is a model and x is a variable, we write µ[x] to denote the value of x in µ, and we extend this notation to terms in the usual way. If X is a set of variables, we denote with X the set obtained by replacing each element x ∈ X with x , and with X i the set obtained by replacing x with x i . If Γ is a set of formulas, we write Γ to denote the formula obtained by taking the conjunction of all its elements. If Γ is unsatisfiable (modulo some theory T ), an unsatisfiable core is a set C ⊆ Γ such that C is still unsatisfiable.
Symbolic Transition Systems. A symbolic transition system S def = X, I, T is a tuple where X is a finite set of (state) variables, I(X) is a formula denoting the initial states of the system, and T (X, X ) is a formula expressing its transition relation. A state s i of S is an assignment to the variables X. A path (execution trace) π = s 0 , s 1 , s 2 , . . . , s k−1 of length k (possibly infinite) for S is a sequence of states such that s 0 |= I and s i ∧ s i+1 {X → X } |= T for all 0 ≤ i < k − 2. We call an unrolling of S of length k the formula Let P (X) be a formula whose assignments represent a property (good states) over the state variables X. The invariant verification problem, denoted with S |= P , is the problem of checking if for all the finite paths s 0 , s 1 , . . . , s k of S, for all i, 0 ≤ i ≤ k, s i |= P . Its dual formulation in terms of reachability of ¬P is the problem of finding a path s 0 , s 1 , . . . , s k of S such that s k |= ¬P . P represents the "good" states, while ¬P represents the "bad" states.  overapproximating it with a formula over the combined theory of linear arithmetic and uninterpreted functions (LRA+EUF). Our main SMT solving procedure follows a classic abstraction refinement loop, in which at each iteration the current overapproximation of the input SMT formula is refined by adding new constraints that rule out one (or possibly more) spurious solutions, until one of the following occurs: (i) the SMT formula becomes unsatisfiable in the LRA+EUF theory; or (ii) the LRA+EUF model for the current overapproximation can be lifted to an NRA model for the original SMT formula; or (iii) the resource budget (e.g. time, memory, number of iterations) is exhausted.
The pseudocode for the top-level algorithm is shown in Fig. 2. We provide more details about its main components in the rest of this section.
Initial Abstraction. The function initial-abstraction takes as input an SMT(NRA) formula ϕ and returns an overapproximation ϕ of it in the LRA+EUF theory.
First, each multiplication expression x * y between two variables 3 occurring in ϕ is replaced by f mul(x, y), where f mul() is a binary uninterpreted function returning a real. We remark that this happens only for non-linear multiplications: expressions like c * x or x * c in which c is a constant are not rewritten.
Then, some simple axioms about multiplication are added to ϕ via static learning. For each f mul(x, y) ∈ ϕ, we add the following axioms: Abstraction Refinement. If the SMT check on the LRA+EUF abstraction returns false (line 7 of Fig. 2), we can conclude that the input formula is unsatisfiable. In this case, Γ contains all the lemmas (discussed later in this section) that were added in the earlier refinements (line 11 of Fig. 2).
Otherwise, we have to check whether the model µ found for ϕ is also a model for the original NRA formula ϕ. Let F muls be the set of all f mul(x, y) terms occurring in ϕ. In its simplest version, the function get-NRA-model checks whether, for all f mul(x, y) . If this is the case, then µ is also a model for the original formula, and get-NRA-model returns true. (We present more sophisticated versions of get-NRA-model below.) Otherwise, let CF muls be the set of all f mul(x, y) terms whose value in µ is different from µ[x] * µ[y]. The function refine generates a set of axioms Γ such that there exists at least one element f mul(x, y) of CF muls such that the formula ϕ ∧ Γ has no model µ that agrees with µ on the values of x, y and f mul( . Intuitively, the axioms Γ block the bad model values for f mul(x, y), making the abstraction more precise by restricting the set of spurious solutions.
In our current implementation, two kinds of lemmas are generated during refinement: tangent lemmas and monotonicity lemmas.
where we can choose a and b as: Basically the equalities in the tangent lemma are providing multiplication lines that enforce the correct value of f mul(x, y) when x = a or y = b. Moreover, the inequalities of the tangent lemma are providing bounds for f mul(x, y) when x and y are not on the multiplication lines.
Monotonicity Lemmas. Let f mul(x, y) and f mul(w, z) be two terms in ϕ, such that . Then, we add the monotonicity lemma (abs(x) ≤ abs(w)∧abs(y) ≤ abs(z)) → abs(f mul(x, y)) ≤ abs(f mul(w, z)), (6) where abs(t) stands for ite(t < 0, −t, t). Finding Models. It is easy to see that our algorithm is expected to perform much better for unsatisfiable instances than for satisfiable ones. The algorithm can return true (meaning that the formula is satisfiable) only if the LRA+EUF solver "guesses" a model that is consistent with all the nonlinear multiplications. In an infinite and dense domain like the reals, the chances that this will happen are close to zero in general. Moreover, our approach is inherently limited, because it can only find models over the rationals. If the input formula is satisfiable, but all its models contain some irrational values, then our algorithm will always abort (or never terminate, if there is no resource budget set). In practice, it is very likely that the same will happen even for formulas admitting a rational solution.
One possibility for addressing this limitation would be to couple our procedure with a complete solver for NRA, to be used for detecting satisfiable cases, in order to implement a more effective version of get-NRA-model. One such possibility is shown in Fig. 3, where we extract the truth assignment ψ induced by the LRA+EUF model µ on the atoms of ϕ: We concretize it by replacing each f mul(x, y) in ψ with x * y, and invoke the complete NRA theory solver on the resulting conjunction of NRA-literals ψ, to check whether it contains at least one solution. Although in general the problem is expected to be simpler than the original input formula because the Boolean structure of ϕ is disregarded, invoking a complete NRA theory solver at each loop iteration of SMT-NRA-check-abstract-ext could be very expensive. Moreover, this would still require a complete NRA theory solver, which might not always be available.
As an alternative, we propose the procedure outlined in Fig. 4, where we extract the truth assignment ψ induced by the LRA+EUF model µ on the atoms of ϕ, and we conjoin to it the multiplication lines: F muls being the usual set of all f mul(x, y) terms occurring in ϕ.
The main idea is to build an LRA+EUF underapproximation ψ * of the NRA formula ψ of Fig. 3, in which all multiplications are forced to be linear. Compared to the previous solution, this has the advantage of requiring a complete SMT(LRA+EUF) solver rather than a (much more expensive) complete NRA solver. Moreover, given the simplicity of the Boolean structure of the underapproximated formula, the check should in general be very cheap. The drawback is that this is (clearly) still an incomplete procedure. However, in our experiments (for which we refer to Sec. 6) we have found it to be surprisingly effective for many problems.
Unlike with the basic implementation of get-NRA-model which considers only one single candidate model at a time, the implementations in Fig.3 and Fig. 4 consider an infinite amount of them, drastically increasing the chances of finding a model.
Correctness and progress. We notice that the procedure in Fig. 2 is correct. In fact, it returns false only if ϕ is NRA-unsatisfiable because by construction ϕ is an overapproximation of ϕ, and all axioms in Γ are valid in any theory interpreting f mul(x, y) as x * y. Also, it returns true only if ϕ is NRA-satisfiable: if get-NRA-model is based only on evaluation, then by construction µ is an LRA+EUF-model for ϕ s.t. each f mul(x, y) equals x * y in µ, so that µ is also a model for ϕ; if get-NRA-model is as in Fig. 3, then µ is an NRA-model of a conjunction of literals ψ which tautologically entails ϕ, so that µ is a model for ϕ; if get-NRA-model is as in Fig. 4, then µ is an LRA+EUF-model of a conjunction of literals ψ * which tautologically entails ϕ and it is s.t. each f mul(x, y) equals x * y in µ, so that µ is a also model for ϕ.
We also notice that the progress of the procedure in Fig. 2 is guaranteed by the refinement step, which rules out significant parts of the search space at every loop by means of the added lemmas.
Important Heuristics for Refinement. The description of refine provided above leaves some flexibility in deciding what axioms to add (and how many of them) at each iteration. It is possible to conceive strategies with an increasing degree of eagerness, from very lazy (e.g. adding only a single axiom per iteration) to more aggressive ones. In our current implementation, we eagerly add all the axioms (2)-(6) that are violated by the current abstract solution µ, leaving the investigation of alternative strategies as future work. However, we found the following two strategies to be crucial for performance.
Tangent lemma frontiers. The tangent lemmas of (2) for a given point (a, b) are based on the fact that the multiplication function x * y is a hyperbolic paraboloid surface, and a tangent plane to such surface cuts the surface into four regions such that in two of the regions the tangent plane is above the surface, whereas in the other two regions the tangent plane is below the surface (see Fig. 1). Each instantiation of (2) for a given point, therefore, can only provide either a lower or an upper bound for a given region. In some cases, this might lead to an infinite refinement loop in which at each iteration the "wrong" bound is refined. In order to address the problem, we use the following strategy. For each f mul(x, y) in the input formula, we maintain a frontier l x , u x , l y , u y with the invariant that whenever x is in the interval [l x , u x ] or y is in the interval [l y , u y ], then f mul(x, y) has both an upper and a lower bound. Initially, the frontiers are set to 0, 0, 0, 0 . Whenever a lemma (2) for f mul(x, y) is instantiated on a point (a, b), we generate further instantiations of (2) and update the frontier as follows: case a < l x and b < l y : instantiate (2) on (a, u y ) and on (u x , b), and set the frontier to a, u x , b, u y ; case a < l x and b > u y : instantiate (2) on (a, l y ) and on (u x , b), and set the frontier to a, u x , l y , b ; case a > u x and b > u y : instantiate (2) on (a, l y ) and on (l x , b), and set the frontier to l x , a, l y , b ; case a > u x and b < l y : instantiate (2) on (a, u y ) and on (l x , b), and set the frontier to l x , a, b, u y . Tangent lemma rounding. The instantiation of a tangent lemma at the point (a, b) has the side-effect of adding the rational constants a, b and a * b to the formula that is solved by the LRA+EUF solver. If such values have large numerators and/or denominators, they might be a source of a significant slow-down for the LRA solver (which works on exact, arbitrary-precision rational arithmetic). We address this issue by observing that, in order to block a bad model µ such that µ[f mul(x, y)] = µ[x] * µ[y], it is sufficient to add one of the two equalities of (2); therefore, instead of instantiating a tangent lemma at (a, b), we can instantiate it at either (a + δ, b) or at (a, b + δ), for any value of δ. In practice, if a (resp. b) is a rational constant with a very large numerator or denominator, instead of instantiating a tangent lemma at (a, b), we instantiate two tangent lemmas at ( a , b) and ( a , b).

From Satisfiability to Verification
We now move from satisfiability checking to verification.
Overview. In principle, the solver described in the previous section could be integrated as a "black box" in any off-the-shelf SMT-based verification algorithm, such as BMC, k-induction, or one of the many extensions of IC3 to the SMT case (e.g. [4,9,18,21]). In practice, however, such black-box integration would hardly be effective, especially in the case of state-of-the-art algorithms like IC3. IC3 requires a very incremental interaction with the underlying SMT engine, which is asked to solve a large number of relatively-cheap queries. The procedure of Sec. 4, however, can be very expensive, especially for satisfiable queries, which are very common in an IC3-like algorithm. 4 Moreover, some of the IC3 extensions mentioned above require the ability of performing (approximated) quantifier eliminations, a functionality not provided by the algorithm of Fig. 2. We propose therefore a white-box integration, in which we lift the abstraction refinement approach of Sec. 4 at the transition system level. We generate an abstract LRA+EUF version of the input NRA transition system, which is then checked with the IC3-based procedure of [9]. In case a counterexample is produced, we use the SMT-NRA-check-abstract-ext algorithm of Fig. 2 to check whether it is spurious. If so, the axioms generated by SMT-NRA-check-abstract-ext are then used to refine the abstraction of the transition system. The pseudo-code of this algorithm is reported in Fig. 6. Similarly to the satisfiability checking case, the initial-abstraction function replaces every non-linear multiplication x * y in the input transition system and property with a f mul(x, y) term, and adds some simple axioms about the behaviour of multiplication to the initial-state and transition-relation formulas of the transition system (see Sec. 4). In the rest of this section, we describe the abstraction refinement algorithm in more detail.
transition system refine-transition-system ( S : transition system, Γ : set of axioms): 1. let X, I, T = S 2. Γ I , Γ T = ∅, ∅ 3. for each γ in Γ : 4. if vars(γ) ⊆ X 0 : 5.  Counterexample Checking and Refinement. When IC3-LRA+EUF-prove returns a counterexample trace π for the abstract system S, we use SMT-NRA-check-abstract-ext to check for its spuriousness. The function get-cex-formula builds a formula ψ to feed to SMT-NRA-check-abstract-ext, whose unsatisfiability implies that π is spurious. The formula ψ is built by unrolling the transition relation of S, and optionally adding constraints that restrict the allowed transitions to be compatible with the states in π. Various heuristics are possible, trading generality for complexity: ψ could be fully constrained by the states in π (thus checking only one abstract counterexample path per iteration); it could be only partially constrained (e.g. by considering only the Boolean variables and/or the state variables occurring only in linear constraints); or it could be left unconstrained, considering only the length of the abstract counterexample. In our current implementation (see Sec. 6), we use the last option, i.e. we only consider the length of π to build a BMC formula that checks for any counterexample of the given length, leaving the investigation of alternative strategies to future work.
If SMT-NRA-check-abstract-ext returns true, the property is violated. In this case, we can use the model found by SMT-NRA-check-abstract-ext to build a counterexample trace for the input system and property.
If SMT-NRA-check-abstract-ext returns false, we use the axioms Γ produced during search to refine the transition system S, using the procedure shown in Fig. 7. Essentially, refine-transition-system translates back the axioms from their unrolled version (i.e. on variables X 0 , X 1 , . . .) to their "single step" version (on variables X and X ), adding each of them either to the initial-states formula or to the transition relation formula. In case an axiom γ spans more than a single transition step (lines 9-10 of Fig. 7), we arbitrarily choose to map the variables with the lowest index as current state variables X, and all the others as next-state variables X . Notice that this might cause some refinement failure, as discussed in the next paragraph.
Reducing the number of axioms to add. In general, not all the axioms generated during a call to SMT-NRA-check-abstract-ext are needed to successfully block a counterexample, especially if eager strategies like those described in Sec. 4 are used. In the long run, having a large number of redundant axioms can be quite harmful for performance. In let C be an unsatisfiable core of ψ 4.
abort # refinement failure Fig. 8. Reducing the axioms needed for refinement. order to mitigate this problem, we apply a filtering strategy (based on unsatisfiable cores) to the set of axioms, before adding them to the transition system. Instead of adding Γ I and Γ T directly to S, we invoke the function shown in Fig. 8. Note that due to the flattening of multi-step axioms described above (lines 9-10 of Fig. 7), the refinement might fail. In this case, our current implementation simply aborts the execution. 5 6 Experimental Analysis Implementation and comparisons.
We have implemented a prototype of the IC3-NRA-prove procedure using the IC3 engine of NUXMV [7] for IC3-LRA+EUF-prove. The code is written in Python, using the PYSMT library [17]. Our implementation, benchmarks, and experimental data are available at https://es-static.fbk.eu/people/griggio/papers/ tacas17-ic3-nra.tar.gz. We have used the following tools for our evaluation. NUXMV-LRA-static: we apply the upfront abstraction of NRA to LRA proposed in [8], running the IC3 engine of NUXMV on the resulting transition system. NRA-BMC-{Z3, DREAL } and NRA-K-induction-{Z3, DREAL }: we have implemented the BMC [3] and k-induction [29] algorithms in Python (using PYSMT), using either Z3 (NRA) or DREAL (NRA) as back-end SMT solver. iSAT3[1e-1] and iSAT3[1e-9]: we have used the latest version of the ISAT3 solver [24], which combines an SMT solver integrating CDCL and interval constraint propagation techniques with an interpolation-based abstaction/refinement algorithm for verification. ISAT3 supports both transition systems and software programs encoded as control flow graphs. Similarly to DREAL, ISAT3 may return a "maybe unsafe" answer and provide a candidate solution identifying the upper and lower bounds on the variables. In the experiments, iSAT3[1e-1] is the configuration suggested by the ISAT3 authors 6 and iSAT3[1e-9] is the same except that the minimum splitting width (msw) parameter is set 5 We remark however that so far we have never observed this behaviour during our experiments. to 10 −9 . We have used a smaller value for the msw to get more precise answers, i.e. "safe" or "unsafe", as suggested in the ISAT3 user manual.
Benchmarks. We have collected a total of 114 NRA benchmarks from various sources.
Handcrafted. This set contains 14 hand-written instances, 13 safe and 1 unsafe.
HYST. This is the biggest set, consisting of 65 benchmarks. These are generated from the Hybrid examples that come with the HYST [2] distribution, by approximating the continuous time by sampling at a fixed time interval. This process is done automatically using an extended version of HYST. Since the generated benchmarks are approximations, we do not know their safety status. The benchmarks contain mostly non-linear behaviour. iSAT3 and iSAT3-CFG. The 11 benchmarks in this set (7 safe, 4 unsafe) are taken from [24] and the ISAT3 examples available online. nuXmv. In this set, we have 2 safe benchmarks which we collected from the nuXmv users' mailing list. These benchmarks have complex boolean structure. SAS13. These 13 benchmarks are generated from the C programs used in [5], but interpreted over NRA instead of the theory of IEEE floating-point numbers. This makes some of the instances unsafe. TCM. We have generated 2 safe benchmarks from the Simulink models (taken from the case study [6]) by first generating the C code using the Embedded Coder 7 and then encoding the program into a symbolic transition system.  (8)  Each column shows a benchmark family, and each entry gives the number of safe/unsafe instances found. For tools working over interval arithmetic, the number of "maybe unsafe" is reported in parentheses. Results. We ran our experiments on a cluster of machines with 2.67GHz Xeon X5650 CPUs and 96GB of RAM, running Scientific Linux 6.7. We used 6GB memory limit and 3600 seconds CPU timeout. The results are summarized in Tables 1 and 2 and in Fig. 9. The plots show the time to solve an instance on the x-axis and the total number of solved instances on the y-axis. Table 1 reports a summary of the solved instances by family, whereas Table 2 shows a comparitive analysis by reporting for each tool the number of uniquely solved instances and the difference of solved instances w.r.t IC3-NRA-prove. We can make the following observations from the experimental results: -IC3-NRA-prove is the best performer overall, and it significantly outperforms all the other approaches on safe instances (where it can solve 9 problems that are out of reach for all the other tools). Interestingly, despite its simplicity, our model finding approach (as outlined in §4) is surprisingly effective, allowing IC3-NRA-prove to find 15 counterexample traces. -The simple abstraction proposed in [8] is quite effective for many families, allowing NUXMV-LRA-static to verify more properties than the approaches based on Kinduction with an NRA solver. However, IC3-NRA-prove results in a clear and very significant improvement, solving more than twice as many instances than NUXMV-LRA-static (and losing only 1). -None of the other tools (with the exception of NUXMV-LRA-static) is able to solve any safe benchmark in the HyComp and nuXmv families. These benchmarks have a non-trivial Boolean structure and a significant linear component. Both IC3-NRA-prove and NUXMV-LRA-static are able to fully exploit the effectiveness of the underlying IC3 engine of NUXMV, outperforming the competitors. However, IC3-NRA-prove is very competitive also on the HYST family, whose instances are mostly non-linear and have very little Boolean structure. -Increasing the default precision of ISAT3 significantly reduces the number of "maybe unsafe" answers, but it doesn't seem to help in solving more benchmarks.
In fact, we remark that even with the increased precision iSAT3[1e-9] classifies 2 safe instances as "maybe unsafe" (whereas in the default configuration, 6 safe instances are classified as "maybe unsafe").

Conclusions and Future Work
We presented a novel abstraction-refinement approach to the verification of transition systems with nonlinear dynamics expressed in the NRA theory. We abstract non-linear multiplication as an uninterpreted function, leveraging efficient invariant checkers for transition systems over LRA and EUF to solve the problem in the abstract space. In case of spurious counterexample, the abstraction of multiplication is incrementally refined by introducing suitable axioms, based on the idea of tangent planes. An extensive experimental evaluation demonstrates that the proposed approach is significantly more effective than approaches directly based on SMT(NRA) solving. This work opens up several important directions. First, we are going to improve the implementation, by integrating all the steps within the NUXMV [7] model checker, and to perform a thorough analysis of the various heuristic choices. Second, we will investigate the potential of the approach for SMT, both for other theories (e.g. NIA) and for extended functionalities (e.g. interpolation). We will also extend the scope of the approach to deal with transcendental functions, look-up tables, and partially axiomatized functions (e.g. gain functions known to be monotonic and of restricted co-domain).
Finally, we are going to investigate the generalization of the approach from transition systems to continuous-time hybrid systems with nonlinear characteristic functions.