Pointwise Regularity in Analysis

• Pointwise regularity is the local study of how a function approximates a simpler model, such as a polynomial or wavelet decay, near a fixed point.
• It applies to various settings including PDEs, stochastic processes, and fractal analysis, and leverages methods like Hölder estimates and Besov decay.
• The framework uses techniques such as compactness, perturbation, and dyadic iteration to precisely capture local behavior in both interior and boundary contexts.

Pointwise regularity is the local study of how a function, a PDE solution, or a stochastic sample path is approximated near a fixed point by a simpler model, typically a polynomial, a homogeneous profile, or a multiscale scaling law. In the PDE literature this is commonly expressed by estimates of the form

$|u(x)-P(x)|\le K|x-x_0|^{k+\alpha},$

while in stochastic and fractal settings it is encoded by pointwise Hölder exponents, Besov-type local oscillation, or wavelet-coefficient decay. Recent work develops this perspective for fully nonlinear elliptic and parabolic equations, divergence-form equations, conical and rough domains, free boundaries, self-affine curves, and non-Gaussian processes, with an emphasis on identifying the precise local model rather than only a global modulus of continuity (Lian et al., 2020, Lian, 2022, Saka, 2015, Daw et al., 2022).

1. Definitions and local models

In the pointwise Schauder framework, a bounded function $f$ is $C^{k,\alpha}$ at $x_0$ if there exists a polynomial $P$ of degree $k$ such that

$|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$

near $x_0$. This definition is used systematically for interior and boundary regularity in fully nonlinear elliptic equations, including the boundary notions $u\in C^{1,\alpha}(x_0)$ and $u\in C^{2,\alpha}(x_0)$, and it is extended in general-form equations to endpoint classes $f$0 and $f$1 corresponding to $f$2 and $f$3 (Lian et al., 2019, Lian et al., 2020).

A boundary variant replaces approximation of the function alone by simultaneous approximation of the boundary geometry. In one standard formulation, $f$4 if, after a rigid motion, the boundary lies between two parallel tubes of thickness $f$5 around the graph of a polynomial $f$6. A further variant appears on cones: for viscosity solutions of $f$7 in $f$8, one says $f$9 if there exist a homogeneous profile $C^{k,\alpha}$0 and a scalar $C^{k,\alpha}$1 such that

$C^{k,\alpha}$2

near the vertex. In the half-space case $C^{k,\alpha}$3, this reduces to the usual $C^{k,\alpha}$4 boundary estimate (Lian, 2022, Lian et al., 2022).

In stochastic analysis, the local model is usually encoded by the pointwise Hölder exponent. For a process $C^{k,\alpha}$5, one says that $C^{k,\alpha}$6 is pointwise Hölder continuous of order $C^{k,\alpha}$7 at $C^{k,\alpha}$8 if there exists a polynomial $C^{k,\alpha}$9 of degree at most $x_0$0 such that

$x_0$1

and the exponent is

$x_0$2

For parameterized affine zipper curves one also distinguishes the liminf exponent

$x_0$3

from the regular exponent

$x_0$4

when the limit exists (Daw et al., 2022, Bárány et al., 2016).

A function-space formulation identifies pointwise regularity with local Besov decay. For self-affine lattice tilings, the pointwise Besov space $x_0$5 is defined by the boundedness of

$x_0$6

and one has the equivalence

$x_0$7

Consequently,

$x_0$8

is exactly the pointwise Hölder exponent (Saka, 2015).

2. Interior pointwise regularity in elliptic and parabolic PDE

A broad interior theory is developed for fully nonlinear elliptic equations in general forms,

$x_0$9

under a structure condition (SC2) that allows quadratic growth in the gradient. The model class includes equations such as

$P$0

and

$P$1

Within this framework, interior pointwise $P$2, $P$3, and $P$4 regularity are obtained, together with endpoint $P$5 and $P$6 results (Lian et al., 2020).

For divergence-form equations, an interior pointwise theory is formulated in Campanato form. In the parabolic setting, $P$7 means that there exists a parabolic polynomial $P$8 of degree $P$9 such that

$k$0

