This paper proposes an accurate and computationally efficient solver of the heat equation on a volumetric domain, through the (r,r)-degree Padé-Chebyshev rational approximation of the ex- ponential representation of the solution. To this end, the heat diffusion problem is converted to a set of r differential equations, which involve only the Laplace-Beltrami operator. The discrete heat equation is equivalent to r sparse, symmetric linear systems and is independent of the volume discretization as a tetrahedral mesh or a regular grid, the evaluation of the Laplacian spectrum, and the selection of a subset of eigenpairs. Our approach has a super-linear computational cost, is free of user-defined parameters, and has an approximation accuracy lower than 10^(-r). Finally, we propose a simple criterion to select the time value that provides the best compromise between approximation accuracy and smoothness of the solution.

Volumetric heat Kernel: Padé-Chebyshev approximation, convergence, and computation

2015

Abstract

This paper proposes an accurate and computationally efficient solver of the heat equation on a volumetric domain, through the (r,r)-degree Padé-Chebyshev rational approximation of the ex- ponential representation of the solution. To this end, the heat diffusion problem is converted to a set of r differential equations, which involve only the Laplace-Beltrami operator. The discrete heat equation is equivalent to r sparse, symmetric linear systems and is independent of the volume discretization as a tetrahedral mesh or a regular grid, the evaluation of the Laplacian spectrum, and the selection of a subset of eigenpairs. Our approach has a super-linear computational cost, is free of user-defined parameters, and has an approximation accuracy lower than 10^(-r). Finally, we propose a simple criterion to select the time value that provides the best compromise between approximation accuracy and smoothness of the solution.
2015
Istituto di Matematica Applicata e Tecnologie Informatiche - IMATI -
Diffusion geometry
Laplacian spectrum
Shape analysis
Volumetric heat kernel
File in questo prodotto:
File Dimensione Formato  
prod_308677-doc_102027.pdf

solo utenti autorizzati

Descrizione: Volumetric heat Kernel: Padé-Chebyshev approximation, convergence, and computation
Tipologia: Versione Editoriale (PDF)
Dimensione 2.77 MB
Formato Adobe PDF
2.77 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.14243/275622
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 4
social impact