2024Activity reportProject-TeamHYCOMES
RNSR: 201321225U- Research center Inria Centre at Rennes University
- Team name: Modélisation hybride & conception par contrats pour les systèmes embarqués multi-physiques
- In collaboration with:Institut de recherche en informatique et systèmes aléatoires (IRISA)
- Domain:Algorithmics, Programming, Software and Architecture
- Theme:Embedded and Real-time Systems
Keywords
Computer Science and Digital Science
- A2.1. Programming Languages
- A2.1.1. Semantics of programming languages
- A2.1.9. Synchronous languages
- A2.1.10. Domain-specific languages
- A2.2. Compilation
- A2.2.1. Static analysis
- A2.2.8. Code generation
- A2.3. Embedded and cyber-physical systems
- A2.3.1. Embedded systems
- A2.3.2. Cyber-physical systems
- A2.3.3. Real-time systems
- A2.4. Formal method for verification, reliability, certification
- A2.4.1. Analysis
- A2.4.3. Proofs
- A6. Modeling, simulation and control
- A6.1. Methods in mathematical modeling
- A6.1.1. Continuous Modeling (PDE, ODE)
- A6.1.5. Multiphysics modeling
- A6.4.3. Observability and Controlability
- A6.4.4. Stability and Stabilization
- A6.4.5. Control of distributed parameter systems
- A6.5. Mathematical modeling for physical sciences
- A7.2.1. Decision procedures
- A7.2.2. Automated Theorem Proving
- A7.2.4. Mechanized Formalization of Mathematics
- A8. Mathematics of computing
- A8.4. Computer Algebra
Other Research Topics and Application Domains
- B4. Energy
- B4.4. Energy delivery
- B4.4.1. Smart grids
- B5. Industry of the future
- B5.1. Factory of the future
- B5.2. Design and manufacturing
- B5.9. Industrial maintenance
- B8. Smart Cities and Territories
- B8.1. Smart building/home
- B8.1.1. Energy for smart buildings
- B8.2. Connected city
- B8.3. Urbanism and urban planning
1 Team members, visitors, external collaborators
Research Scientists
- Benoit Caillaud [Team leader, INRIA, Senior Researcher]
- Albert Benveniste [INRIA, Emeritus]
- Yahao Chen [INRIA, ISFP]
- Khalil Ghorbal [INRIA, Researcher]
PhD Student
- Maxime Bridoux [INRIA]
Technical Staff
- Mathias Malandain [INRIA, Engineer]
Interns and Apprentices
- Morgane Flamant [INRIA, Intern, from May 2024 until Aug 2024]
Administrative Assistant
- Armelle Mozziconacci [CNRS]
2 Overall objectives
Hycomes was created as a local team of the Rennes - Bretagne Atlantique Inria research center in 2013 and has been created as an Inria Project-Team in 2016. The team is focused on two topics in cyber-physical systems design:
- Hybrid systems modeling, with an emphasis on the design of modeling languages in which software systems, in interaction with a complex physical environment, can be modelled, simulated and verified. A special attention is paid to the mathematical rigorous semantics of these languages, and to the correctness (wrt. such semantics) of the simulations and of the static analyses that must be performed during compilation. The Modelica language is the main application field. The team aims at contributing language extensions facilitating the modeling of physical domains which are poorly supported by the Modelica language. The Hycomes team is also designing new structural analysis methods for hybrid (aka. multi-mode) Modelica models. New simulation and verification techniques for large Modelica models are also in the scope of the team.
- Contract-based design and interface theories, with applications to requirements engineering in the context of safety-critical systems design. The objective of our research is to bridge the gap between system-level requirements, often expressed in natural, constrained or semi-formal languages and formal models, that can be simulated and verified.
3 Research program
3.1 Hybrid Systems Modeling
Systems industries today make extensive use of mathematical modeling tools to design computer controlled physical systems. This class of tools addresses the modeling of physical systems with models that are simpler than usual scientific computing problems by using only Ordinary Differential Equations (ODE) and Difference Equations but not Partial Differential Equations (PDE). This family of tools first emerged in the 1980's with SystemBuild by MatrixX (now distributed by National Instruments) followed soon by Simulink by Mathworks, with an impressive subsequent development.
In the early 90's control scientists from the University of Lund (Sweden) realized that the above approach did not support component based modeling of physical systems with reuse 1. For instance, it was not easy to draw an electrical or hydraulic circuit by assembling component models of the various devices. The development of the Omola language by Hilding Elmqvist was a first attempt to bridge this gap by supporting some form of Differential Algebraic Equations (DAE) in the models. Modelica quickly emerged from this first attempt and became in the 2000's a major international concerted effort with the Modelica Consortium. A wider set of tools, both industrial and academic, now exists in this segment 2. In the Electronic Design Automation (EDA) sector, VHDL-AMS was developed as a standard 65 and also enables the use of differential algebraic equations. Several domain-specific languages and tools for mechanical systems or electronic circuits also support some restricted classes of differential algebraic equations. Spice is the historic and most striking instance of these domain-specific languages/tools 3. The main difference is that equations are hidden and the fixed structure of the differential algebraic results from the physical domain covered by these languages.
Despite the fact that these tools are now widely used by a number of engineers, they raise a number of technical difficulties. The meaning of some programs, their mathematical semantics, is indeed ambiguous. A main source of difficulty is the correct simulation of continuous-time dynamics, interacting with discrete-time dynamics: How the propagation of mode switchings should be handled? How to avoid artifacts due to the use of a global ODE solver causing unwanted coupling between seemingly non interacting subsystems? Also, the mixed use of an equational style for the continuous dynamics with an imperative style for the mode changes and resets, is a source of difficulty when handling parallel composition. It is therefore not uncommon that tools return complex warnings for programs with many different suggested hints for fixing them. Yet, these “pathological” programs can still be executed, if wanted so, giving surprising results — See for instance the Simulink examples in 30, 23 and 24.
Indeed this area suffers from the same difficulties that led to the development of the theory of synchronous languages as an effort to fix obscure compilation schemes for discrete time equation based languages in the 1980's. Our vision is that hybrid systems modeling tools deserve similar efforts in theory as synchronous languages did for the programming of embedded systems.
3.2 Background on non-standard analysis
Non-Standard analysis plays a central role in our research on hybrid systems modeling 23, 30, 25, 24, 27, 3. The following text provides a brief summary of this theory and gives some hints on its usefulness in the context of hybrid systems modeling. This presentation is based on our paper 2, a chapter of Simon Bliudze's PhD thesis 37, and a recent presentation of non-standard analysis, not axiomatic in style, due to the mathematician Lindström 90.
Non-standard numbers allowed us to reconsider the semantics of hybrid
systems and propose a radical alternative to the super-dense
time semantics developed by Edward Lee and his team as part of the
Ptolemy II project, where cascades of successive instants can occur in
zero time by using
Non-standard analysis was proposed by Abraham Robinson in the 1960s to allow the explicit manipulation of “infinitesimals” in analysis 84, 56, 52. Robinson's approach is axiomatic; he proposes adding three new axioms to the basic Zermelo-Fraenkel (ZFC) framework. While the need for non-standard analysis (in addition to the usual or standard analysis) has long agitated the mathematical community, it is not our purpose to debate such aspects. The important thing for us is that non-standard analysis allows the use of the non-standard discretization of continuous dynamics “as if” it was operational.
Not surprisingly, such an idea is not novel. Iwasaki et al. 67 first proposed using non-standard analysis to discuss the nature of time in hybrid systems. Bliudze and Krob 36, 37 have also used non-standard analysis as a mathematical support for defining a system theory for hybrid systems. They discuss in detail the notion of “system” and investigate computability issues. The formalization they propose closely follows that of Turing machines, with a memory tape and a control mechanism.
3.3 Structural Analysis of DAE Systems
The Modelica language is based on Differential Algebraic Equations (DAE). The general form of a DAE is given by:
where
Let
The state variables of
A leading variable
DAE are a strict generalization of ordinary differential
equations (ODE), in the sense that it may not be immediate
to rewrite a DAE as an explicit ODE of the form
For a square DAE of dimension
can locally be made explicit, i.e., the Jacobian matrix of
In practice, the problem of automatically finding a minimal
solution
A renowned method for the SA of DAE is the Pantelides method;
however, Pryce's
3.3.1 Pantelides method
In 1988, Pantelides proposed what is probably the most well-known SA method for DAE 81. The main idea of his work is that the structural representation of a DAE can be condensed into a bipartite graph whose left nodes (resp. right nodes) represent the equations (resp. the variables), and in which an edge exists if and only if the variable occurs in the equation.
By detecting specific subsets of the nodes, called Minimally Structurally Singular (MSS) subsets, the Pantelides method iteratively differentiates part of the equations until a perfect matching between the equations and the leading variables is found. One can easily prove that this is a necessary and sufficient condition for the structural nonsingularity of the system.
The main reason why the Pantelides method is not used in our work is that it cannot efficiently be adapted to multimode DAE (mDAE). As a matter of fact, the adjacency graph of a mDAE has both its nodes and edges parametrized by the subset of modes in which they are active; this, in turn, requires that a parametrized Pantelides method must branch every time no mode-independent MSS is found, ultimately resulting, in the worst case, in the enumeration of modes.
3.3.2 Pryce's Sigma-method
Albeit less renowned that the Pantelides method, Pryce's
This matrix is given by:
where
The primal problem consists in finding a maximum-weight perfect matching (MWPM) in the weighted adjacency graph. This is actually an assignment problem for which several standard algorithms exist, such as the push-relabel algorithm 63 or the Edmonds-Karp algorithm 58 to only give a few. However, none of these algorithms are easily parametrizable, even for applications to mDAE systems with a fixed number of variables.
The dual problem consists in finding the component-wise minimal
solution
- Initialize
to the zero vector. - For every
, - For every
, - Repeat Steps 2 and 3 until convergence is reached.
From the results proved by Pryce in 82, it is known
that the above algorithm terminates if and only if it is provided a
MWPM, and that the values it returns are independent of the choice of
a MWPM whenever there exist several such matchings. In particular, a
direct corollary is that the
Another important result is that, if the Pantelides method succeeds
for a given DAE
Working with this method is natural for our works, since the algorithm for solving the dual problem is easily parametrizable for dealing with multimode systems, as shown in our recent paper 43.
3.3.3 Block triangular decomposition
Once structural analysis has been performed, system
- the dependency graph of system
is generated, by taking into account the perfect matching between equations and unknowns ; - the strongly connected components (SCC) in this graph are determined: these will be the equation blocks that have to be solved;
- the block dependency graph is constructed as the
condensation of the dependency graph, from the knowledge of the SCC;
a BTF of system
can be made explicit from this graph.
3.4 Contract-Based Design, Interfaces Theories, and Requirements Engineering
System companies such as automotive and aeronautic companies are facing significant difficulties due to the exponentially raising complexity of their products coupled with increasingly tight demands on functionality, correctness, and time-to-market. The cost of being late to market or of imperfections in the products is staggering as witnessed by the recent recalls and delivery delays that many major car and airplane manufacturers had to bear in the recent years. The root causes of these design problems are complex and relate to a number of issues ranging from design processes and relationships with different departments of the same company and with suppliers, to incomplete requirement specification and testing.
We believe the most promising means to address the challenges in systems engineering is to employ formal design methodologies that seamlessly and coherently combine the various viewpoints of the design space (behavior, time, energy, reliability, ...), that provide the appropriate abstractions to manage the inherent complexity, and that can provide correct-by-construction implementations. The following issues must be addressed when developing new approaches to the design of complex systems:
- The overall design flows for heterogeneous systems and the associated use of models across traditional boundaries are not well developed and understood. Relationships between different teams inside a same company, or between different stake-holders in the supplier chain, are not supported by precise mathematical specifications of the components each party is expected to deliver.
- System requirements capture and analysis is in large part a heuristic process, where informal text and natural language-based techniques in use today are facing significant challenges 70. Formal requirements engineering is in its infancy: mathematical models, formal analysis techniques and links to system implementation must be developed.
- Dealing with variability, uncertainty, and life-cycle issues, such as extensibility of a product family, are not well-addressed using available systems engineering methodologies and tools.
The challenge is to address the entire process and not to consider only local solutions of methodology, tools, and models that ease part of the design.
Contract-based design has been proposed as a new approach to
the system design problem that is rigorous and effective in dealing
with the problems and challenges described before, and that, at the
same time, does not require a radical change in the way industrial
designers carry out their task as it cuts across design flows of
different types.
Indeed, contracts can be used almost everywhere and at nearly all
stages of system design, from early requirements capture, to embedded
computing infrastructure and detailed design involving circuits and
other hardware. Intuitively, a contract captures two properties,
respectively representing the assumptions on the environment and the
guarantees of the system under these assumptions. Hence, a
contract can be defined as a pair
- mathematical foundations of contracts, that enable the design of formal verification frameworks;
- System engineering methodologies and tools, that focus on requirements modeling, contract specification and verification, at multiple abstraction levels.
A detailed bibliography on contract and interface theories for embedded system design can be found in 5. In a nutshell, contract and interface theories fall into two main categories:
-
Assume/guarantee contracts.
By explicitly relying on the notions of assumptions and guarantees, A/G-contracts are intuitive. This makes them appealing for the engineer. In A/G-contracts, assumptions and guarantees are just properties regarding the behavior of a component and of its environment. The typical case is when these properties are formal languages or sets of traces. This includes the class of safety properties 71, 48, 77, 21, 50. Contract theories were initially developed as specification formalisms able to refuse some inputs from the environment 57. A/G-contracts were advocated in 26 and are is still a very active research topic, with several contributions dealing with the timed 34 and probabilistic 40, 41 viewpoints in system design, and even hybrid systems design 80.
-
Automata theoretic interfaces.
Interfaces combine assumptions and guarantees in a single, automata theoretic specification. Most interface theories are based on Lynch's Input/Output Automata 76, 75. Interface Automata 17, 16, 18, 46 focus primarily on parallel composition and compatibility: two interfaces are compatible if there exists at least one environment where they can work together. The idea is that the resulting composition exposes as an interface the needed information to ensure that incompatible pairs of states cannot be reached. This can be achieved by using the possibility, for an Interface Automaton, to refuse some inputs from the environment in a given state. This amounts to the implicit assumption that the environment will never produce any of the refused inputs, when the interface is in this state. Modal Interfaces 83 inherit from both Interface Automata and the originally unrelated notion of Modal Transition System 73, 20, 38, 72. Modal Interfaces are strictly more expressive than Interface Automata by decoupling the I/O orientation of an event and its deontic modalities (mandatory, allowed or forbidden). Informally, a must transition is offered in every component that realizes the modal interface, while a may transition is optional. Research on interface theories is still very active. For instance, timed 19, 31, 33, 54, 53, 32, probabilistic 40, 55 and energy-aware 47 interface theories have been proposed recently.
Requirements Engineering is one of the major concerns in large systems industries today, particularly so in sectors where certification prevails 85. Most requirements engineering tools offer a poor structuring of the requirements and cannot be considered as formal modeling frameworks today. They are nothing less, but nothing more than an informal structured documentation enriched with hyperlinks.
We see Contract-Based Design and Interfaces Theories as innovative tools in support of Requirements Engineering. The Software Engineering community has extensively covered several aspects of Requirements Engineering, in particular:
- the development and use of large and rich ontologies; and
- the use of Model Driven Engineering technology for the structural aspects of requirements and resulting hyperlinks (to tests, documentation, PLM, architecture, and so on).
Behavioral models and properties, however, are not properly encompassed by the above approaches. This is the cause of a remaining gap between this phase of systems design and later phases where formal model based methods involving behavior have become prevalent. We believe that our work on contract-based design and interface theories is best suited to bridge this gap.
3.5 Efficient Symbolic Computation for Sparse Systems
This project consists in exploiting the parsimony of sparse systems to accelerate their symbolic manipulation (quantifiers elimination 51, differential-algebraic reductions 86 etc.).
Let us cite two typical examples as a motivation: Boolean functions (
The current algorithms suffer from a theoretical complexity that is at best exponential (in the size of the input) limiting their use to instances of very modest size. The classic approach to overcome this problem is to develop/use numerical methods (with their limits and intrinsic problems) when possible of course. We aim to explore a different avenue.
In this project, we wish to exploit the structure of sparse systems to push the symbolic approach beyond its theoretical limits. The a priori limited application of our methods for dense systems is compensated by the fact that in practice, the problems are very often structured (in this regard, let us content ourselves with quoting the SAT solvers which successfully tackle industrial instances of a theoretically NP-complete problem).
The idea of exploiting the structure to speed up calculations that are a priori complex is not new. It has notably been developed and successfully used in signal processing via Factor Graphs 74, where one restricts oneself to local propagation of information, guided by an abstract graph which represents the structure of the system overall. Our approach is similar: we basically seek to use expensive algorithms sparingly on only subsystems involving only a small number of variables, thus hoping to reduce the theoretical worst case. One could then legitimately wonder why it is not enough to apply what has already been done on Factor Graphs? The difficulty (and the novelty for that matter) lies in the implementation of this idea for the problems that interest us. Let's start by emphasizing that the propagation of information has a significantly different impact depending on the operator (or quantifier) to be eliminated: a minimization or a summation do not look like a projection at all! This will obviously not prevent us from importing good ideas applicable to our problems and vice versa.
More related to symbolic computation, to our knowledge, at least two recent attempts exist: chordal networks 49 which propose a representation of the ideals of the ring of polynomials (therefore algebraic sets), and triangular block shapes 88, initiated independently and under development in our team and which tackle Boolean functions, or, if you will, the algebraic sets over the field of Booleans. The similarity between the two approaches is striking and suggests that there is a common way of doing things that could be exploited beyond these two examples. It is this unification that interests us in the first place in this project.
We identify three research problems to explore:
-
T1.
Unify several optimization problems on graphs as a single problem parameterized by a cost function, we coin such a problem WAP, for weighted adjacency propagation.
-
T2.
Adapt (and possibly improve) the algorithm of 87 to WAP and consequently to all instances of the single problem detailed in T1.
-
T3.
Propose a unified and modular method consisting of: (1) an elimination algorithm, (2) a data structure and (3) an efficient algorithm to solve the problem (with an adequate cost function).
The work on chordal networks and our work on Boolean functions immediately become special cases. For example, for Boolean functions, one could use Binary Decision Diagrams (BDDs) 39 to represent each piece of the initial system. In fact, the final representation will no longer be a single monolithic BDD as is currently the case, but rather a graph of BDDs. In the same way, an algebraic set will be represented by a graph where each node is a Gröbner basis (or any other data structure used to represent systems of equations).
The structure of the system becomes thus apparent and is exploited to optimize the used representation, opening the way to a better understanding and therefore to a more efficient and better targeted manipulation. Let's remember a simple fact here: symbolic manipulation often solves the problem exactly (without approximation or compromise). Therefore, pushing the limits of applicability of these techniques to scale them can only be appreciated and will undoubtedly have a significant impact on all the areas where they apply and the list is as long as it is varied. (compilation, certification, validation, synthesis, etc.).
4 Application domains
The Hycomes team contributes to the design of mathematical modeling languages and tools, to be used for the design of cyberphysical systems. In a nutshell, two major applications can be clearly identified: (i) our work on the structural analysis of multimode DAE systems has a sizeable impact on the techniques to be used in Modelica tools; (ii) our work on the verification of dynamical systems has an impact on the design methodology for safety-critical cyberphysical systems. These two applications are detailed below.
4.1 Modelica
Mathematical modeling tools are a considerable business, with major actors such as MathWorks, with Matlab/Simulink, or Wolfram, with Mathematica. However, none of these prominent tools are suitable for the engineering of large systems. The Modelica language has been designed with this objective in mind, making the best of the advantages of DAEs to support a component-based approach. Several industries in the energy sector have adopted Modelica as their main systems engineering language.
Although multimode features have been introduced in version 3.3 of the language 59, proper tool support of multimode models is still lagging behind. The reason is not a lack of interest from tool vendors and academia, but rather that multimode DAE systems poses several fundamental difficulties, such as a proper definition of a concept of solutions for multimode DAEs, how to handle mode switchings that trigger a change of system structure, or how impulsive variables should be handled. Our work on multimode DAEs focuses on these crucial issues 28.
Thanks to our IsamDAE software 43, 42, a larger class of Modelica models are expected to be compiled and simulated correctly. This should enable industrial users to have cleaner and simpler multimode Modelica models, with dynamically changing structure of cyberphysical systems. On the longer term, our ambition is to provide efficient code-generation techniques for the Modelica language, supporting, in full generality, multimode DAE systems, with dynamically changing differentiation index, structure and dimension.
The Hycomes team also focuses on scalability problems related to the compilation and simulation of large Modelica models. Digital twins developed by industrial Modelica users in the energy sector tend to be extremely large models, with up to
The Hycomes team is working on a new generation of algorithms for the compilation of the Modelica language, that can scale up to large models. The key contributations are modular index-reduction 4 and block-triangular equation sorting algorithms, that can be applied to incomplete (rectangular) DAE systems.
4.2 Dynamical Systems Verification
In addition to well-defined operational semantics for hybrid systems, one often needs to provide formal guarantees about the behavior of some critical components of the system, or at least its main underlying logic. To do so, we are actively developing new techniques to automatically verify whether a hybrid system complies with its specifications, and/or to infer automatically the envelope within which the system behaves safely. The approaches we developed have been already successfully used to formally verify the intricate logic of the ACAS X, a mid-air collision avoidance system that advises the pilot to go upward or downward to avoid a nearby airplane which requires mixing the continuous motion of the aircraft with the discrete decisions to resolve the potential conflict 68. This challenging example is nothing but an instance of the kind of systems we are targeting: autonomous smart systems that are designed to perform sophisticated tasks with an internal tricky logic. What is even more interesting perhaps is that such techniques can be often "reverted" to actually synthesize missing components so that some property holds, effectively helping the design of such complex systems.
5 Social and environmental responsibility
The expected impact of our research is to allow both better designs and more efficient exploitation of energy production units and distribution networks, enabling large-scale energy savings. At least, this is what we could observe in the context of the FUI ModeliScale collaborative project (2018–2021), focused on electric grids, urban heat networks, and building thermal modeling.
The rationale is as follows: system engineering models are meant to assess the correctness, safety and optimality of a system under design. However, system models are still useful after the system has been put in operation. This is especially true in the energy sector, where systems have an extremely long lifespan (for instance, more than 50 years for some nuclear power plants) and are upgraded periodically, to integrate new technologies. Exactly like in software engineering, where a software and its model co-evolve throughout the lifespan of the software, a co-evolution of the system and its physical models has to be maintained. This is required in order to maintain the safety of the system, but also its optimality.
Moreover, physical models can be instrumental to the optimal exploitation of a system. A typical example are model-predictive control (MPC) techniques, where the model is simulated, during the exploitation of the system, in order to predict system trajectories up to a bounded-time horizon. Optimal control inputs can then be computed by mathematical programming methods, possibly using multiple simulation results. This has been proved to be a practical solution 61, whenever classical optimal control methods are ineffective, for instance, when the system is non-linear or discontinuous. However, this requires the generation of high-performance simulation code, capable of simulating a system much faster than real-time.
The structural analysis techniques implemented in IsamDAE 43 generate a conditional block dependency graph, that can be used to generate high-performance simulation code : static code can be generated for each block of equations, and a scheduling of these blocks can be computed, at runtime, at each mode switching, thanks to an inexpensive topological sort algorithm. In contrast to other approaches (such as 60), no structural analysis, block-triangular decompositions, or automatic differentiation has to be performed at runtime.
6 Highlights of the year
- The design and implementation of an effective procedure to prove the non-existence of particular algebraic invariants, known as Darboux polynomials, in the literature. Such polynomials play an important role in the integrability of dynamical systems and proving their non-existence is often an error-prone painful task (if at all possible). This work was presented at ISSAC'24 10.
- A definition of the concept of solutions for linear piecewise affine differential-algebraic equations (PWA-DAEs). Unlike the conventional perspective that treats jumps as discrete-time dynamics, this work interprets jumps as continuous dynamics, parameterized by a virtual time variable. Moreover, by adapting the hybrid time-domain solution theory for continuous-discrete hybrid systems, we define the concept of jump-flow solutions, for PWA-DAEs leveraging Filippov solutions for differential inclusions. This work was presented at CDC'24 9
7 New software, platforms, open data
7.1 New software
7.1.1 IsamDAE
-
Name:
Implicit Structural Analysis of Multimode DAE systems
-
Keywords:
Structural analysis, Differential algebraic equations, Multimode, Scheduling, Consistent initialization, Code generation
-
Scientific Description:
Modeling languages and tools based on Differential Algebraic Equations (DAE) bring several specific issues that do not exist with modeling languages based on Ordinary Differential Equations. The main problem is the determination of the differentiation index and latent equations. Prior to generating simulation code and calling solvers, the compilation of a model requires a structural analysis step, which reduces the differentiation index to a level acceptable by numerical solvers.
The Modelica language, among others, allows hybrid models with multiple modes, mode-dependent dynamics and state-dependent mode switching. These Multimode DAE (mDAE) systems are much harder to deal with. The main difficulties are (i) the combinatorial explosion of the number of modes, and (ii) the correct handling of mode switchings.
The IsamDAE software aims at providing a compilation chain for mDAE-based modeling languages that make it possible to efficiently generate correct simulation code for multimode models. Novel structural analysis methods for mDAE systems were designed and implemented, based on an implicit representation of the varying structure of such systems. Several standard algorithms, such as J. Pryce's Sigma-method and the Dulmage-Mendelsohn decomposition, were adapted to the multimode case, using Binary Decision Diagrams (BDD) to represent the mode-dependent structure of an mDAE system.
IsamDAE determines, as a function of the mode, the set of latent equations, the leading variables and the state vector. This is then used to compute a conditional dependency graph (CDG) of the system, that can be used to generate simulation code with a mode-dependent scheduling of the blocks of equations. The software is also fit for generating simulation code for the hybrid dynamical system simulation tool Siconos, as well as handling the structural analysis of the multimode consistent initialization problem associated with an mDAE system.
-
Functional Description:
IsamDAE (Implicit Structural Analysis of Multimode DAE systems) is a software library implementing new structural analysis methods for multimode DAE systems, based on an implicit representation of incidence graphs, matchings between equations and variables, and block decompositions. The input of the software is a variable dimension multimode DAE system consisting in a set of guarded equations and guarded variable declarations. It computes a mode-dependent structural index reduction of the multimode system and is able to produce a mode-dependent graph for the scheduling of blocks of equations in long modes, check the structural nonsingularity of the associated consistent initialization problem, or generate simulation code for the nonsmooth dynamical system simulation tool Siconos.
IsamDAE is coded in OCaml, and uses the following packages: GuaCaml by Joan Thibault, MLBDD by Arlen Cox, Menhir by François Pottier and Yann Régis-Gianas, Pprint by François Pottier, Snowflake by Joan Thibault, XML-Light by Nicolas Cannasse and Jacques Garrigue.
-
Release Contributions:
New features:
* XML representations of the structure of a multimode DAE model are accepted as inputs by the IsamDAE tool, in order to enable weak coupling with tools based on existing DAE-based languages. IsamDAE distinguishes between MEL and XML inputs based on the extension of the input file (.mel versus .mdae.xml).
Bug fixes:
* A better handling of the model structure for consistent initialization prevents subtle bugs that were observed for a few models and initial events. Specific error messages are returned when initial equations involve variables that are not active in the corresponding modes.
Performance improvement:
* Better handling of sets of equations/variables labeled with propositional formulas, thanks to an adapted data structure.
Various:
* Verbosity option -v now takes as a parameter an integer ranging from 0 ("quiet") to 5 ("deep debug"). The detailed output of CoSTreD is only available in "deep debug" mode.
- URL:
- Publications:
-
Contact:
Benoit Caillaud
-
Participants:
Benoit Caillaud, Mathias Malandain, Joan Thibault, Alexandre Rocca, Bertrand Provot
7.1.2 modularSigma
-
Name:
A modular Sigma-method for the structural analysis of large DAE systems
-
Keywords:
Differential algebraic equations, Modularity
-
Scientific Description:
A key feature of the Modelica language is its object-oriented nature: components are instances of classes and they can aggregate other components, so that extremely large models can be efficiently designed as "trees of components". However, the structural analysis of Modelica models, a necessary step for generating simulation code, often relies on the flattening of this hierarchical structure, which undermines the scalability of the language and results in widely-used Modelica tools not being able to compile and simulate such large models. This software implements a new algorithm for the modular structural analysis of Modelica models. An adaptation of Pryce's Sigma-method for non-square DAE systems, along with a carefully crafted notion of component interface, make it possible to fully exploit the object tree structure of a model. The structural analysis of a component class can be performed once and for all, only requiring the information provided by the interface of its child components. The resulting method alleviates the exponential computation costs that can be yielded by model flattening, hence, its scalability makes it ideally suited for the modeling and simulation of large cyber-physical systems.
Algorithms implemented in modularSigma are based on the Sigma-method, which reduces the DAE structural index-reduction problem to two complementary linear programs: the primal problem amounts to the computation of a maximal-weight perfect matching of the equation-variable incidence graph of the DAE, while the dual problem consists in the computation of the minimal solution of a difference bound matrix (DBM). Modularity is achieved thanks to a decomposition of both problems, using dynamic programming principles (akin to message passing techniques, that are often used in statistical estimation) and memoization of the intermediate results.
-
Functional Description:
The software performs the index reduction and the bloc-triangular decomposition of large DAE systems, defined as the composition, hiding and renaming of incomplete (rectangular) DAE systems.
-
News of the Year:
The initial purpose of modularSigma has been to benchmark the algorithms detailed in the paper presented at the Modelica'23 conference.
- Publication:
-
Contact:
Benoit Caillaud
-
Participant:
Benoit Caillaud
7.1.3 PosInvSet
-
Name:
Positive Invariant Sets
-
Keywords:
Symbolic computation, Semi-algebraic set, Differential equations
-
Functional Description:
Given a semi-algebraic set S, that is a Boolean combination of equations and inequalities of polynomials, and a polynomial differential equation, we show that an algorithm can effectively decide whether S is a positive invariant set for the considered dynamic, that is, if the initial condition is in S, then the entire trajectory defined by the dynamics belongs to S.
We implemented in Mathematica two different procedures. Both require a backend algorithm for real quantifiers elimination (like the Cylindrical Algebraic Decomposition). One procedure form a monolithic request for the entire problem. The other chop the problem into small pieces following the Boolean structure of the input S.
-
Release Contributions:
Adaptation of the generic procedures to the linear case for scalability. The linear case means linear differential equations and semi-linear sets for the set S.
-
Contact:
Khalil Ghorbal
8 New results
8.1 A Modular Structural Analysis of DAE Systems
Participants: Albert Benveniste, Benoît Caillaud, Mathias Malandain.
In 29, a new modular structural analysis algorithm has been proposed that takes full advantage of the object tree structure of a DAE model. The bedrock of this method is a novel concept of structural analysis-aware interface for components. The essence of a component interface is to capture the necessary information about a Modelica class that needs to be exposed, in order to perform the structural analysis of a component comprising instances of the former class, while hiding away useless information regarding the equations and all protected features it may contain.
In order to compute a component interface, one has to be able to perform the structural analysis of the possibly non-square DAE system that this component encapsulates, and to use the interfaces of the components it aggregates in this analysis.
We base our algorithm on Pryce's
Putting all of this together, it is then possible to perform a modular structural analysis, in which structural analysis is performed at the class level, and the results can then be instantiated for each component of the system model, knowing its context. Hence, structural information at the system level is derived from composing the result of component-level analysis. Modular structural analysis yields huge gains in terms of memory usage and computational costs, as the analysis of a single large-scale DAE is replaced with that of multiple smaller subsystems. Moreover, the analysis is performed at the class level, meaning that a single structural analysis is needed for all system components that are instances of the same class.
In 2024, the modular structural analysis algorihtm has been fully implemented (see Section 7.1.2) and tested against several scalability benchmarks from the literature 45, 78. This experimental work confirms that the algorithm has an empirical complexity that is logarithmic in the size of the model, provided the model has a low tree-width, as it is often the case for energy network infrastructures.
8.2 Fault Diagnosability Analysis of Multi-Mode Systems
Participants: Benoît Caillaud, Mathias Malandain.
This work has been conducted in collaboration with the University of Linköping (Sweden) on the topic of system diagnosis, based on multimode DAE systems.
Fault detection and diagnosis are important for the health monitoring of physical systems. Model-based approaches for single-mode, smooth, systems are a well-established field, supported by a large body of literature covering various approaches like structural methods 35, parity space techniques, and observer-based methods 66.
While single-mode systems are often described using differential algebraic equations (DAEs), the modeling of non-smooth physical systems yields switched DAEs, also known as multimode DAEs (mDAEs), which combine continuous behaviors, defined as solutions of a set of DAE systems, with discrete mode changes 89, 28. Direct application of traditional fault diagnosis methods to all possible configurations of multi-mode systems quickly becomes intractable, as the number of modes tends to be exponential in the size of the system. The method proposed by 69 works around this issue by coupling a mode estimation algorithm with a single-mode diagnosis methodology, akin to just-in-time compilation in computer science. This approach unfortunately puts the burden on solving mode estimation problems, which often turn out to be intractable for the same reason.
Structural fault detectability and isolability are a graph-based method to evaluate diagnosability properties on DAEs 62. It is based on the Dulmage-Mendelsohn decomposition (DM), a building block of the structural analysis of equation systems. In 11, we show how its extension to multimode systems, introduced in 4, can be applied in the context of structural fault detectability and isolability of mmDAEs 64. Building upon our previous research studies, the methods presented in this paper represent advancements in diagnostic methodologies for multi-mode systems, providing novel ways to study the diagnosability of multi-mode systems without enumerating their modes.
The case study used throughout this article is a model of a reconfigurable battery system, in which switching strategies enable to produce an AC output without relying on a central inverter 22. This model is parameterized by the number of battery cells, so that both the inherent complexity associated with the diagnostics of such systems and the scalability of our approaches can be addressed.
8.3 Automated Reasoning For The Existence Of Darboux Polynomials
Participants: Maxime Bridoux, Khalil Ghorbal.
Darboux polynomials are particular algebraic invariants that play an important role in the integrability theory of ordinary differential equations (ODEs). Computation of Darboux polynomials is a central problem in the Prelle-Singer procedure for computing elementary first integrals of planar systems of polynomial ODEs [12], which yields a systematic method for computing elementary closed-form solutions (whenever these exist) to an important class of ordinary differential equations. Owing to this important application, algorithms for generating Darboux polynomials have received considerable attention in computer algebra. More recently, Darboux polynomials have found application in the area of formal safety verification of cyber-physical systems, where the problem of their automatic generation is encountered in the broader context of searching for invariant (and positively invariant) sets.
Darboux generation algorithms are semi-decision
procedures enumerating all Darboux polynomials up to a certain
fixed bound on the total degree. The bound is eventually increased
until finding a (not necessarily irreducible) Darboux polynomial or
reaching memory and/or time limits. Theoretically, the existence
of a bound on the total degree of irreducible Darboux polynomials
is, as of today, an open problem when
Given a polynomial ordinary differential equation (ODE), we devise generic polynomial reduction algorithms to automatically investigate the intertwined relationship between the total degree of (nontrivial) Darboux polynomials and the polynomials defining the ODE. By generic we mean that both the coefficients and the multidegree of the involved polynomials are symbolic. We use Newton polytopes as a light-weight abstraction to select optimal weight monomial orders improving the efficiency of the involved computations. The method works by inferring necessary conditions on both the coefficients and the multidegree for the polynomial to be Darboux. These conditions are then used, via constants' propagation, to restrict the shape of the generic candidate, pinpointing which monomials ought to be preserved by removing the superfluous ones. In some relevant cases, we are able to automatically prove the nonexistence of (nontrivial) Darboux polynomials providing a new toolbox to prove and formally certify that some limit cycles are not algebraic. 10
8.4 Solution concepts for linear piecewise affine or switched differential-algebraic equations
Participants: Yahao Chen.
The standard approach for dealing with (continuous) differential-algebraic equations (DAEs) is to solve the algebraic constraints—either manually or with the aid of computer algebra systems—in terms of a subset of the variables. These solutions are then substituted into the remaining differential equations, which results in a lower-dimensional ordinary differential equation (ODE). However, this approach does not account for the response to inconsistent initial conditions. Additionally, the presence of discontinuities complicates the global solution of the algebraic constraints, as different constraints may be active in different regions of the state-space. This leads to locally defined ODEs with different sets of variables, creating ambiguity in how solutions should be “glued” together when their variables differ. As a result, there is currently no theoretical foundation for studying discontinuous DAEs. Without such a foundation, subsequent studies—such as numerical simulations, stability analyses, and controller designs—lack a sound justification.
As a starting point for investigating discontinuous DAEs, we focus on linear piecewise affine systems (PWA-DAEs), described by
8.5 Pacti: Assume-Guarantee Contracts for Efficient Compositional Analysis and Design
Participants: Albert Benveniste, Benoît Caillaud.
Contract-based design is a method to facilitate modular design of systems. While there has been substantial progress in the theory of contracts, there has been less progress on practical algorithms for the algebraic operations in the theory. In this paper, we present (1) principles to implement a contract-based design tool at scale and (2) Pacti, a tool that can efficiently compute these operations. We illustrate the use of Pacti in a variety of case studies. 8
9 Partnerships and cooperations
9.1 International initiatives
9.1.1 Visits of international scientists
Participants: Khalil Ghorbal.
Andrew Sogokon, Lancaster University (Lecturer in the computer science department), visited our team in June 25-28 2024. Together with Khalil Ghorbal they actively worked on improving the scalability of formal verification techniques to prove the safety of linear dynamical systems towards handling invariants' candidates provided by the recent usage of neural networks for such dynamics. Funding Program iVisit 2024 IRISA.
Participants: Yahao Chen, Benoît Caillaud, Khalil Ghorbal, Albert Benveniste.
Stephan Trenn, University of Groningen (Associate Professor in the Faculty of Science and Engineering), visited our team in January 22-26 2024. A seminar was organized for presenting Trenn's works on the solution and stability analysis of differential-algebraic equations DAEs control systems. With Yahao Chen , he worked on the topic of linear piecewise affine DAEs and state-dependent switched DAEs.
9.2 European initiatives
9.2.1 Other european programs/initiatives
Participants: Benoît Caillaud, Mathias Malandain.
Since 2023 we have been developing an informal, but fruitful collaboration with Erik Frisk and Mattias Krysander (Linköping University, Sweden) on structural methods for the health monitoring of multimode systems. In 2024, we have contributed a decision procedure for the detectability and isolability of faults in multimode systems, that scales up to systems with large mode combinatorics 11 — see Section 8.2 for more details.
9.3 National initiatives
Participants: Benoît Caillaud, Mathias Malandain.
The Hycomes team is currently involved in the development of a digital twin of part of Line B of the Rennes subway. These works are carried out in collaboration with the I4S team (Inria, Rennes and Gustave Eiffel University, Nantes, France).
Focusing on a portion of the viaduct for line B, the considered subsystem consists of the concrete rollways, the electric supply and guiding rails, the defrosting/heating elements embedded in the rollways, and the electric control stations. Its digital twin is designed to predict the heating requirements for defrosting the outdoor subway track and send energy-optimal heating commands to the control stations. It is made of a finite element model of heat balance (diffusion, convection and radiation coupling) in the viaduct, interacting with a 3D environmental model (for drop shadows), measurements (real-time data from thermocouples and weather stations), weather forecasts, and a model of the heating control system. Physical modeling and software architecture are carried out by Benoît Caillaud and Mathias Malandain , as well as Jean Dumoulin and Thibaud Toullier (Gustave Eiffel University), with the involvement of Siemens Mobility (designer and manufacturer) and Keolis (operator).
The purpose of this digital twin, seen as a research artifact, is to be used as a testbed for model-based optimization and health monitoring techniques to be developed jointly between the I4S and Hycomes teams.
10 Dissemination
10.1 Promoting scientific activities
10.1.1 Scientific events: organisation
- Khalil Ghorbal is an organizer of 68NQRT, the local department seminar at Inria Rennes, inviting about 15 speakers per year.
- Benoît Caillaud is associate editor of Research Directions: Cyber-Physical Systems, an open-access journal published by Cambridge, University Press.
- Benoît Caillaud is also serving on the board of the MDPI Computation journal.
10.1.2 Invited talks
- Khalil Ghorbal was invited to the polysys seminar, LIP6, Paris. November 15th 2024. "The Quest For An Algebraic Characterization for Q-matrices".
- Khalil Ghorbal was invited to the Gallinette seminar, Nantes. November 28th 2024. "Automated reasoning for polynomial dynamical systems: characterizing, generating and formally verifying invariant sets".
- February 5-9, 2024. Yahao Chen participated to the E-Pico Winter School at Centrale Nantes: Presentation of one tutorial and Co-Chair for poster session.
- March 4-5, 2024. Yahao Chen participated to the Seminars in Technocentre Renault: Presentation of the works on torque estimations of electric moteur.
10.1.3 Scientific expertise
- Benoît Caillaud has evaluated proposals submitted for funding to the ANR (the French national research funding agency).
10.2 Teaching - Supervision - Juries
10.2.1 Teaching
- Master degree in computer science: Khalil Ghorbal , Category Theory, Monads, and Computation, at ENS Rennes, France.
- Master degree in computer science: Maxime Bridoux , Logic Courses (TD), at ENS Rennes, France.
- Licence degree in computer science: Maxime Bridoux , Imperative Programming (TP Java), at University of Rennes, France.
- Agregation informatique: Khalil Ghorbal and Maxime Bridoux , oral examination and lecture preparation, at ENS Rennes, France.
- Yahao Chen , Advanced control of electric propulsion systems at Central Nantes, France.
- Yahao Chen , Case study application dedicated to electric vehicle topology at Central Nantes, France.
10.2.2 Supervision
- Khalil Ghorbal is supervising the PhD work of Maxime Bridoux , on the broad topic of efficient symbolic computation methods for sparse algebraic systems. Funding Inria (AEx Backbone). Started September 2022. Expected to defend by September 2025.
- Khalil Ghorbal supervised the PhD work of Christelle Kozaily , on Linear Complementarity Systems. Started September 2018. Funding Inria (IPL ModeliScale). Christelle Kozaily defended her PhD thesis, entiteled on the existence of solutions for Linear Complementarity Problems, in October 27th 2024.
- Benoît Caillaud and Khalil Ghorbal supervised the PhD work of Joan Thibault on binary decision diagrams. Started September 2019. Funding University of Rennes. Joan Thibault defended his PhD thesis entitled Introduction to Structure Theory and its Application to Boolean Functions in December 19th 2024.
- Morgane Flamant carried out a 4-months internship (from May to August 2024), under the supervision of Benoît Caillaud and Mathias Malandain , contributing to the collaboration with the I4S project-team on digital twin engineering for structural health monitoring and energy efficiency. She developed numerical models that were used for testing and validating the software developed in the context of this joint project.
10.2.3 Specific official responsibilities in science outreach structures
Albert Benveniste is member of the French National Academy of Technology (he is vice-president at the Digital Node). As such he directed a think tank about Blockchain who published a report entitled “Blockchain, une technologie disruptive avec des enjeux de sûreté, résilience et impact environnemental”.
11 Scientific production
11.1 Major publications
- 1 articleBuilding a Hybrid Systems Modeler on Synchronous Languages Principles.Proceedings of the IEEE1069September 2018, 1568--1592HALDOI
- 2 articleNon-standard semantics of hybrid systems modelers.Journal of Computer and System Sciences783This work was supported by the SYNCHRONICS large scale initiative of INRIA2012, 877-910HALDOIback to text
- 3 articleThe mathematical foundations of physical systems modeling languages.Annual Reviews in Control502020, 72-118HALDOIback to text
- 4 articleAlgorithms for the Structural Analysis of Multimode Modelica Models.Electronics1117September 2022, 1-63HALDOIback to textback to text
- 5 articleContracts for System Design.Foundations and Trends in Electronic Design Automation122-32018, 124-400HALDOIback to text
- 6 articleA Formally Verified Hybrid System for Safe Advisories in the Next-Generation Airborne Collision Avoidance System.International Journal on Software Tools for Technology Transfer196November 2017, 717-741HALDOI
- 7 articleOperational Models for Piecewise-Smooth Systems.ACM Transactions on Embedded Computing Systems (TECS)165sOctober 2017, 185:1--185:19HALDOI
11.2 Publications of the year
International journals
- 8 articlePacti: Assume-Guarantee Contracts for Efficient Compositional Analysis and Design.ACM Transactions on Cyber-Physical SystemsNovember 2024, 1-33HALDOIback to text
International peer-reviewed conferences
- 9 inproceedingsSolution concepts for linear piecewise affine differential-algebraic equations.CDC 2024 - 63rd IEEE Conference on Decision and ControlMilan, Italy2024, 1-6HALback to textback to text
- 10 inproceedingsAutomated Reasoning For The Existence Of Darboux Polynomials.ISSAC 2024 - International Symposium on Symbolic and Algebraic ComputationRaleigh, NC, United StatesACMJuly 2024, 324-333HALDOIback to textback to text
- 11 inproceedingsFault Diagnosability Analysis of Multi-Mode Systems.SAFEPROCESS 2024 - 12th IFAC Symposium on Fault Detection, Supervision and Safety for Technical Processes584Ferrara, ItalyElsevier2024, 210-215HALDOIback to textback to text
Reports & preprints
- 12 miscA Mathematica Package for Certifying the Nonexistence of Darboux Polynomials.July 2024HAL
- 13 miscStability and stabilization of state-dependent switched linear differential-algebraic equations.2024HAL
- 14 miscOn Covering Euclidean Spaces with Q-arrangements of Cones.2024HAL
- 15 miscStrong Structural Controllability Analysis of Structured Networks with Identical Nodes.2024HAL
11.3 Cited publications
- 16 inproceedingsGame Models for Open Systems.Verification: Theory and Practice2772Lecture Notes in Computer ScienceSpringer2003, 269-289DOIback to text
- 17 inproceedingsInterface automata.Proc. of the 9th ACM SIGSOFT International Symposium on Foundations of Software Engineering (FSE'01)ACM Press2001, 109--120DOIback to text
- 18 inproceedingsInterface-based design.In Engineering Theories of Software Intensive Systems, proceedings of the Marktoberdorf Summer SchoolKluwer2004DOIback to text
- 19 inproceedingsTimed Interfaces.Proc. of the 2nd International Workshop on Embedded Software (EMSOFT'02)2491Lecture Notes in Computer ScienceSpringer2002, 108--122DOIback to text
- 20 article20 Years of Modal and Mixed Specifications.Bulletin of European Association of Theoretical Computer Science1942008, URL: https://dblp.org/rec/journals/eatcs/AntonikHLNW08.bibback to text
- 21 bookPrinciples of Model Checking.MIT Press, Cambridge2008, URL: https://mitpress.mit.edu/9780262026499/principles-of-model-checking/back to text
- 22 inproceedingsDesign and Analysis of Battery-Integrated Modular Multilevel Converters for Automotive Powertrain Applications.2021 23rd European Conference on Power Electronics and Applications (EPE'21 ECCE Europe)IEEE2021, P--1back to text
- 23 articleBuilding a hybrid systems modeler on synchronous languages principles.Proceedings of the IEEE1069September 2018, 1568--1592HALDOIback to textback to text
- 24 miscA Type-Based Analysis of Causality Loops In Hybrid Systems Modelers.Deliverable D3.1_1 v 1.0 of the Sys2soft collaborative project ''Physics Aware Software''December 2013HALback to textback to text
- 25 miscSemantics of multi-mode DAE systems.Deliverable D.4.1.1 of the ITEA2 Modrio collaborative projectAugust 2013HALback to text
- 26 inproceedingsMultiple viewpoint contract-based specification and design.Proceedings of the Software Technology Concertation on Formal Methods for Components and Objects (FMCO'07)5382Revised Lectures, Lecture Notes in Computer ScienceAmsterdam, The NetherlandsSpringerOctober 2008DOIback to text
- 27 techreportStructural Analysis of Multimode DAE Systems: summary of results.RR-9387Inria Rennes -- Bretagne AtlantiqueJanuary 2021, 27HALback to text
- 28 articleThe mathematical foundations of physical systems modeling languages.Annual Reviews in Control502020, 72-118HALDOIback to textback to text
- 29 inproceedingsTowards the separate compilation of Modelica: modularity and interfaces for the index reduction of incomplete DAE systems.Linköping Electronic Conference Proceedings204Aachen, GermanyOctober 2023, 10HALDOIback to text
- 30 inproceedingsA type-based analysis of causality loops in hybrid modelers.HSCC '14: International Conference on Hybrid Systems: Computation and ControlProceedings of the 17th international conference on Hybrid systems: computation and control (HSCC '14)Berlin, GermanyACM PressApril 2014, 13HALDOIback to textback to text
- 31 inproceedingsA Compositional Approach on Modal Specifications for Timed Systems..11th International Conference on Formal Engineering Methods (ICFEM'09)5885LNCSRio de Janeiro, BrazilSpringerDecember 2009, 679-697HALback to text
- 32 articleModal event-clock specifications for timed component-based design.Science of Computer Programming772012, 1212-1234HALDOIback to text
- 33 inproceedingsRefinement and Consistency of Timed Modal Specifications..3rd International Conference on Language and Automata Theory and Applications (LATA'09)5457LNCSTarragona, SpainSpringerApril 2009, 152-163HALDOIback to text
- 34 inproceedingsA proposal for real-time interfaces in SPEEDS.Design, Automation and Test in Europe (DATE'10)IEEE2010, 441-446DOIback to text
- 35 inbookDiagnosis and Fault-Tolerant Control.Springer Berlin, Heidelberg09 2006, 109-188DOIback to text
- 36 articleModelling of complex systems: Systems as dataflow machines.Fundamenta Informaticae9122009, 251-274HALDOIback to text
- 37 phdthesisUn cadre formel pour l'étude des systèmes industriels complexes: un exemple basé sur l'infrastructure de l'UMTS.Ecole Polytechnique2006back to textback to text
- 38 articleGraphical versus logical specifications.Theoretical Computer Science10611992, 3-20URL: https://www.sciencedirect.com/science/article/pii/030439759290276LDOIback to text
- 39 articleGraph-Based Algorithms for Boolean Function Manipulation.IEEE Trans. Comput.358August 1986, 677--691URL: http://dx.doi.org/10.1109/TC.1986.1676819DOIback to text
- 40 inproceedingsCompositional design methodology with constraint Markov chains.QEST 2010Williamsburg, Virginia, United StatesSeptember 2010, URL: http://hal.inria.fr/inria-00591578/enDOIback to textback to text
- 41 articleConstraint Markov Chains.Theoretical Computer Science41234May 2011, 4373-4404URL: http://hal.inria.fr/hal-00654003/enDOIback to text
- 42 miscDemo: IsamDAE, an Implicit Structural Analysis Tool for Multimode DAE Systems.PosterApril 2020, 1HALback to text
- 43 inproceedingsImplicit structural analysis of multimode DAE systems.HSCC 2020 - 23rd ACM International Conference on Hybrid Systems: Computation and ControlSydney New South Wales Australia, FranceACMApril 2020, 1-11HALDOIback to textback to textback to text
- 44 articleThe index of general nonlinear DAEs.Numerische Mathematik722dec 1995, 173--196URL: http://dx.doi.org/10.1007/s002110050165DOIback to text
- 45 inproceedingsScalableTestGrids - An Open-Source and Flexible Benchmark Suite to Assess Modelica Tool Performance on Large-Scale Power System Test Cases.Proceedings of the 14th International Modelica ConferenceLinköping Electronic Conference Proceedings181Linköping, SwedenModelica Association and Linköping University Electronic PressSeptember 2021, 351--358DOIback to text
- 46 phdthesisA Framework for Compositional Design and Analysis of Systems.EECS Department, University of California, BerkeleyDec 2007, URL: http://www.eecs.berkeley.edu/Pubs/TechRpts/2007/EECS-2007-174.htmlback to text
- 47 inproceedingsResource Interfaces.Embedded Software, Third International Conference, EMSOFT 2003, Philadelphia, PA, USA, October 13-15, 2003, Proceedings2855Lecture Notes in Computer ScienceSpringer2003, 117-133DOIback to text
- 48 inproceedingsCharacterization of temporal property classes.ICALP623Lecture Notes in Computer ScienceSpringer1992, 474-486DOIback to text
- 49 articleChordal Networks of Polynomial Ideals.SIAM J. Appl. Algebra Geom.112017, 73--110URL: https://doi.org/10.1137/16M106995XDOIback to text
- 50 bookModel Checking.MIT Press1999, URL: https://mitpress.mit.edu/9780262038836/model-checking/back to text
- 51 articlePartial Cylindrical Algebraic Decomposition for Quantifier Elimination.J. Symb. Comput.1231991, 299--328DOIback to text
- 52 bookN. J.N. J. Cutland, eds. Nonstandard Analysis and its Applications.Cambridge Univ. Press1988DOIback to text
- 53 inproceedingsECDAR: An Environment for Compositional Design and Analysis of Real Time Systems.Automated Technology for Verification and Analysis - 8th International Symposium, ATVA 2010, Singapore, September 21-24, 2010. Proceedings2010, 365-370DOIback to text
- 54 inproceedingsTimed I/O automata: a complete specification theory for real-time systems.Proceedings of the 13th ACM International Conference on Hybrid Systems: Computation and Control, HSCC 2010, Stockholm, Sweden, April 12-15, 20102010, 91-100DOIback to text
- 55 inproceedingsAbstract Probabilistic Automata.Verification, Model Checking, and Abstract Interpretation - 12th International Conference, VMCAI 2011, Austin, TX, USA, January 23-25, 2011. Proceedings6538Lecture Notes in Computer Science2011, 324-339DOIback to text
- 56 bookAnalyse non standard.Hermann1989, URL: https://www.editions-hermann.fr/livre/analyse-non-standard-francine-dienerback to text
- 57 bookTrace Theory for Automatic Hierarchical Verification of Speed-Independent Circuits.ACM Distinguished DissertationsMIT Press1989DOIback to text
- 58 articleTheoretical improvements in algorithmic efficiency for network flow problems.Journal of the ACM1921972, 248--264URL: http://dx.doi.org/10.1145/321694.321699DOIback to text
- 59 inproceedingsModelica extensions for Multi-Mode DAE Systems.Proceedings of the 10th International Modelica Conference, March 10-12, 2014, Lund, SwedenLinköping University Electronic Pressmar 2014DOIback to text
- 60 inproceedingsModia-dynamic modeling and simulation with julia.Juliacon'18University College London, UKAugust 2018, URL: https://elib.dlr.de/124133/back to text
- 61 inproceedingsSurvey of industrial applications of embedded model predictive control.2016 European Control Conference (ECC)2016, 601-601DOIback to text
- 62 articleDiagnosability analysis considering causal interpretations for differential constraints.IEEE Transactions on Systems, Man, and Cybernetics-Part A: Systems and Humans4252012, 1216--1229back to text
- 63 inproceedingsA new approach to the maximum flow problem.Proceedings of the eighteenth annual ACM symposium on Theory of computing (STOC'86)1986, URL: http://dx.doi.org/10.1145/12130.12144DOIback to text
- 64 inproceedingsHierarchical Diagnosis Algorithm for Component-Based Multi-Mode Systems.22nd IFAC World CongressIFAC2023back to text
- 65 miscIEEE Standard VHDL Analog and Mixed-Signal Extensions, Std 1076.1-1999.1999, URL: http://dx.doi.org/10.1109/IEEESTD.1999.90578DOIback to text
- 66 bookFault-Diagnosis Systems: An Introduction from Fault Detection to Fault Tolerance.Springer Berlin Heidelberg2006back to text
- 67 inproceedingsModeling time in hybrid Systems: How fast is ``instantaneous''?IJCAI1995, 1773--1781URL: https://www.ijcai.org/Proceedings/95-2/Papers/097.pdfback to text
- 68 inproceedingsFormal verification of ACAS X, an industrial airborne collision avoidance system.2015 International Conference on Embedded Software, EMSOFT 2015, Amsterdam, Netherlands, October 4-9, 2015Amsterdam, NetherlandsIEEE2015, 127--136HALDOIback to text
- 69 articleStructural fault detection and isolation in hybrid systems.IEEE Transactions on Automation Science and Engineering1542017, 1585--1599back to text
- 70 phdthesisPrincipe de transduction sémantique pour l'application de théories d'interfaces sur des documents de spécification.Université Rennes 1 ; Rennes 1April 2021HALback to text
- 71 articleProving the correctness of multiprocess programs.IEEE Trans. Software Eng.321977, 125--143DOIback to text
- 72 inproceedingsOn Modal Refinement and Consistency.Proc. of the 18th International Conference on Concurrency Theory (CONCUR'07)Springer2007, 105--119DOIback to text
- 73 inproceedingsA Modal Process Logic.Proceedings of the Third Annual Symposium on Logic in Computer Science (LICS'88)IEEE1988, 203-210DOIback to text
- 74 articleAn introduction to factor graphs.IEEE Signal Processing Magazine2112004, 28-41DOIback to text
- 75 inproceedingsInput/Output Automata: Basic, Timed, Hybrid, Probabilistic and Dynamic.CONCUR 2003 - Concurrency Theory, 14th International Conference, Marseille, France, September 3-5, 2003, Proceedings2761Lecture Notes in Computer ScienceSpringer2003, 187-188DOIback to text
- 76 articleA Proof of the Kahn Principle for Input/Output Automata.Inf. Comput.8211989, 81-92DOIback to text
- 77 bookTemporal verification of reactive systems: Safety.Springer1995DOIback to text
- 78 inproceedingsAutomated model generation and simplification for district heating and cooling networks.Proceedings of the 13th International Modelica ConferenceLinköping Electronic Conference Proceedings157Regensburg, GermanyModelica Association and Linköping University Electronic PressMarch 2019, 179--186DOIback to text
- 79 articleApplying ``Design by Contract''.Computer2510October 1992, 40--51URL: http://dx.doi.org/10.1109/2.161279DOIback to text
- 80 articleMethodology for the Design of Analog Integrated Interfaces Using Contracts.IEEE Sensors Journal1212Dec. 2012, 3329--3345DOIback to text
- 81 articleThe consistent initialization of differential-algebraic systems.SIAM J. Sci. Stat. Comput.921988, 213--231DOIback to textback to text
- 82 articleA Simple Structural Analysis Method for DAEs.BIT Numerical Mathematics412March 2001, 364--394URL: http://dx.doi.org/10.1023/a:1021998624799DOIback to textback to textback to text
- 83 articleA Modal Interface Theory for Component-based Design.Fundamenta Informaticae1081-22011, 119-149URL: http://hal.inria.fr/inria-00554283/enDOIback to text
- 84 bookNon-Standard Analysis.Princeton Landmarks in Mathematics1996, URL: https://press.princeton.edu/books/paperback/9780691044903/non-standard-analysisback to text
- 85 articleIndustry needs and research directions in requirements engineering for embedded systems.Requirements Engineering172012, 57--78URL: http://link.springer.com/article/10.1007/s00766-011-0144-xDOIback to text
- 86 inproceedingsThe Ritt–Kolchin theory for differential polynomials.Differential Algebra and Related Topics2002, 1-70DOIback to text
- 87 articlePositive-instance driven dynamic programming for treewidth.J. Comb. Optim.3742019, 1283--1311DOIback to text
- 88 articleLeveraging Structural Analysis for Quantified Boolean Formulae.Summer School on Modelling and Verification of Parallel Processes, Grenoble, France6http://khalilghorbal.info/assets/pdf/papers/RBTF_movep.pdf2020back to text
- 89 inbookSwitched Differential Algebraic Equations.Dynamics and Control of Switched Electronic Systems: Advanced Perspectives for Modeling, Simulation and Control of Power ConvertersF.Francesco Vasca and L.Luigi Iannelli, eds. LondonSpringer London2012, 189--216DOIback to textback to text
- 90 incollectionAn invitation to nonstandard analysis.Nonstandard Analysis and its ApplicationsCambridge Univ. Press1988, 1--105DOIback to text