-
Dual Dictionaries in Linear Programming
Authors:
Patrick T. Perkins,
Xiang Gao
Abstract:
In order to use the Dual Simplex Method, one needs to prove a certain bijection between the dictionaries associated with the primal problem and those associated with its dual. We give a short conceptual proof of why this bijection exists.
In order to use the Dual Simplex Method, one needs to prove a certain bijection between the dictionaries associated with the primal problem and those associated with its dual. We give a short conceptual proof of why this bijection exists.
△ Less
Submitted 19 September, 2023;
originally announced October 2023.
-
Experimental and computational comparison of freeze-thaw induced pressure generation in red and sugar maple
Authors:
Maryam Zarrinderakht,
Isabell Konrad,
Timothy R. Wilmot,
Timothy D. Perkins,
Abby K. van den Berg,
John M. Stockie
Abstract:
Sap exudation is the process whereby trees such as sugar (Acer saccharum) and red maple (Acer rubrum) generate unusually high positive stem pressure in response to repeated cycles of freeze and thaw. This elevated xylem pressure permits the sap to be harvested over a period of several weeks and hence is a major factor in the viability of the maple syrup industry. The extensive literature on sap ex…
▽ More
Sap exudation is the process whereby trees such as sugar (Acer saccharum) and red maple (Acer rubrum) generate unusually high positive stem pressure in response to repeated cycles of freeze and thaw. This elevated xylem pressure permits the sap to be harvested over a period of several weeks and hence is a major factor in the viability of the maple syrup industry. The extensive literature on sap exudation documents competing hypotheses regarding the physical and biological mechanisms that drive positive pressure generation in maple, but to date relatively little effort has been expended on devising mathematical models for the exudation process. In this paper, we utilize an existing model of Graf et al. [J. Roy. Soc. Interface 12:20150665, 2015] that describes heat and mass transport within the multiphase gas-liquid-ice mixture in the porous xylem tissue. The model captures the inherent multiscale nature of xylem transport by including phase change and osmotic transport in wood cells on the microscale, which is coupled to heat transport through the tree stem on the macroscale. A parametric study based on simulations with synthetic temperature data identifies the model parameters that have greatest impact on stem pressure build-up. Measured daily temperature fluctuations are then used as model inputs and the resulting simulated pressures are compared directly with experimental measurements taken from mature red and sugar maple stems during the sap harvest season. The results demonstrate that our multiscale freeze-thaw model reproduces realistic exudation behavior, thereby providing novel insights into the specific physical mechanisms that dominate positive pressure generation in maple trees.
△ Less
Submitted 3 January, 2024; v1 submitted 5 June, 2021;
originally announced June 2021.
-
Reinforcement Learning via Parametric Cost Function Approximation for Multistage Stochastic Programming
Authors:
Saeed Ghadimi,
Raymond T. Perkins,
Warren B. Powell
Abstract:
The most common approaches for solving stochastic resource allocation problems in the research literature is to either use value functions ("dynamic programming") or scenario trees ("stochastic programming") to approximate the impact of a decision now on the future. By contrast, common industry practice is to use a deterministic approximation of the future which is easier to understand and solve,…
▽ More
The most common approaches for solving stochastic resource allocation problems in the research literature is to either use value functions ("dynamic programming") or scenario trees ("stochastic programming") to approximate the impact of a decision now on the future. By contrast, common industry practice is to use a deterministic approximation of the future which is easier to understand and solve, but which is criticized for ignoring uncertainty. We show that a parameterized version of a deterministic lookahead can be an effective way of handling uncertainty, while enjoying the computational simplicity of a deterministic lookahead. We present the parameterized lookahead model as a form of policy for solving a stochastic base model, which is used as the basis for optimizing the parameterized policy. This approach can handle complex, high-dimensional state variables, and avoids the usual approximations associated with scenario trees. We formalize this approach and demonstrate its use in the context of a complex, nonstationary energy storage problem.
△ Less
Submitted 2 January, 2020;
originally announced January 2020.
-
Stochastic Optimization with Parametric Cost Function Approximations
Authors:
Raymond T. Perkins III,
Warren B. Powell
Abstract:
A widely used heuristic for solving stochastic optimization problems is to use a deterministic rolling horizon procedure, which has been modified to handle uncertainty (e.g. buffer stocks, schedule slack). This approach has been criticized for its use of a deterministic approximation of a stochastic problem, which is the major motivation for stochastic programming. We recast this debate by identif…
▽ More
A widely used heuristic for solving stochastic optimization problems is to use a deterministic rolling horizon procedure, which has been modified to handle uncertainty (e.g. buffer stocks, schedule slack). This approach has been criticized for its use of a deterministic approximation of a stochastic problem, which is the major motivation for stochastic programming. We recast this debate by identifying both deterministic and stochastic approaches as policies for solving a stochastic base model, which may be a simulator or the real world. Stochastic lookahead models (stochastic programming) require a range of approximations to keep the problem tractable. By contrast, so-called deterministic models are actually parametrically modified cost function approximations which use parametric adjustments to the objective function and/or the constraints. These parameters are then optimized in a stochastic base model which does not require making any of the types of simplifications required by stochastic programming. We formalize this strategy and describe a gradient-based stochastic search strategy to optimize the parameters.
△ Less
Submitted 14 March, 2017;
originally announced March 2017.
-
Dirichlet problems on graphs with ends
Authors:
Tony Perkins
Abstract:
In classical potential theory, one can solve the Dirichlet problem on unbounded domains such as the upper half plane. These domains have two types of boundary points; the usual finite boundary points and another point at infinity. W. Woess has solved a discrete version of the Dirichlet problem on the ends of graphs analogous to having multiple points at infinity and no finite boundary. Whereas C.…
▽ More
In classical potential theory, one can solve the Dirichlet problem on unbounded domains such as the upper half plane. These domains have two types of boundary points; the usual finite boundary points and another point at infinity. W. Woess has solved a discrete version of the Dirichlet problem on the ends of graphs analogous to having multiple points at infinity and no finite boundary. Whereas C. Kiselman has solved a similar version of the Dirichlet problem on graphs analogous to bounded domains. In this work, we combine the two ideas to solve a version of the Dirichlet problem on graphs with finitely many ends and boundary points of the Kiselman type.
△ Less
Submitted 24 June, 2014;
originally announced June 2014.
-
Scaling Properties of Paths on Graphs
Authors:
R. Edwards,
E. Foxall,
T. J. Perkins
Abstract:
Let $G$ be a directed graph on finitely many vertices and edges, and assign a positive weight to each edge on $G$. Fix vertices $u$ and $v$ and consider the set of paths that start at $u$ and end at $v$, self-intersecting in any number of places along the way. For each path, sum the weights of its edges, and then list the path weights in increasing order. The asymptotic behaviour of this sequence…
▽ More
Let $G$ be a directed graph on finitely many vertices and edges, and assign a positive weight to each edge on $G$. Fix vertices $u$ and $v$ and consider the set of paths that start at $u$ and end at $v$, self-intersecting in any number of places along the way. For each path, sum the weights of its edges, and then list the path weights in increasing order. The asymptotic behaviour of this sequence is described, in terms of the structure and type of strongly connected components on the graph. As a special case, for a Markov chain the asymptotic probability of paths obeys either a power law scaling or a weaker type of scaling, depending on the structure of the transition matrix. This generalizes previous work by Mandelbrot and others, who established asymptotic power law scaling for special classes of Markov chains.
△ Less
Submitted 12 June, 2013;
originally announced June 2013.
-
Convex and subharmonic functions on graphs
Authors:
Matthew Burke,
Tony Perkins
Abstract:
We explore the relationship between convex and subharmonic functions on discrete sets. Our principal concern is to determine the setting in which a convex function is necessarily subharmonic. We initially consider the primary notions of convexity on graphs and show that more structure is needed to establish the desired result. To that end, we consider a notion of convexity defined on lattice-like…
▽ More
We explore the relationship between convex and subharmonic functions on discrete sets. Our principal concern is to determine the setting in which a convex function is necessarily subharmonic. We initially consider the primary notions of convexity on graphs and show that more structure is needed to establish the desired result. To that end, we consider a notion of convexity defined on lattice-like graphs generated by normed abelian groups. For this class of graphs, we are able to prove that all convex functions are subharmonic.
△ Less
Submitted 1 April, 2013;
originally announced April 2013.
-
Green functions, the fine topology and restoring coverings
Authors:
Tony L. Perkins
Abstract:
There are several equivalent ways to define continuous harmonic functions $H(K)$ on a compact set $K$ in $\mathbb R^n$. One may let $H(K)$ be the unform closures of all functions in $C(K)$ which are restrictions of harmonic functions on a neighborhood of $K$, or take $H(K)$ as the subspace of $C(K)$ consisting of functions which are finely harmonic on the fine interior of $K$. In \cite{DG74} it wa…
▽ More
There are several equivalent ways to define continuous harmonic functions $H(K)$ on a compact set $K$ in $\mathbb R^n$. One may let $H(K)$ be the unform closures of all functions in $C(K)$ which are restrictions of harmonic functions on a neighborhood of $K$, or take $H(K)$ as the subspace of $C(K)$ consisting of functions which are finely harmonic on the fine interior of $K$. In \cite{DG74} it was shown that these definitions are equivalent. Using a localization result of \cite{BH78} one sees that a function $h\in H(K)$ if and only if it is continuous and finely harmonic on on every fine connected component of the fine interior of $K$. Such collection of sets are usually called {\it restoring}.
Another equivalent definition of $H(K)$ was introduced in \cite{P97} using the notion of Jensen measures which leads another restoring collection of sets. The main goal of this paper is to reconcile the results in \cite{DG74} and \cite{P97}.
To study these spaces, two notions of Green functions have previously been introduced. One by \cite{P97} as the limit of Green functions on domains $D_j$ where the domains $D_j$ are decreasing to $K$, and alternatively following \cite{F72, F75} one has the fine Green function on the fine interior of $K$. Our Theorem \ref{T:green_equiv} shows that these are equivalent notions.
In Section \ref{S:Jensen} a careful study of the set of Jensen measures on $K$, leads to an interesting extension result (Corollary \ref{C:extend}) for superharmonic functions. This has a number of applications. In particular we show that the two restoring coverings are the same. We are also able to extend some results of \cite{GL83} and \cite{P97} to higher dimensions.
△ Less
Submitted 17 December, 2010;
originally announced December 2010.
-
Non-compact versions of Edwards' Theorem
Authors:
Nihat G. Gogus,
Tony L. Perkins,
Evgeny A. Poletsky
Abstract:
Edwards' Theorem establishes duality between a convex cone in the space of continuous functions on a compact space and the set of representing or Jensen measures for this cone. In this paper we prove non-compact versions of this theorem.
Edwards' Theorem establishes duality between a convex cone in the space of continuous functions on a compact space and the set of representing or Jensen measures for this cone. In this paper we prove non-compact versions of this theorem.
△ Less
Submitted 23 December, 2010; v1 submitted 17 December, 2010;
originally announced December 2010.
-
The Dirichlet Problem for Harmonic Functions on Compact Sets
Authors:
Tony Perkins
Abstract:
For any compact set $K\subset \mathbb{R}^n$ we develop the theory of Jensen measures and subharmonic peak points, which form the set $\mathcal{O}_K$, to study the Dirichlet problem on $K$. Initially we consider the space $h(K)$ of functions on $K$ which can be uniformly approximated by functions harmonic in a neighborhood of $K$ as possible solutions. As in the classical theory, our Theorem 8.1 sh…
▽ More
For any compact set $K\subset \mathbb{R}^n$ we develop the theory of Jensen measures and subharmonic peak points, which form the set $\mathcal{O}_K$, to study the Dirichlet problem on $K$. Initially we consider the space $h(K)$ of functions on $K$ which can be uniformly approximated by functions harmonic in a neighborhood of $K$ as possible solutions. As in the classical theory, our Theorem 8.1 shows $C(\mathcal{O}_K)\cong h(K)$ for compact sets with $\mathcal{O}_K$ closed. However, in general a continuous solution cannot be expected even for continuous data on $\rO_K$ as illustrated by Theorem 8.1. Consequently, we show that the solution can be found in a class of finely harmonic functions. Moreover by Theorem 8.7, in complete analogy with the classical situation, this class is isometrically isomorphic to $C_b(\mathcal{O}_K)$ for all compact sets $K$.
△ Less
Submitted 21 January, 2011; v1 submitted 30 April, 2010;
originally announced April 2010.