{"title":"Meshless weighting coefficients for arbitrary nodes: The efficient computation to machine precision using hyper-dual numbers","authors":"Jason L. Roberts","doi":"10.1016/j.advengsoft.2024.103753","DOIUrl":null,"url":null,"abstract":"<div><p>A computationally efficient algorithm to calculate the weighting coefficients required to evaluate derivatives for arbitrary multi-dimensional distributions of points is presented. The iterative algorithm guarantees IEEE 754 64-bit precision (at least 15 significant decimal digits) for the weighting coefficients. Convergence acceleration is achieved through the use of a Taylor series of up to third order, and hyper-dual numbers to obtain the derivatives required for the Taylor series. The method is applied as part of a finite point solution for three test examples, a Poisson equation, creeping flow around a cylinder, and heat conduction in a triangular annulus. The open source FORTRAN-90 implementation has been optimised for random distributions of points in 1 to 3 dimensions.</p></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"197 ","pages":"Article 103753"},"PeriodicalIF":4.0000,"publicationDate":"2024-08-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Advances in Engineering Software","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0965997824001601","RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
A computationally efficient algorithm to calculate the weighting coefficients required to evaluate derivatives for arbitrary multi-dimensional distributions of points is presented. The iterative algorithm guarantees IEEE 754 64-bit precision (at least 15 significant decimal digits) for the weighting coefficients. Convergence acceleration is achieved through the use of a Taylor series of up to third order, and hyper-dual numbers to obtain the derivatives required for the Taylor series. The method is applied as part of a finite point solution for three test examples, a Poisson equation, creeping flow around a cylinder, and heat conduction in a triangular annulus. The open source FORTRAN-90 implementation has been optimised for random distributions of points in 1 to 3 dimensions.
期刊介绍:
The objective of this journal is to communicate recent and projected advances in computer-based engineering techniques. The fields covered include mechanical, aerospace, civil and environmental engineering, with an emphasis on research and development leading to practical problem-solving.
The scope of the journal includes:
• Innovative computational strategies and numerical algorithms for large-scale engineering problems
• Analysis and simulation techniques and systems
• Model and mesh generation
• Control of the accuracy, stability and efficiency of computational process
• Exploitation of new computing environments (eg distributed hetergeneous and collaborative computing)
• Advanced visualization techniques, virtual environments and prototyping
• Applications of AI, knowledge-based systems, computational intelligence, including fuzzy logic, neural networks and evolutionary computations
• Application of object-oriented technology to engineering problems
• Intelligent human computer interfaces
• Design automation, multidisciplinary design and optimization
• CAD, CAE and integrated process and product development systems
• Quality and reliability.