Two to Six Dimensional Numerical Solutions to the Poisson Eigenvalue Partial Differential Equation Using Generalized Multiquadrics
- 1 Convergent Solutions, Livermore, CA, USA
- 2 Advanpix, LLC, Yokohama, Japan
Abstract
We solve numerically an eigenvalue elliptic partial differential equation (PDE) ranging from two to six dimensions using the generalized multiquadric (GMQ) radial basis functions (RBFs). Two discretization methods are employed. The first method is similar to the classic mesh-based discretization method requiring n centers per dimension or a total n d points. The second method is based upon n randomly generated points in ℜ d requiring far fewer data centers than the classic mesh method. Instead of having a crisp boundary, we form a “fuzzy” boundary. Using the analytic or numerical inverse interior and boundary operators, we find the local and global minima and maxima to cull discretization points. We also find that the GMQ-RBF “flatness” can be controlled by increasing the GMQ exponential, β . We perform a search to find the smallest root mean squared error (RMSE) by varying the exponent, the maximum, the minimum range of the GMQ shape parameter, and boundary influence, with the exponential having the most influence. Because the GMQ-RBFs are essentially nonlinear, it follows that the starting point of the simple search influences the end result. The optimal algorithm for high dimensional PDEs is still the subject of much research and may wait for the common place availability of massively parallel quantum computers for even higher dimensional PDEs and integral equations (IEs).
- Bellman, R.E. and Rand Corporation (1957) Dynamic Programming. Princeton University Press.
- Naito, R. and Yamada, T. (2023) Deep High-Order Splitting Method for Semilinear Degenerate PDEs and Application to High-Dimensional Nonlinear Pricing Models. Digital Finance , 6, 693-725. https://doi.org/10.1007/s42521-023-00091-z
- Yu, W. and Mascagni, M. (2023) Monte Carlo Methods for Partial Differential Equations with Applications to Electronic Design Automation. Department of Computer Science, Florida State University. https://doi.org/10.1007/978-981-19-3250-2.
- Kuo, F.Y. and Nuyens, D. (2016) Application of Quasi-Monte Carlo Methods to Elliptic PDEs with Random Diffusion Coefficients: A Survey of Analysis and Implementation. Foundations of Computational Mathematics , 16, 1631-1696. https://doi.org/10.1007/s10208-016-9329-5
- Kim, H.H. and Yang, H.J. (2024) Domain Decomposition Algorithms for Neural Network Approximation of Partial Differential Equations. In: Dostál, Z., et al ., Eds., Domain Decomposition Methods in Science and Engineering XXVII , Springer, 27-37. https://doi.org/10.1007/978-3-031-50769-4_3
- Hu, Q., Basir, S. and Senocak, I. (2025) Non-Overlapping, Schwarz-Type Domain Decomposition Method for Physics and Equality Constrained Artificial Neural Networks. Computer Methods in Applied Mechanics and Engineering , 436, Article ID: 117706. https://doi.org/10.1016/j.cma.2024.117706
- Huang, J., Guo, W. and Cheng, Y. (2023) Adaptive Sparse Grid Discontinuous Galerkin Method: Review and Software Implementation. Communications on Applied Mathematics and Computation , 6, 501-532. https://doi.org/10.1007/s42967-023-00268-8
- Chen, W., Hwang, H. and Tsai, T. (2012) Maxima-Finding Algorithms for Multidimensional Samples: A Two-Phase Approach. Computational Geometry , 45, 33-53. https://doi.org/10.1016/j.comgeo.2011.08.001
- Hu, Z., Shukla, K., Karniaakis, G.E. and Kawaguchi, K. (2023) Tackling the Curse of Dimensionality with Physics-Informed Neural Networks. Neural Networks , 176, Article ID: 106368.
- Cohen, S.N., Jiang, D. and Sirignano, J. (2023) Curse of Dimensionality with Physics-Informed Neural Networks. Journal of Machine Learning Research , 24, 1-49.
- Sarra, S.A. and Kansa, E.J. (2009) Multiquadric Radial Basis Function Approximation Methods for the Numerical Solution of Partial Differential Equations. Computational Mechanics, Tech Science Press.