首页 > 最新文献

Monte Carlo Methods and Applications最新文献

英文 中文
Frontmatter 头版头条
Q4 Mathematics Pub Date : 2023-03-01 DOI: 10.1515/mcma-2023-frontmatter1
{"title":"Frontmatter","authors":"","doi":"10.1515/mcma-2023-frontmatter1","DOIUrl":"https://doi.org/10.1515/mcma-2023-frontmatter1","url":null,"abstract":"","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136270828","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Linking the Monte Carlo radiative transfer algorithm to the radiative transfer equation 将蒙特卡罗辐射传递算法与辐射传递方程联系起来
IF 0.9 Q4 Mathematics Pub Date : 2023-01-27 DOI: 10.1515/mcma-2023-2001
P. J. Valades-Pelayo, M. Ramirez-Cabrera, A. Balbuena-Ortega
Abstract This manuscript presents a short route to justify the widely used Monte Carlo Radiative Transfer (MCRT) algorithm straight from the Radiative Transfer Equation (RTE). In this regard, this paper starts deriving a probability measure obtained from the integral formulation of the RTE under a unidirectional point source in an infinite domain. This derivation only requires the analytical integration of the first two terms of a perturbation expansion. Although derivations have been devised to clarify the relationship between the MCRT and the RTE, they tend to be rather long and elaborate. Considering how simple it is to justify the MCRT from a loose probabilistic interpretation of the photon’s physical propagation process, the decay in popularity of former approaches relating MCRT to the RTE is entirely understandable. Unfortunately, all of this has given the false impression that MCRT and the RTE are not that closely related, to the point that recent works have explicitly stated that no direct link exists between them. This work presents a simpler route demonstrating how the MCRT algorithm emerges to statistically sample the RTE explicitly through Markov chains, further clarifying the method’s foundations. Although compact, the derivation proposed in this work does not skip any fundamental step, preserving mathematical rigor while giving specific expressions and functions. Thus, this derivation can help devise efficient ways to statistically sample the RTE for different scenarios or when coupling the MCRT method with other methods traditionally grounded in the RTE, such as the Spherical Harmonics and Discrete Ordinates methods.
摘要本文直接从辐射传递方程(RTE)出发,提出了一条简单的途径来证明广泛使用的蒙特卡罗辐射传递(MCRT)算法。在这方面,本文开始推导一个概率测度,该测度是从无限域中单向点源下RTE的积分公式中获得的。这种推导只需要对扰动展开的前两项进行分析积分。尽管推导是为了澄清MCRT和RTE之间的关系而设计的,但它们往往相当冗长和详细。考虑到从光子物理传播过程的松散概率解释来证明MCRT是多么简单,以前将MCRT与RTE相关的方法的流行程度下降是完全可以理解的。不幸的是,所有这些都给人一种错误的印象,即MCRT和RTE没有那么紧密的联系,以至于最近的作品明确表示它们之间不存在直接联系。这项工作提供了一个更简单的途径,展示了MCRT算法是如何通过马尔可夫链显式地对RTE进行统计采样的,进一步阐明了该方法的基础。尽管紧凑,但这项工作中提出的推导并没有跳过任何基本步骤,在给出特定表达式和函数的同时保持了数学的严谨性。因此,这种推导可以帮助设计有效的方法来对不同场景的RTE进行统计采样,或者当将MCRT方法与传统上基于RTE的其他方法相结合时,例如球面谐波和离散坐标方法。
{"title":"Linking the Monte Carlo radiative transfer algorithm to the radiative transfer equation","authors":"P. J. Valades-Pelayo, M. Ramirez-Cabrera, A. Balbuena-Ortega","doi":"10.1515/mcma-2023-2001","DOIUrl":"https://doi.org/10.1515/mcma-2023-2001","url":null,"abstract":"Abstract This manuscript presents a short route to justify the widely used Monte Carlo Radiative Transfer (MCRT) algorithm straight from the Radiative Transfer Equation (RTE). In this regard, this paper starts deriving a probability measure obtained from the integral formulation of the RTE under a unidirectional point source in an infinite domain. This derivation only requires the analytical integration of the first two terms of a perturbation expansion. Although derivations have been devised to clarify the relationship between the MCRT and the RTE, they tend to be rather long and elaborate. Considering how simple it is to justify the MCRT from a loose probabilistic interpretation of the photon’s physical propagation process, the decay in popularity of former approaches relating MCRT to the RTE is entirely understandable. Unfortunately, all of this has given the false impression that MCRT and the RTE are not that closely related, to the point that recent works have explicitly stated that no direct link exists between them. This work presents a simpler route demonstrating how the MCRT algorithm emerges to statistically sample the RTE explicitly through Markov chains, further clarifying the method’s foundations. Although compact, the derivation proposed in this work does not skip any fundamental step, preserving mathematical rigor while giving specific expressions and functions. Thus, this derivation can help devise efficient ways to statistically sample the RTE for different scenarios or when coupling the MCRT method with other methods traditionally grounded in the RTE, such as the Spherical Harmonics and Discrete Ordinates methods.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2023-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43267975","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Methodology for nonparametric bias reduction in kernel regression estimation 核回归估计中的非参数偏差约简方法
IF 0.9 Q4 Mathematics Pub Date : 2023-01-10 DOI: 10.1515/mcma-2022-2130
Y. Slaoui
Abstract In this paper, we propose and investigate two new kernel regression estimators based on a bias reduction transformation technique. We study the properties of these estimators and compare them with Nadaraya–Watson’s regression estimator and Slaoui’s (2016) regression estimator. It turns out that, with an adequate choice of the parameters of the two proposed estimators, the rate of convergence of two estimators will be faster than the two classical estimators, and the asymptotic MISE (mean integrated squared error) will be smaller than the two classical estimators. We corroborate these theoretical results through simulations and a real Malaria dataset.
摘要在本文中,我们提出并研究了两种新的基于偏差减少变换技术的核回归估计量。我们研究了这些估计量的性质,并将其与Nadaraya–Watson回归估计量和Slaoui(2016)回归估计量进行了比较。结果表明,在适当选择两个估计量的参数的情况下,两个估计的收敛速度将快于两个经典估计量,并且渐近均方误差将小于两个经典估算量。我们通过模拟和真实的疟疾数据集证实了这些理论结果。
{"title":"Methodology for nonparametric bias reduction in kernel regression estimation","authors":"Y. Slaoui","doi":"10.1515/mcma-2022-2130","DOIUrl":"https://doi.org/10.1515/mcma-2022-2130","url":null,"abstract":"Abstract In this paper, we propose and investigate two new kernel regression estimators based on a bias reduction transformation technique. We study the properties of these estimators and compare them with Nadaraya–Watson’s regression estimator and Slaoui’s (2016) regression estimator. It turns out that, with an adequate choice of the parameters of the two proposed estimators, the rate of convergence of two estimators will be faster than the two classical estimators, and the asymptotic MISE (mean integrated squared error) will be smaller than the two classical estimators. We corroborate these theoretical results through simulations and a real Malaria dataset.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2023-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49604358","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Development and implementation of branching random walk on spheres algorithms for solving the 2D elastostatics Lamé equation 求解二维弹性静力学lam<s:1>方程的球上分支随机行走算法的开发与实现
IF 0.9 Q4 Mathematics Pub Date : 2023-01-10 DOI: 10.1515/mcma-2022-2131
I. Shalimova, K. Sabelfeld
Abstract In this paper, we address a long-standing open problem in stochastic simulation: construction of a random walk on spheres (RWS) algorithm for solving a system of elasticity equations, known as the Lamé equation. Many attempts to generalize the classic probabilistic representations like the Kac formula for parabolic and scalar elliptic equations failed. A different approach based on a branching random walk on spheres (BRWS) introduced in our paper of 1995 [K. K. Sabelfeld and D. Talay, Integral formulation of the boundary value problems and the method of random walk on spheres, Monte Carlo Methods Appl. 1 1995, 1, 1–34] made little progress in solving this problem. In the present study, we further improve the BRWS algorithm by a special implementation of a branching anisotropic random walk on spheres process.
摘要在本文中,我们解决了随机模拟中一个长期存在的开放问题:构造求解弹性方程组(即Lamé方程)的随机球上行走(RWS)算法。许多推广经典概率表示的尝试,如抛物型和标量椭圆方程的Kac公式,都失败了。我们在1995年的论文[K.K.Sabelfeld和D.Talay,边值问题的积分公式和球上随机行走方法,蒙特卡罗方法应用1 1995,1,1–34]中介绍了一种基于分支随机行走(BRWS)的不同方法,在解决这个问题方面进展甚微。在本研究中,我们通过一个特殊的分支各向异性球体随机行走过程来进一步改进BRWS算法。
{"title":"Development and implementation of branching random walk on spheres algorithms for solving the 2D elastostatics Lamé equation","authors":"I. Shalimova, K. Sabelfeld","doi":"10.1515/mcma-2022-2131","DOIUrl":"https://doi.org/10.1515/mcma-2022-2131","url":null,"abstract":"Abstract In this paper, we address a long-standing open problem in stochastic simulation: construction of a random walk on spheres (RWS) algorithm for solving a system of elasticity equations, known as the Lamé equation. Many attempts to generalize the classic probabilistic representations like the Kac formula for parabolic and scalar elliptic equations failed. A different approach based on a branching random walk on spheres (BRWS) introduced in our paper of 1995 [K. K. Sabelfeld and D. Talay, Integral formulation of the boundary value problems and the method of random walk on spheres, Monte Carlo Methods Appl. 1 1995, 1, 1–34] made little progress in solving this problem. In the present study, we further improve the BRWS algorithm by a special implementation of a branching anisotropic random walk on spheres process.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2023-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46049317","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Parameter least-squares estimation for time-inhomogeneous Ornstein–Uhlenbeck process 时间非均匀Ornstein-Uhlenbeck过程的参数最小二乘估计
IF 0.9 Q4 Mathematics Pub Date : 2022-11-02 DOI: 10.1515/mcma-2022-2127
G. Pramesti
Abstract We address the least-squares estimation of the drift coefficient parameter θ = ( λ , A , B , ω p ) theta=(lambda,A,B,omega_{p}) of a time-inhomogeneous Ornstein–Uhlenbeck process that is observed at high frequency, in which the discretized step size ℎ satisfies h → 0 hto 0 . In this paper, under the conditions n ⁢ h → ∞ nhtoinfty and n ⁢ h 2 → 0 nh^{2}to 0 , we prove the consistency and the asymptotic normality of the estimators. We obtain the convergence of the parameters at rate n ⁢ h sqrt{nh} , except for ω p omega_{p} at n 3 ⁢ h 3 sqrt{n^{3}h^{3}} . To verify our theoretical findings, we do a simulation study. We then illustrate the use of the proposed model in fitting the energy use of light fixtures in one Belgium household and the stock exchange.
摘要我们讨论了在高频下观测到的时间非均匀Ornstein–Uhlenbeck过程的漂移系数参数θ=(λ,A,B,ωp)θ=(lambda,A,B,omega_{p})的最小二乘估计,其中离散化的步长ℎ 满足h→ 0小时到0。在本文中,在条件n h→ ∞ nhtoinfty和n h 2→ 0nh^{2}到0,我们证明了估计量的一致性和渐近正态性。我们得到了在速率n h sqrt{nh}下参数的收敛性,除了ω^{3}h^{3} }。为了验证我们的理论发现,我们进行了一项模拟研究。然后,我们说明了所提出的模型在拟合比利时一个家庭和证券交易所的灯具能源使用方面的用途。
{"title":"Parameter least-squares estimation for time-inhomogeneous Ornstein–Uhlenbeck process","authors":"G. Pramesti","doi":"10.1515/mcma-2022-2127","DOIUrl":"https://doi.org/10.1515/mcma-2022-2127","url":null,"abstract":"Abstract We address the least-squares estimation of the drift coefficient parameter θ = ( λ , A , B , ω p ) theta=(lambda,A,B,omega_{p}) of a time-inhomogeneous Ornstein–Uhlenbeck process that is observed at high frequency, in which the discretized step size ℎ satisfies h → 0 hto 0 . In this paper, under the conditions n ⁢ h → ∞ nhtoinfty and n ⁢ h 2 → 0 nh^{2}to 0 , we prove the consistency and the asymptotic normality of the estimators. We obtain the convergence of the parameters at rate n ⁢ h sqrt{nh} , except for ω p omega_{p} at n 3 ⁢ h 3 sqrt{n^{3}h^{3}} . To verify our theoretical findings, we do a simulation study. We then illustrate the use of the proposed model in fitting the energy use of light fixtures in one Belgium household and the stock exchange.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2022-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44053654","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Global random walk on grid algorithm for solving Navier–Stokes and Burgers equations 求解Navier-Stokes和Burgers方程的全局随机网格行走算法
IF 0.9 Q4 Mathematics Pub Date : 2022-10-28 DOI: 10.1515/mcma-2022-2126
K. Sabelfeld, Oleg Bukhasheev
Abstract The global random walk on grid method (GRWG) is developed for solving two-dimensional nonlinear systems of equations, the Navier–Stokes and Burgers equations. This study extends the GRWG which we have earlier developed for solving the nonlinear drift-diffusion-Poisson equation of semiconductors (Physica A 556 (2020), Article ID 124800). The Burgers equation is solved by a direct iteration of a system of linear drift-diffusion equations, while the Navier–Stokes equation is solved in the stream function-vorticity formulation.
摘要提出了求解二维非线性方程组Navier–Stokes和Burgers方程的全局随机网格行走方法。本研究扩展了我们早期开发的用于求解半导体非线性漂移扩散泊松方程的GRWG(Physica A 556(2020),文章ID 124800)。Burgers方程是通过线性漂移扩散方程组的直接迭代求解的,而Navier–Stokes方程是通过流函数涡度公式求解的。
{"title":"Global random walk on grid algorithm for solving Navier–Stokes and Burgers equations","authors":"K. Sabelfeld, Oleg Bukhasheev","doi":"10.1515/mcma-2022-2126","DOIUrl":"https://doi.org/10.1515/mcma-2022-2126","url":null,"abstract":"Abstract The global random walk on grid method (GRWG) is developed for solving two-dimensional nonlinear systems of equations, the Navier–Stokes and Burgers equations. This study extends the GRWG which we have earlier developed for solving the nonlinear drift-diffusion-Poisson equation of semiconductors (Physica A 556 (2020), Article ID 124800). The Burgers equation is solved by a direct iteration of a system of linear drift-diffusion equations, while the Navier–Stokes equation is solved in the stream function-vorticity formulation.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2022-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43562871","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Simulation of drift-diffusion process at high Péclet numbers by the random walk on spheres method 用球上随机游走法模拟高psamclet数下的漂移扩散过程
IF 0.9 Q4 Mathematics Pub Date : 2022-10-28 DOI: 10.1515/mcma-2022-2128
K. Sabelfeld, Ivan Aksyuk
Abstract In this paper, we address the problem of flow simulation at high Péclet numbers by the random walk on spheres (RWS) method. Conventional deterministic methods here face difficulties related to high solution gradients near the boundary in the region known as the boundary layer. In the finite-difference methods, this leads to introduction of very fine meshes which in turn causes problems of stability and high dimensions. The RWS algorithm is mesh free, but the high Péclet number flows should probably also affect the efficiency of simulations. However, it turns out that the RWS algorithm can be well adapted to this case. We present an analysis of the RWS algorithm for different examples of flows with high Péclet number. Simulations are carried out for different boundary conditions and for two-layered material with different diffusion coefficients of exciton’s mobility.
摘要本文用球上随机游走(RWS)方法解决了高psamclet数下的流动模拟问题。传统的确定性方法在这里面临着与边界层附近区域的高解梯度有关的困难。在有限差分方法中,这会导致引入非常精细的网格,从而导致稳定性和高维问题。RWS算法是无网格的,但过高的psamclet数流可能也会影响模拟的效率。然而,事实证明RWS算法可以很好地适应这种情况。本文对具有高psamclet数的不同流例的RWS算法进行了分析。在不同边界条件下,对具有不同激子迁移率扩散系数的双层材料进行了模拟。
{"title":"Simulation of drift-diffusion process at high Péclet numbers by the random walk on spheres method","authors":"K. Sabelfeld, Ivan Aksyuk","doi":"10.1515/mcma-2022-2128","DOIUrl":"https://doi.org/10.1515/mcma-2022-2128","url":null,"abstract":"Abstract In this paper, we address the problem of flow simulation at high Péclet numbers by the random walk on spheres (RWS) method. Conventional deterministic methods here face difficulties related to high solution gradients near the boundary in the region known as the boundary layer. In the finite-difference methods, this leads to introduction of very fine meshes which in turn causes problems of stability and high dimensions. The RWS algorithm is mesh free, but the high Péclet number flows should probably also affect the efficiency of simulations. However, it turns out that the RWS algorithm can be well adapted to this case. We present an analysis of the RWS algorithm for different examples of flows with high Péclet number. Simulations are carried out for different boundary conditions and for two-layered material with different diffusion coefficients of exciton’s mobility.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2022-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42307771","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Monte Carlo method for parabolic equations involving fractional Laplacian 包含分数阶拉普拉斯式的抛物方程的蒙特卡罗方法
IF 0.9 Q4 Mathematics Pub Date : 2022-10-27 DOI: 10.1515/mcma-2022-2129
Caiyu Jiao, Changpin Li
Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.
摘要应用蒙特卡罗方法求解了具有分数阶拉普拉斯式的线性抛物型方程的Dirichlet问题。该方法利用了由具有跳跃的对称稳定lsamvy过程驱动的相关随机微分方程的弱逼近思想。我们利用跳跃适应方案来近似lsamvy过程,给出了精确的边界退出时间。当解的正则性较低时,通过去掉lsamvy过程的小跳变,建立了数值格式,并给出了收敛阶。当解具有较高的正则性时,我们用一个简单的过程代替小的跳跃,建立一个高阶的数值格式,然后显示更高的收敛阶。最后,给出了十维和一百维情况下的数值实验,验证了理论估计,并证明了所提格式对高维抛物方程的数值效率。
{"title":"Monte Carlo method for parabolic equations involving fractional Laplacian","authors":"Caiyu Jiao, Changpin Li","doi":"10.1515/mcma-2022-2129","DOIUrl":"https://doi.org/10.1515/mcma-2022-2129","url":null,"abstract":"Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2022-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44779092","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A random walk algorithm to estimate a lower bound of the star discrepancy 一种估计星差下界的随机游动算法
IF 0.9 Q4 Mathematics Pub Date : 2022-10-21 DOI: 10.1515/mcma-2022-2125
Maryam Alsolami, M. Mascagni
Abstract In many Monte Carlo applications, one can substitute the use of pseudorandom numbers with quasirandom numbers and achieve improved convergence. This is because quasirandom numbers are more uniform that pseudorandom numbers. The most common measure of that uniformity is the star discrepancy. Moreover, the main error bound in quasi-Monte Carlo methods, called the Koksma–Hlawka inequality, has the star discrepancy in the formulation. A difficulty with this bound is that computing the star discrepancy is very costly. The star discrepancy can be computed by evaluating a function called the local discrepancy at a number of points. The supremum of these local discrepancy values is the star discrepancy. If we have a point set in [ 0 , 1 ] s {[0,1]^{s}} with N members, we need to compute the local discrepancy at N s {N^{s}} points. In fact, computing star discrepancy is NP-hard. In this paper, we will consider an approximate algorithm for a lower bound on the star discrepancy based on using a random walk through some of the N s {N^{s}} points. This approximation is much less expensive that computing the star discrepancy, but still accurate enough to provide information on convergence. Our numerical results show that the random walk algorithm has the same convergence rate as the Monte Carlo method, which is O ( N - 1 2 {O(N^{-frac{1}{2}}} ).
摘要在许多蒙特卡罗应用中,可以用拟随机数代替伪随机数的使用,从而达到提高收敛性的目的。这是因为准随机数比伪随机数更均匀。这种均匀性最常见的测量方法是恒星差异。此外,拟蒙特卡罗方法的主要误差界,称为Koksma-Hlawka不等式,在公式中具有星形差异。这个界限的一个困难是,计算恒星差异的成本非常高。星形差异可以通过在若干点上计算一个称为局部差异的函数来计算。这些局部差值的最大值是星形差值。如果我们有一个在[0,1]s {[0,1]^{s}}中有N个成员的点集,我们需要计算N s {N^{s}}点上的局部差异。事实上,计算恒星差异是np困难的。在本文中,我们将考虑一种基于随机遍历一些N s {N^{s}}点的星差下界的近似算法。这种近似方法比计算恒星差异要便宜得多,但仍然足够精确,可以提供关于收敛的信息。我们的数值结果表明,随机漫步算法具有与蒙特卡罗方法相同的收敛速度,即O(N - 1 2 {O(N^{-frac{1}{2}}})。
{"title":"A random walk algorithm to estimate a lower bound of the star discrepancy","authors":"Maryam Alsolami, M. Mascagni","doi":"10.1515/mcma-2022-2125","DOIUrl":"https://doi.org/10.1515/mcma-2022-2125","url":null,"abstract":"Abstract In many Monte Carlo applications, one can substitute the use of pseudorandom numbers with quasirandom numbers and achieve improved convergence. This is because quasirandom numbers are more uniform that pseudorandom numbers. The most common measure of that uniformity is the star discrepancy. Moreover, the main error bound in quasi-Monte Carlo methods, called the Koksma–Hlawka inequality, has the star discrepancy in the formulation. A difficulty with this bound is that computing the star discrepancy is very costly. The star discrepancy can be computed by evaluating a function called the local discrepancy at a number of points. The supremum of these local discrepancy values is the star discrepancy. If we have a point set in [ 0 , 1 ] s {[0,1]^{s}} with N members, we need to compute the local discrepancy at N s {N^{s}} points. In fact, computing star discrepancy is NP-hard. In this paper, we will consider an approximate algorithm for a lower bound on the star discrepancy based on using a random walk through some of the N s {N^{s}} points. This approximation is much less expensive that computing the star discrepancy, but still accurate enough to provide information on convergence. Our numerical results show that the random walk algorithm has the same convergence rate as the Monte Carlo method, which is O ( N - 1 2 {O(N^{-frac{1}{2}}} ).","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2022-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41818243","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Superposition of forward and backward motion 前后运动叠加
IF 0.9 Q4 Mathematics Pub Date : 2022-10-19 DOI: 10.1515/mcma-2022-2124
Manfred Harringer
Abstract I consider black body radiation. The wall of the black body exchanges photons with the radiation field in equilibrium, therefore with a common temperature in Planck’s radiation law. The underlying process of radiation consists of creation and annihilation of photons. I want to present an alternate model of motions, where the process of radiation consists of small steps in positive and negative direction, not zero in mean. The detection of radiation consists of storing and restoring of packages of energy. I get an analogue of Planck’s radiation law, where the common temperature emerges from the underlying common model of small steps. The object of the law is not the radiation, but a storage of packages of energy, which belongs to the wall of the black body.
摘要我认为黑体辐射。黑体壁在平衡状态下与辐射场交换光子,因此具有普朗克辐射定律中的共同温度。辐射的基本过程包括光子的产生和湮灭。我想提出一个交替的运动模型,其中辐射过程由正负方向的小步组成,而不是平均值为零。辐射探测包括能量包的储存和恢复。我得到了普朗克辐射定律的类似物,其中共同的温度来自于小台阶的基本共同模型。法律的对象不是辐射,而是能量包的储存,它属于黑体的壁。
{"title":"Superposition of forward and backward motion","authors":"Manfred Harringer","doi":"10.1515/mcma-2022-2124","DOIUrl":"https://doi.org/10.1515/mcma-2022-2124","url":null,"abstract":"Abstract I consider black body radiation. The wall of the black body exchanges photons with the radiation field in equilibrium, therefore with a common temperature in Planck’s radiation law. The underlying process of radiation consists of creation and annihilation of photons. I want to present an alternate model of motions, where the process of radiation consists of small steps in positive and negative direction, not zero in mean. The detection of radiation consists of storing and restoring of packages of energy. I get an analogue of Planck’s radiation law, where the common temperature emerges from the underlying common model of small steps. The object of the law is not the radiation, but a storage of packages of energy, which belongs to the wall of the black body.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2022-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45461671","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
期刊
Monte Carlo Methods and Applications
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1