What is the method of interpolation in the function pdeinterpolant?

4 views (last 30 days)
What is the numerical method being used for interpolation in the function "pdeinterpolant"? I couldn't find any info on the actual algorithm being used (e.g. using the command doc pdeinterpolant). In particular I'd like to know whether I can expect C0 or C1 continuity in the interpolant when used on 2D quadratic triangulations.

Answers (1)

Abhas
Abhas on 30 Sep 2024
Hi Alex,
The "pdeInterpolant" function in MATLAB is used to interpolate solutions of partial differential equations (PDEs) at arbitrary points within the geometry. However, the specific numerical method used for interpolation in pdeInterpolant is not explicitly documented in the MATLAB documentation.
From general principles of finite element methods (FEM) and interpolation, we can infer some details:
  • Continuity: For 2D quadratic triangulations, you can typically expect "C0 continuity". This means the interpolant is continuous, but its first derivative may not be continuous across element boundaries. While "C1 continuity" requires both the function and its first derivative to be continuous across element boundaries. Achieving "C1 continuity" typically involves using more complex elements or additional constraints, which is uncommon for standard quadratic elements in finite element methods.
  • Interpolation Method: The interpolation is likely based on shape functions used in "FEM". For quadratic elements, these shape functions is usually done using second-degree polynomials, which provide a higher degree of accuracy compared to linear elements.
You may refer to the below MathWorks documentation links if you need more precise control over the interpolation method or continuity:
  1. https://www.mathworks.com/help/pde/ug/pde.stationaryresults.interpolatesolution.html
  2. https://www.mathworks.com/help/pde/ug/pde.stationaryresults.evaluategradient.html

Categories

Find more on Genomics and Next Generation Sequencing in Help Center and File Exchange

Products


Release

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!