Under small-BMO assumptions on the leading coefficients and appropriate pointwise regularity of lower-order terms and data, weak solutions satisfy pointwise $k$1 and $k$2 estimates; if $k$3 vanishes to order $k$4, the required smoothness of the coefficients can be weakened by $k$5 or $k$6, depending on the term. The same framework yields a characterization of nodal sets: $k$7 with each stratum contained in a finite union of $k$8-dimensional $k$9 manifolds (Lian, 2024).

For elliptic and parabolic equations with divergence-free drifts,

$|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$0

interior pointwise $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$1 regularity for any $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$2 is proved under one of three scale-invariant smallness conditions on the drift: a Morrey-type $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$3 smallness, an $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$4-energy smallness, or a $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$5 smallness. The proof is based on the energy inequality and the perturbation technique, combined with compactness and Campanato iteration (Lian, 2024).

A distinct line treats locally uniformly elliptic equations, where uniform ellipticity holds only on bounded regions of $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$6-space. For equations

$|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$7

with $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$8 only $|f(x)-P(x)|\le K|x-x_0|^{k+\alpha}$9-uniformly elliptic, pointwise $x_0$0, $x_0$1, and $x_0$2 regularity are obtained under smallness assumptions on $x_0$3, on $x_0$4, or on the deviation from a model operator $x_0$5. The applications include the prescribed mean curvature equation, the Monge–Ampère equation, the $x_0$6-Hessian equations, the $x_0$7-Hessian quotient equations, and the Lagrangian mean curvature equation; the accompanying remarks state that the smallness assumptions are necessary in most cases (Lian et al., 2024).

3. Boundary, singular geometry, and free boundaries

Boundary pointwise regularity sharpens classical boundary Schauder theory by requiring only pointwise geometric control of the boundary at a single point. For fully nonlinear elliptic equations,

$x_0$8

if $x_0$9, $u\in C^{1,\alpha}(x_0)$0, and $u\in C^{1,\alpha}(x_0)$1 satisfies a Morrey-type bound, then $u\in C^{1,\alpha}(x_0)$2; if $u\in C^{1,\alpha}(x_0)$3, $u\in C^{1,\alpha}(x_0)$4, and $u\in C^{1,\alpha}(x_0)$5, then $u\in C^{1,\alpha}(x_0)$6 (Lian et al., 2019). A parabolic analogue for

$u\in C^{1,\alpha}(x_0)$7

establishes boundary pointwise $u\in C^{1,\alpha}(x_0)$8 regularity for any $u\in C^{1,\alpha}(x_0)$9, using flat-boundary model problems, compactness, and scaling (Lian et al., 2022).

For Dirichlet and oblique derivative problems with pointwise regular data, a higher-order boundary theory proves that if $u\in C^{2,\alpha}(x_0)$0, $u\in C^{2,\alpha}(x_0)$1, $u\in C^{2,\alpha}(x_0)$2, and $u\in C^{2,\alpha}(x_0)$3, then $u\in C^{2,\alpha}(x_0)$4. In the oblique case, if

$u\in C^{2,\alpha}(x_0)$5

then $u\in C^{2,\alpha}(x_0)$6, $u\in C^{2,\alpha}(x_0)$7, and $u\in C^{2,\alpha}(x_0)$8 imply $u\in C^{2,\alpha}(x_0)$9. These results are then applied to obstacle-type and one-phase problems to obtain $f$00 regularity of free boundaries from $f$01 regularity (Lian et al., 2022).

The conical setting replaces Taylor polynomials by the homogeneous profile dictated by the geometry. Let $f$02 be a $f$03-smooth open cap and $f$04 the cone based on $f$05. For bounded viscosity solutions of

$f$06

with $f$07 uniformly elliptic and positively $f$08-homogeneous, there exist $f$09, $f$10, and $f$11 such that

$f$12

where $f$13 is the known homogeneous solution vanishing on $f$14. This estimate leads directly to Liouville theorems on cones: if $f$15 at infinity, then $f$16 is a multiple of $f$17; if $f$18, then $f$19 (Lian, 2022).

