-
A Unified View of Polarity for Functions
Authors:
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
We propose a unified view of the polarity of functions, that encompasses all specific definitions, generalizes several well-known properties and provides new results. We show that bipolar sets and bipolar functions are isomorphic lattices. Also, we explore three possible notions of polar subdifferential associated with a nonnegative function, and we make the connection with the notion of alignemen…
▽ More
We propose a unified view of the polarity of functions, that encompasses all specific definitions, generalizes several well-known properties and provides new results. We show that bipolar sets and bipolar functions are isomorphic lattices. Also, we explore three possible notions of polar subdifferential associated with a nonnegative function, and we make the connection with the notion of alignement of vectors.
△ Less
Submitted 22 October, 2024;
originally announced October 2024.
-
A Two-Timescale Decision-Hazard-Decision Formulation for Storage Usage Values Calculation
Authors:
Camila Martinez Parra,
Michel de Lara,
Jean-Philippe Chancelier,
Pierre Carpentier,
Jean-Marc Janin,
Manuel Ruiz
Abstract:
The penetration of renewable energies requires additional storages to deal with intermittency. Accordingly, there is growing interest in evaluating the opportunity cost (usage value) associated with stored energy in large storages, a cost obtained by solving a multistage stochastic optimization problem. Today, to compute usage values under uncertainties, an adequacy resource problem is solved usin…
▽ More
The penetration of renewable energies requires additional storages to deal with intermittency. Accordingly, there is growing interest in evaluating the opportunity cost (usage value) associated with stored energy in large storages, a cost obtained by solving a multistage stochastic optimization problem. Today, to compute usage values under uncertainties, an adequacy resource problem is solved using stochastic dynamic programming assuming a hazard-decision information structure. This modelling assumes complete knowledge of the coming week uncertainties, which is not adapted to the system operation as the intermittency occurs at smaller timescale. We equip the twotimescale problem with a new information structure considering planning and recourse decisions: decision-hazard-decision. This structure is used to decompose the multistage decision-making process into a nonanticipative planning step in which the on/off decisions for the thermal units are made, and a recourse step in which the power modulation decisions are made once the uncertainties have been disclosed. In a numerical case, we illustrate how usage values are sensitive as how the disclosure of information is modelled.
△ Less
Submitted 30 August, 2024;
originally announced August 2024.
-
Multistage stochastic optimization of a mono-site hydrogen infrastructure by decomposition techniques
Authors:
Raian Lefgoum,
Sezin Afsar,
Pierre Carpentier,
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
The development of hydrogen infrastructures requires to reduce their costs. In this paper, we develop a multistage stochastic optimization model for the management of a hydrogen infrastructure which consists of an electrolyser, a compressor and a storage to serve a transportation demand. This infrastructure is powered by three different sources: on-site photovoltaic panels (PV), renewable ene…
▽ More
The development of hydrogen infrastructures requires to reduce their costs. In this paper, we develop a multistage stochastic optimization model for the management of a hydrogen infrastructure which consists of an electrolyser, a compressor and a storage to serve a transportation demand. This infrastructure is powered by three different sources: on-site photovoltaic panels (PV), renewable energy through a power purchase agreement (PPA) and the power grid. We consider uncertainties affecting on-site photovoltaic production and hydrogen demand. Renewable energy sources are emphasized in the hydrogen production process to ensure eligibility for a subsidy, which is awarded if the proportion of nonrenewable electricity usage stays under a predetermined threshold. We solve the multistage stochastic optimization problem using a decomposition method based on Lagrange duality. The numerical results indicate that the solution to this problem, formulated as a policy, achieves a small duality gap, thus proving the effectiveness of this approach.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Learning with Fitzpatrick Losses
Authors:
Seta Rakotomandimby,
Jean-Philippe Chancelier,
Michel de Lara,
Mathieu Blondel
Abstract:
Fenchel-Young losses are a family of convex loss functions, encompassing the squared, logistic and sparsemax losses, among others. Each Fenchel-Young loss is implicitly associated with a link function, for mapping model outputs to predictions. For instance, the logistic loss is associated with the soft argmax link function. Can we build new loss functions associated with the same link function as…
▽ More
Fenchel-Young losses are a family of convex loss functions, encompassing the squared, logistic and sparsemax losses, among others. Each Fenchel-Young loss is implicitly associated with a link function, for mapping model outputs to predictions. For instance, the logistic loss is associated with the soft argmax link function. Can we build new loss functions associated with the same link function as Fenchel-Young losses? In this paper, we introduce Fitzpatrick losses, a new family of convex loss functions based on the Fitzpatrick function. A well-known theoretical tool in maximal monotone operator theory, the Fitzpatrick function naturally leads to a refined Fenchel-Young inequality, making Fitzpatrick losses tighter than Fenchel-Young losses, while maintaining the same link function for prediction. As an example, we introduce the Fitzpatrick logistic loss and the Fitzpatrick sparsemax loss, counterparts of the logistic and the sparsemax losses. This yields two new tighter losses associated with the soft argmax and the sparse argmax, two of the most ubiquitous output layers used in machine learning. We study in details the properties of Fitzpatrick losses and in particular, we show that they can be seen as Fenchel-Young losses using a modified, target-dependent generating function. We demonstrate the effectiveness of Fitzpatrick losses for label proportion estimation.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Orbital perturbation coupling of primary oblateness and solar radiation pressure
Authors:
Martin Lara,
Elena Fantino,
Roberto Flores
Abstract:
Solar radiation pressure can have a substantial long-term effect on the orbits of high area-to-mass ratio spacecraft, such as solar sails. We present a study of the coupling between radiation pressure and the gravitational perturbation due to polar flattening. Removing the short-period terms via perturbation theory yields a time-dependent two-degree-of-freedom Hamiltonian, depending on one physica…
▽ More
Solar radiation pressure can have a substantial long-term effect on the orbits of high area-to-mass ratio spacecraft, such as solar sails. We present a study of the coupling between radiation pressure and the gravitational perturbation due to polar flattening. Removing the short-period terms via perturbation theory yields a time-dependent two-degree-of-freedom Hamiltonian, depending on one physical and one dynamical parameter. While the reduced model is non-integrable in general, assuming coplanar orbits (i.e., both Spacecraft and Sun on the equator) results in an integrable invariant manifold. We discuss the qualitative features of the coplanar dynamics, and find three regions of the parameters space characterized by different regimes of the reduced flow. For each regime, we identify the fixed points and their character. The fixed points represent frozen orbits, configurations for which the long-term perturbations cancel out to the order of the theory. They are advantageous from the point of view of station keeping, allowing the orbit to be maintained with minimal propellant consumption. We complement existing studies of the coplanar dynamics with a more rigorous treatment, deriving the generating function of the canonical transformation that underpins the use of averaged equations. Furthermore, we obtain an analytical expression for the bifurcation lines that separate the regions with different qualitative flow.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Higher-order composition of short- and long-period effects for improving analytical ephemeris computation
Authors:
Martin Lara,
Elena Fantino,
Hadi Susanto,
Roberto Flores
Abstract:
The construction of an analytic orbit theory that takes into account the main effects of the Geopotential is notably simplified when splitting the removal of periodic effects in several stages. Conversely, this splitting of the analytical solution into several transformations reduces the evaluation efficiency for dense ephemeris output. However, the advantage is twofold when the different parts of…
▽ More
The construction of an analytic orbit theory that takes into account the main effects of the Geopotential is notably simplified when splitting the removal of periodic effects in several stages. Conversely, this splitting of the analytical solution into several transformations reduces the evaluation efficiency for dense ephemeris output. However, the advantage is twofold when the different parts of the mean-to-osculating transformation are composed into a single transformation. To show that, Brouwer's solution is extended to the second order of the zonal harmonic of the second degree by the sequential elimination of short- and long-period terms. Then, the generating functions of the different transformations are composed into a single one, from which a single mean-to-osculating transformation is derived. The new, unique transformation notably speeds up the evaluation process, commonly improving evaluation efficiency by at least one third with respect to the customary decomposition of the analytical solution into three different parts.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
Room-temperature polariton repulsion and ultra-strong coupling for a non-trivial topological one-dimensional tunable Fibonacci-conjugated porous-Silicon photonic quasi-crystal showing quasi bound-states-in-the-continuum
Authors:
Atzin David Ruiz Pérez,
Salvador Escobar Guerrero,
Ma. Del Rocío Nava Lara,
Jorge-Alejandro Reyes-Esqueda
Abstract:
Room temperature strong coupling from CdSeS/Zn quantum-dots embedded into a tunable porous-silicon Fibonacci-conjugated array could be observed when exciton's energy was tuned either to the photonic-edge or the defect in the middle of the pseudo-bandgap region of the 1D cavity. Both, the photonic-edge and the defect could be identified as topological edge modes and quasi-bound-states-in-the-contin…
▽ More
Room temperature strong coupling from CdSeS/Zn quantum-dots embedded into a tunable porous-silicon Fibonacci-conjugated array could be observed when exciton's energy was tuned either to the photonic-edge or the defect in the middle of the pseudo-bandgap region of the 1D cavity. Both, the photonic-edge and the defect could be identified as topological edge modes and quasi-bound-states-in-the-continuum, where large density of states and field localization over a wider bandwidth produce a broadband Purcell enhancement, helping to optimize the coupling among the exciton and the 1D photonic quasi-crystal despite the natural difficulty to make the quantum dots to penetrate the cavity pores. A clear repulsion among polaritons, amounting to almost 8 meV for in-plane k values when the cavity energy is larger than the exciton one (blue k-detuning), was measured when increasing the incident light fluence, marking the potential of this non-trivial topological array for achieving polariton quantum blockade. Evidence for ultra-strong coupling, where a shift as large as 20 meV, could be found when the defect of the pseudo-bandgap region of the cavity was tuned to the exciton.
△ Less
Submitted 29 May, 2023; v1 submitted 17 May, 2023;
originally announced May 2023.
-
On mean elements in artificial satellite theory
Authors:
Martin Lara
Abstract:
The merits of a perturbation theory based on a mean to osculating transformation that is pure periodic in the fast angle are investigated. The exact separation of the purely short-period effects of the perturbed Keplerian dynamics from the long-period mean frequencies is achieved by a non-canonical transformation, which, therefore, cannot be computed by Hamiltonian methods. For this case, the evol…
▽ More
The merits of a perturbation theory based on a mean to osculating transformation that is pure periodic in the fast angle are investigated. The exact separation of the purely short-period effects of the perturbed Keplerian dynamics from the long-period mean frequencies is achieved by a non-canonical transformation, which, therefore, cannot be computed by Hamiltonian methods. For this case, the evolution of the mean elements strictly adheres to the average behavior of the osculating orbit. However, due to the inescapable truncation of perturbation solutions, the fact that this theory confines the long-period variations of the semimajor axis into the mean variation equations, how tiny they may be, can have adverse effects in the accuracy of long-term semi-analytic propagations based on it
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
Towards unraveling calibration biases in medical image analysis
Authors:
María Agustina Ricci Lara,
Candelaria Mosquera,
Enzo Ferrante,
Rodrigo Echeveste
Abstract:
In recent years the development of artificial intelligence (AI) systems for automated medical image analysis has gained enormous momentum. At the same time, a large body of work has shown that AI systems can systematically and unfairly discriminate against certain populations in various application scenarios. These two facts have motivated the emergence of algorithmic fairness studies in this fiel…
▽ More
In recent years the development of artificial intelligence (AI) systems for automated medical image analysis has gained enormous momentum. At the same time, a large body of work has shown that AI systems can systematically and unfairly discriminate against certain populations in various application scenarios. These two facts have motivated the emergence of algorithmic fairness studies in this field. Most research on healthcare algorithmic fairness to date has focused on the assessment of biases in terms of classical discrimination metrics such as AUC and accuracy. Potential biases in terms of model calibration, however, have only recently begun to be evaluated. This is especially important when working with clinical decision support systems, as predictive uncertainty is key for health professionals to optimally evaluate and combine multiple sources of information. In this work we study discrimination and calibration biases in models trained for automatic detection of malignant dermatological conditions from skin lesions images. Importantly, we show how several typically employed calibration metrics are systematically biased with respect to sample sizes, and how this can lead to erroneous fairness analysis if not taken into consideration. This is of particular relevance to fairness studies, where data imbalance results in drastic sample size differences between demographic sub-groups, which, if not taken into account, can act as confounders.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
Jupiter Science Enabled by ESA's Jupiter Icy Moons Explorer
Authors:
Leigh N. Fletcher,
Thibault Cavalié,
Davide Grassi,
Ricardo Hueso,
Luisa M. Lara,
Yohai Kaspi,
Eli Galanti,
Thomas K. Greathouse,
Philippa M. Molyneux,
Marina Galand,
Claire Vallat,
Olivier Witasse,
Rosario Lorente,
Paul Hartogh,
François Poulet,
Yves Langevin,
Pasquale Palumbo,
G. Randall Gladstone,
Kurt D. Retherford,
Michele K. Dougherty,
Jan-Erik Wahlund,
Stas Barabash,
Luciano Iess,
Lorenzo Bruzzone,
Hauke Hussmann
, et al. (25 additional authors not shown)
Abstract:
ESA's Jupiter Icy Moons Explorer (JUICE) will provide a detailed investigation of the Jovian system in the 2030s, combining a suite of state-of-the-art instruments with an orbital tour tailored to maximise observing opportunities. We review the Jupiter science enabled by the JUICE mission, building on the legacy of discoveries from the Galileo, Cassini, and Juno missions, alongside ground- and spa…
▽ More
ESA's Jupiter Icy Moons Explorer (JUICE) will provide a detailed investigation of the Jovian system in the 2030s, combining a suite of state-of-the-art instruments with an orbital tour tailored to maximise observing opportunities. We review the Jupiter science enabled by the JUICE mission, building on the legacy of discoveries from the Galileo, Cassini, and Juno missions, alongside ground- and space-based observatories. We focus on remote sensing of the climate, meteorology, and chemistry of the atmosphere and auroras from the cloud-forming weather layer, through the upper troposphere, into the stratosphere and ionosphere. The Jupiter orbital tour provides a wealth of opportunities for atmospheric and auroral science: global perspectives with its near-equatorial and inclined phases, sampling all phase angles from dayside to nightside, and investigating phenomena evolving on timescales from minutes to months. The remote sensing payload spans far-UV spectroscopy (50-210 nm), visible imaging (340-1080 nm), visible/near-infrared spectroscopy (0.49-5.56 $μ$m), and sub-millimetre sounding (near 530-625\,GHz and 1067-1275\,GHz). This is coupled to radio, stellar, and solar occultation opportunities to explore the atmosphere at high vertical resolution; and radio and plasma wave measurements of electric discharges in the Jovian atmosphere and auroras. Cross-disciplinary scientific investigations enable JUICE to explore coupling processes in giant planet atmospheres, to show how the atmosphere is connected to (i) the deep circulation and composition of the hydrogen-dominated interior; and (ii) to the currents and charged particle environments of the external magnetosphere. JUICE will provide a comprehensive characterisation of the atmosphere and auroras of this archetypal giant planet.
△ Less
Submitted 26 October, 2023; v1 submitted 20 April, 2023;
originally announced April 2023.
-
Fundamental groups of proper varieties are finitely presented
Authors:
Marcin Lara,
Vasudevan Srinivas,
Jakob Stix
Abstract:
It was recently proven by Esnault, Shusterman and the second named author, that the étale fundamental group of a connected smooth projective variety over an algebraically closed field $k$ is finitely presented. In this note, we extend this result to all connected proper schemes over $k$.
It was recently proven by Esnault, Shusterman and the second named author, that the étale fundamental group of a connected smooth projective variety over an algebraically closed field $k$ is finitely presented. In this note, we extend this result to all connected proper schemes over $k$.
△ Less
Submitted 26 May, 2023; v1 submitted 16 March, 2023;
originally announced March 2023.
-
Decomposition Methods for Dynamically Monotone Two-Time-Scale Stochastic Optimization Problems
Authors:
Tristan Rigaut,
Pierre Carpentier,
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
In energy management, it is common that strategic investment decisions (storage capacity, production units) are made at a slow time scale, whereas operational decisions (storage, production) are made at a fast time scale: for such problems, the total number of decision stages may be huge. In this paper, we consider multistage stochastic optimization problems with two time-scales, and we propose a…
▽ More
In energy management, it is common that strategic investment decisions (storage capacity, production units) are made at a slow time scale, whereas operational decisions (storage, production) are made at a fast time scale: for such problems, the total number of decision stages may be huge. In this paper, we consider multistage stochastic optimization problems with two time-scales, and we propose a time block decomposition scheme to address them numerically. More precisely, our approach relies on two assumptions. On the one hand, we suppose slow time scale stagewise independence of the noise process: the random variables that occur during a slow time scale interval are independent of those at another slow time scale interval. This makes it possible to use Dynamic Programming at the slow time scale. On the other hand, we suppose a dynamically monotone property for the problem under consideration, which makes it possible to obtain bounds. Then, we present two algorithmic methods to compute upper and lower bounds for slow time scale Bellman value functions. Both methods rely respectively on primal and dual decomposition of the Bellman equation applied at the slow time scale. We assess the methods tractability and validate their efficiency by solving a battery management problem where the fast time scale operational decisions have an impact on the storage current capacity, hence on the strategic decisions to renew the battery at the slow time scale.
△ Less
Submitted 7 March, 2023;
originally announced March 2023.
-
Universal behavior in complex-mediated reactions: Dynamics of S(1D)+ o-D2 --> D + SD at low collision energies
Authors:
Manuel Lara,
P. G. Jambrina,
F. J. Aoiz
Abstract:
Reactive and elastic cross-sections, and rate coefficients, have been calculated for the S(1D)+ D2 (v=0, j=0) reaction using a modified hyperspherical quantum reactive scattering method. The considered collision energy ranges from the ultracold regime, where only one partial wave is open, up to the Langevin regime, where many of them contribute. This work presents the extension of the quantum calc…
▽ More
Reactive and elastic cross-sections, and rate coefficients, have been calculated for the S(1D)+ D2 (v=0, j=0) reaction using a modified hyperspherical quantum reactive scattering method. The considered collision energy ranges from the ultracold regime, where only one partial wave is open, up to the Langevin regime, where many of them contribute. This work presents the extension of the quantum calculations, which were compared with the experimental results in a previous work, down to energies in the cold and ultracold domains. Results are analyzed and compared with the universal case of the quantum defect theory by Jachymski et al. [Phys. Rev. Lett. 110, 213202 (2013)]. State-to-state integral and differential cross sections are also shown covering the ranges of low-thermal, cold and ultracold collision energy regimes. It is found that at E/k_B T < 1 K there are substantial departures from the expected statistical behavior, and that dynamical features become increasingly important with decreasing collision energy, leading to vibrational excitation.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
A Model for Understanding and Reducing Developer Burnout
Authors:
Bianca Trinkenreich,
Klaas-Jan Stol,
Igor Steinmacher,
Marco Gerosa,
Anita Sarma,
Marcelo Lara,
Michael Feathers,
Nicholas Ross,
Kevin Bishop
Abstract:
Job burnout is a type of work-related stress associated with a state of physical or emotional exhaustion that also involves a sense of reduced accomplishment and loss of personal identity. Burnt out can affect one's physical and mental health and has become a leading industry concern and can result in high workforce turnover. Through an empirical study at Globant, a large multi-national company, w…
▽ More
Job burnout is a type of work-related stress associated with a state of physical or emotional exhaustion that also involves a sense of reduced accomplishment and loss of personal identity. Burnt out can affect one's physical and mental health and has become a leading industry concern and can result in high workforce turnover. Through an empirical study at Globant, a large multi-national company, we created a theoretical model to evaluate the complex interplay among organizational culture, work satisfaction, and team climate, and how they impact developer burnout. We conducted a survey of developers in software delivery teams (n=3281) to test our model and analyzed the data using structural equation modeling, moderation, and multi-group analysis. Our results show that Organizational Culture, Climate for Learning, Sense of Belonging, and Inclusiveness are positively associated with Work Satisfaction, which in turn is associated with Reduced Burnout. Our model generated through a large-scale survey can guide organizations in how to reduce workforce burnout by creating a climate for learning, inclusiveness in teams, and a generative organizational culture where new ideas are welcome, information is actively sought and bad news can be shared without fear.
△ Less
Submitted 24 January, 2023; v1 submitted 22 January, 2023;
originally announced January 2023.
-
Contributions on complexity bounds for Deterministic Partially Observed Markov Decision Process
Authors:
Cyrille Vessaire,
Jean-Philippe Chancelier,
Michel de Lara,
Pierre Carpentier,
Alejandro Rodríguez-Martínez
Abstract:
Markov Decision Processes (Mdps) form a versatile framework used to model a wide range of optimization problems. The Mdp model consists of sets of states, actions, time steps, rewards, and probability transitions. When in a given state and at a given time, the decision maker's action generates a reward and determines the state at the next time step according to the probability transition function.…
▽ More
Markov Decision Processes (Mdps) form a versatile framework used to model a wide range of optimization problems. The Mdp model consists of sets of states, actions, time steps, rewards, and probability transitions. When in a given state and at a given time, the decision maker's action generates a reward and determines the state at the next time step according to the probability transition function. However, Mdps assume that the decision maker knows the state of the controlled dynamical system. Hence, when one needs to optimize controlled dynamical systems under partial observation, one often turns toward the formalism of Partially Observed Markov Decision Processes (Pomdp). Pomdps are often untractable in the general case as Dynamic Programming suffers from the curse of dimensionality. Instead of focusing on the general Pomdps, we present a subclass where transitions and observations mappings are deterministic: Deterministic Partially Observed Markov Decision Processes (Det-Pomdp). That subclass of problems has been studied by (Littman, 1996) and (Bonet, 2009). It was first considered as a limit case of Pomdps by Littman, mainly used to illustrate the complexity of Pomdps when considering as few sources of uncertainties as possible. In this paper, we improve on Littman's complexity bounds. We then introduce and study an even simpler class: Separated Det-Pomdps and give some new complexity bounds for this class. This new class of problems uses a property of the dynamics and observation to push back the curse of dimensionality.
△ Less
Submitted 20 January, 2023;
originally announced January 2023.
-
Differentiability and Regularization of Parametric Convex Value Functions in Stochastic Multistage Optimization
Authors:
Adrien Le Franc,
Jean-Philippe Chancelier,
Pierre Carpentier,
Michel de Lara
Abstract:
In multistage decision problems, it is often the case that an initial strategic decision (such as investment) is followed by many operational ones (operating the investment). Such initial strategic decision can be seen as a parameter affecting a multistage decision problem. More generally, we study in this paper a standard multistage stochastic optimization problem depending on a parameter. When t…
▽ More
In multistage decision problems, it is often the case that an initial strategic decision (such as investment) is followed by many operational ones (operating the investment). Such initial strategic decision can be seen as a parameter affecting a multistage decision problem. More generally, we study in this paper a standard multistage stochastic optimization problem depending on a parameter. When the parameter is fixed, Stochastic Dynamic Programming provides a way to compute the optimal value of the problem. Thus, the value function depends both on the state (as usual) and on the parameter. Our aim is to investigate on the possibility to efficiently compute gradients of the value function with respect to the parameter, when these objects exist. When nondifferentiable, we propose a regularization method based on the Moreau-Yosida envelope. We present a numerical test case from day-ahead power scheduling.
△ Less
Submitted 6 February, 2023; v1 submitted 20 December, 2022;
originally announced December 2022.
-
Duality Between Lagrangians and Rockafellians
Authors:
Michel de Lara
Abstract:
In his monograph \emph{Conjugate Duality and Optimization}, Rockafellar puts forward a ``perturbation + duality'' method to obtain a dual problem for an original minimization problem. First, one embeds the minimization problem into a family of perturbed problems (thus giving a so-called perturbation function); the perturbation of the original function to be minimized has recently been called…
▽ More
In his monograph \emph{Conjugate Duality and Optimization}, Rockafellar puts forward a ``perturbation + duality'' method to obtain a dual problem for an original minimization problem. First, one embeds the minimization problem into a family of perturbed problems (thus giving a so-called perturbation function); the perturbation of the original function to be minimized has recently been called a Rockafellian. Second, when the perturbation variable belongs to a primal vector space paired, by a bilinear form, with a dual vector space, one builds a Lagrangian from a Rockafellian; one also obtains a so-called dual function (and a dual problem). The method has been extended from Fenchel duality to generalized convexity: when the perturbation belongs to a primal set paired, by a coupling function, with a dual set, one also builds a Rockafellian from a Lagrangian. Following these paths, we highlight a duality between Lagrangians and Rockafellians. Where the material mentioned above mostly focuses on moving from Rockafellian to Lagrangian, we treat them equally and display formulas that go both ways. We propose a definition of Lagrangian-Rockafellian couples. We characterize these latter as dual functions, with respect to a coupling, and also in terms of generalized convex functions. The duality between perturbation and dual functions is not as clear cut.
△ Less
Submitted 10 March, 2023; v1 submitted 23 November, 2022;
originally announced November 2022.
-
Time Consistency for Multistage Stochastic Optimization Problems under Constraints in Expectation
Authors:
Pierre Carpentier,
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
We consider sequences-indexed by time (discrete stages)-of families of multistage stochastic optimization problems. At each time, the optimization problems in a family are parameterized by some quantities (initial states, constraint levels.. .). In this framework, we introduce an adapted notion of time consistent optimal solutions, that is, solutions that remain optimal after truncation of the pas…
▽ More
We consider sequences-indexed by time (discrete stages)-of families of multistage stochastic optimization problems. At each time, the optimization problems in a family are parameterized by some quantities (initial states, constraint levels.. .). In this framework, we introduce an adapted notion of time consistent optimal solutions, that is, solutions that remain optimal after truncation of the past and that are optimal for any values of the parameters. We link this time consistency notion with the concept of state variable in Markov Decision Processes for a class of multistage stochastic optimization problems incorporating state constraints at the final time, either formulated in expectation or in probability. For such problems, when the primitive noise random process is stagewise independent and takes a finite number of values, we show that time consistent solutions can be obtained by considering a finite dimensional state variable. We illustrate our results on a simple dam management problem.
△ Less
Submitted 29 August, 2022;
originally announced August 2022.
-
A torsion-based solution to the hyperbolic regime of the J2-problem
Authors:
Martin Lara,
Alessandro Masat,
Camilla Colombo
Abstract:
A popular intermediary in the theory of artificial satellites is obtained after the elimination of parallactic terms from the J2-problem Hamiltonian. The resulting quasi-Keplerian system is in turn converted into the Kepler problem by a torsion. When this reduction process is applied to unbounded orbits the solution is made of Keplerian hyperbolae. For this last case, we show that the torsion-base…
▽ More
A popular intermediary in the theory of artificial satellites is obtained after the elimination of parallactic terms from the J2-problem Hamiltonian. The resulting quasi-Keplerian system is in turn converted into the Kepler problem by a torsion. When this reduction process is applied to unbounded orbits the solution is made of Keplerian hyperbolae. For this last case, we show that the torsion-based solution provides an effective alternative to the Keplerian approximation customarily used in flyby computations. Also, we check that the extension of the torsion-based solution to higher orders of the oblateness coefficient yields the expected convergence of asymptotic solutions to the true orbit.
△ Less
Submitted 24 July, 2022;
originally announced July 2022.
-
Accelerating high order discontinuous Galerkin solvers using neural networks: 3D compressible Navier-Stokes equations
Authors:
Fernando Manrique de Lara,
Esteban Ferrer
Abstract:
We propose to accelerate a high order discontinuous Galerkin solver using neural networks. We include a corrective forcing to a low polynomial order simulation to enhance its accuracy. The forcing is obtained by training a deep fully connected neural network, using a high polynomial order simulation but only for a short time frame. With this corrective forcing, we can run the low polynomial order…
▽ More
We propose to accelerate a high order discontinuous Galerkin solver using neural networks. We include a corrective forcing to a low polynomial order simulation to enhance its accuracy. The forcing is obtained by training a deep fully connected neural network, using a high polynomial order simulation but only for a short time frame. With this corrective forcing, we can run the low polynomial order simulation faster (with large time steps and low cost per time step) while improving its accuracy.
We explored this idea for a 1D Burgers' equation in (Marique and Ferrer, CAF 2022), and we have extended this work to the 3D Navier-Stokes equations, with and without a Large Eddy Simulation closure model. We test the methodology with the turbulent Taylor Green Vortex case and for various Reynolds numbers (30, 200 and 1600). In addition, the Taylor Green Vortex evolves with time and covers laminar, transitional, and turbulent regimes, as time progresses.
The proposed methodology proves to be applicable to a variety of flows and regimes. The results show that the corrective forcing is effective in all Reynolds numbers and time frames (excluding the initial flow development). We can train the corrective forcing with a polynomial order of 8, to increase the accuracy of simulations from a polynomial order 3 to 6, when correcting outside the training time frame. The low order correct solution is 4 to 5 times faster than a simulation with comparable accuracy (polynomial order 6).
Additionally, we explore changes in the hyperparameters and use transfer learning to speed up the training. We observe that it is not useful to train a corrective forcing using a different flow condition. However, an already trained corrective forcing can be used to initialise a new training (at the correct flow conditions) to obtain an effective forcing with only a few training iterations.
△ Less
Submitted 23 July, 2022;
originally announced July 2022.
-
HORSES3D: a high-order discontinuous Galerkin solver for flow simulations and multi-physics applications
Authors:
E. Ferrer,
G. Rubio,
G. Ntoukas,
W. Laskowski,
O. A. Mariño,
S. Colombo,
A. Mateo-Gabín,
F. Manrique de Lara,
D. Huergo,
J. Manzanero,
A. M. Rueda-Ramírez,
D. A. Kopriva,
E. Valero
Abstract:
We present the latest developments of our High-Order Spectral Element Solver (HORSES3D), an open source high-order discontinuous Galerkin framework, capable of solving a variety of flow applications, including compressible flows (with or without shocks), incompressible flows, various RANS and LES turbulence models, particle dynamics, multiphase flows, and aeroacoustics. We provide an overview of t…
▽ More
We present the latest developments of our High-Order Spectral Element Solver (HORSES3D), an open source high-order discontinuous Galerkin framework, capable of solving a variety of flow applications, including compressible flows (with or without shocks), incompressible flows, various RANS and LES turbulence models, particle dynamics, multiphase flows, and aeroacoustics. We provide an overview of the high-order spatial discretisation (including energy/entropy stable schemes) and anisotropic p-adaptation capabilities. The solver is parallelised using MPI and OpenMP showing good scalability for up to 1000 processors. Temporal discretisations include explicit, implicit, multigrid, and dual time-stepping schemes with efficient preconditioners. Additionally, we facilitate meshing and simulating complex geometries through a mesh-free immersed boundary technique. We detail the available documentation and the test cases included in the GitHub repository.
△ Less
Submitted 20 June, 2022;
originally announced June 2022.
-
Optimization of a domestic microgrid equipped with solar panel and battery: Model Predictive Control and Stochastic Dual Dynamic Programming approaches
Authors:
François Pacaud,
Pierre Carpentier,
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
In this study, a microgrid with storage (battery, hot water tank) and solar panel is considered. We benchmark two algorithms, MPC and SDDP, that yield online policies to manage the microgrid, and compare them with a rule based policy. Model Predictive Control (MPC) is a well-known algorithm which models the future uncertainties with a deterministic forecast. By contrast, Stochastic Dual Dynamic Pr…
▽ More
In this study, a microgrid with storage (battery, hot water tank) and solar panel is considered. We benchmark two algorithms, MPC and SDDP, that yield online policies to manage the microgrid, and compare them with a rule based policy. Model Predictive Control (MPC) is a well-known algorithm which models the future uncertainties with a deterministic forecast. By contrast, Stochastic Dual Dynamic Programming (SDDP) models the future uncertainties as stagewise independent random variables with known probability distributions. We present a scheme, based on out-of-sample validation, to fairly compare the two online policies yielded by MPC and SDDP. Our numerical studies put to light that MPC and SDDP achieve significant gains compared to the rule based policy, and that SDDP overperforms MPC not only on average but on most of the out-of-sample assessment scenarios.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
Earth satellite dynamics by Picard iterations
Authors:
Martin Lara
Abstract:
The main effects of the Earth's oblateness on the motion of artificial satellites are usually derived from the variation of parameters equations of an average representation of the oblateness disturbing function. Rather, we approach their solution under the strict mathematical assumptions of Picard's iterative method. Our approach recovers the known linear trends of the right ascension of the asce…
▽ More
The main effects of the Earth's oblateness on the motion of artificial satellites are usually derived from the variation of parameters equations of an average representation of the oblateness disturbing function. Rather, we approach their solution under the strict mathematical assumptions of Picard's iterative method. Our approach recovers the known linear trends of the right ascension of the ascending node and the argument of the perigee, but differs from the accepted solution in the value of the mean motion. This amended rate radically improves the in-track errors of typical orbit propagations. In addition, our truncation of the Picard iterations solution to its secular terms includes the corrections that must be applied to the osculating initial conditions in the right propagation of the mean dynamics.
△ Less
Submitted 5 May, 2022;
originally announced May 2022.
-
Variants of the de Jong fundamental group
Authors:
Piotr Achinger,
Marcin Lara,
Alex Youcis
Abstract:
For a rigid space $X$, we answer two questions of de Jong about the category $\mathbf{Cov}^\mathrm{adm}_X$ of coverings which are locally in the admissible topology on $X$ the disjoint union of finite etale coverings: we show that this class is different from the one used by de Jong, but still gives a tame infinite Galois category. In addition, we prove that the objects of…
▽ More
For a rigid space $X$, we answer two questions of de Jong about the category $\mathbf{Cov}^\mathrm{adm}_X$ of coverings which are locally in the admissible topology on $X$ the disjoint union of finite etale coverings: we show that this class is different from the one used by de Jong, but still gives a tame infinite Galois category. In addition, we prove that the objects of $\mathbf{Cov}^\mathrm{et}_X$ (with the analogous definition) correspond precisely to locally constant sheaves for the pro-etale topology defined by Scholze.
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
Detection of iron emission lines and a temperature inversion on the dayside of the ultra-hot Jupiter KELT-20b
Authors:
F. Yan,
A. Reiners,
E. Pallé,
D. Shulyak,
M. Stangret,
K. Molaverdikhani,
L. Nortmann,
P. Mollière,
Th. Henning,
N. Casasayas-Barris,
D. Cont,
G. Chen,
S. Czesla,
A. Sánchez-López,
M. López-Puertas,
I. Ribas,
A. Quirrenbach,
J. A. Caballero,
P. J. Amado,
D. Galadí-Enríquez,
S. Khalafinejad,
L. M. Lara,
D. Montes,
G. Morello,
E. Nagel
, et al. (3 additional authors not shown)
Abstract:
Ultra-hot Jupiters (UHJs) are gas giants with very high equilibrium temperatures. In recent years, multiple chemical species, including various atoms and ions, have been discovered in their atmospheres. Most of these observations have been performed with transmission spectroscopy, although UHJs are also ideal targets for emission spectroscopy due to their strong thermal radiation. We present high-…
▽ More
Ultra-hot Jupiters (UHJs) are gas giants with very high equilibrium temperatures. In recent years, multiple chemical species, including various atoms and ions, have been discovered in their atmospheres. Most of these observations have been performed with transmission spectroscopy, although UHJs are also ideal targets for emission spectroscopy due to their strong thermal radiation. We present high-resolution thermal emission spectroscopy of the transiting UHJ KELT-20b/MASCARA-2b. The observation was performed with the CARMENES spectrograph at orbital phases before and after the secondary eclipse. We detected atomic Fe using the cross-correlation technique. The detected Fe lines are in emission, which unambiguously indicates a temperature inversion on the dayside hemisphere. We furthermore retrieved the temperature structure with the detected Fe lines. The result shows that the atmosphere has a strong temperature inversion with a temperature of $4900\pm{700}$ K and a pressure of $10^{-4.8_{-1.1}^{+1.0}}$ bar at the upper layer of the inversion. A joint retrieval of the CARMENES data and the TESS secondary eclipse data returns a temperature of $2550_{-250}^{+150}$ K and a pressure of $10^{-1.5_{-0.6}^{+0.7}}$ bar at the lower layer of the temperature inversion. The detection of such a strong temperature inversion is consistent with theoretical simulations that predict an inversion layer on the dayside of UHJs. The joint retrieval of the CARMENES and TESS data demonstrates the power of combing high-resolution emission spectroscopy with secondary eclipse photometry in characterizing atmospheric temperature structures.
△ Less
Submitted 21 January, 2022;
originally announced January 2022.
-
Multistage Optimization of a Petroleum Production System with Material Balance Model
Authors:
Cyrille Vessaire,
Jean-Philippe Chancelier,
Michel de Lara,
Pierre Carpentier,
Alejandro Rodríguez-Martínez,
Anna Roberts
Abstract:
In this paper, we propose a mathematical formulation for the management of an oil production network as a multistage optimization problem. The reservoir is modeled as a controlled dynamical system by using material balance equations. We use a dynamic programming algorithm to solve the optimization problem. Two numerical applications illustrate our work: the first one consists in optimizing the pro…
▽ More
In this paper, we propose a mathematical formulation for the management of an oil production network as a multistage optimization problem. The reservoir is modeled as a controlled dynamical system by using material balance equations. We use a dynamic programming algorithm to solve the optimization problem. Two numerical applications illustrate our work: the first one consists in optimizing the production of a gas reservoir, whereas the second one tackles an oil reservoir with water injection.
△ Less
Submitted 20 September, 2022; v1 submitted 4 January, 2022;
originally announced January 2022.
-
The Capra-subdifferential of the l0 pseudonorm
Authors:
Adrien Le Franc,
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
The l0 pseudonorm counts the nonzero coordinates of a vector. It is often used in optimization problems to enforce the sparsity of the solution. However, this function is nonconvex and noncontinuous, and optimization problems formulated with l0 in the objective function or in the constraints are hard to solve in general. Recently, a new family of coupling functions - called Capra (constant along p…
▽ More
The l0 pseudonorm counts the nonzero coordinates of a vector. It is often used in optimization problems to enforce the sparsity of the solution. However, this function is nonconvex and noncontinuous, and optimization problems formulated with l0 in the objective function or in the constraints are hard to solve in general. Recently, a new family of coupling functions - called Capra (constant along primal rays) - has proved to induce relevant generalized Fenchel-Moreau conjugacies to handle the l0 pseudonorm. In particular, under a suitable choice of source norm on the Euclidean space used in the definition of the Capra coupling - the function l0 is Capra-subdifferentiable, hence is Capra-convex. In this article, we give explicit formulations for the Capra subdifferential of l0, when the source norm is a lp norm with p larger that 1. We illustrate our results with graphical visualizations of the Capra subdifferential of l0 for the Euclidean source norm.
△ Less
Submitted 18 August, 2022; v1 submitted 31 December, 2021;
originally announced December 2021.
-
Note on the analytical integration of circumterrestrial orbits
Authors:
Martin Lara
Abstract:
The acclaimed merits of analytical solutions based on a fictitious time developed in the 1970's were partially overvalued due to a common misuse of classical analytical solutions based on the physical time that were taken as reference. With the main problem of the artificial satellite theory as a model, we carry out a more objective comparison of both kinds of theories. We find that the proper ini…
▽ More
The acclaimed merits of analytical solutions based on a fictitious time developed in the 1970's were partially overvalued due to a common misuse of classical analytical solutions based on the physical time that were taken as reference. With the main problem of the artificial satellite theory as a model, we carry out a more objective comparison of both kinds of theories. We find that the proper initialization of classical solutions notably balances the performance of the two distinct approaches in what respects to accuracy. Besides, extension of both kinds of satellite theories to higher orders show additional pros and cons of each different perturbation approach, thus providing complementary information to prospective users on which kind of analytical solution may better support their needs.
△ Less
Submitted 17 October, 2021;
originally announced October 2021.
-
Causal Inference Theory with Information Dependency Models
Authors:
Benjamin Heymann,
Michel de Lara,
Jean-Philippe Chancelier
Abstract:
Inferring the potential consequences of an unobserved event is a fundamental scientific question. To this end, Pearl's celebrated do-calculus provides a set of inference rules to derive an interventional probability from an observational one. In this framework, the primitive causal relations are encoded as functional dependencies in a Structural Causal Model (SCM), which are generally mapped into…
▽ More
Inferring the potential consequences of an unobserved event is a fundamental scientific question. To this end, Pearl's celebrated do-calculus provides a set of inference rules to derive an interventional probability from an observational one. In this framework, the primitive causal relations are encoded as functional dependencies in a Structural Causal Model (SCM), which are generally mapped into a Directed Acyclic Graph (DAG) in the absence of cycles. In this paper, by contrast, we capture causality without reference to graphs or functional dependencies, but with information fields and Witsenhausen's intrinsic model. The three rules of do-calculus reduce to a unique sufficient condition for conditional independence, the topological separation, which presents interesting theoretical and practical advantages over the d-separation. With this unique rule, we can deal with systems that cannot be represented with DAGs, for instance systems with cycles and/or 'spurious' edges. We treat an example that cannot be handled-to the extent of our knowledge-with the tools of the current literature. We also explain why, in the presence of cycles, the theory of causal inference might require different tools, depending on whether the random variables are discrete or continuous.
△ Less
Submitted 9 August, 2021; v1 submitted 6 August, 2021;
originally announced August 2021.
-
Topological Conditional Separation
Authors:
Michel de Lara,
Jean-Philippe Chancelier,
Benjamin Heymann
Abstract:
Pearl's d-separation is a foundational notion to study conditional independence between random variables. We define the topological conditional separation and we show that it is equivalent to the d-separation, extended beyond acyclic graphs, be they finite or infinite.
Pearl's d-separation is a foundational notion to study conditional independence between random variables. We define the topological conditional separation and we show that it is equivalent to the d-separation, extended beyond acyclic graphs, be they finite or infinite.
△ Less
Submitted 6 August, 2021;
originally announced August 2021.
-
Conditional Separation as a Binary Relation. A Coq Assisted Proof
Authors:
Jean-Philippe Chancelier,
Michel de Lara,
Benjamin Heymann
Abstract:
The concept of d-separation holds a pivotal role in causality theory, serving as a fundamental tool for deriving conditional independence properties from causal graphs. Pearl defined the d-separation of two subsets conditionally on a third one. In this study, we present a novel perspective by showing i) how the d-separation can be extended beyond acyclic graphs, possibly infinite, and ii) how…
▽ More
The concept of d-separation holds a pivotal role in causality theory, serving as a fundamental tool for deriving conditional independence properties from causal graphs. Pearl defined the d-separation of two subsets conditionally on a third one. In this study, we present a novel perspective by showing i) how the d-separation can be extended beyond acyclic graphs, possibly infinite, and ii) how it can be expressed and characterized as a binary relation between vertices. Compared to the typical perspectives in causality theory, our equivalence opens the door to more compact and computational proofing techniques, because the language of binary relations is well adapted to equational reasoning. Additionally, and of independent interest, the proofs of the results presented in this paper are checked with the Coq proof assistant.
△ Less
Submitted 2 April, 2024; v1 submitted 6 August, 2021;
originally announced August 2021.
-
Specialization for the pro-étale fundamental group
Authors:
Piotr Achinger,
Marcin Lara,
Alex Youcis
Abstract:
For a formal scheme $\mathfrak{X}$ of finite type over a complete rank one valuation ring, we construct a specialization morphism \[ π^{\rm dJ}_1(\mathfrak{X}_η) \to π^{\rm proet}_1(\mathfrak{X}_k) \] from the de Jong fundamental group of the rigid generic fiber to the Bhatt-Scholze pro-étale fundamental group of the special fiber. The construction relies on an interplay between admissible blowups…
▽ More
For a formal scheme $\mathfrak{X}$ of finite type over a complete rank one valuation ring, we construct a specialization morphism \[ π^{\rm dJ}_1(\mathfrak{X}_η) \to π^{\rm proet}_1(\mathfrak{X}_k) \] from the de Jong fundamental group of the rigid generic fiber to the Bhatt-Scholze pro-étale fundamental group of the special fiber. The construction relies on an interplay between admissible blowups of $\mathfrak{X}$ and normalizations of the irreducible components of $\mathfrak{X}_k$, and employs the Berthelot tubes of these irreducible components in an essential way. Using related techniques, we show that under certain smoothness and semistability assumptions, covering spaces in the sense of de Jong of a smooth rigid space which are tame satisfy étale descent.
△ Less
Submitted 14 July, 2021;
originally announced July 2021.
-
Minimization Interchange Theorem on Posets
Authors:
Jean-Philippe Chancelier,
Michel de Lara,
Benoît Tran
Abstract:
Interchange theorems between minimization and integration are useful in optimization, especially in optimal control and in stochastic optimization. In this article, we establish a generalized minimization interchange theorem, where integration is replaced by a monotone mapping between posets (partially ordered sets). As an application, we recover, and slightly extend, classical results from the li…
▽ More
Interchange theorems between minimization and integration are useful in optimization, especially in optimal control and in stochastic optimization. In this article, we establish a generalized minimization interchange theorem, where integration is replaced by a monotone mapping between posets (partially ordered sets). As an application, we recover, and slightly extend, classical results from the literature, and we tackle the case of the Choquet integral. Our result provides insight on the mechanisms behind existing interchange results.
△ Less
Submitted 13 July, 2021;
originally announced July 2021.
-
Decentralized Multistage Optimization of Large-Scale Microgrids under Stochasticity
Authors:
François Pacaud,
Michel de Lara,
Jean-Philippe Chancelier,
Pierre Carpentier
Abstract:
Microgrids are recognized as a relevant tool to absorb decentralized renewable energies in the energy mix. However, the sequential handling of multiple stochastic productions and demands, and of storage, make their management a delicate issue. We add another layer of complexity by considering microgrids where different buildings stand at the nodes of a network and are connected by the arcs; some b…
▽ More
Microgrids are recognized as a relevant tool to absorb decentralized renewable energies in the energy mix. However, the sequential handling of multiple stochastic productions and demands, and of storage, make their management a delicate issue. We add another layer of complexity by considering microgrids where different buildings stand at the nodes of a network and are connected by the arcs; some buildings host local production and storage capabilities, and can exchange with others their energy surplus. We formulate the problem as a multistage stochastic optimization problem, corresponding to the minimization of the expected temporal sum of operational costs, while satisfying the energy demand at each node, for all time. The resulting mathematical problem has a large-scale nature, exhibiting both spatial and temporal couplings. However, the problem displays a network structure that makes it amenable to a mix of spatial decomposition-coordination with temporal decomposition methods. We conduct numerical simulations on microgrids of different sizes and topologies, with up to 48 nodes and 64 state variables. Decomposition methods are faster and provide more efficient policies than a state-of-the-art Stochastic Dual Dynamic Programming algorithm. Moreover, they scale almost linearly with the state dimension, making them a promising tool to address more complex microgrid optimal management problems.
△ Less
Submitted 8 June, 2021;
originally announced June 2021.
-
Decomposition-Coordination Method for Finite Horizon Bandit Problems
Authors:
Michel de Lara,
Benjamin Heymann,
Jean-Philippe Chancelier
Abstract:
Optimally solving a multi-armed bandit problem suffers the curse of dimensionality. Indeed, resorting to dynamic programming leads to an exponential growth of computing time, as the number of arms and the horizon increase. We introduce a decompositioncoordination heuristic, DeCo, that turns the initial problem into parallelly coordinated one-armed bandit problems. As a consequence, we obtain a com…
▽ More
Optimally solving a multi-armed bandit problem suffers the curse of dimensionality. Indeed, resorting to dynamic programming leads to an exponential growth of computing time, as the number of arms and the horizon increase. We introduce a decompositioncoordination heuristic, DeCo, that turns the initial problem into parallelly coordinated one-armed bandit problems. As a consequence, we obtain a computing time which is essentially linear in the number of arms. In addition, the decomposition provides a theoretical lower bound on the regret. For the two-armed bandit case, dynamic programming provides the exact solution, which is almost matched by the DeCo heuristic. Moreover, in numerical simulations with up to 100 rounds and 20 arms, DeCo outperforms classic algorithms (Thompson sampling and Kullback-Leibler upper-confidence bound) and almost matches the theoretical lower bound on the regret for 20 arms.
△ Less
Submitted 21 May, 2024; v1 submitted 2 June, 2021;
originally announced June 2021.
-
Best Convex Lower Approximations of the l 0 Pseudonorm on Unit Balls
Authors:
Thomas Bittar,
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
Whereas the norm of a vector measures amplitude (and is a 1-homogeneous function), sparsity is measured by the 0-homogeneous l0 pseudonorm, which counts the number of nonzero components. We propose a family of conjugacies suitable for the analysis of 0-homogeneous functions. These conjugacies are derived from couplings between vectors, given by their scalar product divided by a 1-homogeneous norma…
▽ More
Whereas the norm of a vector measures amplitude (and is a 1-homogeneous function), sparsity is measured by the 0-homogeneous l0 pseudonorm, which counts the number of nonzero components. We propose a family of conjugacies suitable for the analysis of 0-homogeneous functions. These conjugacies are derived from couplings between vectors, given by their scalar product divided by a 1-homogeneous normalizing factor. With this, we characterize the best convex lower approximation of a 0-homogeneous function on the unit ''ball'' of a normalization function (i.e. a norm without the requirement of subadditivity). We do the same with the best convex and 1-homogeneous lower approximation. In particular, we provide expressions for the tightest convex lower approximation of the l0 pseudonorm on any unit ball, and we show that the tightest norm which minorizes the l0 pseudonorm on the unit ball of any lp-norm is the l1-norm. We also provide the tightest convex lower convex approximation of the l0 pseudonorm on the unit ball of any norm.
△ Less
Submitted 31 May, 2021;
originally announced May 2021.
-
Rank-Based Norms, Capra-Conjugacies and the Rank Function
Authors:
Paul Barbier,
Jean-Philippe Chancelier,
Michel de Lara,
Valentin Paravy
Abstract:
We consider the space of matrices, with given number of rows and of columns, equipped with the classic trace scalar product. With any matrix (source) norm, we associate a coupling, called Capra, between the space of matrices and itself. Then, we compute the Capra conjugate and biconjugate of the rank function. They are expressed in function of a sequence of rank-based norms, more precisely general…
▽ More
We consider the space of matrices, with given number of rows and of columns, equipped with the classic trace scalar product. With any matrix (source) norm, we associate a coupling, called Capra, between the space of matrices and itself. Then, we compute the Capra conjugate and biconjugate of the rank function. They are expressed in function of a sequence of rank-based norms, more precisely generalized r-rank and dual r-rank matrix norms associated with the matrix source norm. We deduce a lower bound of the rank function given by a variational formula which involves the generalized r-rank norms. In the case of the Frobenius norm, we show that the rank function is equal to the variational formula.
△ Less
Submitted 6 February, 2023; v1 submitted 31 May, 2021;
originally announced May 2021.
-
A new radial, natural, higher order intermediary of the main problem four decades after the elimination of the parallax
Authors:
Martin Lara
Abstract:
Simplifications in dealing with the equation of the center when the short-period effects are removed from a zonal Hamiltonian are commonly attributed to the elimination of parallactic terms. But this interpretation is incorrect, and the simplifications rather stem from the removal of concomitant long-period terms, an outcome that can also be achieved without need of eliminating the parallax. To sh…
▽ More
Simplifications in dealing with the equation of the center when the short-period effects are removed from a zonal Hamiltonian are commonly attributed to the elimination of parallactic terms. But this interpretation is incorrect, and the simplifications rather stem from the removal of concomitant long-period terms, an outcome that can also be achieved without need of eliminating the parallax. To show that, a Lie transforms simplification is invented that augments the exponents of the inverse of the radius and still achieves analogous simplifications in handling the equation of the center to those provided by the classical elimination of the parallax simplification. The particular case in which the new transformation does not modify the exponents of the parallactic terms of the original problem, leads to a new intermediary of the main problem that, while keeping higher order effects of J2, is formally analogous to Cid's first order radial intermediary.
△ Less
Submitted 30 May, 2021;
originally announced May 2021.
-
Design of distant retrograde orbits based on a higher order analytical solution
Authors:
Martin Lara
Abstract:
Because the distant retrograde orbits dynamics inherently depends on special functions, approximate analytical solutions in the literature are commonly constrained to providing rough approximations of the qualitative behavior. We rely on perturbation methods and succeed in computing a higher order analytical solution that extends the range of applicability of previous solutions to the problem of r…
▽ More
Because the distant retrograde orbits dynamics inherently depends on special functions, approximate analytical solutions in the literature are commonly constrained to providing rough approximations of the qualitative behavior. We rely on perturbation methods and succeed in computing a higher order analytical solution that extends the range of applicability of previous solutions to the problem of relative, quasi-satellite orbits with large librations. Besides, the analytical solution provides two design parameters that are effectively used in the computation of periodic, distant retrograde orbits without constraining to the typical 1:1 resonance.
△ Less
Submitted 30 May, 2021;
originally announced May 2021.
-
Geometric arcs and fundamental groups of rigid spaces
Authors:
Piotr Achinger,
Marcin Lara,
Alex Youcis
Abstract:
We develop the notion of a geometric covering of a rigid space X, which yields a much larger class of covering spaces than that studied previously by de Jong. Geometric coverings of X are closed under disjoint unions and are étale local on X. If X is connected, its geometric coverings form a tame infinite Galois category, and hence are classified by a topological group. The definition is based on…
▽ More
We develop the notion of a geometric covering of a rigid space X, which yields a much larger class of covering spaces than that studied previously by de Jong. Geometric coverings of X are closed under disjoint unions and are étale local on X. If X is connected, its geometric coverings form a tame infinite Galois category, and hence are classified by a topological group. The definition is based on the property of lifting of "geometric arcs," making it similar to geometric coverings of schemes studied by Bhatt and Scholze as well as semicoverings of topological spaces introduced by Brazas.
△ Less
Submitted 23 March, 2022; v1 submitted 11 May, 2021;
originally announced May 2021.
-
A vehicle routing problem for biological sample transportation in healthcare: mathematical formulations and a metaheuristic approach
Authors:
Mario Benini,
Paolo Detti,
Garazi Zabalo Manrique de Lara
Abstract:
In this paper, a real-world transportation problem is addressed, concerning the collection and the transportation of biological sample tubes from sampling points to a main hospital. Blood and other biological samples are collected in different centers during morning hours. Then, the samples are transported to the main hospital, for their analysis, by a fleet of vehicles located in geographically d…
▽ More
In this paper, a real-world transportation problem is addressed, concerning the collection and the transportation of biological sample tubes from sampling points to a main hospital. Blood and other biological samples are collected in different centers during morning hours. Then, the samples are transported to the main hospital, for their analysis, by a fleet of vehicles located in geographically distributed depots. Each sample has a limited lifetime and must arrive to the main hospital within that time. If a sample cannot arrive to the hospital within the lifetime, either is discarded or must be processed in dedicated facilities called Spoke Centers.Two Mixed Integer Linear Programming formulations and an Adaptive Large Neighborhood Search (ALNS) metaheuristic algorithm have been developed for the problem. Computational experiments on different sets of instances based on real-life data provided by the Local Healthcare Authority of Bologna, Italy, are presented. A comparison on small instances with the optimal solutions obtained by the formulations shows the effectiveness of the proposed ALNS algorithm. On real-life instances, different batching policies of the samples are evaluated. The results show that the ALNS algorithm is able to find solutions in which all the samples are delivered on time, while in the real case about the 40% [5] of the samples is delivered late.
△ Less
Submitted 12 March, 2021;
originally announced April 2021.
-
Evidence of energy-, recombination-, and photon-limited escape regimes in giant planet H/He atmospheres
Authors:
M. Lampón,
M. López-Puertas,
S. Czesla,
A. Sánchez-López,
L. M. Lara,
M. Salz,
J. Sanz-Forcada,
K. Molaverdikhani,
A. Quirrenbach,
E. Pallé,
J. A. Caballero,
Th. Henning,
L. Nortmann,
P. J. Amado,
D. Montes,
A. Reiners,
I. Ribas
Abstract:
Hydrodynamic escape is the most efficient atmospheric mechanism of planetary mass loss and has a large impact on planetary evolution. Three hydrodynamic escape regimes have been identified theoretically: energy-limited, recombination-limited, and photon-limited. However, no evidence of these regimes had been reported until now. Here, we report evidence of these three regimes via an analysis of hel…
▽ More
Hydrodynamic escape is the most efficient atmospheric mechanism of planetary mass loss and has a large impact on planetary evolution. Three hydrodynamic escape regimes have been identified theoretically: energy-limited, recombination-limited, and photon-limited. However, no evidence of these regimes had been reported until now. Here, we report evidence of these three regimes via an analysis of helium I triplet at 10830 angstroms and Ly-$α$ absorption involving a 1D hydrodynamic model that allows us to estimate hydrogen recombination and advection rates. In particular, we show that HD 209458 b is in the energy-limited regime, HD 189733 b is in the recombination-limited regime, and GJ 3470 b is in the photon-limited regime. These exoplanets can be considered as benchmark cases for their respective regimes.
△ Less
Submitted 18 April, 2021;
originally announced April 2021.
-
Conditional Infimum and Hidden Convexity in Optimization
Authors:
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
Detecting hidden convexity is one of the tools to address nonconvex minimization problems. After giving a formal definition of hidden convexity, we introduce the notion of conditional infimum, as it will prove instrumental in detecting hidden convexity. We develop the theory of the conditional infimum, and we establish a tower property, relevant for minimization problems. Thus equipped, we provide…
▽ More
Detecting hidden convexity is one of the tools to address nonconvex minimization problems. After giving a formal definition of hidden convexity, we introduce the notion of conditional infimum, as it will prove instrumental in detecting hidden convexity. We develop the theory of the conditional infimum, and we establish a tower property, relevant for minimization problems. Thus equipped, we provide a sufficient condition for hidden convexity in nonconvex minimization problems. We illustrate our result on nonconvex quadratic minimization problems. We conclude with perspectives for using the conditional infimum in relation to the so-called S-procedure, to couplings and conjugacies, and to lower bound convex programs.
△ Less
Submitted 12 April, 2021;
originally announced April 2021.
-
Kuhn's Equivalence Theorem for Games in Product Form
Authors:
Benjamin Heymann,
Michel de Lara,
Jean-Philippe Chancelier
Abstract:
We propose an alternative to the tree representation of extensive form games. Games in product form represent information with $σ$-fields over a product set, and do not require an explicit description of the play temporality, as opposed to extensive form games on trees. This representation encompasses games with a continuum of actions, randomness and players, as well as games for which the play or…
▽ More
We propose an alternative to the tree representation of extensive form games. Games in product form represent information with $σ$-fields over a product set, and do not require an explicit description of the play temporality, as opposed to extensive form games on trees. This representation encompasses games with a continuum of actions, randomness and players, as well as games for which the play order cannot be determined in advance. We adapt and prove Kuhn's theorem-regarding equivalence between mixed and behavioral strategies under perfect recall-for games in product form with continuous action sets.
△ Less
Submitted 13 July, 2022; v1 submitted 12 April, 2021;
originally announced April 2021.
-
A theorem on meromorphic descent and the specialization of the pro-étale fundamental group
Authors:
Marcin Lara,
Jiu-Kang Yu,
Lei Zhang
Abstract:
Given a Noetherian formal scheme $\hat X$ over ${\rm Spf}(R)$, where $R$ is a complete DVR, we first prove a theorem of meromorphic descent along a possibly infinite cover of $\hat{X}$. Using this we construct a specialization functor from the category of continuous representations of the pro-étale fundamental group of the special fiber to the category of $F$-divided sheaves on the generic fiber.…
▽ More
Given a Noetherian formal scheme $\hat X$ over ${\rm Spf}(R)$, where $R$ is a complete DVR, we first prove a theorem of meromorphic descent along a possibly infinite cover of $\hat{X}$. Using this we construct a specialization functor from the category of continuous representations of the pro-étale fundamental group of the special fiber to the category of $F$-divided sheaves on the generic fiber. This specialization functor partially recovers the specialization functor of the étale fundamental groups. We also express the pro-étale fundamental group of a connected scheme $X$ of finite type over a field as coproducts and quotients of the free group and the étale fundamental groups of the normalizations of the irreducible components of $X$ and those of its singular loci.
△ Less
Submitted 6 February, 2022; v1 submitted 21 March, 2021;
originally announced March 2021.
-
A Fresh Geometrical Look at the General S-Procedure
Authors:
Michel de Lara,
Jean-Baptiste Hiriart-Urruty
Abstract:
We revisit the S-procedure for general functions with "geometrical glasses". We thus delineate a necessary condition, and almost a sufficient condition, to have the S-procedure valid. Everything is expressed in terms of convexity of augmented sets (convex hulls, conical hulls) of images built from the data functions.
We revisit the S-procedure for general functions with "geometrical glasses". We thus delineate a necessary condition, and almost a sufficient condition, to have the S-procedure valid. Everything is expressed in terms of convexity of augmented sets (convex hulls, conical hulls) of images built from the data functions.
△ Less
Submitted 12 May, 2021; v1 submitted 11 February, 2021;
originally announced February 2021.
-
Optimal Joint Allocation of Efforts in Inclusive Fitness by Related Individuals
Authors:
Michel de Lara
Abstract:
Families are places of affection and cooperation, but also of conflict. In his famous paper Parent-Offspring Conflict, Robert L. Trivers builds upon W. D. Hamilton's concept of inclusive fitness to argue for genetic conflict in parent-offspring relationships, and to derive numerical predictions on the intensity of the conflict. We propose a mathematical model of game theory that depicts how each m…
▽ More
Families are places of affection and cooperation, but also of conflict. In his famous paper Parent-Offspring Conflict, Robert L. Trivers builds upon W. D. Hamilton's concept of inclusive fitness to argue for genetic conflict in parent-offspring relationships, and to derive numerical predictions on the intensity of the conflict. We propose a mathematical model of game theory that depicts how each member of a family allocates her resource budget to maximize her inclusive fitness; this latter is made of the sum of personal fitness plus the sum of relatives fitnesses weighted by Wright's coefficients of relationship. We define an optimal allocation profile as a Nash equilibrium, and we characterize the solutions in function of resource budgets, coefficients of relationship and derivatives of personal fitnesses.
△ Less
Submitted 17 November, 2020;
originally announced November 2020.
-
Constant Along Primal Rays Conjugacies and Generalized Convexity for Functions of the Support
Authors:
Jean-Philippe Chancelier,
Michel de Lara
Abstract:
The support of a vector in R d is the set of indices with nonzero entries. Functions of the support have the property to be 0-homogeneous and, because of that, the Fenchel conjugacy fails to provide relevant analysis. In this paper, we define the coupling Capra between R d and itself by dividing the classic Fenchel scalar product coupling by a given (source) norm on R d. Our main result is that, w…
▽ More
The support of a vector in R d is the set of indices with nonzero entries. Functions of the support have the property to be 0-homogeneous and, because of that, the Fenchel conjugacy fails to provide relevant analysis. In this paper, we define the coupling Capra between R d and itself by dividing the classic Fenchel scalar product coupling by a given (source) norm on R d. Our main result is that, when both the source norm and its dual norm are orthant-strictly monotonic, any nondecreasing finite-valued function of the support mapping is Capra-convex, that is, is equal to its Capra-biconjugate (generalized convexity). We also establish that any such function is the composition of a proper convex lower semi continuous function on R d with the normalization mapping on the unit sphere (hidden convexity), and that, when normalized, it admits a variational formulation, which involves a family of generalized local-K-support dual norms.
△ Less
Submitted 23 October, 2020;
originally announced October 2020.
-
A data-driven approach to constraining the atmospheric temperature structure of KELT-9b
Authors:
L. Fossati,
D. Shulyak,
A. G. Sreejith,
T. Koskinen,
M. E. Young,
P. E. Cubillos,
L. M. Lara,
K. France,
M. Rengel,
P. W. Cauley,
J. D. Turner,
A. Wyttenbach,
F. Yan
Abstract:
Context. Observationally constraining the atmospheric temperature-pressure (TP) profile of exoplanets is an important step forward for improving planetary atmosphere models, further enabling one to place the detection of spectral features and the measurement of atomic and molecular abundances through transmission and emission spectroscopy on solid ground. Aims. The aim is to constrain the TP profi…
▽ More
Context. Observationally constraining the atmospheric temperature-pressure (TP) profile of exoplanets is an important step forward for improving planetary atmosphere models, further enabling one to place the detection of spectral features and the measurement of atomic and molecular abundances through transmission and emission spectroscopy on solid ground. Aims. The aim is to constrain the TP profile of the ultra-hot Jupiter KELT-9b by fitting synthetic spectra to the observed H$α$ and H$β$ lines and identify why self-consistent planetary TP models are unable to fit the observations. Methods. We construct 126 one-dimensional TP profiles varying the lower and upper atmospheric temperatures, as well as the location and gradient of the temperature rise. For each TP profile, we compute transmission spectra of the H$α$ and H$β$ lines employing the Cloudy radiative transfer code, which self-consistently accounts for non-local thermodynamic equilibrium (NLTE) effects. Results. The TP profiles leading to best fit the observations are characterised by an upper atmospheric temperature of 10000-11000 K and by an inverted temperature profile at pressures higher than 10$^{-4}$ bar. We find that the assumption of local thermodynamic equilibrium (LTE) leads to overestimate the level population of excited hydrogen by several orders of magnitude, and hence to significantly overestimate the strength of the Balmer lines. The chemical composition of the best fitting models indicate that the high upper atmospheric temperature is most likely driven by metal photoionisation and that FeII and FeIII have comparable abundances at pressures lower than 10$^{-6}$ bar, possibly making the latter detectable. Conclusions. Modelling the atmospheres of ultra-hot Jupiters requires one to account for metal photoionisation. [abridged]
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
Brouwer's satellite solution redux
Authors:
Martin Lara
Abstract:
Brouwer's solution to the artificial satellite problem is revisited to show that the complete Hamiltonian reduction is rather achieved in the plain Poincaré's style, through a single canonical transformation, than using a sequence of partial reductions based on von Zeipel's alternative for dealing with perturbed degenerate Hamiltonian systems.
Brouwer's solution to the artificial satellite problem is revisited to show that the complete Hamiltonian reduction is rather achieved in the plain Poincaré's style, through a single canonical transformation, than using a sequence of partial reductions based on von Zeipel's alternative for dealing with perturbed degenerate Hamiltonian systems.
△ Less
Submitted 1 November, 2020; v1 submitted 22 September, 2020;
originally announced September 2020.