2023Activity reportProject-TeamHEPHAISTOS

7.1.1 ALIAS

• Name:
  Algorithms Library of Interval Analysis for Systems
• Keyword:
  Interval analysis
• Functional Description:

  The ALIAS library whose development started in 1998, is a collection of procedures based on interval analysis for systems solving and optimization.

  ALIAS is made of two parts:

  ALIAS-C++ : the C++ library (87 000 code lines) which is the core of the algorithms

  ALIAS-Maple : the Maple interface for ALIAS-C++ (55 000 code lines). This interface allows one to specify a solving problem within Maple and get the results within the same Maple session. The role of this interface is not only to generate the C++ code automatically, but also to perform an analysis of the problem in order to improve the efficiency of the solver. Furthermore, a distributed implementation of the algorithms is available directly within the interface.

• URL:
  http://­www-sop.­inria.­fr/­hephaistos/­developpements/­main.­html
• Contact:
  Jean-Pierre Merlet
• Participants:
  Jean-Pierre Merlet, Odile Pourtallier

7.1.2 CAVATOI

• Name:
  Smart digital contact-book
• Keywords:
  Web Application, Home care, Data analysis, Machine learning, Anomaly detection, Data visualization
• Functional Description:
  This piece of software is a web application designed for collection, analysis and visualisation of activities data of an elderly people based on the testimonies of helpers/visitors. It includes functionalities for detecting and analysing slight changes, and based on AI algorithms deciding wether or not alerting authorised family or helpers.
• Release Contributions:
  Fully functional proof of concept
• News of the Year:
  Initially developed at Inria Start-up Studio as the technological part of a transfer process, this software will be used and enhanced in our scientific activities.
• Authors:
  Pierre Pigeon, Yves Papegay
• Contact:
  Yves Papegay

Cable-driven parallel robot

A completely new version of our old modular MARIONET-CRANE cable-driven parallel robot prototype is now installed in the robotic hall, for experiments in the field of walking assistance and health monitoring for frail people.

This installation takes benefit and improves the two previous ones, installed for art performances in exhibition center in Amilly (45) and in Mouans-Sartoux (06) - see below 7.3.

Software environment for solving parametric equations with IA

Our work in robotics has led us to investigate the use of IA for finding the real root(s) of parametric non linear square system of equations i.e. systems which have as many unknowns as equations but with equation coefficients that are functions of parameters that are assumed to be bounded. Such a system has usually several solutions and their number depends on the parameters values and cannot be determined in advance. An example of such a system is presented in section 8.1.1 together with the principle of the method. Just as an outline the method requires to have determined the solutions (not necessarily all of them) for a small set $𝒮$ of specific systems. We are then able to create from $𝒮$ different learning sets i.e. a set of samples that are described by a pair ($P_i,S_i$) where $P_i$ is a valued parameters vector and $S_i$ is a solution of the system for the given $P_i$. These learning sets are obtained by using a structure rule based on the concept of linear aspect: two samples ($P_1,S_1$), ($P_2,S_2$) belong to the same aspect if when following a linear path between $P_1$ and $P_2$ we don't cross a singular system and the solution path leads from $S_1$ to $S_2$. These learning sets are used to create neural networks using a specific training process based not only on the decrease of a loss function but also on an index that quantifies how many exact solutions can be obtained over all samples of the training set. Creating these networks is relatively computer intensive but then getting the exact solutions from the networks predictions for a specific system is very fast. Consequently using this approach is appropriate when a sufficient number of system instances have to be solved, this number depending upon the computation time of alternate methods. Furthermore the process have numerous steps that can benefit from a distributed implementation both for creating the networks but also when solving a given system and we are currently investigating the use of low power consumption AI processors to speed up the calculation.

We cannot guarantee to obtain all solutions (although our extensive tests have shown that in general we will miss very few solutions) but we have developed a self-learning process that may allow to reduce the number of missed solution(s) as soon as new system instances are solved . The method is generic and we are currently developing a software framework that takes as input $𝒮$ and produces an AI-based solver with minimal manual intervention.

8.1.1 Kinematics and AI

