Related papers: Numerical implementation of some reweighted path i…
The spectral density of various ensembles of sparse symmetric random matrices is analyzed using the cavity method. We consider two cases: matrices whose associated graphs are locally tree-like, and sparse covariance matrices. We derive a…
The treatment of two-dimensional random walks in the quarter plane leads to Markov processes which involve semi-infinite matrices having Toeplitz or block Toeplitz structure plus a low-rank correction. Finding the steady state probability…
We introduce an extension to the Weighted Ensemble (WE) path sampling method to restrict sampling to a one dimensional path through a high dimensional phase space. Our method, which is based on the finite-temperature string method, permits…
Full-waveform inversion (FWI) is a powerful geophysical imaging technique that infers high-resolution subsurface physical parameters by solving a non-convex optimization problem. However, due to limitations in observation, e.g., limited…
Feynman's path integral is herein generalized to the nonextensive canonical density matrix based on Tsallis entropy. This generalization is done in two ways by using unnormalized and normalized constraints. Firstly, we consider the path…
We construct in the small-time setting the upper and lower estimates for the transition probability density of a L\'evy process in $\rn$. Our approach relies on the complex analysis technique and the asymptotic analysis of the inverse…
In prior work \cite{AD} of Lars Andersson and Bruce K. Driver, the path space with finite interval over a compact Riemannian manifold is approximated by finite dimensional manifolds $H_{x,\P} (M)$ consisting of piecewise geodesic paths…
We develop a computationally efficient and robust algorithm for generating pseudo-random samples from a broad class of smooth probability distributions in one and two dimensions. The algorithm is based on inverse transform sampling with a…
We consider a method for the approximation of iterated stochastic integrals of arbitrary multiplicity $k$ $(k\in \mathbb{N})$ with respect to the infinite-dimensional $Q$-Wiener process using the mean-square approximation method of iterated…
The Chiral Random Matrix Model or the Gaussian Penner Model (generalized Laguerre ensemble) is re-examined in the light of the results which have been found in double well matrix models [D97,BD99] and subtleties discovered in the single…
We study the finite convergence of iterative methods for solving convex feasibility problems. Our key assumptions are that the interior of the solution set is nonempty and that certain overrelaxation parameters converge to zero, but with a…
The restricted Feynman path integrals (RFPIs) have been proposed to study continuous quantum measurements in physics. The RFPIs are heuristically determined in terms of the usual probability amplitude multiplied by weight for each path,…
Bayesian methods are a popular choice for statistical inference in small-data regimes due to the regularization effect induced by the prior. In the context of density estimation, the standard nonparametric Bayesian approach is to target the…
We present an extension of the Levenberg-Marquardt algorithm for fitting multichannel nuclear cross section data. Our approach offers a practical and robust alternative to conventional trust-region methods for analyzing experimental data.…
When sampling multi-modal probability distributions, correctly estimating the relative probability of each mode, even when the modes have been discovered and locally sampled, remains challenging. We test a simple reweighting scheme designed…
We provide dual algorithms for sampling the space of abstract simplicial complexes on a fixed number of vertices. We develop a generative and descriptive sampler designed with heuristics to help balance the combinatorial multiplicities of…
In the framework of solid mechanics, the task of deriving material parameters from experimental data has recently re-emerged with the progress in full-field measurement capabilities and the renewed advances of machine learning. In this…
Path integral techniques for the density matrix of a one-dimensional statistical system near a boundary previously employed in black-hole physics are applied to providing a new interpretation of the density matrix renormalization group: its…
We develop randomized matrix-free algorithms for estimating partial traces, a generalization of the trace arising in quantum physics and chemistry. Our algorithm improves on the typicality-based approach used in [T. Chen and Y-C. Cheng,…
The presence of erratic or unstable paths in standard kinetic Monte Carlo simulations significantly undermines the accurate simulation and sampling of transition pathways. While typically reliable methods, such as the Gillespie algorithm,…