-
Exact discovery is polynomial for sparse causal Bayesian networks
Authors:
Felix L. Rios,
Giusi Moffa,
Jack Kuipers
Abstract:
Causal Bayesian networks are widely used tools for summarising the dependencies between variables and elucidating their putative causal relationships. Learning networks from data is computationally hard in general. The current state-of-the-art approaches for exact causal discovery are integer linear programming over the underlying space of directed acyclic graphs, dynamic programming and shortest-…
▽ More
Causal Bayesian networks are widely used tools for summarising the dependencies between variables and elucidating their putative causal relationships. Learning networks from data is computationally hard in general. The current state-of-the-art approaches for exact causal discovery are integer linear programming over the underlying space of directed acyclic graphs, dynamic programming and shortest-path searches over the space of topological orders, and constraint programming combining both. For dynamic programming over orders, the computational complexity is known to be exponential base 2 in the number of variables in the network. We demonstrate how to use properties of Bayesian networks to prune the search space and lower the computational cost, while still guaranteeing exact discovery. When including new path-search and divide-and-conquer criteria, we prove optimality in quadratic time for matchings, and polynomial time for any network class with logarithmically-bound largest connected components. In simulation studies we observe the polynomial dependence for sparse networks and that, beyond some critical value, the logarithm of the base grows with the network density. Our approach then out-competes the state-of-the-art at lower densities. These results therefore pave the way for faster exact causal discovery in larger and sparser networks.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
An upgraded 0.4-meter telescope fleet for Las Cumbres Observatory's Educational and Science Programs
Authors:
Daniel-Rolf Harbeck,
Brook Taylor,
Annie Kirby,
Mark Bowman,
Steve Foale,
Kal Kadlec,
Curtis McCully,
Matthew Daily,
Jon DeVera,
Dave Douglass,
Mark Willis,
Ian Baker,
Nikolaus Volgenau,
Patrick Conway,
Brian Haworth,
Jesus Estrada,
Edward Gomez,
Sandy Seale,
Alice Hopkinson,
Fernando Rios,
Prerana Kotapali,
Lisa Storrie-Lombardi,
Wayne Rosing
Abstract:
Las Cumbres Observatory (LCOGT) operates a global network of robotic 0.4, 1.0, and 2.0-meter telescopes to facilitate scientific research and education in time-domain astronomy. LCOGT's flagship educational program, Global Sky Partners (GSP), awards up to 1500 hours per year of telescope time to individuals and organizations that run their own, fully supported, educational programs. The GSP has a…
▽ More
Las Cumbres Observatory (LCOGT) operates a global network of robotic 0.4, 1.0, and 2.0-meter telescopes to facilitate scientific research and education in time-domain astronomy. LCOGT's flagship educational program, Global Sky Partners (GSP), awards up to 1500 hours per year of telescope time to individuals and organizations that run their own, fully supported, educational programs. The GSP has a presence in 40 countries and 45% of the Partners target under-served, under-represented, and developing world audiences. The degradation and obsolescence of the original 0.4-meter telescope network prompted LCOGT to update the fleet of 10 telescopes to a new system consisting of predominantly off-the-shelf products. New PlaneWave DeltaRho 350 telescopes with Gemini Focuser/Rotators, LCOGT filter wheels, and QHY600 CMOS cameras, complement the original, custom-built mount. The deployment of all ten telescopes was completed in March 2024. We describe the design and performance of this new system and its components. We comment on modifications made to the QHY600 cameras, as well as on the treatment of random telegraph noise of its CMOS detectors within our data processing system BANZAI. The new telescope network supports the GSP program as well as multiple key science projects, including follow-up observations for the TESS satellite mission.
△ Less
Submitted 16 May, 2024;
originally announced May 2024.
-
On a fibrational construction for optics, lenses, and Dialectica categories
Authors:
Matteo Capucci,
Bruno Gavranović,
Abdullah Malik,
Francisco Rios,
Jonathan Weinberger
Abstract:
Categories of lenses/optics and Dialectica categories are both comprised of bidirectional morphisms of basically the same form. In this work we show how they can be considered a special case of an overarching fibrational construction, generalizing Hofstra's construction of Dialectica fibrations and Spivak's construction of generalized lenses. This construction turns a tower of Grothendieck fibrati…
▽ More
Categories of lenses/optics and Dialectica categories are both comprised of bidirectional morphisms of basically the same form. In this work we show how they can be considered a special case of an overarching fibrational construction, generalizing Hofstra's construction of Dialectica fibrations and Spivak's construction of generalized lenses. This construction turns a tower of Grothendieck fibrations into another tower of fibrations by iteratively twisting each of the components, using the opposite fibration construction.
△ Less
Submitted 12 June, 2024; v1 submitted 24 March, 2024;
originally announced March 2024.
-
Scalable Structure Learning for Sparse Context-Specific Causal Systems
Authors:
Felix Leopoldo Rios,
Alex Markham,
Liam Solus
Abstract:
Several approaches to graphically representing context-specific relations among jointly distributed categorical variables have been proposed, along with structure learning algorithms. While existing optimization-based methods have limited scalability due to the large number of context-specific models, the constraint-based methods are more prone to error than even constraint-based DAG learning algo…
▽ More
Several approaches to graphically representing context-specific relations among jointly distributed categorical variables have been proposed, along with structure learning algorithms. While existing optimization-based methods have limited scalability due to the large number of context-specific models, the constraint-based methods are more prone to error than even constraint-based DAG learning algorithms since more relations must be tested. We present a hybrid algorithm for learning context-specific models that scales to hundreds of variables while testing no more constraints than standard DAG learning algorithms. Scalable learning is achieved through a combination of an order-based MCMC algorithm and sparsity assumptions analogous to those typically invoked for DAG models. To implement the method, we solve a special case of an open problem recently posed by Alon and Balogh. The method is shown to perform well on synthetic data and real world examples, in terms of both accuracy and scalability.
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
Optimal boundary regularity and a Hopf-type lemma for Dirichlet problems involving the logarithmic Laplacian
Authors:
Víctor Hernández-Santamaría,
Luis Fernando López Ríos,
Alberto Saldaña
Abstract:
We study the optimal boundary regularity of solutions to Dirichlet problems involving the logarithmic Laplacian. Our proofs are based on the construction of suitable barriers via the Kelvin transform and direct computations. As applications of our results, we show a Hopf-type Lemma for nonnegative weak solutions and the uniqueness of solutions to some nonlinear problems.
We study the optimal boundary regularity of solutions to Dirichlet problems involving the logarithmic Laplacian. Our proofs are based on the construction of suitable barriers via the Kelvin transform and direct computations. As applications of our results, we show a Hopf-type Lemma for nonnegative weak solutions and the uniqueness of solutions to some nonlinear problems.
△ Less
Submitted 3 July, 2024; v1 submitted 31 January, 2024;
originally announced January 2024.
-
Adiabatic Quantum Support Vector Machines
Authors:
Prasanna Date,
Dong Jun Woun,
Kathleen Hamilton,
Eduardo A. Coello Perez,
Mayanka Chandra Shekhar,
Francisco Rios,
John Gounley,
In-Saeng Suh,
Travis Humble,
Georgia Tourassi
Abstract:
Adiabatic quantum computers can solve difficult optimization problems (e.g., the quadratic unconstrained binary optimization problem), and they seem well suited to train machine learning models. In this paper, we describe an adiabatic quantum approach for training support vector machines. We show that the time complexity of our quantum approach is an order of magnitude better than the classical ap…
▽ More
Adiabatic quantum computers can solve difficult optimization problems (e.g., the quadratic unconstrained binary optimization problem), and they seem well suited to train machine learning models. In this paper, we describe an adiabatic quantum approach for training support vector machines. We show that the time complexity of our quantum approach is an order of magnitude better than the classical approach. Next, we compare the test accuracy of our quantum approach against a classical approach that uses the Scikit-learn library in Python across five benchmark datasets (Iris, Wisconsin Breast Cancer (WBC), Wine, Digits, and Lambeq). We show that our quantum approach obtains accuracies on par with the classical approach. Finally, we perform a scalability study in which we compute the total training times of the quantum approach and the classical approach with increasing number of features and number of data points in the training dataset. Our scalability results show that the quantum approach obtains a 3.5--4.5 times speedup over the classical approach on datasets with many (millions of) features.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
Identifying overparameterization in Quantum Circuit Born Machines
Authors:
Andrea Delgado,
Francisco Rios,
Kathleen E. Hamilton
Abstract:
In machine learning, overparameterization is associated with qualitative changes in the empirical risk landscape, which can lead to more efficient training dynamics. For many parameterized models used in statistical learning, there exists a critical number of parameters, or model size, above which the model is constructed and trained in the overparameterized regime. There are many characteristics…
▽ More
In machine learning, overparameterization is associated with qualitative changes in the empirical risk landscape, which can lead to more efficient training dynamics. For many parameterized models used in statistical learning, there exists a critical number of parameters, or model size, above which the model is constructed and trained in the overparameterized regime. There are many characteristics of overparameterized loss landscapes. The most significant is the convergence of standard gradient descent to global or local minima of low loss. In this work, we study the onset of overparameterization transitions for quantum circuit Born machines, generative models that are trained using non-adversarial gradient-based methods. We observe that bounds based on numerical analysis are in general good lower bounds on the overparameterization transition. However, bounds based on the quantum circuit's algebraic structure are very loose upper bounds. Our results indicate that fully understanding the trainability of these models remains an open question.
△ Less
Submitted 10 July, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
On a generalized Cahn-Hilliard model with p-Laplacian
Authors:
Raffaele Folino,
Luis Fernando Lopez Rios,
Marta Strani
Abstract:
A generalized Cahn-Hilliard model in a bounded interval of the real line with no-flux boundary conditions is considered. The label "generalized" refers to the fact that we consider a concentration dependent mobility, the $p$-Laplace operator with $p>1$ and a double well potential of the form $F(u)=\frac{1}{2θ}|1-u^2|^θ$, with $θ>1$; these terms replace, respectively, the constant mobility, the lin…
▽ More
A generalized Cahn-Hilliard model in a bounded interval of the real line with no-flux boundary conditions is considered. The label "generalized" refers to the fact that we consider a concentration dependent mobility, the $p$-Laplace operator with $p>1$ and a double well potential of the form $F(u)=\frac{1}{2θ}|1-u^2|^θ$, with $θ>1$; these terms replace, respectively, the constant mobility, the linear Laplace operator and the $C^2$ potential satisfying $F"(\pm1)>0$, which are typical of the standard Cahn-Hilliard model. After investigating the associated stationary problem and highlighting the differences with the standard results, we focus the attention on the long time dynamics of solutions when $θ\geq p>1$. In the $critical$ $θ=p>1$, we prove $exponentially$ $slow$ $motion$ of profiles with a transition layer structure, thus extending the well know results of the standard model, where $θ=p=2$; conversely, in the $supercritical$ case $θ>p>1$, we prove $algebraic$ $slow$ $motion$ of layered profiles.
△ Less
Submitted 22 February, 2022;
originally announced February 2022.
-
Variational formulation for fractional hyperbolic problems in the theory of viscoelasticity
Authors:
Luis Fernando López Ríos,
Julián Bravo-Castillero
Abstract:
In this article a theoretical framework for problems involving fractional equations of hyperbolic type arising in the theory of viscoelasticity is presented. Based on the Galerkin method, a variational problem of the fractionary viscoelasticity is studied. An appropriate functional setting is introduced in order to establish the existence, uniqueness and a priori estimates for weak solutions. This…
▽ More
In this article a theoretical framework for problems involving fractional equations of hyperbolic type arising in the theory of viscoelasticity is presented. Based on the Galerkin method, a variational problem of the fractionary viscoelasticity is studied. An appropriate functional setting is introduced in order to establish the existence, uniqueness and a priori estimates for weak solutions. This framework is developed in close concordance with important physical quantities of the theory of viscoelasticity
△ Less
Submitted 19 August, 2021;
originally announced August 2021.
-
Benchpress: A Scalable and Versatile Workflow for Benchmarking Structure Learning Algorithms
Authors:
Felix L. Rios,
Giusi Moffa,
Jack Kuipers
Abstract:
Describing the relationship between the variables in a study domain and modelling the data generating mechanism is a fundamental problem in many empirical sciences. Probabilistic graphical models are one common approach to tackle the problem. Learning the graphical structure for such models is computationally challenging and a fervent area of current research with a plethora of algorithms being de…
▽ More
Describing the relationship between the variables in a study domain and modelling the data generating mechanism is a fundamental problem in many empirical sciences. Probabilistic graphical models are one common approach to tackle the problem. Learning the graphical structure for such models is computationally challenging and a fervent area of current research with a plethora of algorithms being developed. To facilitate the benchmarking of different methods, we present a novel Snakemake workflow, called Benchpress for producing scalable, reproducible, and platform-independent benchmarks of structure learning algorithms for probabilistic graphical models. Benchpress is interfaced via a simple JSON-file, which makes it accessible for all users, while the code is designed in a fully modular fashion to enable researchers to contribute additional methodologies. Benchpress currently provides an interface to a large number of state-of-the-art algorithms from libraries such as BDgraph, BiDAG, bnlearn, causal-learn, gCastle, GOBNILP, pcalg, r.blip, scikit-learn, TETRAD, and trilearn as well as a variety of methods for data generating models and performance evaluation. Alongside user-defined models and randomly generated datasets, the workflow also includes a number of standard datasets and graphical models from the literature, which may be included in a benchmarking study. We demonstrate the applicability of this workflow for learning Bayesian networks in five typical data scenarios. The source code and documentation is publicly available from http://benchpressdocs.readthedocs.io.
△ Less
Submitted 4 December, 2023; v1 submitted 8 July, 2021;
originally announced July 2021.
-
Long-time behavior of solutions to the generalized Allen-Cahn model with degenerate diffusivity
Authors:
Raffaele Folino,
Luis F. López Ríos,
Ramón G. Plaza
Abstract:
The generalized Allen-Cahn equation,
\[
u_t=\varepsilon^2(D(u)u_x)_x-\frac{\varepsilon^2}2D'(u)u_x^2-F'(u),
\]
with nonlinear diffusion, $D = D(u)$, and potential, $F = F(u)$, of the form
\[
D(u) = |1-u^2|^{m}, \quad \text{or} \quad D(u) = |1-u|^{m}, \quad m >1,
\]
and
\[
F(u)=\frac{1}{2n}|1-u^2|^{n}, \qquad n\geq2,
\]
respectively, is studied. These choices correspond to a…
▽ More
The generalized Allen-Cahn equation,
\[
u_t=\varepsilon^2(D(u)u_x)_x-\frac{\varepsilon^2}2D'(u)u_x^2-F'(u),
\]
with nonlinear diffusion, $D = D(u)$, and potential, $F = F(u)$, of the form
\[
D(u) = |1-u^2|^{m}, \quad \text{or} \quad D(u) = |1-u|^{m}, \quad m >1,
\]
and
\[
F(u)=\frac{1}{2n}|1-u^2|^{n}, \qquad n\geq2,
\]
respectively, is studied. These choices correspond to a reaction function that can be derived from a double well potential, and to a generalized degenerate diffusivity coefficient depending on the density $u$ that vanishes at one or at the two wells, $u = \pm 1$. It is shown that interface layer solutions that are equal to $\pm 1$ except at a finite number of thin transitions of width $\varepsilon$ persist for an either exponentially or algebraically long time, depending upon the interplay between the exponents $n$ and $m$. For that purpose, energy bounds for a renormalized effective energy potential of Ginzburg-Landau type are derived.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
Spectral stability of monotone traveling fronts for reaction diffusion-degenerate Nagumo equations
Authors:
J. Francisco Leyva,
Luis F. López Ríos,
Ramón G. Plaza
Abstract:
This paper establishes the spectral stability of monotone traveling front solutions for reaction-diffusion equations where the reaction function is of Nagumo (or bistable) type and with diffusivities which are density dependent and degenerate at zero (one of the equilibrium points of the reaction). Spectral stability is understood as the property that the spectrum of the linearized operator around…
▽ More
This paper establishes the spectral stability of monotone traveling front solutions for reaction-diffusion equations where the reaction function is of Nagumo (or bistable) type and with diffusivities which are density dependent and degenerate at zero (one of the equilibrium points of the reaction). Spectral stability is understood as the property that the spectrum of the linearized operator around the wave, acting on an exponentially weighted space, is contained in the complex half plane with non-positive real part. Three different types of monotone waves are studied: (i) stationary diffusion-degenerate fronts, connecting the two stable equilibria of the reaction; (ii) traveling diffusion-degenerate fronts connecting zero with the unstable equilibrium; and, (iii) non-degenerate fronts. In the first two cases, the degeneracy is responsible of the loss of hyperbolicity of the asymptotic coefficient matrices of the spectral problem at one of the end points, precluding the application of standard techniques to locate the essential spectrum. This difficulty is overcome with a suitable partition of the spectrum, a generalized convergence of operators technique, the analysis of singular (or Weyl) sequences and the use of energy estimates. The monotonicity of the fronts, as well as detailed descriptions of the decay structure of eigenfunctions on a case by case basis, are key ingredients to show that all traveling fronts under consideration are spectrally stable in a suitably chosen exponentially weighted $L^2$ energy space.
△ Less
Submitted 7 December, 2022; v1 submitted 16 September, 2019;
originally announced September 2019.
-
Holographic Code Rate
Authors:
Noah Bray-Ali,
David Chester,
Dugan Hammock,
Marcelo M. Amaral,
Klee Irwin,
Michael F. Rios
Abstract:
Holographic codes grown with perfect tensors on regular hyperbolic tessellations using an inflation rule protect quantum information stored in the bulk from errors on the boundary provided the code rate is less than one. Hyperbolic geometry bounds the holographic code rate and guarantees quantum error correction for codes grown with any inflation rule on all regular hyperbolic tessellations in a c…
▽ More
Holographic codes grown with perfect tensors on regular hyperbolic tessellations using an inflation rule protect quantum information stored in the bulk from errors on the boundary provided the code rate is less than one. Hyperbolic geometry bounds the holographic code rate and guarantees quantum error correction for codes grown with any inflation rule on all regular hyperbolic tessellations in a class whose size grows exponentially with the rank of the perfect tensors for rank five and higher. For the tile completion inflation rule, holographic triangle codes have code rate more than one but all others perform quantum error correction.
△ Less
Submitted 25 August, 2019;
originally announced August 2019.
-
Software Citation Implementation Challenges
Authors:
Daniel S. Katz,
Daina Bouquin,
Neil P. Chue Hong,
Jessica Hausman,
Catherine Jones,
Daniel Chivvis,
Tim Clark,
Mercè Crosas,
Stephan Druskat,
Martin Fenner,
Tom Gillespie,
Alejandra Gonzalez-Beltran,
Morane Gruenpeter,
Ted Habermann,
Robert Haines,
Melissa Harrison,
Edwin Henneken,
Lorraine Hwang,
Matthew B. Jones,
Alastair A. Kelly,
David N. Kennedy,
Katrin Leinweber,
Fernando Rios,
Carly B. Robinson,
Ilian Todorov
, et al. (2 additional authors not shown)
Abstract:
The main output of the FORCE11 Software Citation working group (https://www.force11.org/group/software-citation-working-group) was a paper on software citation principles (https://doi.org/10.7717/peerj-cs.86) published in September 2016. This paper laid out a set of six high-level principles for software citation (importance, credit and attribution, unique identification, persistence, accessibilit…
▽ More
The main output of the FORCE11 Software Citation working group (https://www.force11.org/group/software-citation-working-group) was a paper on software citation principles (https://doi.org/10.7717/peerj-cs.86) published in September 2016. This paper laid out a set of six high-level principles for software citation (importance, credit and attribution, unique identification, persistence, accessibility, and specificity) and discussed how they could be used to implement software citation in the scholarly community. In a series of talks and other activities, we have promoted software citation using these increasingly accepted principles. At the time the initial paper was published, we also provided guidance and examples on how to make software citable, though we now realize there are unresolved problems with that guidance. The purpose of this document is to provide an explanation of current issues impacting scholarly attribution of research software, organize updated implementation guidance, and identify where best practices and solutions are still needed.
△ Less
Submitted 21 May, 2019;
originally announced May 2019.
-
Sequential sampling of junction trees for decomposable graphs
Authors:
Jimmy Olsson,
Tetyana Pavlenko,
Felix L. Rios
Abstract:
The junction-tree representation provides an attractive structural property for organizing a decomposable graph. In this study, we present two novel stochastic algorithms, which we call the junction-tree expander and junction-tree collapser for sequential sampling of junction trees for decomposable graphs. We show that recursive application of the junction-tree expander, expanding incrementally th…
▽ More
The junction-tree representation provides an attractive structural property for organizing a decomposable graph. In this study, we present two novel stochastic algorithms, which we call the junction-tree expander and junction-tree collapser for sequential sampling of junction trees for decomposable graphs. We show that recursive application of the junction-tree expander, expanding incrementally the underlying graph with one vertex at a time, has full support on the space of junction trees with any given number of underlying vertices. On the other hand, the junction-tree collapser provides a complementary operation for removing vertices in the underlying decomposable graph of a junction tree, while maintaining the junction tree property. A direct application of our suggested algorithms is demonstrated in a sequential-Monte-Carlo setting designed for sampling from distributions on spaces of decomposable graphs. Numerical studies illustrate the utility of the proposed algorithms for combinatorial computations on decomposable graphs and junction trees. All the methods proposed in the paper are implemented in the Python library trilearn.
△ Less
Submitted 15 February, 2021; v1 submitted 2 June, 2018;
originally announced June 2018.
-
Bayesian learning of weakly structural Markov graph laws using sequential Monte Carlo methods
Authors:
Jimmy Olsson,
Tetyana Pavlenko,
Felix L. Rios
Abstract:
We present a sequential sampling methodology for weakly structural Markov laws, arising naturally in a Bayesian structure learning context for decomposable graphical models. As a key component of our suggested approach, we show that the problem of graph estimation, which in general lacks natural sequential interpretation, can be recast into a sequential setting by proposing a recursive Feynman-Kac…
▽ More
We present a sequential sampling methodology for weakly structural Markov laws, arising naturally in a Bayesian structure learning context for decomposable graphical models. As a key component of our suggested approach, we show that the problem of graph estimation, which in general lacks natural sequential interpretation, can be recast into a sequential setting by proposing a recursive Feynman-Kac model that generates a flow of junction tree distributions over a space of increasing dimensions. We focus on particle McMC methods to provide samples on this space, in particular on particle Gibbs (PG), as it allows for generating McMC chains with global moves on an underlying space of decomposable graphs. To further improve the PG mixing properties, we incorporate a systematic refreshment step implemented through direct sampling from a backward kernel. The theoretical properties of the algorithm are investigated, showing that the proposed refreshment step improves the performance in terms of asymptotic variance of the estimated distribution. The suggested sampling methodology is illustrated through a collection of numerical examples demonstrating high accuracy in Bayesian graph structure learning in both discrete and continuous graphical models.
△ Less
Submitted 23 July, 2019; v1 submitted 31 May, 2018;
originally announced May 2018.
-
Graphical posterior predictive classifier: Bayesian model averaging with particle Gibbs
Authors:
Tatjana Pavlenko,
Felix Leopoldo Rios
Abstract:
In this study, we present a multi-class graphical Bayesian predictive classifier that incorporates the uncertainty in the model selection into the standard Bayesian formalism. For each class, the dependence structure underlying the observed features is represented by a set of decomposable Gaussian graphical models. Emphasis is then placed on the Bayesian model averaging which takes full account of…
▽ More
In this study, we present a multi-class graphical Bayesian predictive classifier that incorporates the uncertainty in the model selection into the standard Bayesian formalism. For each class, the dependence structure underlying the observed features is represented by a set of decomposable Gaussian graphical models. Emphasis is then placed on the Bayesian model averaging which takes full account of the class-specific model uncertainty by averaging over the posterior graph model probabilities. An explicit evaluation of the model probabilities is well known to be infeasible. To address this issue, we consider the particle Gibbs strategy of Olsson et al. (2018b) for posterior sampling from decomposable graphical models which utilizes the Christmas tree algorithm of Olsson et al. (2018a) as proposal kernel. We also derive a strong hyper Markov law which we call the hyper normal Wishart law that allow to perform the resultant Bayesian calculations locally. The proposed predictive graphical classifier reveals superior performance compared to the ordinary Bayesian predictive rule that does not account for the model uncertainty, as well as to a number of out-of-the-box classifiers.
△ Less
Submitted 6 June, 2018; v1 submitted 21 July, 2017;
originally announced July 2017.
-
A Categorical Model for a Quantum Circuit Description Language (Extended Abstract)
Authors:
Francisco Rios,
Peter Selinger
Abstract:
Quipper is a practical programming language for describing families of quantum circuits. In this paper, we formalize a small, but useful fragment of Quipper called Proto-Quipper-M. Unlike its parent Quipper, this language is type-safe and has a formal denotational and operational semantics. Proto-Quipper-M is also more general than Quipper, in that it can describe families of morphisms in any symm…
▽ More
Quipper is a practical programming language for describing families of quantum circuits. In this paper, we formalize a small, but useful fragment of Quipper called Proto-Quipper-M. Unlike its parent Quipper, this language is type-safe and has a formal denotational and operational semantics. Proto-Quipper-M is also more general than Quipper, in that it can describe families of morphisms in any symmetric monoidal category, of which quantum circuits are but one example. We design Proto-Quipper-M from the ground up, by first giving a general categorical model of parameters and state. The distinction between parameters and state is also known from hardware description languages. A parameter is a value that is known at circuit generation time, whereas a state is a value that is known at circuit execution time. After finding some interesting categorical structures in the model, we then define the programming language to fit the model. We cement the connection between the language and the model by proving type safety, soundness, and adequacy properties.
△ Less
Submitted 1 March, 2018; v1 submitted 8 June, 2017;
originally announced June 2017.
-
A Prior Distribution over Directed Acyclic Graphs for Sparse Bayesian Networks
Authors:
Felix L. Rios,
John M. Noble,
Timo J. T. Koski
Abstract:
The main contribution of this article is a new prior distribution over directed acyclic graphs, which gives larger weight to sparse graphs. This distribution is intended for structured Bayesian networks, where the structure is given by an ordered block model. That is, the nodes of the graph are objects which fall into categories (or blocks); the blocks have a natural ordering. The presence of a re…
▽ More
The main contribution of this article is a new prior distribution over directed acyclic graphs, which gives larger weight to sparse graphs. This distribution is intended for structured Bayesian networks, where the structure is given by an ordered block model. That is, the nodes of the graph are objects which fall into categories (or blocks); the blocks have a natural ordering. The presence of a relationship between two objects is denoted by an arrow, from the object of lower category to the object of higher category. The models considered here were introduced in Kemp et al. (2004) for relational data and extended to multivariate data in Mansinghka et al. (2006). The prior over graph structures presented here has an explicit formula. The number of nodes in each layer of the graph follow a Hoppe Ewens urn model.
We consider the situation where the nodes of the graph represent random variables, whose joint probability distribution factorises along the DAG. We describe Monte Carlo schemes for finding the optimal aposteriori structure given a data matrix and compare the performance with Mansinghka et al. (2006) and also with the uniform prior.
△ Less
Submitted 25 April, 2015;
originally announced April 2015.
-
Quantum Communication with Continuum Single-Photon Pulses
Authors:
F. F. S. Rios,
R. V. Ramos
Abstract:
In this work, we analyze the behavior of continuum single-photon pulses in some quantum communication schemes. In particular, we consider the single-photon interference in a Mach-Zenhder interferometer, the HOM interference and the quantum bit commitment protocol.
In this work, we analyze the behavior of continuum single-photon pulses in some quantum communication schemes. In particular, we consider the single-photon interference in a Mach-Zenhder interferometer, the HOM interference and the quantum bit commitment protocol.
△ Less
Submitted 15 May, 2014;
originally announced May 2014.
-
Estado Actual de la Práctica de la Ingeniería de Software en México
Authors:
Reyes Juárez-Ramírez,
Karen Cortés Verdín,
Beatriz Angélica Toscano de la Torre,
Hanna Oktaba,
Carlos Alberto Fernández-y-Fernández,
Brenda Leticia Flores Ríos,
Fabiola Angulo Molina
Abstract:
The software engineering is a relatively new discipline compared to other sciences, since the origins of the term itself dates back to the years 1968 and 1969. At present, the market and the software industry have a significant relevance in several countries of the world; however, although Mexico is immersed in this race, has not even reached the level of success achieved in other countries in thi…
▽ More
The software engineering is a relatively new discipline compared to other sciences, since the origins of the term itself dates back to the years 1968 and 1969. At present, the market and the software industry have a significant relevance in several countries of the world; however, although Mexico is immersed in this race, has not even reached the level of success achieved in other countries in this sector. This paper presents an overview of the situation that keeps the practice of software engineering in Mexico, with emphasis on the academic realm. It shows a compilation of scientific research activity carried out in universities, as well as a brief analysis of undergraduate educational programs including the software engineering discipline . At the end, future work to be done is proposed in order to find a point of convergence between academia and industry, and also to support the flourishing of this business which somehow will have a positive impact on the economy of our country.
△ Less
Submitted 24 October, 2013;
originally announced October 2013.
-
Observability of Nonlinear Complex Networks in the Presence of Symmetries: A Graphical Approach
Authors:
Afroza Shirin,
Dionicio F. Rios,
Francesco Sorrentino
Abstract:
Reconstructing the states of the nodes of a dynamical network is a problem of fundamental importance in the study of neuronal and genetic networks. An underlying related problem is that of observability, i.e., identifying the conditions under which such a reconstruction is possible. In this paper we study observability of complex dynamical networks, where we consider the effects of network symmetr…
▽ More
Reconstructing the states of the nodes of a dynamical network is a problem of fundamental importance in the study of neuronal and genetic networks. An underlying related problem is that of observability, i.e., identifying the conditions under which such a reconstruction is possible. In this paper we study observability of complex dynamical networks, where we consider the effects of network symmetries on observability. We present an efficient algorithm that returns a minimal set of necessary sensor nodes for observability in the presence of symmetries.
△ Less
Submitted 30 March, 2017; v1 submitted 23 August, 2013;
originally announced August 2013.