As mentioned last year we have started to work on the direct kinematics (DK) of cable-driven parallel robot (CDPR) having sagging cables (i.e. being elastic and having their own mass). The direct kinematics consists in determining the pose(s) of the platform for given cable length $L_0$. Sagging cable may be modeled by the Irvine textbook planar model where the upper attachment point $A$ of the cable is supposed to be grounded and be the origin of a planar frame with horizontal and vertical axis. The model provides the coordinates of the lowest attachment point $B$ of the cable as functions of the cable length $L_0$ at rest and of the horizontal and vertical components $F_x,F_z$ of the force applied at $B$:

where $E$ is the Young modulus of the cable material (which characterizes its elasticity), $\mu$ its linear density and $A_0$ the area of a cross-section of the cable that are all supposed to be known. Note that if $x_b,z_b$ and $L_0$ are known the 2 unknowns $F_x,F_z$ of these equations may be found using neural networks 17.

We have proposed last year a preliminary version of a DK solver using AI and we have finalized this year a very efficient DK solver 16.

The pose of a CDPR is defined by a vector $𝐗$ which has 6 components: three coordinates of the center of mass of the platform and 3 angles which describes the platform orientation. For a CDPR with $n$ cables we have $2n+6$ unknowns (the $2n$$F_x,F_z$ and the components of $𝐗$). As the $x_b,z_b$ of each cable can be expressed as functions of $𝐗$ we have $2n$ Irvine equations. In a static case we must ensure the mechanical equilibrium of the platform, which provide 6 equations in $𝐗$ and in the $F_x,F_z$. Hence solving the DK consists in solving a square system $ℱ$ of $2n+6$ equations, that has usually several solutions whose number cannot be determined in advance. For solving the DK we assume that we have been able to determine solutions of $ℱ$ (not necessarily all of them) for a few set of $n{L}_0$. Hence we have a set of $(L_0^j,S_1^j,S_2^j,...S_m^j)$ where $S_k^j$ are solutions of $ℱ$ for the set of lengths $L_0^j$. These data can be represented graphically: an horizontal line, called level, represents a set of lengths $L_0^j$ and points, called nodes, on this line represents the $S_k^j$. A node is there represented by a pair ($L_0^j,S_k^j$). Consider two nodes at different levels ($L_0^j,S_k^j$), ($L_0^u,S_l^u$). We define a linear path in the $L_0$ space as $L_0=L_0^j+\lambda (L_0^u-L_0^j)$ and we use a continuation scheme with $S_k^j$ as starting point and $\lambda$ as continuation parameter for determining a solution of $ℱ$ for $L_0=L_0^u$. Starting from $\lambda =0,S=S_k^j$ we incrementally increase $\lambda$ and use Newton with $S$ as guess to determine a solution $S$ of $ℱ$. This scheme stops if either $\lambda =1$ or if Newton encounters a singularity. Note that if $\lambda =1$, then $S$ may be a solution that was unknown for $L_0^u$.

If two nodes ($L_0^j,S_k^j$), ($L_0^u,S_l^u$) are connected with this scheme we establish an edge between them, allowing to establish a directed graph called the solution graph. In this graph we may identify the non redundant circuits i.e. a list of successively connected nodes that contains at most one node at a given level. In robotics terms these circuits define a path in the $L_0$ space and an associated singularity-free path for the platform that allow to connect the pose of the first node of the circuit to the pose of the last one and all poses of the nodes of the circuit belong to the same variety called an aspect. Poses belonging to the same aspect have the property that any two poses in the same aspect can be joined by following a singularity-free path in the $L_0$ space. The non redundant circuits will be used to obtain learning sets for the neural networks that we will create: when creating the edges with the continuation scheme we may store sets of samples ($L_0,S$) using storage rules such that

• the maximal change between the current $L_0$ and the previously stored sample is larger than a given threshold
• the distance between the current location of $G$ and the previously stored one is larger than a given threshold

If not enough samples have been obtained for a given circuit we may add samples just by selecting a stored sample ($L_0^u,S^u$), choosing an arbitrary unit vector $𝐯$ in the $L_0$ space and applying a continuation scheme $L_0=L_0^u+\lambda 𝐯$ starting at $\lambda =0$, solution $S^u$ while storing samples along this path by using the same storage rules.

