-
Mode-based estimation of the center of symmetry
Authors:
José E. Chacón,
Javier Fernández Serrano
Abstract:
In the mean-median-mode triad of univariate centrality measures, the mode has been overlooked for estimating the center of symmetry in continuous and unimodal settings. This paper expands on the connection between kernel mode estimators and M-estimators for location, bridging the gap between the nonparametrics and robust statistics communities. The variance of modal estimators is studied in terms…
▽ More
In the mean-median-mode triad of univariate centrality measures, the mode has been overlooked for estimating the center of symmetry in continuous and unimodal settings. This paper expands on the connection between kernel mode estimators and M-estimators for location, bridging the gap between the nonparametrics and robust statistics communities. The variance of modal estimators is studied in terms of a bandwidth parameter, establishing conditions for an optimal solution that outperforms the household sample mean. A purely nonparametric approach is adopted, modeling heavy-tailedness through regular variation. The results lead to an estimator proposal that includes a novel one-parameter family of kernels with compact support, offering extra robustness and efficiency. The effectiveness and versatility of the new method are demonstrated in a real-world case study and a thorough simulation study, comparing favorably to traditional and more competitive alternatives. Several myths about the mode are clarified along the way, reopening the quest for flexible and efficient nonparametric estimators.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
The SCIP Optimization Suite 9.0
Authors:
Suresh Bolusani,
Mathieu Besançon,
Ksenia Bestuzheva,
Antonia Chmiela,
João Dionísio,
Tim Donkiewicz,
Jasper van Doornmalen,
Leon Eifler,
Mohammed Ghannam,
Ambros Gleixner,
Christoph Graczyk,
Katrin Halbig,
Ivo Hedtke,
Alexander Hoen,
Christopher Hojny,
Rolf van der Hulst,
Dominik Kamp,
Thorsten Koch,
Kevin Kofler,
Jurgen Lentz,
Julian Manns,
Gioni Mexi,
Erik Mühmer,
Marc E. Pfetsch,
Franziska Schlösser
, et al. (6 additional authors not shown)
Abstract:
The SCIP Optimization Suite provides a collection of software packages for mathematical optimization, centered around the constraint integer programming (CIP) framework SCIP. This report discusses the enhancements and extensions included in the SCIP Optimization Suite 9.0. The updates in SCIP 9.0 include improved symmetry handling, additions and improvements of nonlinear handlers and primal heuris…
▽ More
The SCIP Optimization Suite provides a collection of software packages for mathematical optimization, centered around the constraint integer programming (CIP) framework SCIP. This report discusses the enhancements and extensions included in the SCIP Optimization Suite 9.0. The updates in SCIP 9.0 include improved symmetry handling, additions and improvements of nonlinear handlers and primal heuristics, a new cut generator and two new cut selection schemes, a new branching rule, a new LP interface, and several bug fixes. The SCIP Optimization Suite 9.0 also features new Rust and C++ interfaces for SCIP, new Python interface for SoPlex, along with enhancements to existing interfaces. The SCIP Optimization Suite 9.0 also includes new and improved features in the LP solver SoPlex, the presolving library PaPILO, the parallel framework UG, the decomposition framework GCG, and the SCIP extension SCIP-SDP. These additions and enhancements have resulted in an overall performance improvement of SCIP in terms of solving time, number of nodes in the branch-and-bound tree, as well as the reliability of the solver.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Bayesian taut splines for estimating the number of modes
Authors:
José E. Chacón,
Javier Fernández Serrano
Abstract:
The number of modes in a probability density function is representative of the complexity of a model and can also be viewed as the number of subpopulations. Despite its relevance, there has been limited research in this area. A novel approach to estimating the number of modes in the univariate setting is presented, focusing on prediction accuracy and inspired by some overlooked aspects of the prob…
▽ More
The number of modes in a probability density function is representative of the complexity of a model and can also be viewed as the number of subpopulations. Despite its relevance, there has been limited research in this area. A novel approach to estimating the number of modes in the univariate setting is presented, focusing on prediction accuracy and inspired by some overlooked aspects of the problem: the need for structure in the solutions, the subjective and uncertain nature of modes, and the convenience of a holistic view that blends local and global density properties. The technique combines flexible kernel estimators and parsimonious compositional splines in the Bayesian inference paradigm, providing soft solutions and incorporating expert judgment. The procedure includes feature exploration, model selection, and mode testing, illustrated in a sports analytics case study showcasing multiple companion visualisation tools. A thorough simulation study also demonstrates that traditional modality-driven approaches paradoxically struggle to provide accurate results. In this context, the new method emerges as a top-tier alternative, offering innovative solutions for analysts.
△ Less
Submitted 8 May, 2024; v1 submitted 11 July, 2023;
originally announced July 2023.
-
Primordial black hole superradiance and evaporation in the string axiverse
Authors:
Marco Calzà,
João G. Rosa,
Filipe Serrano
Abstract:
In the string axiverse scenario, light primordial black holes may spin up due to the Hawking emission of a large number of light (sub-MeV) axions. We show that this may trigger superradiant instabilities associated with a heavier axion during the black holes' evolution, and study the coupled dynamics of superradiance and evaporation. We find, in particular, that the present black hole mass-spin di…
▽ More
In the string axiverse scenario, light primordial black holes may spin up due to the Hawking emission of a large number of light (sub-MeV) axions. We show that this may trigger superradiant instabilities associated with a heavier axion during the black holes' evolution, and study the coupled dynamics of superradiance and evaporation. We find, in particular, that the present black hole mass-spin distribution should follow the superradiance threshold condition for black hole masses below the value at which the superradiant cloud forms, for a given heavy axion mass. Furthermore, we show that the decay of the heavy axions within the superradiant cloud into photon pairs may lead to a distinctive line in the black hole's emission spectrum, superimposed on its electromagnetic Hawking emission.
△ Less
Submitted 15 June, 2023;
originally announced June 2023.
-
Enabling Research through the SCIP Optimization Suite 8.0
Authors:
Ksenia Bestuzheva,
Mathieu Besançon,
Wei-Kun Chen,
Antonia Chmiela,
Tim Donkiewicz,
Jasper van Doornmalen,
Leon Eifler,
Oliver Gaul,
Gerald Gamrath,
Ambros Gleixner,
Leona Gottwald,
Christoph Graczyk,
Katrin Halbig,
Alexander Hoen,
Christopher Hojny,
Rolf van der Hulst,
Thorsten Koch,
Marco Lübbecke,
Stephen J. Maher,
Frederic Matter,
Erik Mühmer,
Benjamin Müller,
Marc E. Pfetsch,
Daniel Rehfeldt,
Steffan Schlein
, et al. (10 additional authors not shown)
Abstract:
The SCIP Optimization Suite provides a collection of software packages for mathematical optimization centered around the constraint integer programming framework SCIP. The focus of this paper is on the role of the SCIP Optimization Suite in supporting research. SCIP's main design principles are discussed, followed by a presentation of the latest performance improvements and developments in version…
▽ More
The SCIP Optimization Suite provides a collection of software packages for mathematical optimization centered around the constraint integer programming framework SCIP. The focus of this paper is on the role of the SCIP Optimization Suite in supporting research. SCIP's main design principles are discussed, followed by a presentation of the latest performance improvements and developments in version 8.0, which serve both as examples of SCIP's application as a research tool and as a platform for further developments. Further, the paper gives an overview of interfaces to other programming and modeling languages, new features that expand the possibilities for user interaction with the framework, and the latest developments in several extensions built upon SCIP.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
Global Optimization of Mixed-Integer Nonlinear Programs with SCIP 8
Authors:
Ksenia Bestuzheva,
Antonia Chmiela,
Benjamin Müller,
Felipe Serrano,
Stefan Vigerske,
Fabian Wegscheider
Abstract:
For over ten years, the constraint integer programming framework SCIP has been extended by capabilities for the solution of convex and nonconvex mixed-integer nonlinear programs (MINLPs). With the recently published version 8.0, these capabilities have been largely reworked and extended. This paper discusses the motivations for recent changes and provides an overview of features that are particula…
▽ More
For over ten years, the constraint integer programming framework SCIP has been extended by capabilities for the solution of convex and nonconvex mixed-integer nonlinear programs (MINLPs). With the recently published version 8.0, these capabilities have been largely reworked and extended. This paper discusses the motivations for recent changes and provides an overview of features that are particular to MINLP solving in SCIP. Further, difficulties in benchmarking global MINLP solvers are discussed and a comparison with several state-of-the-art global MINLP solvers is provided.
△ Less
Submitted 2 January, 2023;
originally announced January 2023.
-
A characterization of maximal homogeneous-quadratic-free sets
Authors:
Gonzalo Muñoz,
Joseph Paat,
Felipe Serrano
Abstract:
The intersection cut framework was introduced by Balas in 1971 as a method for generating cutting planes in integer optimization. In this framework, one uses a full-dimensional convex $S$-free set, where $S$ is the feasible region of the integer program, to derive a cut separating $S$ from a non-integral vertex of a linear relaxation of $S$. Among all $S$-free sets, it is the inclusion-wise maxima…
▽ More
The intersection cut framework was introduced by Balas in 1971 as a method for generating cutting planes in integer optimization. In this framework, one uses a full-dimensional convex $S$-free set, where $S$ is the feasible region of the integer program, to derive a cut separating $S$ from a non-integral vertex of a linear relaxation of $S$. Among all $S$-free sets, it is the inclusion-wise maximal ones that yield the strongest cuts. Recently, this framework has been extended beyond the integer case in order to obtain cutting planes in non-linear settings. In this work, we consider the specific setting when $S$ is defined by a homogeneous quadratic inequality. In this 'quadratic-free' setting, every function $Γ: D^m \to D^n$, where $D^k$ is the unit disk in $\mathbb{R}^k$, generates a representation of a quadratic-free set. While not every $Γ$ generates a maximal quadratic free set, it is the case that every full-dimensional maximal quadratic free set is generated by some $Γ$. Our main result shows that the corresponding quadratic-free set is full-dimensional and maximal if and only if $Γ$ is non-expansive and satisfies a technical condition. This result yields a broader class of maximal $S$-free sets than previously known. Our result stems from a new characterization of maximal $S$-free sets (for general $S$ beyond the quadratic setting) based on sequences that 'expose' inequalities defining the $S$-free set.
△ Less
Submitted 26 September, 2023; v1 submitted 9 November, 2022;
originally announced November 2022.
-
Bump hunting through density curvature features
Authors:
José E. Chacón,
Javier Fernández Serrano
Abstract:
Bump hunting deals with finding in sample spaces meaningful data subsets known as bumps. These have traditionally been conceived as modal or concave regions in the graph of the underlying density function. We define an abstract bump construct based on curvature functionals of the probability density. Then, we explore several alternative characterizations involving derivatives up to second order. I…
▽ More
Bump hunting deals with finding in sample spaces meaningful data subsets known as bumps. These have traditionally been conceived as modal or concave regions in the graph of the underlying density function. We define an abstract bump construct based on curvature functionals of the probability density. Then, we explore several alternative characterizations involving derivatives up to second order. In particular, a suitable implementation of Good and Gaskins' original concave bumps is proposed in the multivariate case. Moreover, we bring to exploratory data analysis concepts like the mean curvature and the Laplacian that have produced good results in applied domains. Our methodology addresses the approximation of the curvature functional with a plug-in kernel density estimator. We provide theoretical results that assure the asymptotic consistency of bump boundaries in the Hausdorff distance with affordable convergence rates. We also present asymptotically valid and consistent confidence regions bounding curvature bumps. The theory is illustrated through several use cases in sports analytics with datasets from the NBA, MLB and NFL. We conclude that the different curvature instances effectively combine to generate insightful visualizations.
△ Less
Submitted 8 May, 2024; v1 submitted 30 July, 2022;
originally announced August 2022.
-
Adaptive Cut Selection in Mixed-Integer Linear Programming
Authors:
Mark Turner,
Thorsten Koch,
Felipe Serrano,
Michael Winkler
Abstract:
Cutting plane selection is a subroutine used in all modern mixed-integer linear programming solvers with the goal of selecting a subset of generated cuts that induce optimal solver performance. These solvers have millions of parameter combinations, and so are excellent candidates for parameter tuning. Cut selection scoring rules are usually weighted sums of different measurements, where the weight…
▽ More
Cutting plane selection is a subroutine used in all modern mixed-integer linear programming solvers with the goal of selecting a subset of generated cuts that induce optimal solver performance. These solvers have millions of parameter combinations, and so are excellent candidates for parameter tuning. Cut selection scoring rules are usually weighted sums of different measurements, where the weights are parameters. We present a parametric family of mixed-integer linear programs together with infinitely many family-wide valid cuts. Some of these cuts can induce integer optimal solutions directly after being applied, while others fail to do so even if an infinite amount are applied. We show for a specific cut selection rule, that any finite grid search of the parameter space will always miss all parameter values, which select integer optimal inducing cuts in an infinite amount of our problems. We propose a variation on the design of existing graph convolutional neural networks, adapting them to learn cut selection rule parameters. We present a reinforcement learning framework for selecting cuts, and train our design using said framework over MIPLIB 2017 and a neural network verification data set. Our framework and design show that adaptive cut selection does substantially improve performance over a diverse set of instances, but that finding a single function describing such a rule is difficult. Code for reproducing all experiments is available at https://github.com/Opt-Mucca/Adaptive-Cutsel-MILP.
△ Less
Submitted 23 February, 2023; v1 submitted 22 February, 2022;
originally announced February 2022.
-
The SCIP Optimization Suite 8.0
Authors:
Ksenia Bestuzheva,
Mathieu Besançon,
Wei-Kun Chen,
Antonia Chmiela,
Tim Donkiewicz,
Jasper van Doornmalen,
Leon Eifler,
Oliver Gaul,
Gerald Gamrath,
Ambros Gleixner,
Leona Gottwald,
Christoph Graczyk,
Katrin Halbig,
Alexander Hoen,
Christopher Hojny,
Rolf van der Hulst,
Thorsten Koch,
Marco Lübbecke,
Stephen J. Maher,
Frederic Matter,
Erik Mühmer,
Benjamin Müller,
Marc E. Pfetsch,
Daniel Rehfeldt,
Steffan Schlein
, et al. (10 additional authors not shown)
Abstract:
The SCIP Optimization Suite provides a collection of software packages for mathematical optimization centered around the constraint integer programming framework SCIP. This paper discusses enhancements and extensions contained in version 8.0 of the SCIP Optimization Suite. Major updates in SCIP include improvements in symmetry handling and decomposition algorithms, new cutting planes, a new plugin…
▽ More
The SCIP Optimization Suite provides a collection of software packages for mathematical optimization centered around the constraint integer programming framework SCIP. This paper discusses enhancements and extensions contained in version 8.0 of the SCIP Optimization Suite. Major updates in SCIP include improvements in symmetry handling and decomposition algorithms, new cutting planes, a new plugin type for cut selection, and a complete rework of the way nonlinear constraints are handled. Additionally, SCIP 8.0 now supports interfaces for Julia as well as Matlab. Further, UG now includes a unified framework to parallelize all solvers, a utility to analyze computational experiments has been added to GCG, dual solutions can be postsolved by PaPILO, new heuristics and presolving methods were added to SCIP-SDP, and additional problem classes and major performance improvements are available in SCIP-Jack.
△ Less
Submitted 16 December, 2021;
originally announced December 2021.
-
Semiparametric bivariate extreme-value copulas
Authors:
Javier Fernández Serrano
Abstract:
Extreme-value copulas arise as the limiting dependence structure of component-wise maxima. Defined in terms of a functional parameter, they are one of the most widespread copula families due to their flexibility and ability to capture asymmetry. Despite this, meeting the complex analytical properties of this parameter in an unconstrained setting remains a challenge, restricting most uses to models…
▽ More
Extreme-value copulas arise as the limiting dependence structure of component-wise maxima. Defined in terms of a functional parameter, they are one of the most widespread copula families due to their flexibility and ability to capture asymmetry. Despite this, meeting the complex analytical properties of this parameter in an unconstrained setting remains a challenge, restricting most uses to models with very few parameters or nonparametric models. Focusing on the bivariate case, we propose a novel semiparametric approach. Our procedure relies on a series of transformations, including Williamson's transform and starting from a zero-integral spline. Without further constraints, wholly compliant solutions can be efficiently obtained through maximum likelihood estimation, leveraging gradient optimization. We successfully conducted several experiments on simulated and real-world data. Our method outperforms another well-known nonparametric technique over small and medium-sized samples in various settings. Its expressiveness is illustrated with precious data gathered by the gravitational wave detection LIGO and Virgo collaborations.
△ Less
Submitted 24 March, 2022; v1 submitted 23 September, 2021;
originally announced September 2021.
-
On obtaining the convex hull of quadratic inequalities via aggregations
Authors:
Santanu S. Dey,
Gonzalo Munoz,
Felipe Serrano
Abstract:
A classical approach for obtaining valid inequalities for a set involves weighted aggregations of the inequalities that describe such set. When the set is described by linear inequalities, thanks to the Farkas lemma, we know that every valid inequality can be obtained using aggregations. When the inequalities describing the set are two quadratics, Yildiran showed that the convex hull of the set is…
▽ More
A classical approach for obtaining valid inequalities for a set involves weighted aggregations of the inequalities that describe such set. When the set is described by linear inequalities, thanks to the Farkas lemma, we know that every valid inequality can be obtained using aggregations. When the inequalities describing the set are two quadratics, Yildiran showed that the convex hull of the set is given by at most two aggregated inequalities. In this work, we study the case of a set described by three or more quadratic inequalities. We show that, under technical assumptions, the convex hull of a set described by three quadratic inequalities can be obtained via (potentially infinitely many) aggregated inequalities. We also show, through counterexamples, that it is unlikely to have a similar result if either the technical conditions are relaxed, or if we consider four or more inequalities.
△ Less
Submitted 23 June, 2021;
originally announced June 2021.
-
On Generalized Surrogate Duality in Mixed-Integer Nonlinear Programming
Authors:
Benjamin Müller,
Gonzalo Muñoz,
Maxime Gasse,
Ambros Gleixner,
Andrea Lodi,
Felipe Serrano
Abstract:
The most important ingredient for solving mixed-integer nonlinear programs (MINLPs) to global epsilon-optimality with spatial branch and bound is a tight, computationally tractable relaxation. Due to both theoretical and practical considerations, relaxations of MINLPs are usually required to be convex. Nonetheless, current optimization solver can often successfully handle a moderate presence of no…
▽ More
The most important ingredient for solving mixed-integer nonlinear programs (MINLPs) to global epsilon-optimality with spatial branch and bound is a tight, computationally tractable relaxation. Due to both theoretical and practical considerations, relaxations of MINLPs are usually required to be convex. Nonetheless, current optimization solver can often successfully handle a moderate presence of nonconvexities, which opens the door for the use of potentially tighter nonconvex relaxations. In this work, we exploit this fact and make use of a nonconvex relaxation obtained via aggregation of constraints: a surrogate relaxation. These relaxations were actively studied for linear integer programs in the 70s and 80s, but they have been scarcely considered since. We revisit these relaxations in an MINLP setting and show the computational benefits and challenges they can have. Additionally, we study a generalization of such relaxation that allows for multiple aggregations simultaneously and present the first algorithm that is capable of computing the best set of aggregations. We propose a multitude of computational enhancements for improving its practical performance and evaluate the algorithm's ability to generate strong dual bounds through extensive computational experiments.
△ Less
Submitted 1 December, 2019;
originally announced December 2019.
-
Maximal quadratic-free sets
Authors:
Gonzalo Muñoz,
Felipe Serrano
Abstract:
The intersection cut paradigm is a powerful framework that facilitates the generation of valid linear inequalities, or cutting planes, for a potentially complex set S. The key ingredients in this construction are a simplicial conic relaxation of S and an S-free set: a convex zone whose interior does not intersect S. Ideally, such S-free set would be maximal inclusion-wise, as it would generate a d…
▽ More
The intersection cut paradigm is a powerful framework that facilitates the generation of valid linear inequalities, or cutting planes, for a potentially complex set S. The key ingredients in this construction are a simplicial conic relaxation of S and an S-free set: a convex zone whose interior does not intersect S. Ideally, such S-free set would be maximal inclusion-wise, as it would generate a deeper cutting plane. However, maximality can be a challenging goal in general. In this work, we show how to construct maximal S-free sets when S is defined as a general quadratic inequality. Our maximal S-free sets are such that efficient separation of a vertex in LP-based approaches to quadratically constrained problems is guaranteed. To the best of our knowledge, this work is the first to provide maximal quadratic-free sets.
△ Less
Submitted 29 November, 2019; v1 submitted 27 November, 2019;
originally announced November 2019.
-
Visible points, the separation problem, and applications to MINLP
Authors:
Felipe Serrano
Abstract:
In this paper we introduce a technique to produce tighter cutting planes for mixed-integer non-linear programs. Usually, a cutting plane is generated to cut off a specific infeasible point. The underlying idea is to use the infeasible point to restrict the feasible region in order to obtain a tighter domain. To ensure validity, we require that every valid cut separating the infeasible point from t…
▽ More
In this paper we introduce a technique to produce tighter cutting planes for mixed-integer non-linear programs. Usually, a cutting plane is generated to cut off a specific infeasible point. The underlying idea is to use the infeasible point to restrict the feasible region in order to obtain a tighter domain. To ensure validity, we require that every valid cut separating the infeasible point from the restricted feasible region is still valid for the original feasible region. We translate this requirement in terms of the separation problem and the reverse polar. In particular, if the reverse polar of the restricted feasible region is the same as the reverse polar of the feasible region, then any cut valid for the restricted feasible region that separates the infeasible point, is valid for the feasible region. We show that the reverse polar of the visible points of the feasible region from the infeasible point coincides with the reverse polar of the feasible region. In the special where the feasible region is described by a single non-convex constraint intersected with a convex set we provide a characterization of the visible points. Furthermore, when the non-convex constraint is quadratic the characterization is particularly simple. We also provide an extended formulation for a relaxation of the visible points when the non-convex constraint is a general polynomial. Finally, we give some conditions under which for a given set there is an inclusion-wise smallest set, in some predefined family of sets, whose reverse polars coincide.
△ Less
Submitted 18 July, 2019;
originally announced July 2019.
-
On the Relation between the Extended Supporting Hyperplane Algorithm and Kelley's Cutting Plane Algorithm
Authors:
Felipe Serrano,
Robert Schwarz,
Ambros Gleixner
Abstract:
Recently, Kronqvist et al.~\cite{KronqvistLundellWesterlund2016} rediscovered the supporting hyperplane algorithm of Veinott~\cite{Veinott1967} and demonstrated its computational benefits for solving convex mixed-integer nonlinear programs.
In this paper we derive the algorithm from a geometric point of view.
This enables us to show that the supporting hyperplane algorithm is equivalent to Kel…
▽ More
Recently, Kronqvist et al.~\cite{KronqvistLundellWesterlund2016} rediscovered the supporting hyperplane algorithm of Veinott~\cite{Veinott1967} and demonstrated its computational benefits for solving convex mixed-integer nonlinear programs.
In this paper we derive the algorithm from a geometric point of view.
This enables us to show that the supporting hyperplane algorithm is equivalent to Kelley's cutting plane algorithm~\cite{J.E.Kelley1960} applied to a particular reformulation of the problem.
As a result, we extend the applicability of the supporting hyperplane algorithm to convex problems represented by general, not necessarily convex, differentiable functions that satisfy a mild condition.
△ Less
Submitted 20 May, 2019;
originally announced May 2019.
-
Using two-dimensional Projections for Stronger Separation and Propagation of Bilinear Terms
Authors:
Benjamin Müller,
Felipe Serrano,
Ambros Gleixner
Abstract:
One of the most fundamental ingredients in mixed-integer nonlinear programming solvers is the well-known McCormick relaxation for a product of two variables x and y over a box-constrained domain. The starting point of this paper is the fact that the convex hull of the graph of xy can be much tighter when computed over a strict, non-rectangular subset of the box. In order to exploit this in practic…
▽ More
One of the most fundamental ingredients in mixed-integer nonlinear programming solvers is the well-known McCormick relaxation for a product of two variables x and y over a box-constrained domain. The starting point of this paper is the fact that the convex hull of the graph of xy can be much tighter when computed over a strict, non-rectangular subset of the box. In order to exploit this in practice, we propose to compute valid linear inequalities for the projection of the feasible region onto the x-y-space by solving a sequence of linear programs akin to optimization-based bound tightening. These valid inequalities allow us to employ results from the literature to strengthen the classical McCormick relaxation. As a consequence, we obtain a stronger convexification procedure that exploits problem structure and can benefit from supplementary information obtained during the branch-and-bound algorithm such as an objective cutoff. We complement this by a new bound tightening procedure that efficiently computes the best possible bounds for x, y, and xy over the available projections. Our computational evaluation using the academic solver SCIP exhibit that the proposed methods are applicable to a large portion of the public test library MINLPLib and help to improve performance significantly.
△ Less
Submitted 10 January, 2020; v1 submitted 13 March, 2019;
originally announced March 2019.
-
Intersection cuts for factorable MINLP
Authors:
Felipe Serrano
Abstract:
Given a factorable function f, we propose a procedure that constructs a concave underestimator of f that is tight at a given point. These underestimators can be used to generate intersection cuts. A peculiarity of these underestimators is that they do not rely on a bounded domain. We propose a strengthening procedure for the intersection cuts that exploits the bounds of the domain. Finally, we pro…
▽ More
Given a factorable function f, we propose a procedure that constructs a concave underestimator of f that is tight at a given point. These underestimators can be used to generate intersection cuts. A peculiarity of these underestimators is that they do not rely on a bounded domain. We propose a strengthening procedure for the intersection cuts that exploits the bounds of the domain. Finally, we propose an extension of monoidal strengthening to take advantage of the integrality of the non-basic variables.
△ Less
Submitted 7 December, 2018;
originally announced December 2018.
-
Semi-exact solutions of the Razavy potential
Authors:
Qian Dong,
F. A. Serrano,
Guo-Hua Sun,
Jian Jing,
Shi-Hai Dong
Abstract:
In this work we study the quantum system with the symmetric Razavy potential and show how to find its exact solutions. We find that the solutions are given by the confluent Heun functions. The eigenvalues have to be calculated numerically. The properties of the wave functions depending on $m$ are illustrated graphically for a given potential parameter $ξ$. We find that the even and odd wave functi…
▽ More
In this work we study the quantum system with the symmetric Razavy potential and show how to find its exact solutions. We find that the solutions are given by the confluent Heun functions. The eigenvalues have to be calculated numerically. The properties of the wave functions depending on $m$ are illustrated graphically for a given potential parameter $ξ$. We find that the even and odd wave functions with definite parity are changed to odd and even wave functions when the potential parameter $m$ increases. This arises from the fact that the parity, which is a defined symmetry for very small $m$, is completely violated for large $m$. We also notice that the energy levels $ε_{i}$ decrease with the increasing potential parameter $m$.
△ Less
Submitted 16 May, 2018;
originally announced May 2018.
-
Fisher's information for the position-dependent mass Schrödinger system
Authors:
B. J. Falaye,
F. A. Serrano,
Shi-Hai Dong
Abstract:
This study presents the Fisher information for the position-dependent mass Schrödinger equation with hyperbolical potential {$V(x)=-V_0{\rm csch}^2(ax)$}. The analysis of the quantum-mechanical probability for the ground and exited states $(n=0, 1, 2)$ has been obtained via the Fisher's information. This controls both chemical and physical properties of some molecular systems. The Fisher informati…
▽ More
This study presents the Fisher information for the position-dependent mass Schrödinger equation with hyperbolical potential {$V(x)=-V_0{\rm csch}^2(ax)$}. The analysis of the quantum-mechanical probability for the ground and exited states $(n=0, 1, 2)$ has been obtained via the Fisher's information. This controls both chemical and physical properties of some molecular systems. The Fisher information is considered only for $x>0$ due to the singular point at $x=0$. We found that Fisher-information-based uncertainty relation and the Cramer-Rao inequality holds. Some relevant numerical results are presented. The results presented shows that the Cramer-Rao and the Heisenberg products in both spaces provide a natural measure for anharmonicity of {$-V_0{\rm csch}^2(ax)$}
△ Less
Submitted 29 September, 2015;
originally announced September 2015.
-
Structural and Dynamical Patterns on Online Social Networks: the Spanish May 15th Movement as a case study
Authors:
Javier Borge-Holthoefer,
Alejandro Rivero,
Iñigo García,
Elisa Cauhé,
Alfredo Ferrer,
Darío Ferrer,
David Francos,
David Iñiguez,
María Pilar Pérez,
Gonzalo Ruiz,
Francisco Sanz,
Fermín Serrano,
Cristina Viñas,
Alfonso Tarancón,
Yamir Moreno
Abstract:
The number of people using online social networks in their everyday life is continuously growing at a pace never saw before. This new kind of communication has an enormous impact on opinions, cultural trends, information spreading and even in the commercial success of new products. More importantly, social online networks have revealed as a fundamental organizing mechanism in recent country-wide s…
▽ More
The number of people using online social networks in their everyday life is continuously growing at a pace never saw before. This new kind of communication has an enormous impact on opinions, cultural trends, information spreading and even in the commercial success of new products. More importantly, social online networks have revealed as a fundamental organizing mechanism in recent country-wide social movements. In this paper, we provide a quantitative analysis of the structural and dynamical patterns emerging from the activity of an online social network around the ongoing May 15th (15M) movement in Spain. Our network is made up by users that exchanged tweets in a time period of one month, which includes the birth and stabilization of the 15M movement. We characterize in depth the growth of such dynamical network and find that it is scale-free with communities at the mesoscale. We also find that its dynamics exhibits typical features of critical systems such as robustness and power-law distributions for several quantities. Remarkably, we report that the patterns characterizing the spreading dynamics are asymmetric, giving rise to a clear distinction between information sources and sinks. Our study represent a first step towards the use of data from online social media to comprehend modern societal dynamics.
△ Less
Submitted 8 July, 2011;
originally announced July 2011.
-
Threefolds with nef anticanonical bundles
Authors:
Thomas Peternell,
Fernando Serrano
Abstract:
We prove that the Albanese map of a smooth projective threefold, whose anticanonical bundle is nef, is a surjective submersion. We also investigate morphisms of threefolds to curves and surfaces whose relative anticanonical bundle are nef.
We prove that the Albanese map of a smooth projective threefold, whose anticanonical bundle is nef, is a surjective submersion. We also investigate morphisms of threefolds to curves and surfaces whose relative anticanonical bundle are nef.
△ Less
Submitted 24 February, 1998; v1 submitted 5 February, 1998;
originally announced February 1998.
-
Development of a Spanish Version of the Xerox Tagger
Authors:
Fernando Sánchez León,
Amalio F. Nieto Serrano
Abstract:
This paper describes work performed withing the CRATER ({\em C}orpus {\em R}esources {\em A}nd {\em T}erminology {\em E}xt{\em R}action, MLAP-93/20) project, funded by the Commission of the European Communities. In particular, it addresses the issue of adapting the Xerox Tagger to Spanish in order to tag the Spanish version of the ITU (International Telecommunications Union) corpus. The model im…
▽ More
This paper describes work performed withing the CRATER ({\em C}orpus {\em R}esources {\em A}nd {\em T}erminology {\em E}xt{\em R}action, MLAP-93/20) project, funded by the Commission of the European Communities. In particular, it addresses the issue of adapting the Xerox Tagger to Spanish in order to tag the Spanish version of the ITU (International Telecommunications Union) corpus. The model implemented by this tagger is briefly presented along with some modifications performed on it in order to use some parameters not probabilistically estimated. Initial decisions, like the tagset, the lexicon and the training corpus are also discussed. Finally, results are presented and the benefits of the {\em mixed model} justified.
△ Less
Submitted 19 May, 1995;
originally announced May 1995.