Timbwaoga A. J. Ouermi, Robert M Kirby, Martin Berzins
{"title":"算法xxxx: HiPPIS一种结构化网格高阶保正映射软件","authors":"Timbwaoga A. J. Ouermi, Robert M Kirby, Martin Berzins","doi":"arxiv-2310.08818","DOIUrl":null,"url":null,"abstract":"Polynomial interpolation is an important component of many computational\nproblems. In several of these computational problems, failure to preserve\npositivity when using polynomials to approximate or map data values between\nmeshes can lead to negative unphysical quantities. Currently, most\npolynomial-based methods for enforcing positivity are based on splines and\npolynomial rescaling. The spline-based approaches build interpolants that are\npositive over the intervals in which they are defined and may require solving a\nminimization problem and/or system of equations. The linear polynomial\nrescaling methods allow for high-degree polynomials but enforce positivity only\nat limited locations (e.g., quadrature nodes). This work introduces open-source\nsoftware (HiPPIS) for high-order data-bounded interpolation (DBI) and\npositivity-preserving interpolation (PPI) that addresses the limitations of\nboth the spline and polynomial rescaling methods. HiPPIS is suitable for\napproximating and mapping physical quantities such as mass, density, and\nconcentration between meshes while preserving positivity. This work provides\nFortran and Matlab implementations of the DBI and PPI methods, presents an\nanalysis of the mapping error in the context of PDEs, and uses several 1D and\n2D numerical examples to demonstrate the benefits and limitations of HiPPIS.","PeriodicalId":501256,"journal":{"name":"arXiv - CS - Mathematical Software","volume":"19 2","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2023-10-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Algorithm xxxx: HiPPIS A High-Order Positivity-Preserving Mapping Software for Structured Meshes\",\"authors\":\"Timbwaoga A. J. Ouermi, Robert M Kirby, Martin Berzins\",\"doi\":\"arxiv-2310.08818\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"Polynomial interpolation is an important component of many computational\\nproblems. In several of these computational problems, failure to preserve\\npositivity when using polynomials to approximate or map data values between\\nmeshes can lead to negative unphysical quantities. Currently, most\\npolynomial-based methods for enforcing positivity are based on splines and\\npolynomial rescaling. The spline-based approaches build interpolants that are\\npositive over the intervals in which they are defined and may require solving a\\nminimization problem and/or system of equations. The linear polynomial\\nrescaling methods allow for high-degree polynomials but enforce positivity only\\nat limited locations (e.g., quadrature nodes). This work introduces open-source\\nsoftware (HiPPIS) for high-order data-bounded interpolation (DBI) and\\npositivity-preserving interpolation (PPI) that addresses the limitations of\\nboth the spline and polynomial rescaling methods. HiPPIS is suitable for\\napproximating and mapping physical quantities such as mass, density, and\\nconcentration between meshes while preserving positivity. This work provides\\nFortran and Matlab implementations of the DBI and PPI methods, presents an\\nanalysis of the mapping error in the context of PDEs, and uses several 1D and\\n2D numerical examples to demonstrate the benefits and limitations of HiPPIS.\",\"PeriodicalId\":501256,\"journal\":{\"name\":\"arXiv - CS - Mathematical Software\",\"volume\":\"19 2\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2023-10-13\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - CS - Mathematical Software\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2310.08818\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - CS - Mathematical Software","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2310.08818","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Algorithm xxxx: HiPPIS A High-Order Positivity-Preserving Mapping Software for Structured Meshes
Polynomial interpolation is an important component of many computational
problems. In several of these computational problems, failure to preserve
positivity when using polynomials to approximate or map data values between
meshes can lead to negative unphysical quantities. Currently, most
polynomial-based methods for enforcing positivity are based on splines and
polynomial rescaling. The spline-based approaches build interpolants that are
positive over the intervals in which they are defined and may require solving a
minimization problem and/or system of equations. The linear polynomial
rescaling methods allow for high-degree polynomials but enforce positivity only
at limited locations (e.g., quadrature nodes). This work introduces open-source
software (HiPPIS) for high-order data-bounded interpolation (DBI) and
positivity-preserving interpolation (PPI) that addresses the limitations of
both the spline and polynomial rescaling methods. HiPPIS is suitable for
approximating and mapping physical quantities such as mass, density, and
concentration between meshes while preserving positivity. This work provides
Fortran and Matlab implementations of the DBI and PPI methods, presents an
analysis of the mapping error in the context of PDEs, and uses several 1D and
2D numerical examples to demonstrate the benefits and limitations of HiPPIS.