Regarding the neural networks we use multi-layer perceptron (MLP) with a specific training procedure. It is partly based on a decrease of the MSE loss but when a new MLP is obtained we run a verification procedure to obtain the success rate of the MLP. For calculating this rate the verification procedure runs the MLP on each sample of the training set and use the solution prediction of the MLP as guess for a Newton scheme. The success rate is then calculated as the percentage of samples such that the verification procedure leads to the expected solution of $ℱ$ (e.g if the success rate is 100, then the verification procedure obtains the expected solution for all samples in the training set). Clearly if during the training we find a MLP that has a 100% success rate, then it is not necessary to continue the training and furthermore given the size of our problem the training time is low (a few minutes). Consequently as we have no a-priori knowledge on the number of layers, number of neurons in the layers and activation function that may lead to a MLP that has a 100% success rate we may perform a systematic approach by testing all combinations of MLP that have from 2 to 10 layers, between 10 to 200 neurons with a step size of 20 and activation function among a set of 3 classical one (ReLU, LeakyReLU, CELU). We also use an adaptive learning rate usually starting at ${10}^{-3}$ and decreasing by 2 the learning rate if after 15 optimizer step the loss has not decreased and stopping the learning if the learning rate is lower than ${10}^{-6}$. For a given training set we stop the systematic MLP exploration as soon as a MLP with 100% success rate is obtained. If a 100% success rate is not reached we retain the MLP having the highest success rate and we consider in sequence the sample(s) whose solution has not been found and use them as starting point for creating a new learning set using random unit vectors in the $L_0$ space and a continuation process to obtain new samples. This set is used to train a new MLP that will allow to find the missed solution and possibly other missed solution(s). New MLPs are then added until all solutions of the training set are found. After having processed all the training sets we obtain a set of MLPs such that the the MLP prediction used as a guess for Newton (a procedure called local solver) will provide the solution for all the samples of the training sets.

For example for a CDPR having 8 cables (and consequently the DK has 22 equations) we have used 11 initial sets of $L_0$ with between 8 and 24 DK solutions to establish 1133 MLPs that cover all solutions of the training sets. However some of these MLPs are redundant (they provide solutions that will also be discovered by other MLPs) and it is interesting to reduce the number of necessary MLPs while still getting all solutions. For that purpose let $𝒰$ be the set of all samples over all training sets. We run in sequence all MLPs on $𝒰$ and determine the MLP that provides the largest number of solutions ${𝒱}_1$ over $𝒰$ and store it in a list. We repeat this procedure for the remaining MLPs over the set of samples in $𝒰-{𝒱}_1$ until this list of samples is empty. In our example the final number of MLPs is 138 and we are able to establish a DK solver by using the local solver for each MLP.

We may then assess the quality of this DK solver on verification sets that provide all DK solutions for a given set of $L_0$. These sets are established as the training sets by selecting random unit vectors in the $L_0$ space but at the opposite of the training set the continuation scheme follows all solutions from the starting point, which implies that the number of DK solutions for each sample is the same than the number of solution of the starting point. We also ensure that the $L_0$ of the samples in the verification set are all different from the one used in the training sets. The result are very good with about 99.9% of the samples for which all solutions are found. But two cases may occur:

• the DK solver find solution(s) for some samples that were not present in the verification set: we just update the verification set. Note that this update may also be used to update the solution graph by discovering new nodes and possibly new non redundant circuits which will provide new MLPs
• the DK solver miss solution(s) for some samples: we just use the missed solution to create new learning sets and their corresponding MLPs

After this update all solutions are found for all samples of the verification set. We may then repeat the procedure with a new verification set. However we cannot guarantee that the current DK solver will find all solutions in all case. For example we have found examples where we have assessed 10 verification sets with 100% success rate but for the 11-th one solution has been missed, imposing the creation of a new MLP.

The computation time for this DK solver may be decomposed into the time for establishing the necessary MLPs and the exploitation time for solving one instance of the DK. In a sequential process about 100 hours are necessary to establish the solver while solving one instance of the DK takes about 1 second (to be compared to the 20+ hours that are required when using interval analysis or continuation). However these times may be improved as there are numerous steps for establishing the MLPs that may benefit from a parallel implementation while the local solvers of DK solver may be run in parallel. This year we have investigated the use of the IA processor JETSON Nano but the result was somewhat disappointing because of the limited floating point power of the GPU. Anyhow this AI solver is intended to be used when numerous DK instances have to be solved and here as soon as at least 5 DK instances have to be solved the AI solver is faster than the alternate methods.

