A fast multipole method for stellar dynamics

IF 16.281
Walter Dehnen
{"title":"A fast multipole method for stellar dynamics","authors":"Walter Dehnen","doi":"10.1186/s40668-014-0001-7","DOIUrl":null,"url":null,"abstract":"<p>The approximate computation of all gravitational forces between <i>N</i> interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than <math><mi>O</mi><mo>(</mo><mi>N</mi><mo>)</mo></math> operations. FMM groups particles into spatially bounded cells and uses cell-cell interactions to approximate the force at <i>any</i> position within the sink cell by a Taylor expansion obtained from the multipole expansion of the source cell. By employing a novel estimate for the errors incurred in this process, I minimise the computational effort required for a given accuracy and obtain a well-behaved distribution of force errors. For relative force errors of ~10<sup>?7</sup>, the computational costs exhibit an empirical scaling of <math><mo>∝</mo><msup>\n <mi>N</mi>\n <mn>0.87</mn>\n </msup></math>. My implementation (running on a 16 core node) out-performs a GPU-based direct summation with comparable force errors for <math><mi>N</mi><mo>?</mo><msup>\n <mrow>\n <mn>10</mn>\n </mrow>\n <mrow>\n <mn>5</mn>\n </mrow>\n </msup></math>.</p>","PeriodicalId":523,"journal":{"name":"Computational Astrophysics and Cosmology","volume":null,"pages":null},"PeriodicalIF":16.2810,"publicationDate":"2014-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1186/s40668-014-0001-7","citationCount":"53","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computational Astrophysics and Cosmology","FirstCategoryId":"4","ListUrlMain":"https://link.springer.com/article/10.1186/s40668-014-0001-7","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 53

Abstract

The approximate computation of all gravitational forces between N interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than O(N) operations. FMM groups particles into spatially bounded cells and uses cell-cell interactions to approximate the force at any position within the sink cell by a Taylor expansion obtained from the multipole expansion of the source cell. By employing a novel estimate for the errors incurred in this process, I minimise the computational effort required for a given accuracy and obtain a well-behaved distribution of force errors. For relative force errors of ~10?7, the computational costs exhibit an empirical scaling of N 0.87 . My implementation (running on a 16 core node) out-performs a GPU-based direct summation with comparable force errors for N? 10 5 .

Abstract Image

恒星动力学的快速多极方法
用快速多极法(FMM)近似计算N个相互作用粒子之间的所有引力,可以达到直接求和的精度,但需要的操作少于O(N)次。FMM将粒子分组到空间有界的单元中,并利用单元间的相互作用,通过从源单元的多极展开得到的泰勒展开来近似吸收单元内任何位置的力。通过对这一过程中产生的误差采用新的估计,我最大限度地减少了给定精度所需的计算工作量,并获得了良好的力误差分布。相对力误差为~10?7、计算成本的经验标度为∝N 0.87。我的实现(在16核节点上运行)优于基于gpu的直接求和,对于N?5 . b。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
自引率
0.00%
发文量
0
期刊介绍: Computational Astrophysics and Cosmology (CompAC) is now closed and no longer accepting submissions. However, we would like to assure you that Springer will maintain an archive of all articles published in CompAC, ensuring their accessibility through SpringerLink's comprehensive search functionality.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信