Boundary pointwise regularity has also been extended beyond smooth domains. On $f$20-uniform domains, a notion of weak solution with nonzero boundary data is introduced through the auxiliary functions

$f$21

together with subsolution inequalities. Under an admissibility condition and pointwise regularity assumptions $f$22, $f$23, boundary pointwise $f$24 regularity follows, and linearity with respect to harmonic functions yields $f$25 and $f$26 under stronger assumptions on the boundary and data (Guan et al., 22 Sep 2025). For divergence-form equations with distributional coefficients, if a boundary point satisfies a measure condition, the multiplier $f$27 is form-bounded with sufficiently small constant, and the nonhomogeneous terms satisfy Dini decay, then the weak solution is continuous there in the $f$28-sense, with an $f$29-mean oscillation estimate implying pointwise Hölder regularity when the Dini remainder is integrable (Jingqi et al., 2024).

Free-boundary regularity is another setting where pointwise expansions are decisive. In the parabolic obstacle problem

$f$30

Lindgren and Monneau derive a Weiss-type monotonicity formula, a singular-point monotonicity formula, and a second-order Taylor expansion at singular free boundary points under pointwise Dini and double-Dini conditions. Under Dini continuity of $f$31, the regular set is locally a parabolic $f$32-surface, while the singular set is locally contained in a union of parabolic $f$33 manifolds (Lindgren et al., 2013).

4. Wavelet, Besov, and multifractal formulations

In self-affine and multiresolution settings, pointwise regularity is encoded by decay of local oscillation or wavelet coefficients. For self-affine lattice tilings associated with an integer dilation matrix $f$34, local oscillation is measured by

$f$35

and the corresponding Besov norm is

$f$36

Under multiresolution analysis, one also has a wavelet characterization, and pointwise regularity at $f$37 is equivalent to decay of nearby wavelet coefficients: $f$38 whenever $f$39 (Saka, 2015).

For parameterized affine zipper fractal curves, dominated splitting of index-$f$40 provides a dynamical mechanism for well-defined local anisotropy. The pressure function $f$41, defined as the unique real root of

$f$42

is continuous, strictly increasing, concave, and $f$43. Under SOSC, non-degeneracy, and dominated splitting, the Hausdorff dimensions of level sets of the pointwise Hölder exponent are given, on the stated ranges, by the Legendre-transform formula

$f$44

Under Assumption A, the same formula extends to the full spectrum for the regular exponent $f$45, and Assumption A is equivalent to the existence of $f$46 for Lebesgue-a.e. point. The de Rham curve is a concrete application: for $f$47, Assumption A holds, the curve is differentiable Lebesgue-a.e. with derivative zero, and the nondifferentiability set has dimension $f$48, where $f$49 (Bárány et al., 2016).

Wavelet methods can also be used algorithmically. The Iterated Amplitude Adjusted Wavelet Transform preserves the data histogram and the wavelet-coefficient magnitudes at each scale and location. Because the estimate

$f$50

links the pointwise Hölder exponent to the modulus of wavelet coefficients, preserving those magnitudes preserves the pointwise Hölder regularity structure. The method is designed to preserve the multifractal properties of the data while randomizing the spatial arrangement of singularities, and it is used for testing oscillating singularities and velocity–intermittency coupling in turbulence (Keylock, 2017).

5. Stochastic processes and refined pointwise oscillation

For stochastic processes, pointwise regularity is often finer than the sole value of the Hölder exponent. In the generalized Rosenblatt process $f$51, with self-similarity exponent $f$52, wavelet methods identify three types of points on any interval. Ordinary points satisfy a modulus involving $f$53, rapid points a stronger modulus involving $f$54, and slow points satisfy a modulus with no logarithm. The wavelet-leader characterization is

$f$55

and the proof uses a Meyer-wavelet expansion, kernel cancellation for $f$56, and tail bounds in the second Wiener chaos. Compared with fractional Brownian motion, the logarithmic corrections differ by square-roots, reflecting the replacement of Gaussian first-chaos tails by exponential second-chaos tails (Daw et al., 2022).

