skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, January 16 until 2:00 AM ET on Friday, January 17 due to maintenance. We apologize for the inconvenience.


Title: A Symmetric Interior Penalty Method for an Elliptic Distributed Optimal Control Problem with Pointwise State Constraints
Abstract

We construct a symmetric interior penalty method for an elliptic distributed optimal control problem with pointwise state constraints on general polygonal domains.The resulting discrete problems are quadratic programs with simple box constraints that can be solved efficiently by a primal-dual active set algorithm.Both theoretical analysis and corroborating numerical results are presented.

 
more » « less
Award ID(s):
2208404
PAR ID:
10508358
Author(s) / Creator(s):
; ;
Publisher / Repository:
DeGruyter
Date Published:
Journal Name:
Computational Methods in Applied Mathematics
Volume:
23
Issue:
3
ISSN:
1609-4840
Page Range / eLocation ID:
565 to 589
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Background

    Genome-scale metabolic network models and constraint-based modeling techniques have become important tools for analyzing cellular metabolism. Thermodynamically infeasible cycles (TICs) causing unbounded metabolic flux ranges are often encountered. TICs satisfy the mass balance and directionality constraints but violate the second law of thermodynamics. Current practices involve implementing additional constraints to ensure not only optimal but also loopless flux distributions. However, the mixed integer linear programming problems required to solve become computationally intractable for genome-scale metabolic models.

    Results

    We aimed to identify the fewest needed constraints sufficient for optimality under the loopless requirement. We found that loopless constraints are required only for the reactions that share elementary flux modes representing TICs with reactions that are part of the objective function. We put forth the concept of localized loopless constraints (LLCs) to enforce this minimal required set of loopless constraints. By combining with a novel procedure for minimal null-space calculation, the computational time for loopless flux variability analysis (ll-FVA) is reduced by a factor of 10–150 compared to the original loopless constraints and by 4–20 times compared to the current fastest method Fast-SNP with the percent improvement increasing with model size. Importantly, LLCs offer a scalable strategy for loopless flux calculations for multi-compartment/multi-organism models of large sizes, for example, shortening the CPU time for ll-FVA from 35 h to less than 2 h for a model with more than104 reactions.

    Availability and implementation

    Matlab functions are available in the Supplementary Material or at https://github.com/maranasgroup/lll-FVA

    Supplementary information

    Supplementary data are available at Bioinformatics online.

     
    more » « less
  2. Abstract Motivation

    A chronogram is a dated phylogenetic tree whose branch lengths have been scaled to represent time. Such chronograms are computed based on available date estimates (e.g. from dated fossils), which provide absolute time constraints for one or more nodes of an input undated phylogeny, coupled with an appropriate underlying model for evolutionary rates variation along the branches of the phylogeny. However, traditional methods for phylogenetic dating cannot take into account relative time constraints, such as those provided by inferred horizontal transfer events. In many cases, chronograms computed using only absolute time constraints are inconsistent with known relative time constraints.

    Results

    In this work, we introduce a new approach, Dating Trees using Relative constraints (DaTeR), for phylogenetic dating that can take into account both absolute and relative time constraints. The key idea is to use existing Bayesian approaches for phylogenetic dating to sample posterior chronograms satisfying desired absolute time constraints, minimally adjust or ‘error-correct’ these sampled chronograms to satisfy all given relative time constraints, and aggregate across all error-corrected chronograms. DaTeR uses a constrained optimization framework for the error-correction step, finding minimal deviations from previously assigned dates or branch lengths. We applied DaTeR to a biological dataset of 170 Cyanobacterial taxa and a reliable set of 24 transfer-based relative constraints, under six different molecular dating models. Our extensive analysis of this dataset demonstrates that DaTeR is both highly effective and scalable and that its application can significantly improve estimated chronograms.

    Availability and implementation

    Freely available from https://compbio.engr.uconn.edu/software/dater/

    Supplementary information

    Supplementary data are available at Bioinformatics online.

     
    more » « less
  3. Abstract Research Highlights

    This research addresses the apparent contradiction that children are sensitive to constraints in experimental paradigms but are ofteninsensitiveto constraints in the real world.

    One explanation for this discrepancy is that constraints in the real world tend to be stable over time and may lose their salience.

    When probabilistic constraints (i.e., when a second option is available but hard to access) are stable, children become de‐sensitized to constraints across trials.

    First observing contexts with greater choice increases children's sensitivity to both probabilistic and deterministic constraints.

     
    more » « less
  4. Abstract

    This paper examines constraints and their role in scientific explanation. Common views in the philosophical literature suggest that constraints are non-causal and that they provide non-causal explanations. While much of this work focuses on examples from physics, this paper explores constraints from other fields, including neuroscience, physiology, and the social sciences. I argue that these cases involve constraints that are causal and that provide a unique type of causal explanation. This paper clarifies what it means for a factor to be a constraint, when such constraints are causal, and how they figure in scientific explanation.

     
    more » « less
  5. Abstract Motivation

    Simulation is an essential technique for generating biomolecular data with a ‘known’ history for use in validating phylogenetic inference and other evolutionary methods. On longer time scales, simulation supports investigations of equilibrium behavior and provides a formal framework for testing competing evolutionary hypotheses. Twenty years of molecular evolution research have produced a rich repertoire of simulation methods. However, current models do not capture the stringent constraints acting on the domain insertions, duplications, and deletions by which multidomain architectures evolve. Although these processes have the potential to generate any combination of domains, only a tiny fraction of possible domain combinations are observed in nature. Modeling these stringent constraints on domain order and co-occurrence is a fundamental challenge in domain architecture simulation that does not arise with sequence and gene family simulation.

    Results

    Here, we introduce a stochastic model of domain architecture evolution to simulate evolutionary trajectories that reflect the constraints on domain order and co-occurrence observed in nature. This framework is implemented in a novel domain architecture simulator, DomArchov, using the Metropolis–Hastings algorithm with data-driven transition probabilities. The use of a data-driven event module enables quick and easy redeployment of the simulator for use in different taxonomic and protein function contexts. Using empirical evaluation with metazoan datasets, we demonstrate that domain architectures simulated by DomArchov recapitulate properties of genuine domain architectures that reflect the constraints on domain order and adjacency seen in nature. This work expands the realm of evolutionary processes that are amenable to simulation.

    Availability and implementation

    DomArchov is written in Python 3 and is available at http://www.cs.cmu.edu/~durand/DomArchov. The data underlying this article are available via the same link.

    Supplementary information

    Supplementary data are available at Bioinformatics online.

     
    more » « less