Finally as mentioned in section 7.2 the principle of the DK solver is generic and may be used for any parametric square equations system.

8.1.2 Preventive maintenance and AI

The kinematics equations governing the behavior of CDPR with sagging cables have been presented presented in section 8.1.1. They are used for control purpose as they establish the relationship between the actuated variables (the cable lengths $L_0$) and the pose of the CDPR platform under the assumptions that the $L_0$ can be measured accurately. Special winches with cable guide have been designed for that purpose but have the drawback of restricting the CDPR workspace and of not being completely reliable as the guide is working only under the assumption that there is a sufficient tension in the cable, a condition that cannot always be enforced. At the opposite simple multi-layered winch offers a compact solution while providing a very large workspace but with the drawback that the cable length deduced from the drum rotation is more approximate. We have proposed to have redundant sensing i.e. to complement the cable length measurements with additional sensors:

• having planar lidars that allow to estimate the location of the platform in the 3D space
• having angular sensors on the cables, located at a known distance from the attachment point of the cables on the platform, that measure the slope of the cable at the sensor location. It can be shown that the cable slope at this point is a function of $F_x,F_z,L_0,\mu ,A_0$ and $E$.

These two sensing methods have been implemented on one of our prototype and have shown to be very effective. Here we investigate if this sensing redundancy may also allow to assess the wear of the cables which will decrease the $E$ of the cables with the objective of improving the accuracy of the CDPR by using more accurate $E$ and of preventive maintenance e.g. changing cables before their breakdown.

When controlling the cable lengths to reach a given pose ${𝐗}_{𝐢}$ we assume that the $E$ of all cables are identical and equal to $E_i$. If this assumption is wrong the CDPR will reach a pose ${𝐗}_{𝐫}\ne {𝐗}_{𝐢}$ and the lidars will allow to calculate the distance $Q_d$ between ${𝐗}_{𝐫},{𝐗}_{𝐢}$ while the angular sensors will measure a change $Q_a$ compared to the inclination of the cable at ${𝐗}_{𝐢}$. These measurements have an error margin and therefore can be exploited if the changes are larger than given thresholds $\Delta Q_d,\Delta Q_a$.

First we have to perform a sensitivity analysis i.e to determine how much change in the $E$ is required to obtain $Q_d\ge \Delta Q_d$ and/or $Q_a\ge \Delta Q_a$. Indeed if measurable changes in the $Q$ are obtained only when the $E$ becomes very low our objectives cannot be reached. Clearly the changes in $Q_a,Q_d$ are dependent upon the pose and upon the way the $E$ of the cables are changing. We have therefore selected a set of representative cable lengths $L_0$ and solved the DK for these lengths which provides the different possible poses of the platform. We have then defined a variation law for the $E$ which uses a unit vector $𝐯$ in the $E$ space (with only negative of null components) such that the vector of the $E$ is equal to $𝐄={𝐄}_{𝐢}+\lambda 𝐯$. Consequently $Q_d$ and the $Q_a$ of the cables are functions of $\lambda$ and a continuation scheme may be used to follow the curves $Q_d\left(\lambda \right),Q_a\left(\lambda \right)$ for determining the first ${\lambda }_d$ such that $Q_d\left({\lambda }_d\right)\ge \Delta Q_d$ and ${\lambda }_a$ such $Q_a\left({\lambda }_a\right)\ge \Delta Q_a$. Clearly the ${\lambda }_a,{\lambda }_d$ are functions of the vector $𝐯$: we have decided to test the wear of a single cable (one component of $𝐯$ equal to -1, the others to 0) and random $𝐯$. We have also discarded cases where the ${\lambda }_d,{\lambda }_a$ correspond to less than 30% of the $E_i$. The analysis is not yet complete but we have found cases where we have measurable changes in $Q$ for moderate changes in the $E$ so that it appears that changes in the $E$ are observable.

After considering the sensitivity we will have to solve the inverse problem, i.e. determining the current value of all the $E$ of the cables. We will have to determine how many $Q$ are required to safely determine these values. Fortunately change in the $E$ is a slow process so that this estimation may be based on records of the $Q$ obtained during trajectories. Still the inverse problem implies solving a set of non-linear equations and our approach will be based on AI. Indeed we may obtain training set in simulation and we believe that, as for the DK, we will have to structure the data and to train several networks for obtaining reliable estimations.