For the multifractional Brownian motion $f$57, the pointwise exponent is governed exactly by the Hurst function. Under Condition 1.3, namely that for each $f$58 there exist $f$59 and constants $f$60 such that

$f$61

one has almost surely

$f$62

Moreover, every non-empty open interval contains at least one slow point $f$63, meaning that

$f$64

The same paper shows that this property persists for a non self-similar process and that the assumption on $f$65 can be weakened to the single-point logarithmic modulus

$f$66

Under this weaker condition, multifractional wavelet series and locally-deforming wavelet expansions still satisfy $f$67 almost surely and exhibit the same slow/ordinary/rapid trichotomy (Esser et al., 2023).

A common misconception is that the pointwise exponent alone fully describes local behavior. The Rosenblatt and multifractional Brownian results show that this is false: points with the same exponent may still differ by logarithmic corrections in their modulus of continuity. The distinction between ordinary, rapid, and slow points is precisely a distinction at fixed exponent (Daw et al., 2022, Esser et al., 2023).

6. Methods, structural consequences, and neighboring notions

Across these works, the main proofs are scale-by-scale. In PDE, the dominant pattern is compactness plus perturbation: one subtracts an approximating polynomial or profile, rescales, proves precompactness, passes to a model limit equation, and iterates. This appears in fully nonlinear elliptic and parabolic settings, in divergence form, and in boundary problems, often with Campanato decay as the output norm (Lian, 2024, Lian et al., 2022, Lian et al., 2020). On cones, the same philosophy is combined with a Hopf lemma, a boundary Lipschitz estimate, an improvement-of-oscillation lemma in rings, and dyadic iteration (Lian, 2022). For fully fractional parabolic equations $f$68, the proof instead uses an integral representation, a decomposition $f$69, a directional-average perturbation estimate for the fractional heat kernel, and new equivalent definitions of pointwise function spaces, yielding pointwise $f$70 or logarithmic regularity according to whether $f$71 or $f$72 (Guo et al., 8 Mar 2026).

Pointwise estimates can also enforce geometric structure. In generated Jacobian equations, Guillen and Kitagawa prove an Aleksandrov estimate and a Sharp-Growth estimate under the analogue of the Ma–Trudinger–Wang condition, denoted $f$73. These estimates imply strict $f$74-convexity and, when $f$75 is $f$76 in the $f$77-variable, interior

$f$78

The theory applies to the near-field reflector problem, where the reciprocal radial function $f$79 is represented through a generating function built from ellipsoids, and the resulting reflector surface is $f$80 (Guillen et al., 2015).

Several papers explicitly show that pointwise regularity is stronger than coarse Hölder control. In the conical fully nonlinear theory, the estimate $f$81 is described as strictly stronger than global $f$82 or even local $f$83 bounds because it identifies the exact homogeneous blow-up limit at the vertex (Lian, 2022). Likewise, smallness assumptions in locally uniformly elliptic problems are not treated as purely technical: the stated remarks assert that they are necessary in most cases, since otherwise the blow-up operators may leave the uniformly elliptic regime (Lian et al., 2024).

A neighboring but distinct theme is pointwise convergence rather than pointwise regularity. For dispersive equations on compact symmetric spaces $f$84 of rank $f$85 or $f$86, the Sobolev threshold $f$87 is sufficient for almost-everywhere convergence of

$f$88

to $f$89 as $f$90. In the $f$91-biinvariant rank-$f$92 case, the sufficiency improves to $f$93 for the Schrödinger, Boussinesq, and Beam equations, while no maximal-function estimate can hold for $f$94 (Dewan et al., 10 Dec 2025). This does not define pointwise regularity in the usual Hölder or polynomial-approximation sense, but it shows how local convergence questions interact with Sobolev smoothness thresholds.

Taken together, these results suggest that pointwise regularity is less a single theorem than a family of local asymptotic principles. Depending on the problem, the asymptotic model may be a Taylor polynomial, a homogeneous boundary profile, a wavelet-decay law, a Besov oscillation rate, or a modulus corrected by logarithms. What remains constant is the objective: to identify the exact local scale at which a function begins to resemble its canonical model.