ISSN 2738-0971 | eISSN 2738-1013

PERFORMANCE, ACCURACY, AND PREDICTIVE CHUNK OPTIMIZATION OF PARALLEL SIMPSON AND TRAPEZOIDAL INTEGRATION METHODS

Authors

Sanja Đurović ORCID 0000-0001-7847-9953
Faculty of Sciences and Mathematics, University of Priština in Kosovska Mitrovica, Kosovska Mitrovica, Serbia
Marko Smilić ORCID 0000-0003-2227-2536
Faculty of Sciences and Mathematics, University of Priština in Kosovska Mitrovica, Kosovska Mitrovica, Serbia
Časlav Stefanović ORCID 0000-0003-4385-3356
Faculty of Sciences and Mathematics, University of Priština in Kosovska Mitrovica, Kosovska Mitrovica, Serbia
Stefan Panić ORCID 0000-0002-5868-1764
Faculty of Sciences and Mathematics, University of Priština in Kosovska Mitrovica, Kosovska Mitrovica, Serbia
Milan Savić ORCID 0000-0002-7692-001X
Faculty of Sciences and Mathematics, University of Priština in Kosovska Mitrovica, Kosovska Mitrovica, Serbia

Keywords

Parallel numerical integration, Simpson’s rule, Trapezoidal rule, Adaptive quadrature, Chunk optimization

Abstract

Parallel numerical integration is a fundamental component of many scientific and engineering applications, yet its performance depends strongly on workload partitioning strategies. In shared-memory environments, the choice of chunk size directly affects load balancing, scheduling overhead, and cache utilization, particularly when adaptive numerical methods are employed. This paper presents an analysis of the performance, accuracy, and scalability of four numerical integration techniques: Simpson’s rule, Adaptive Simpson’s rule, the Trapezoidal rule, and the Adaptive Trapezoidal rule, implemented in a multi-threaded execution model. Execution-time distributions are analyzed across a wide range of chunk sizes and thread counts, revealing substantial performance variability and method-dependent optimal chunk configurations. In addition to performance evaluation, numerical accuracy is assessed to highlight trade-offs between computational efficiency and error reduction. Based on empirical observations, optimal chunk ranges are identified and used to construct a predictive model that estimates near-optimal chunk configurations as a function of the number of threads. The proposed approach enables efficient parameter selection without exhaustive tuning and demonstrates the feasibility of predictive chunk optimization for parallel numerical integration.

Published
2026/04/21
Issue
Vol. 16 No. 1 (2026)
Section
Original Scientific Paper

Cite this article

Checking

Citation styles are being prepared.