8.1.3 Green robotics

A large number of industrial robots are used mostly only for simple pick-and-place operation, taking an object at a given position and moving it to another one. In most cases this task is performed by serial robot which are not energy efficient (beside the load the robot actuation has to move part of the robot own weight) while using a computer and various electronic boards that are largely under exploited which constitutes a waste of energy and resources. We intend to investigate the use for this task of cable-driven parallel robot, which are 25% more energy efficient, and to design an electro-mechanical mechanism for executing pick-and place trajectories while getting rid of the computer and of most electronic resources.

8.2.1 Human activity recognition

Participants: Jean-Pierre Merlet, Yves Papegay, Odile Pourtallier, Eric Wajnberg.

Human activity recognition (HAR) is a major topic in the team. We are focusing on monitoring mobility and displacements (we are not yet interested in recognizing the individual action of our subject) using a sensor-based approach, excluding vision which is intrusive and even prohibited in some places for legal reasons. For that purpose we have in the previous years developed a smart barrier combining redundant passive infrared motion detectors and infrared distance sensors. Smart barriers have been implemented in Ehpad Valrose, a new retirement house in which a specific infrastructure has been put in place to accommodate research works and in Institut Claude Pompidou, an Alzheimer day care hospital from 2019 to 2020. These two long term experiments have allowed us to determine that essential points in HAR are to determine what is possible to measure, the sensor types, how to retrieve and process sensor data, how effective are the quality of measurements on a long term basis and the level of monitoring that is acceptable for frail peoples and their helpers while providing significant and reliable data for the medical community in spite of the uncertainties both in the measurements and in the system modeling. These samples of questions will become central in our work.

We are interested in and currently investigating the use of another type of barriers, based on Lidar's. Measurements produced by lidars are richers and allow a better reliability on the number and the directions of simultaneous crossings.

8.3.1 Optimized flower visiting strategy of bees

Three years ago, through an international scientific cooperation with Israeli scientists located at the Ben-Gurion University of the Negev, and as explained in an activity report provided previously (2020), we developed a probabilistic model whose aim was to understand a strange - and up to now not understandable - reproductive behavior of the potter wasp insect Delta dimidiatipenne (Hymenoptera). Females of this insect lay their eggs in mud chambers provisioned with caterpillars they capture to feed their young. Experimental observations indicate that many of the caterpillars collected by this wasp are actually parasitized by a small gregarious parasitoid wasp and are providing a lower amount of food for the wasp progeny. As a result, all players in the interaction perish - the young potter wasp cannot fully exploit the caterpillars and presumably starve to death; and the small parasitoids complete their development, but cannot break out of the mud and remain trapped in the sealed pot. Following such observation, we developed a probabilistic model trying to understand under what environmental conditions such a striking phenomenon (i.e., bringing back to the nest parasitized caterpillars), can be maintained, despite the high cost to all players. After several problems to develop this modeling work (some of them were due to the current political situation in Israel), we were able to submit a manuscript to a good international journal (Behavioural Processes). By the end September, we received a decision (minor revision) and a revised version of the text was resubmitted by mid-November. We are now waiting on the editorial decision from the journal.

8.3.2 Other research activities developed this year

Several other activities were developed this (and the previous) years, all lead to publications in international journals:

• A work 12 has been published with Odile Pourtallier , originating from an international scientific cooperation with Israeli scientists located at the University of Haifa to understand the optimal foraging decision of bees foraging for nectar. For this, we developed an optimization deterministic model trying to understand what should be the optimized flower visiting strategy, taking into account the ability of the foraging animals to learn the quality of the different flowers they are visiting.
• It is difficult to optimize the treatment of crops with insecticides again pests, since these chemical compounds also kill the natural enemies of the pests that are present, hence hampering our ability to protect the crop efficiently. With colleagues from different countries (USA, Israel, Finland and India), we developed a theoretical model to optimize the use of pesticides to control crop pests, taking also into account the incomes for the farmers. 13
• With an international team, we also published a review articles in a top-level international journal 14, discussing how a trait-based approach can be used – mainly based on international databases – to identify which natural enemy species can be used to develop efficient biological control programs against crop pests. Such an approach poses several problems that remain to be solved, and we proposed different way to do this.
• Finally, with Brazilian colleagues (from the University of São Paulo), we published 15 a theoretical model to optimize the used of genetically modified corn plants to control corn pests in the field. The model is based on an individual-based spatially explicit Monte Carlo simulation.

