Explicit solutions of Genz test integrals

IF 2.8 2区 数学 Q1 MATHEMATICS, APPLIED Applied Mathematics Letters Pub Date : 2025-04-01 Epub Date: 2024-12-26 DOI:10.1016/j.aml.2024.109444
Vesa Kaarnioja
{"title":"Explicit solutions of Genz test integrals","authors":"Vesa Kaarnioja","doi":"10.1016/j.aml.2024.109444","DOIUrl":null,"url":null,"abstract":"<div><div>A collection of test integrals introduced by Genz (1984) has remained popular to this day for assessing the robustness of high-dimensional numerical integration algorithms. However, the explicit solutions to these integrals do not appear to be readily available in the existing literature: typically the true values of the test integrals are simply approximated using “overkill” numerical solutions. In this paper, analytic solutions are presented for the Genz test integrals <span><span><span><math><mrow><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>⋯</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>cos</mo><mrow><mo>(</mo><mrow><mn>2</mn><mi>π</mi><msub><mrow><mi>w</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>+</mo><munderover><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mo>)</mo></mrow><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>⋯</mo><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><msup><mrow><mn>2</mn></mrow><mrow><mi>d</mi></mrow></msup><mo>cos</mo><mrow><mo>(</mo><mrow><mn>2</mn><mi>π</mi><msub><mrow><mi>w</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>+</mo><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac><munderover><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mo>)</mo></mrow><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><mfrac><mrow><mo>sin</mo><mrow><mo>(</mo><mfrac><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mn>2</mn></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><mo>,</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>⋯</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><mfrac><mrow><mn>1</mn></mrow><mrow><msubsup><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo>−</mo><mn>2</mn></mrow></msubsup><mo>+</mo><msup><mrow><mrow><mo>(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>−</mo><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mn>2</mn></mrow></msup></mrow></mfrac><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>⋯</mo><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mrow><mo>(</mo><mrow><mo>arctan</mo><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><mo>arctan</mo><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>−</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>,</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>⋯</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><msup><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><munderover><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mrow><mo>(</mo><mi>d</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>⋯</mo><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn>1</mn></mrow><mrow><mi>d</mi><mo>!</mo><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><munder><mrow><mo>∑</mo></mrow><mrow><mi>u</mi><mo>⊆</mo><mrow><mo>{</mo><msub><mrow><mi>c</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>}</mo></mrow></mrow></munder><mfrac><mrow><msup><mrow><mrow><mo>(</mo><mo>−</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>#</mi><mi>u</mi></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><munder><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>∈</mo><mi>u</mi></mrow></munder><mi>i</mi></mrow></mfrac><mo>,</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>⋯</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>exp</mo><mrow><mo>(</mo><mrow><mo>−</mo><munderover><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msubsup><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn>2</mn></mrow></msubsup><msup><mrow><mrow><mo>(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>−</mo><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mn>2</mn></mrow></msup></mrow><mo>)</mo></mrow><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>⋯</mo><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msup><mrow><mi>π</mi></mrow><mrow><mi>d</mi><mo>/</mo><mn>2</mn></mrow></msup></mrow><mrow><msup><mrow><mn>2</mn></mrow><mrow><mi>d</mi></mrow></msup></mrow></mfrac><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><mfrac><mrow><mi>erf</mi><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><mi>erf</mi><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>−</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><mo>,</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>⋯</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>exp</mo><mrow><mo>(</mo><mrow><mo>−</mo><munderover><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mrow><mo>|</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>−</mo><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>|</mo></mrow></mrow><mo>)</mo></mrow><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>⋯</mo><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><mfrac><mrow><mo>exp</mo><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>−</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>−</mo><mo>exp</mo><mrow><mo>(</mo><mo>−</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><mo>,</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><msub><mrow><mi>w</mi></mrow><mrow><mn>1</mn></mrow></msub></mrow></munderover><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><msub><mrow><mi>w</mi></mrow><mrow><mn>2</mn></mrow></msub></mrow></munderover><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>⋯</mo><munderover><mrow><mo>∫</mo></mrow><mrow><mn>0</mn></mrow><mrow><mn>1</mn></mrow></munderover><mo>exp</mo><mrow><mo>(</mo><mrow><munderover><mrow><mo>∑</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mo>)</mo></mrow><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mi>d</mi></mrow></msub><mo>⋯</mo><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>3</mn></mrow></msub><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>2</mn></mrow></msub><mspace></mspace><mi>d</mi><msub><mrow><mi>x</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn></mrow></munderover><mrow><mo>(</mo><mo>exp</mo><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>−</mo><mn>1</mn><mo>)</mo></mrow><mfrac><mrow><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>3</mn></mrow><mrow><mi>d</mi></mrow></munderover><mrow><mo>(</mo><mo>exp</mo><mrow><mo>(</mo><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>−</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><munderover><mrow><mo>∏</mo></mrow><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>d</mi></mrow></munderover><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><mo>,</mo></mrow></math></span></span></span> where <span><math><mrow><mi>d</mi><mo>∈</mo><msub><mrow><mi>Z</mi></mrow><mrow><mo>+</mo></mrow></msub></mrow></math></span>, <span><math><mrow><mn>0</mn><mo>&lt;</mo><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>&lt;</mo><mn>1</mn></mrow></math></span>, and <span><math><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>∈</mo><msub><mrow><mi>R</mi></mrow><mrow><mo>+</mo></mrow></msub></mrow></math></span> for all <span><math><mrow><mi>i</mi><mo>∈</mo><mrow><mo>{</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mi>d</mi><mo>}</mo></mrow></mrow></math></span>.</div></div>","PeriodicalId":55497,"journal":{"name":"Applied Mathematics Letters","volume":"163 ","pages":"Article 109444"},"PeriodicalIF":2.8000,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Applied Mathematics Letters","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0893965924004646","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2024/12/26 0:00:00","PubModel":"Epub","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0

Abstract

A collection of test integrals introduced by Genz (1984) has remained popular to this day for assessing the robustness of high-dimensional numerical integration algorithms. However, the explicit solutions to these integrals do not appear to be readily available in the existing literature: typically the true values of the test integrals are simply approximated using “overkill” numerical solutions. In this paper, analytic solutions are presented for the Genz test integrals 0101cos(2πw1+i=1dcixi)dxddx1=2dcos(2πw1+12i=1dci)i=1dsin(ci2)ci,0101i=1d1ci2+(xiwi)2dxddx1=i=1dci(arctan(ciwi)+arctan(ciciwi)),0101(1+i=1dcixi)(d+1)dxddx1=1d!i=1dciu{c1,,cd}(1)#u1+iui,0101exp(i=1dci2(xiwi)2)dxddx1=πd/22di=1derf(ciwi)+erf(ciciwi)ci,0101exp(i=1dci|xiwi|)dxddx1=i=1dexp(ciwici)exp(ciwi)ci,0w10w20101exp(i=1dcixi)dxddx3dx2dx1=i=12(exp(ciwi)1)i=3d(exp(ci)1)i=1dci, where dZ+, 0<wi<1, and ciR+ for all i{1,,d}.
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
Genz测试积分的显式解
Genz(1984)引入的一组测试积分至今仍很流行,用于评估高维数值积分算法的鲁棒性。然而,在现有的文献中,这些积分的显式解似乎并不容易得到:通常,测试积分的真实值只是简单地用“过度”数值解来近似。在本文中,分析解决方案提出了贞测试积分∫01⋯∫01 cos(2πw1 +∑i = 1 dcixi) dxd⋯dx1 = 2 dco(2πw1 + 12∑我dci = 1)∏i = 1 dsin ci2 ci,∫01⋯∫01∏i = 1 d1ci−2 +(ξ−wi) 2 dxd⋯dx1 dci =∏i = 1(反正切(ciwi) +反正切(ci−ciwi)), 01⋯∫∫01(1 +∑i = 1 dcixi)−(d + 1) dxd⋯dx1 = 1 d !∏i = 1 dci∑u⊆{c1,…,cd}(−1)# u1 +∑我∈ui, 01⋯∫∫01 exp(−∑i = 1 dci2(ξ−wi) 2) dxd⋯dx1 22 =πd / d∏i = 1 derf (ciwi) +小块土地(ci−ciwi) ci, 01⋯∫∫01 exp(−∑i = 1 dci | xi−wi |) dxd⋯dx1 =∏i = 1 dexp (ciwi−ci)−exp(−ciwi) ci, w1 0 w2∫∫∫0 01⋯∫01 exp(∑i = 1 dcixi) dxd⋯dx3dx2dx1 =∏i = 12 (exp (ciwi)−1)∏i = 3 d (exp (ci)−1)∏我dci = 1,d∈Z +, 0 & lt; wi< ci∈R + 1,我∈{1,…,d}。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 去求助
来源期刊
Applied Mathematics Letters
Applied Mathematics Letters 数学-应用数学
CiteScore
7.70
自引率
5.40%
发文量
347
审稿时长
10 days
期刊介绍: The purpose of Applied Mathematics Letters is to provide a means of rapid publication for important but brief applied mathematical papers. The brief descriptions of any work involving a novel application or utilization of mathematics, or a development in the methodology of applied mathematics is a potential contribution for this journal. This journal''s focus is on applied mathematics topics based on differential equations and linear algebra. Priority will be given to submissions that are likely to appeal to a wide audience.
期刊最新文献
Inverse spectral problem for glassy state relaxation approximated by Prony series Global weak solutions of a macroscopic model of traffic flow with a source A note on the averaging principle for ordinary differential equations depending on the slow time On the minimal wave speed for a multi-group infection-age-structured epidemic model incorporating vaccination and spatial diffusion Blow-up of solutions of the hyperbolic Prandtl system
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1