9.1.1 Visits to international teams

9.2.1 Other european programs/initiatives

• Hephaistos is part of the euROBIN, the Network of Excellence on AI and robotics that was launched in 2021

10.1.1 Scientific events: organisation

10.1.2 Journal

10.1.3 Invited talks

• Jean-Pierre Merlet has given a talk about AI and robotics and one about interval analysis at the joint ACCENTAURI-HEPHAISTOS seminar at Inria.
• Jean-Pierre Merlet presented the HEPHAISTOS activities on handicap during INRIA Scientific days, during the Handicap week and to the "Conseil de l'Âge".
• Jean-Pierre Merlet talked about the past of robotics activities during the 40th years of INRIA Sophia-Antipolis Méditerrannée research center.
• Eric Wajnberg has given an invited talk entitled "A modeling approach to improve the efficiency of augmentative biological control with arthropod natural enemies" during the Second International Latin-American Conference on Biological Control (Juazeiro, Brazil).
• Eric Wajnberg has given an invited talk entitled "Stochastic Individual-Based Model (IBM) in parasitoid ecology and Genetic Algorithm" during the Third Symposium on ecology, evolution and diversity (Federal University of São Paulo, Brazil).

10.1.4 Leadership within the scientific community

• Jean-Pierre Merlet is a member of the IFToMM (International Federation for the Promotion of Mechanism and Machine Science) technical Committees on History and on Computational Kinematics. He is a member of the IFToMM Executive Council Publication Advisory Board and an IEEE Fellow.
• Jean-Pierre Merlet has got an Emeritus research director position at INRIA and an emeritus senior chair of 3IA Côte d’Azur since December 2023.
• Jean-Pierre Merlet was a member of the scientific committee of the CNRS GDR robotique until September 2023 and is now member of the advisory "Conseil des Sages".

10.1.5 Scientific expertise

• Jean-Pierre Merlet is a Nominator for the Japan's Prize. He is the head of the robotics GDR Publication Committee that has produced the 2023 report on "recommended supports for publication" for journals and conferences that does not provide a ranking but advice according to robotics topics.
• Jean-Pierre Merlet is active in the "Robotics in Provence" industrial cluster.
• Yves Papegay is a member of the OpenMath Society, building an extensible standard for representing the semantics of mathematical objects.
• Eric Wajnberg is an appointed member of the Academic Committee of the Hebrew University of Jerusalem, since June 2022, for four years. He has been an appointed member of the International Advisory Board of the “International Center for Excellence in Biological Control”, from August 2018 to August 2023.

10.1.6 Research administration

• Yves Papegay is the head of local CUMI (Committe of users of the numerical resources and tools).
• Jean-Pierre Merlet has been an elected member of INRIA Scientific Committee until mid-January.
• Jean-Pierre Merlet has been a corresponding member of INRIA ethical committee (COERLE) and a member of the Ethical Committee of Université Côte d'Azur (CER) until September.
• Odile Pourtallier has been member of the local researcher recruitment jury.
• Odile Pourtallier is a member of the local Transform committee.

10.2.1 Teaching

• Jean-Pierre Merlet has taught 15 hours on parallel robots to Master ISC (M2) at University of Toulon.
• Yves Papegay has taught the "Robotics and AI" course of L3IA at Université Côte d'Azur.
• Yves Papegay has taught the "Oriented-Object Programming" Course of L2Info at University of French Polynesia.

10.2.2 Supervision

• Jean-Pierre Merlet is supervisor of the PhD of Romain Tissot . Together with Yves Papegay , he is supervising the PhD of Clara Thomas .
• Jean-Pierre Merlet has been the supervisor of the internship of Julien Bondyfalat whose purpose was to investigate the use of an IA parallel processor for interval analysis algorithms and deep learning for problems with low training time but requiring multiple neural networks and training sets.
• Yves Papegay is supervising an industrial project of engineering school students, about design and realisation of a modular and multi-scale cable-driven parallel cranes.

10.2.3 Juries

• Jean-Pierre Merlet has been member of the PhD defense committte of M. Metillon at LS2N Nantes,
• Eric Wajnberg has been member of 10 PhD defense committees at the University of Palermo (Barbaccia P., Consentino B.B., Ponte M., Alinc T., Auteri N., Bambina P., Borgia D., Funsten C.C., Giuga L., Ippolito M., Roma E., Sofia S.).

10.3.1 Education

• Jean-Pierre Merlet and Yves Papegay are active members of the dissemination initiative Terra Numerica.

10.3.2 Interventions

• Jean-Pierre Merlet has given a general talk about ethics for Unica HR4S program.
• Eric Wajnberg has given two public talks in the scope of Science pour Tous 06, one about "Pourquoi la sexualité ? Le regard de la biologie" in the city of Roussillon, the other one about "Lutter contre les ravageurs de culture - De l'usage des pesticides vers des alternatives plus respectueuses de l'environnement" in Aspremont.

International journals

• 11 articleS.Sébastien Briot and J.-P.Jean-Pierre Merlet. Direct Kinematic Singularities and Stability Analysis of Sagging Cable-driven Parallel Robots.IEEE Transactions on Robotics393June 2023, 2240-2254HAL
• 12 articleT.Tamar Keasar, O.Odile Pourtallier and E.Eric Wajnberg. Can sociality facilitate learning of complex tasks? Lessons from bees and flowers.Philosophical Transactions of the Royal Society B: Biological Sciences37820210402January 2023HALDOIback to text
• 13 articleT.Tamar Keasar, E.Eric Wajnberg, G.George Heimpel, I.Ian Hardy, L. S.Liora Shaltiel Harpaz, D.Daphna Gottlieb and S.Saskya van Nouhuys. Dynamic Economic Thresholds for Insecticide Applications Against Agricultural Pests: Importance of Pest and Natural Enemy Migration.Journal of Economic Entomology1162April 2023, 321-330HALDOIback to text
• 14 articleM.Michal Segoli, P.Paul Abram, J.Jacintha Ellers, G.Gili Greenbaum, I. C.Ian C.W. Hardy, G.George Heimpel, T.Tamar Keasar, P.Paul Ode, A.Asaf Sadeh and E.Eric Wajnberg. Trait-based approaches to predicting biological control success: challenges and prospects.Trends in Ecology and Evolution389September 2023, 802-811HALDOIback to text
• 15 articleM.Maysa Tomé, I.Igor Weber, A.Adriano Garcia, J.Josemeri. Jamielniak, E.Eric Wajnberg, M.Mirian Hay-Roe and W.Wesley Godoy. Modeling fall armyworm resistance in Bt-maize areas during crop and off-seasons.Journal of Pest Science964September 2023, 1539-1550HALDOIback to text

International peer-reviewed conferences

• 16 inproceedingsJ.-P.Jean-Pierre Merlet. Advances in the use of neural network for solving the direct kinematics of CDPR with sagging cables.CabeleCon - 6th International conference on calble-driven parallel robotsNantes, FranceSpringerJune 2023HALDOIback to text
• 17 inproceedingsJ.-P.Jean-Pierre Merlet. Irvine cable equations and neural networks.IFToMM WC 2023 - 16th World Congress of the International Federation for the Promotion of Mechanism and Machine ScienceTokyo, JapanNovember 2023HALback to text
• 18 inproceedingsJ.-P.J-P Merlet and Y.Y Papegay. The new exhibition Blind machines, a large 3D printing machine.2023 IEEE International Conference on Robotics and Automation (ICRA)ICRA 2023 - International Conference on Robotics and AutomationLondres (Grande Bretagne), United KingdomMay 2023HALDOI

Scientific books

• 19 bookI. C.Ian C.W. Hardy and E.Eric Wajnberg. Jervis's Insects as Natural Enemies: Practical Perspectives.Springer International Publishing2023, XVIII, 770 p.HALDOI

Scientific book chapters

• 20 inbookM.Mark Fellowes, J.Jacques van Alphen, K.K. Shameer, I.Ian Hardy, E.Eric Wajnberg and M.Mark Jervis. Foraging Behaviour.Jervis's Insects as Natural Enemies: Practical PerspectivesSpringer International PublishingNovember 2023, 1-104HALDOI