Error Analysis of Lagrange Interpolation on Tetrahedrons
This paper describes the analysis of Lagrange interpolation errors on tetrahedrons. In many textbooks, the error analysis of Lagrange interpolation is conducted under geometric assumptions such as shape regularity or the (generalized) maximum angle condition. In this paper, we present a new estimation in which the error is bounded in terms of the diameter and projected circumradius of the tetrahedron. Because we do not impose any geometric restrictions on the tetrahedron itself, our error estimation may be applied to any tetrahedralizations of domains including very thin tetrahedrons.
\section{Introduction}\label{intro}
Lagrange interpolation on tetrahedrons and the associated error
estimates are important subjects in numerical analysis. In particular,
they are crucial in the error analysis of finite element methods.
Throughout this paper, $K \subset \mathbb R^3$ denotes a tetrahedron
with vertices $\mathbf{x}_i$, $i=1,\cdots,4$, and tetrahedrons are assumed to
be closed sets. Let $\lambda_i$ be the barycentric coordinates of a
tetrahedron with respect to $\mathbf{x}_i$. By definition,
we have $0 \le \lambda_i \le 1$, $\sum_{i=1}^{4} \lambda_i =1$.  Let $\mathbb N_0$
be the set of nonnegative integers, and
$\gamma = (a_1,\cdots,a_{4}) \in \mathbb N_0^{4}$ be a multi-index.
Let $k$ be a positive integer. If $|\gamma| := \sum_{i=1}^{4}a_i = k$, then
${\gamma}/{k} := \left({a_1}/{k}, \cdots, {a_{4}}/{k}\right)$
can be regarded as a barycentric coordinate in $K$.
The set $\Sigma^k(K)$ of points on $K$ is defined by
\begin{equation*}
   \Sigma^k(K) := \left\{
   \frac{\gamma}{k} \in K \Bigm| |\gamma| = k, \;
     \gamma \in \mathbb N_0^{4} \right\}.
\end{equation*}
For $k$, $\mathcal{P}_k$ is the set of all polynomials of three variables
whose degree is at most $k$.  For a continuous function $v \in C^0(K)$, 
the Lagrange interpolation $\mathcal{I}_K^k v \in \mathcal{P}_k$ of degree $k$ is defined as
\[
   v(\mathbf{x}) = (\mathcal{I}_K^k v)(\mathbf{x}), \quad \forall \mathbf{x} \in \Sigma^k(K).
\]
We attempt to obtain an upper bound of the error
$|v - \mathcal{I}_K^k v|_{m,p,K}$ for integers $0 \le m \le k$, where
$|\cdot|_{m,p,K}$ is the usual Sobolev semi-norm.

For now, let $K$ be a $d$-simplicial element ($d = 2, 3$).  Let
$h_K := \mathrm{diam}(K)$ and $\rho_K$ be the diameter of its inscribed
sphere.  For the error analysis of Lagrange interpolation on simplicial
elements, many textbooks on finite element methods, such as those by
Ciarlet \cite{Ciarlet}, Brenner and Scott \cite{BrennerScott}, and
Ern and Guermond \cite{ErnGuermond}, explain the following theorem.

\begin{theorem}[Shape-regularity]\label{thm:S-regularity}
Let $\sigma > 0$ be a constant.
If $h_K/\rho_K \le \sigma$, then there
 exists a constant $C=C(\sigma)$ independent of $h_K$ such that
\begin{equation}
  |v - \mathcal{I}_K^1 v|_{1,2,K} \le C  h_K |v|_{2,2,K},
   \qquad  \forall v \in H^2(K).
  \label{eq:standard-est}
\end{equation}
\end{theorem}
Note that \eqref{eq:standard-est} does not hold if $d \ge 4$.

For the case of triangles, shape-regularity is equivalent to
the \textit{minimum angle condition} \cite{Zlamal}.  It is known
that shape-regularity is not an optimal condition on the geometry
of triangles.  If the maximum angle of a triangle is less than a fixed
constant $\theta_1 < \pi$, then the estimation \eqref{eq:standard-est}
holds with a constant $C = C(\theta_1)$.  This condition is known as
the \textit{maximum angle condition} \cite{BabuskaAziz, Jamet}. 

The present authors recently reported
an error estimation in terms of the circumradius of a triangle
\cite{KobayashiTsuchiya1, KobayashiTsuchiya3, KobayashiTsuchiya4}.
Let $R_K$ be the circumradius of a triangle $K$.
\begin{theorem}[Circumradius estimates]\label{thm:circumradius-est}
Let $K$ be an arbitrary triangle.  Let
$1 \le p \le \infty$, and $k$, $m$ be integers such that $k \ge 1$ and
$0 \le  m \le k$.  Then, for the $k$th-order Lagrange interpolation
 $\mathcal{I}_K^k$ on $K$, the estimation
\begin{equation} \label{eq:circumradius-est}
   |v - \mathcal{I}_K^k v|_{m,p,K} \le C
   \left(\frac{R_K}{h_K}\right)^m h_K^{k+1-m} |v|_{k+1,p,K}
  =  C R_K^m h_K^{k+1-2m} |v|_{k+1,p,K}
\end{equation}
holds for any $v \in W^{k+1,p}(K)$,
where the constant $C=C(k,m,p)$ is independent of the geometry of $K$.
\end{theorem}
Note that the circumradius estimation 
\eqref{eq:circumradius-est} is closely related to the definition of
surface area \cite{KobayashiTsuchiya2}.

The aim of this paper is to derive a similar error estimation
to Theorem~\ref{thm:circumradius-est} for tetrahedrons under no specific
geometric restrictions.

To extend the circumradius estimation \eqref{eq:circumradius-est}
to tetrahedrons, an immediate idea is to replace the circumradius
of a triangle with the radius of the circumsphere of a tetrahedron.
However, this idea can be immediately rejected by considering 
the tetrahedron $K$ with vertices
$\mathbf{x}_1 = (h,0,0)^\top$, $\mathbf{x}_2 = (-h,0,0)^\top$,
$\mathbf{x}_3 = (0,-h,h^\alpha)^\top$, and $\mathbf{x}_4 = (0,h,h^\alpha)^\top$
with $h > 0$ and $\alpha > 0$.
Setting $v_1(x,y,z) := x^2 - h^2 + h^{2-\alpha}z$,
we have that $\mathcal{I}_K^1 v_1 \equiv 0$, and a simple computation yields
$|v_1 - \mathcal{I}_K^1 v_1|_{1,\infty,K} = |v_1|_{1,\infty,K} \ge h^{2-\alpha}$
and $|v_1|_{2,\infty,K} = 2$.  Hence, if $\alpha > 2$,
an inequality such as \eqref{eq:circumradius-est} does
not hold for tetrahedrons, although the radius of the circumsphere of the above $K$
converges to $0$ as $h \to 0$.  Tetrahedrons such as the above $K$ are
called \textit{slivers} \cite{Sliver}.

Thus, we introduce the projected circumradius, also denoted by
$R_K$, of a tetrahedron $K$ in Section~\ref{section:Circumradius}.
We then obtain an error estimation for tetrahedrons that is
fundamentally similar to \eqref{eq:circumradius-est}
in Theorem~\ref{thm:main}.

In Section~\ref{Sect:Prei}, we state some preliminary results used in this paper.
In particular, we define the standard position and the projected
circumradius for tetrahedrons.  In Section~\ref{Sect:Squeezing}, we
recall definitions related to the quotient differences of functions
with multi-variables. We then reconfirm the Squeezing Theorem for
tetrahedrons.  In Section~\ref{sect:singular}, we obtain the error
estimation of Lagrange interpolation in terms of the singular values of
a linear transformation.  In Section~\ref{sec:geometric}, we present
a geometric interpretation of the singular values of the linear
transformation, and finally obtain the main theorem.


\section{Preliminaries}\label{Sect:Prei}
\subsection{Notation}
Let $d \ge 1$ be a positive integer and $\mathbb R^d$ be the $d$-dimensional
Euclidean space.  We denote the Euclidean norm of
$\mathbf{x} = (x_1,\cdots,x_d)^\top \in \mathbb R^d$ by $|\mathbf{x}|$. 
We always regard $\mathbf{x} \in \mathbb R^d$ as a column vector.
For a matrix $A$ and $\mathbf{x} \in \mathbb R^d$, $A^\top$ and $\mathbf{x}^\top$
denote their transpositions.

For $\delta = (\delta_1,...,\delta_d) \in \mathbb N_{0}^d$,
the multi-index $\partial^\delta$ of partial differentiation 
(in the sense of the distribution) is defined by
\[
    \partial^\delta = \partial_\mathbf{x}^\delta
    := \frac{\partial^{|\delta|}\ }
   {\partial x_1^{\delta_1}...\partial x_d^{\delta_d}}, \qquad
   |\delta| := \delta_1 + ... + \delta_d.
\]

Let $\Omega \subset \mathbb R^d$ be a (bounded) domain.  The usual Lebesgue
space is denoted by $L^p(\Omega)$ for $1 \le p \le \infty$.
For a positive integer $k$, the Sobolev space $W^{k,p}(\Omega)$ is
defined by
$\displaystyle
  W^{k,p}(\Omega) := 
  \left\{v \in L^p(\Omega) \, | \, \partial^\delta v \in L^p(\Omega), \,
   |\delta| \le k\right\}$.
For $1 \le p < \infty$, the norm and semi-norm of $W^{k,p}(\Omega)$ are defined
by
\begin{gather*}
  |v|_{k,p,\Omega} := 
  \biggl(\sum_{|\delta|=k} |\partial^\delta v|_{0,p,\Omega}^p
   \biggr)^{1/p}, \quad   \|v\|_{k,p,\Omega} := 
  \biggl(\sum_{0 \le m \le k} |v|_{m,p,\Omega}^p \biggr)^{1/p}
\end{gather*}
and $\displaystyle   |v|_{k,\infty,\Omega} := 
  \max_{|\delta|=k} \left\{\mathrm{ess}
   \sup_{\hspace{-5mm}\mathbf{x} \in\Omega}|\partial^\delta v(\mathbf{x})|\right\}$,
 $\displaystyle   \|v\|_{k,\infty,\Omega} := 
  \max_{0 \le m  \le k} \left\{|v|_{m,\infty,\Omega}\right\}$.


\subsection{The imbedding theorem}
Let $1 < p \le \infty$. From Sobolev's imbedding theorem and Morry's
inequality, we have the continuous imbeddings
\begin{gather*}
   W^{2,p}(K) \subset C^{1,1-3/p}(K), \quad p > 3, \\
   W^{2,3}(K) \subset W^{1,q}(K) \subset C^{0,1-3/q}(K),
      \quad \forall q > 3, \\
  W^{2,p}(K) \subset W^{1,3p/(3-p)}(K) \subset C^{0,2-3/p}(K),
      \quad \frac{3}{2} < p < 3, \\
  W^{3,3/2}(K) \subset W^{2,3}(K) \subset W^{1,q}(K)
  \subset C^{0,1-3/q}(K), \quad \forall q > 3, \\
  W^{3,p}(K) \subset W^{2,3p/(3-p)}(K) \subset W^{1,3p/(3-2p)}(K)
  \subset C^{0,3-3/p}(K), \quad 1 < p < \frac{3}{2}.
\end{gather*}
Although Morry's inequality may not be applied, the continuous
imbedding $W^{3,1}(K)$ $\subset C^{0}(K)$ still holds.
For the imbedding theorem, see \cite{AdamsFournier},
\cite{Brezis}, and \cite{KJF}.  In the following, we assume that $p$
is such that the imbedding $W^{k+1,p}(K) \subset C^{0}(K)$
holds, that is,
\begin{align*}
     1 \le p \le \infty, \quad  \text{ if } k+1 \ge 3
      \quad \text{ and } \quad
    \frac{3}{2} < p \le \infty, \quad \text{ if } k+1 =2.
\end{align*}

\subsection{The reference tetrahedrons}
Let $\widehat{K}$ and $\widetilde{K}$ be tetrahedrons that have the 
following vertices (see Figure~\ref{ref_tetra}):
\begin{align*}
  \widehat{K} \text{ has the vertices } \; (0,0,0)^\top, \; (1,0,0)^\top, \;
  (0,1,0)^\top, \; (0,0,1)^\top, \\
  \widetilde{K} \text{ has the vertices } \; (0,0,0)^\top, \; (1,0,0)^\top, \;
  (1,1,0)^\top, \; (0,0,1)^\top.
\end{align*}
\begin{figure}[tbhp]
\centering
\includegraphics[width=9cm]{reference-tetra-crop.eps}
\caption{The reference tetrahedrons $\widehat{K}$ and $\widetilde{K}$.}
\label{ref_tetra}
\end{figure}
These tetrahedrons are called the \textbf{reference tetrahedrons}.
In this paper, we denote the reference tetrahedrons by 
$\mathbf{K}$, that is, $\mathbf{K}$ is either of $\{\widehat{K},\widetilde{K}\}$.


\subsection{Standard position of tetrahedrons}
When considering tetrahedrons, it is convenient to define their ``standard
coordinates.''  Take any tetrahedron $K$ with vertices
$\mathbf{x}_i$, $i = 1, \dots, 4$.  The facet $B$ with vertices $\mathbf{x}_1$,
$\mathbf{x}_2$, $\mathbf{x}_3$ is regarded as the base of $K$. Let $\alpha$ and
$\beta$, $0 < \beta \le \alpha$, be the longest and shortest
lengths of the edges of $B$.   Without loss of generality, we assume
that $\mathbf{x}_1\mathbf{x}_2$ is the longest edge of $B$;
$|\mathbf{x}_1 - \mathbf{x}_2| = \alpha$.     Consider cutting $\mathbb R^3$
with the plane that contains the midpoint of the edge $\mathbf{x}_1\mathbf{x}_2$ and
is perpendicular to the vector $\mathbf{x}_1 - \mathbf{x}_2$.  Then, there exist
two cases: (i) $\mathbf{x}_3$ and $\mathbf{x}_4$ belong to the same half-space,
or (ii) $\mathbf{x}_3$ and $\mathbf{x}_4$ belong to different half-spaces.
Let $\gamma := |\mathbf{x}_1 - \mathbf{x}_4|$.
Under appropriate rotation, translation, and reflection operations,
these situations can be written using the parameters
\begin{equation}
 \begin{cases}
  0 < \beta \le \alpha, \quad 0 < \gamma, \quad s_1^2 + t_1^2  = 1, \;
 s_1 > 0, \; t_1 > 0, \quad \beta s_1 \le \frac{\alpha}{2}, \\
 s_{21}^2 + s_{22}^2 + t_2^2 = 1, \; t_2 > 0, \quad 
  \gamma s_{21} \le \frac{\alpha}{2},
 \end{cases}
 \label{eq:tetra-param}
\end{equation}
as
\begin{subequations} \label{eq:position}
\begin{gather}
 \mathbf{x}_1 = (0,0,0)^\top, \; \mathbf{x}_2 = (\alpha,0,0)^\top, \;
  \mathbf{x}_4 = (\gamma s_{21}, \gamma s_{22}, \gamma t_2)^\top,\\
  \begin{cases}
   \mathbf{x}_3 = (\beta s_1, \beta t_1,0)^\top & \text{for the case (i)} \\
   \mathbf{x}_3 = (\alpha - \beta s_1, \beta t_1,0)^\top & \text{for the case (ii)}
  \end{cases}.
\end{gather}
\end{subequations}
We refer to the coordinates in \eqref{eq:position}
as the \textit{standard position} of a tetrahedron.  In the following,
we sometimes write $h_B := \alpha$.  Let $R_B$ be the circumradius
of $B$.


\subsection{The projected circurmradius of tetrahedrons}
\label{section:Circumradius}
Suppose that a tetrahedron $K$ is at the standard position.
Let $\theta \in \mathbb R$ be such that
$-\frac{\pi}{2} \le \theta \le \frac{\pi}{2}$.
Let $\delta_\theta$ be the linear transformation defined by
the matrix
\[
   \begin{pmatrix}
      \cos\theta & -\sin\theta & 0 \\
          0 & 0 & 0 \\
          0 & 0 & 1
   \end{pmatrix}
   = \begin{pmatrix}
      1 & 0 & 0 \\
      0 & 0 & 0 \\
      0 & 0 & 1
   \end{pmatrix}
   \begin{pmatrix}
      \cos\theta & -\sin\theta & 0 \\
      \sin\theta & \cos\theta & 0 \\
          0 & 0 & 1
   \end{pmatrix}.
\]
That is, $\delta_\theta$ is a composite transformation of rotation
about the $z$-axis with angle $\theta$ and projection to the $xz$-plane.
The image $\delta_\theta(K)$ is a triangle on the $xz$-plane.  Let
$R_\theta$ be the circumradius of $\delta_\theta(K)$.  Define
\[
   R_P := \max_{\theta \in [-\pi/2,\pi/2]} R_\theta.
\]
The \textit{projected circumradius} $R_K$ of a tetrahedron $K$ is
defined by
\begin{equation}\label{def:RK}
   R_K := \min_{B} \frac{R_BR_P}{h_B},
\end{equation}
where the minimum is taken over all the facets of $K$.

\subsection{The setting of error estimation}
We define the set $\mathcal{T}_p^{k}(K) \subset W^{k+1,p}(K)$ by
\begin{gather*}
   \mathcal{T}_p^{k}(K) := \left\{v \in W^{k+1,p}(K) \Bigm|
     v(\mathbf{x}) = 0, \; \forall \mathbf{x} \in \Sigma^k(K) \right\}.
\end{gather*}
For Lagrange interpolation $\mathcal{I}_K^k(v)$, it is clear from the definition that
\[
     v - \mathcal{I}_K^k v \in \mathcal{T}_p^{k}(K), \quad \forall 
     v \in W^{k+1,p}(K).
\]
For an integer $m$ such that $0 \le m \le k$, $B_p^{m,k}(K)$
is defined by
\begin{gather*}
    B_p^{m,k}(K) := \sup_{v \in \mathcal{T}_p^{k}(K)}
    \frac{|v|_{m,p,K}}{|v|_{k+1,p,K}}.
\end{gather*}
Note that we have
\begin{gather*}
 B_p^{m,k}(K) = \inf\left\{C ; 
   |v - \mathcal{I}_K^k v|_{m,p,K} \le C |v|_{k+1,p,K}, \
   \forall v \in W^{k+1,p}(K)\right\},
\end{gather*}
that is, $B_p^{m,k}(K)$ is the \textit{best} constant $C$ for the
error estimation
\[
   |v - \mathcal{I}_K^k v|_{m,p,K} \le C |v|_{k+1,p,K}, \qquad 
    \forall v \in W^{k+1,p}(K).
\]
Therefore, we try to obtain an upper bound of $B_p^{m,k}(K)$ in terms of
geometric quantities of $K$.


\section{The Squeezing Theorem}\label{Sect:Squeezing}
Let $a$, $b \in\mathbb R$ be such that $0 < a \le 1$ and $0 < b$.
We then define the \textit{squeezing map} $sq_1^{ab}:R^3 \to \mathbb R^3$ by
\[
  sq_1^{ab}(x,y,z) := (x, a y, b z)^\top,
   \quad (x,y,z)^\top \in \mathbb R^3.
\]
Let $K_{ab} := sq_1^{ab}(\mathbf{K})$.  In this section, we will prove the following theorem.

\begin{theorem}\label{squeezingtheorem}
Let $k$ be a positive integer and $1 \le p \le \infty$.  There exists a
constant $C_{k,m,p}$ independent of $a$ and $b$ such that, for $m = 0, \cdots, k$,
\begin{gather*}
   B_p^{m,k}(K_{ab}) := \hspace{-0.2cm}
   \sup_{v \in \mathcal{T}_p^{k}(K_{ab})}
  \frac{|v|_{m,p,K_{ab}}}{|v|_{k+1,p,K_{ab}}}
   \le \max\{1,b^{k+1-m}\}C_{k,m,p},
\end{gather*}
where $p$ is taken such that
\begin{gather}
  \begin{cases}
      2 < p \le \infty & \text{ if } k - m = 0, \\ 
      \frac{3}{2} < p \le \infty & \text{ if } k = 1, \; m = 0,\\
      1 \le p \le \infty & \text{ if } k \ge 2 \text{ and } \; k-m \ge 1.
   \end{cases}
   \label{p-cond}
\end{gather}
\end{theorem}
Although the proof of Theorem~\ref{squeezingtheorem} is very similar to that of
\cite[Theorem~1.3]{KobayashiTsuchiya4}, we provide a sketch here for
readers' convenience.  Note that the restriction 
$2 < p \le \infty$ for the case $k=m$ comes from the continuity of
the trace operator $t:W^{1,p}(\mathbf{K}) \ni v \mapsto v|_S \in L^1(S)$, where
$S \subset \mathbf{K}$ is a non-degenerate segment (see
\cite[Section~3]{KobayashiTsuchiya4}).
Using the counter-example given by Shenk \cite{Shenk}, we find that
this restriction cannot be improved.

To prove Theorem~\ref{squeezingtheorem}, we recall the definitions of
difference quotients of multi-variable functions, the rectangular
parallelepiped $\square_\gamma^\delta$ defined by the lattice points
$\mathbf{x}_\gamma$, $\Delta^\delta\mathbf{x}_\gamma$ in $\mathbf{K}$, and the integral
$\int_{\square_\gamma^\delta} v$ introduced in
\cite[Section~2]{KobayashiTsuchiya4}.  The set
$\Xi_p^{\delta,k} \subset W^{k+1-|\delta|,p}(\mathbf{K})$ is then defined by
\begin{align*}
  \Xi_p^{\delta,k} & := \left\{ v \in W^{k+1-|\delta|,p}(\mathbf{K}) \Bigm| 
    \int_{\square_{\gamma}^{\delta}} v = 0,\quad \forall
   \square_{\gamma}^{\delta} \subset \mathbf{K} \right\}.
\end{align*}
Note that $u \in \mathcal{T}_p^{k}(\mathbf{K})$ implies
$\partial^\delta u \in \Xi_p^{\delta,k}$.

\begin{lemma}\label{lem32} 
We have $\Xi_p^{\delta,k} \cap \mathcal{P}_{k-|\delta|} = \{0\}$.
That is, if $q \in \mathcal{P}_{k-|\delta|}$ belongs to $\Xi_p^{\delta,k}$,
then $q=0$.
\end{lemma}

\textit{Proof.}
 Note that
$\mathrm{dim} \mathcal{P}_{k-|\delta|} = \#\{\square_{lp}^\delta \subset \mathbf{K} \}$.
For example, if $k=4$ and $|\delta|=3$, then
$\mathrm{dim}\mathcal{P}_1 = 4$.  This corresponds to the fact that,
in $\mathbf{K}$, there are four cubes of size $1/4$  for $\delta=(1,1,1)$ and
there are four rectangles of size $1/2 \times 1/4$ for $\delta=(1,2,0)$.
All their vertices (corners) belong to $\Sigma^4(\mathbf{K})$
(see Figure~2). Now, suppose that $v \in \mathcal{P}_{k-|\delta|}$ satisfies
$\int_{\square_{lp}^{\delta}} q = 0$ for all
$\square_{lp}^{\delta} \subset \mathbf{K}$.  These conditions are linearly
independent and determine $q = 0$ uniquely. \hfill $\square$

\begin{figure}[htbp]\label{cubes_tetra}
\centering
\subfloat{\includegraphics[width=4.5cm]{kyoko.eps}} \qquad
\subfloat{\includegraphics[width=3cm]{takako.eps}} \qquad
\subfloat{\includegraphics[width=3cm]{michiko.eps}}
\caption{The four cubes and four rectangles in $\mathbf{K}$.}
\end{figure}

\vspace{0.3cm}
The constant $A_p^{\delta,k}$ is defined by
\begin{align*}
   A_p^{\delta,k} := \sup_{v \in \Xi_p^{\delta,k}} \frac{|v|_{0,p,\mathbf{K}}}
  {|v|_{k+1-|\delta|,p,\mathbf{K}}}.
\end{align*}
The following lemma is an extension of \cite[Lemma~2.1]{BabuskaAziz}.

\begin{lemma}\label{lem33}
Let $p$ be such that $2 < p \le \infty$ if $k+1-|\delta|=1$ or
$1 \le p \le \infty$ if $k+1-|\delta|\ge 2$.  We then have
$A_p^{\delta,k} < \infty$.
\end{lemma}

\textit{Proof.}
See the proof of \cite[Lemma~3.3]{KobayashiTsuchiya4}.
\hfill $\square$

\vspace{0.3cm}
\textit{Proof of} Theorem~\ref{squeezingtheorem}.
First, let $1 \le p < \infty$ and $1 \le m \le k$.
For a multi-index $\gamma = (n_1, n_2, n_3) \in \mathbb N_0^3$ and a real $t \neq 0$,
set $(a, b)^{\gamma t} := a^{n_2 t}b^{n_3 t}$.  Take an arbitrary 
$v \in \mathcal{T}_p^k(K_{ab})$ and pull it back to
$u := v\circ sq_1^{ab} \in \mathcal{T}_p^k(\mathbf{K})$.   We then have
{\allowdisplaybreaks
\begin{align*}
  \frac{|v|_{m,p,K_{ab}}^p}{|v|_{k+1,p,K_{ab}}^p}
   & = \frac{\sum_{|\gamma| = m} \frac{m!}{\gamma!}
      (a,b)^{-\gamma p}\left|\partial^{\gamma}u\right|_{0,p,\mathbf{K}}^p}
      {\sum_{|\delta|= k+1} \frac{(k+1)!}{\delta!}
     (a,b)^{- \delta p} \left|\partial^{\delta} u \right|_{0,p,\mathbf{K}}^p} \\
  & = \frac{\sum_{|\gamma| = m}\frac{m!}{\gamma!}(a,b)^{-\gamma p}
         \left|\partial^{\gamma}u \right|_{0,p,\mathbf{K}}^p }
       {\sum_{|\gamma| = m} \frac{m!}{\gamma!}(a,b)^{-\gamma p}
       \left(\sum_{|\eta| = k+1-m} \frac{(k+1-m)!}
             {\eta!(a,b)^{\eta p} }\left|\partial^{\eta}
                 (\partial^{\gamma}u)\right|_{0,p,\mathbf{K}}^p\right)} \\
  & \le \frac{\max\{1,b^{(k+1-m)p}\}
   \sum_{|\gamma| = m}\frac{m!}{\gamma!}(a,b)^{-\gamma p}
         \left|\partial^{\gamma}u \right|_{0,p,\mathbf{K}}^p }
       {\sum_{|\gamma| = m} \frac{m!}{\gamma!}(a,b)^{-\gamma p}
       \left(\sum_{|\eta| = k+1-m} \frac{(k+1-m)!}
             {\eta! }\left|\partial^{\eta}
                 (\partial^{\gamma}u)\right|_{0,p,\mathbf{K}}^p\right)} \\
   & = \max\{1,b^{(k+1-m)p}\}
\frac{\sum_{|\gamma| = m} \frac{m!}{\gamma!}(a,b)^{-\gamma p}
         \left|\partial^{\gamma}u \right|_{0,p,\mathbf{K}}^p}
       {\sum_{|\gamma| = m} \frac{m!}{\gamma!}
      (a,b)^{-\gamma p}
        \left|\partial^{\gamma}u \right|_{k+1-m,p,\mathbf{K}}^p} \\
    & \le \max\{1,b^{(k+1-m)p}\}\max_{|\gamma|=m} 
    \left(A_p^{\gamma,k}\right)^p
   = \left(\max\{1,b^{k+1-m}\}C_{k,m,p}\right)^p,
\end{align*}
}
where $C_{k,m,p} := \max_{|\gamma|=m} A_p^{\gamma,k}$.
Proofs for the cases $1 \le p < \infty$ with $m=0$ and
$p = \infty$ with $0 \le m \le k$ are very similar.  See
the proof of \cite[Theorem~1.3]{KobayashiTsuchiya4}.  
\hfill $\square$


\vspace{0.3cm}
We now generalize the squeezing map.
Let $\alpha$, $\beta$, and $\gamma \in\mathbb R$ be such that
$0 < \beta \le \alpha$ and $0 < \gamma$.
We then define the \textit{squeezing map}
$sq_2^{\alpha\beta\gamma}:\mathbb R^3 \to \mathbb R^3$ by
\[
  sq_2^{\alpha\beta\gamma}(x,y,z) := (\alpha x, \beta y, \gamma z)^\top,
   \quad (x,y,z)^\top \in \mathbb R^3.
\]
Let $K_{\alpha\beta\gamma} := sq_2^{\alpha\beta\gamma}(\mathbf{K})$.  Note that
$G_{\alpha} := sq_2^{\alpha\alpha\alpha}$ is a similar transformation.
Let $\Omega \subset \mathbb R^3$ be a domain and an arbitrary function 
$v \in W^{k,p}(\Omega)$ be pulled-back to
$u := v \circ G_\alpha \in W^{k,p}(\Omega_\alpha)$ with
$\Omega_\alpha := G_{1/\alpha}(\Omega)$.
It is straightforward to check that 
\[
  |v|_{k,p,\Omega} = \alpha^{3/p-k}|u|_{k,p,\Omega_\alpha}.
\]

Now, take an arbitrary function $v \in \mathcal{T}_p^k(K_{\alpha\beta\gamma})$ and define
$u := v\circ sq_2^{\alpha\beta\gamma}$.  Then,
$u = u_1 \circ sq_1^{\frac{\beta}{\alpha},\frac{\gamma}{\alpha}}$ with
$u_1 := v \circ G_\alpha$, because
$sq_2^{\alpha\beta\gamma} = G_\alpha\circ sq_1^{\frac{\beta}{\alpha},\frac{\gamma}{\alpha}}$.
Therefore, it follows from Theorem~\ref{squeezingtheorem} that
\begin{align*}
   \frac{|v|_{m,p,K_{\alpha\beta\gamma}}}{|v|_{k+1,p,K_{\alpha\beta\gamma}}} = \alpha^{k+1-m}
    \frac{|u_1|_{m,p,K_{\frac{\beta}{\alpha},\frac{\gamma}{\alpha}}}}
   {|u_1|_{k+1,p,K_{\frac{\beta}{\alpha},\frac{\gamma}{\alpha}}}}
  & \le \alpha^{k+1-m} \max \left\{1,
   \left(\frac{\gamma}{\alpha}\right)^{k+1-m}\right\} C_{k,m,p} \\
  & = \left(\max\{\alpha,\gamma\}\right)^{k+1-m}C_{k,m,p}.
\end{align*}
Hence, we have derived the following corollary.

\begin{corollary}\label{cor:squeezing}
 Let $\alpha$, $\beta$, and $\gamma \in\mathbb R$ be such that
$0 < \beta \le \alpha$ and $0 < \gamma$.  
Let $K_{\alpha\beta\gamma} := sq_2^{\alpha\beta\gamma}(\mathbf{K})$. 
Assume that $k \ge 1$, $0 \le m \le k$, and $p$ is taken as
\eqref{p-cond}.   We then have
\begin{align*}
  B_p^{m,k}(K_{\alpha\beta\gamma}) := \hspace{-0.2cm} \sup_{v \in \mathcal{T}_p^{k}(K_{\alpha\beta\gamma})}
  \frac{|v|_{m,p,K_{\alpha\beta\gamma}}}{|v|_{k+1,p,K_{\alpha\beta\gamma}}}
  \le \left(\max\{\alpha,\gamma\}\right)^{k+1-m}C_{k,m,p}. \label{extension-3d}
\end{align*}
\end{corollary}


\section{Error estimates of Lagrange interpolation on general
 tetrahedrons}
\label{sect:singular}
In this section, we obtain an error estimation for Lagrange interpolation
on general tetrahedrons.  To this end, we apply the method developed in
\cite{KobayashiTsuchiya3}.

Recall that an arbitrary tetrahedron $K$ is written as
\eqref{eq:position} with parameters \eqref{eq:tetra-param}.  First, we
confirm that $K$ is obtained from the reference tetrahedron $\mathbf{K}$ by an
affine linear transformation.  Define the matrices $\widehat{A}$,
$\widetilde{A}$, $G \in GL(3,\mathbb R)$ by
\begin{align*}
   \widehat{A} := \begin{pmatrix}
        1  & s_{1} & s_{21} \\
        0  & t_{1} & s_{22} \\
        0  & 0     & t_2    
       \end{pmatrix}, \quad
   \widetilde{A} := \begin{pmatrix}
        1  & - s_{1} & s_{21} \\
        0  & t_{1} & s_{22} \\
        0  & 0     & t_2    
       \end{pmatrix}, \quad
  G := \begin{pmatrix}
        \alpha  & 0 & 0 \\
        0  & \beta & 0 \\
        0  & 0     & \gamma    
       \end{pmatrix}.
\end{align*}
We then have $K = \widehat{A}G(\widehat{K})$ for case (i) or
$K = \widetilde{A}G(\widetilde{K})$ for case (ii), that is,
$K = \widehat{A}(K_{\alpha\beta\gamma})$ or $K = \widetilde{A}(K_{\alpha\beta\gamma})$.
Note that $\widehat{A}$ and $\widetilde{A}$ can be decomposed as
$\widehat{A} = X\widehat{Y}$ and $\widetilde{A} = X\widetilde{Y}$
with
\begin{align*}
  X := \begin{pmatrix}
        1 & 0 & s_{21} \\
	0 & 1 & s_{22} \\
        0 & 0 & t_2
       \end{pmatrix}, \qquad
  \widehat{Y} := \begin{pmatrix}
        1 & s_1 & 0 \\
	0 & t_1 & 0 \\
        0 & 0 & 1
       \end{pmatrix}, \qquad
  \widetilde{Y} := \begin{pmatrix}
        1 & -s_1 & 0 \\
	0 & t_1 & 0 \\
        0 & 0 & 1
       \end{pmatrix},
\end{align*}
respectively.  We consider the singular values of $\widehat{A}$,
$\widetilde{A}$, $X$, $\widehat{Y}$, and $\widetilde{Y}$.
A straightforward computation yields 
\begin{gather*}
  \det\left(X^\top X - \mu I\right) 
    = (1 - \mu)\left(\mu^2 - 2\mu + t_2^2\right), \\
  \det\left(\widehat{Y}^\top \widehat{Y} - \mu I\right) 
  = \det\left(\widetilde{Y}^\top \widetilde{Y} - \mu I\right) 
    = (1 - \mu)\left(\mu^2 - 2\mu + t_1^2\right),
\end{gather*}
the eigenvalues of which are
$\mu = 1$, $1 \pm \sqrt{1 - t_i^2} = 1 \pm \mathbf{s}_i$, $i = 1, 2$, where
$\mathbf{s}_1 := |s_1|$ and $\mathbf{s}_2 := \sqrt{s_{21}^2 + s_{22}^2}$.
Therefore, for $\mathbf{a} \in \mathbb R^3$, we have
\begin{gather*}
  (1 - \mathbf{s}_2)|\mathbf{a}|^2 \le |X\mathbf{a}|^2 \le (1 + \mathbf{s}_2)^2|\mathbf{a}|^2, \\
  (1 - \mathbf{s}_1)|\mathbf{a}|^2 \le |Z\mathbf{a}|^2 \le (1 + \mathbf{s}_1)|\mathbf{a}|^2,
  \qquad Z = \widehat{Y} \text{ or } Z = \widetilde{Y},  \\
 \prod_{i=1}^2(1 - \mathbf{s}_i)|\mathbf{a}|^2 \le |V\mathbf{a}|^2 \le 
 \prod_{i=1}^2(1 + \mathbf{s}_i)|\mathbf{a}|^2, \qquad
 V = \widehat{A} \text{ or } V = \widetilde{A}.
\end{gather*}

Let $K := V(K_{\alpha\beta\gamma})$, where $V = \widehat{A}$ or $V = \widetilde{A}$.
A function $v\in W^{r,p}(K)$ is pulled-back
to a function $u \in W^{r,p}(K_{\alpha\beta\gamma})$ by $u(\mathbf{x}) := v(V\mathbf{x})$.
Using inequality (2.1) in \cite{KobayashiTsuchiya3}, we have
\begin{gather*}
  \frac{\prod_{i=1}^2(1-\mathbf{s}_i)^{r}}{t_1^{2r}t_2^{2r}}
  \sum_{|\gamma|=r} (\partial_\mathbf{x}^\gamma u)^2
 \le \sum_{|\gamma|=r} (\partial_\mathbf{y}^\gamma v)^2 \le
\frac{\prod_{i=1}^2(1+\mathbf{s}_i)^{r}}{t_1^{2r}t_2^{2r}}
  \sum_{|\gamma|=r} (\partial_\mathbf{x}^\gamma u)^r,
\end{gather*}
where $\mathbf{y} := V\mathbf{x}$.  The fact that
$\mathrm{det}A = t_1t_2$ and inequalities (2.2), (2.3) in
\cite{KobayashiTsuchiya3} then give
\begin{align*}
   |v|_{m,p,K} & \le 3^{m\tau(p)} 
 \frac{\prod_{i=1}^2(1+\mathbf{s}_i)^{m/2}}
       {(t_1 t_2)^{m-1}}  |u|_{m,p,K_{\alpha\beta\gamma}}, \\
  |v|_{k+1,p,K} & \ge 3^{-(k+1)\tau(p)}   
\frac{\prod_{i=1}^2(1-\mathbf{s}_i)^{(k+1)/2}}{(t_1t_2)^{k}}
  |u|_{k+1,p,K_{\alpha\beta\gamma}},   \\
   \frac{|v|_{m,p,K}}{|v|_{k+1,p,K}} & \le 3^{(k+1+m)\tau(p)}
   \frac{(t_1t_2)^{k+1-m}\prod_{i=1}^2(1+\mathbf{s}_i)^{m/2}|u|_{m,p,K_{\alpha\beta\gamma}}}
     {\prod_{i=1}^2(1-\mathbf{s}_i)^{(k+1)/2}|u|_{k+1,p,K_{\alpha\beta\gamma}}} \\
 & = 3^{(k+1+m)\tau(p)} \frac{\prod_{i=1}^2(1+\mathbf{s}_i)^{(k+1)/2}}
  {\prod_{i=1}^2(1-\mathbf{s}_i)^{m/2}}
  \frac{|u|_{m,p,K_{\alpha\beta\gamma}}}{|u|_{k+1,p,K_{\alpha\beta\gamma}}} \\
  & \le C  \frac{\left(\max\{\alpha,\gamma\}\right)^{k+1-m}}
  {\prod_{i=1}^2(1-\mathbf{s}_i)^{m/2}},
\end{align*}
where
\begin{align*}
  \tau(p) := \begin{cases}
    \frac{1}{p} - \frac{1}{2}, & 1 \le p \le 2 \\
    \frac{1}{2} - \frac{1}{p} , & 2 \le p \le \infty
      \end{cases}, \qquad
  C := 3^{(k+1+m)\tau(p)} 2^{k+1} C_{k,m,p}.
\end{align*}
(See the inequalities in \cite[p.496]{KobayashiTsuchiya3}.)
Note that the constant $C$ depends only on $k$, $m$, and $p$.
Hence, we have derived the following theorem:
\begin{theorem}\label{thmB}
Let $K$ be an arbitrary tetrahedron with vertices given by
\eqref{eq:position} with the parameters in \eqref{eq:tetra-param}.
Let $k$, $m$ be integers with $k \ge 1$ and $0 \le m \le k$.
Let $p$ be taken as in Theorem~\ref{squeezingtheorem} according to $k$ and $m$.
Then, we have
\[
 B_p^{m,k}(K) := \sup_{v \in \mathcal{T}_p^k}\frac{|v|_{m,p,K}}{|v|_{k+1,p,K}}
  \le C  \frac{\left(\max\{\alpha,\gamma\}\right)^{k+1-m}}
  {\prod_{i=1}^2(1-\mathbf{s}_i)^{m/2}},
\]
where $C = C(k,m,p)$ is a constant independent of $K$.
\end{theorem}


\section{A geometric interpretation and the main theorem}
\label{sec:geometric}
In this section, we consider the geometric meaning of the quantity
$\prod_{i=1}^2(1-\mathbf{s}_i)^{-1/2}$ that appeared in Theorem~\ref{thmB}.
Recall that $K$ is a tetrahedron with vertices given by
\eqref{eq:position} and the parameters in \eqref{eq:tetra-param}.
Recall also that $B$ is the base of $K$ with vertices
$\mathbf{x}_1$, $\mathbf{x}_2$, $\mathbf{x}_3$.  The circumradius $R_B$
can then be written as
\[
   R_B = \frac{\sqrt{\alpha^2 - 2\alpha\beta s_1 + \beta^2}}{2t_1}.
\]
Because
\begin{align*}
   \alpha^2 - 2\alpha\beta s_1 + \beta^2 & =
  \frac{\alpha^2(2 - s_1^2)}{4} + \alpha
   \left(\frac{\alpha}{2} - \beta s_1\right) 
   + \left(\frac{\alpha s_1}{2} - \beta\right)^2 \\
   & \ge \frac{\alpha^2(2 - s_1^2)}{4} \ge \frac{\alpha^2}{4},
\end{align*}
we have
\begin{align}
   R_B \ge \frac{\alpha}{4 t_1} = \frac{\alpha}{4 \sqrt{1 - \mathbf{s}_1^2}} 
    \ge \frac{h_B}{4\sqrt{2} \sqrt{1 - \mathbf{s}_1}},
   \label{RB-est}
\end{align}
where we have used the definition $h_B = \alpha$.

Recall that $R_P$ was defined in Section~\ref{section:Circumradius}.
We will show that there exists a constant $C$ independent of $K$
such that
\begin{align}
    R_P \ge C\frac{\max\{\alpha,\gamma\}}{\sqrt{1 - \mathbf{s}_2}}.
   \label{RP-est}
\end{align}

Let $\mathbf{x}_3 = (\eta,\xi,0)^\top$, that is, $\xi = \beta t_1$ and
$\eta = \beta s_1$ or $\eta = \alpha - \beta s_1$.  From the
assumption, we have $0 < \eta < \alpha$ and $\xi > 0$.  Note that 
\begin{align*}
   \delta_\theta(\mathbf{x}_1) = (0,0,0)^\top, \quad
  & \delta_\theta(\mathbf{x}_2) = (\alpha\cos\theta,0,0)^\top, \\
   \delta_\theta(\mathbf{x}_3) = (\eta\cos\theta-\xi\sin\theta,0,0)^\top, 
  \quad &
   \delta_\theta(\mathbf{x}_4) = \gamma
  (s_{21}\cos\theta-s_{22}\sin\theta,0,t_2)^\top.
\end{align*}

Let $\underbar{x} < \overline{\mathrm{x}}$ be the $x$-coordinates of the end points of the base of
$\delta_\theta(K)$.  The assumptions given in \eqref{eq:tetra-param}
yield
\begin{align*}
   \underbar{x} \le 0, \quad \overline{\mathrm{x}} = \alpha\cos\theta, \quad
   0 \le \theta \le \frac{\pi}{2}, \qquad
   \underbar{x} =0, \quad \overline{\mathrm{x}} \ge \alpha\cos\theta, \quad
   -\frac{\pi}{2} \le \theta \le 0. 
\end{align*}
Defining
$w = w(\theta) := s_{21}\cos\theta - s_{22}\sin\theta$, $R_\theta$ is
written as
\[
   R_\theta = \frac{1}{2\gamma t_2}
   \left((\overline{\mathrm{x}} - \gamma w)^2 + \gamma^2t_2^2\right)^{1/2}
   \left((\underbar{x} - \gamma w)^2 + \gamma^2t_2^2\right)^{1/2}.
\]
Take an arbitrary
$\theta \in \left[-\frac{\pi}{3},\frac{\pi}{3}\right]$.
Suppose that the inequality
\begin{align}
  \gamma w \le \frac{\underbar{x} + \overline{\mathrm{x}}}{2}
  \label{key-inequality1}
\end{align}
holds and there exists a constant $C_1$ independent of
$\theta$ such that
\begin{align}
 |\underbar{x} - \gamma w| \ge C_1 \gamma \mathbf{s}_2.  \label{key-inequality2}
\end{align}
We then have
\begin{align*}
   (\underbar{x} - \gamma w)^2 + \gamma^2 t_2^2 & \ge C_1^2
   (1 - t_2^2) + \gamma^2 t_2^2 \ge
   \min\{1, C_1^2\} \gamma^2, \\
   (\overline{\mathrm{x}} - \gamma w)^2 + \gamma^2 t_2^2 & \ge \max\left\{
      (\underbar{x} - \gamma w)^2 + \gamma^2 t_2^2, \;
      \frac{\alpha^2\cos^2\theta}{4} + \gamma^2 t_2^2 \right\} \\
  & \ge \max\left\{\min\{1,C_1^2\} \gamma^2,\, C_2^2\alpha^2\right\} \\
  & \ge \min\{C_1^2,C_2^2\}\max\{\alpha^2,\gamma^2\},
\end{align*}
where $C_2 = C_2(\theta) := \frac{\cos\theta}{2}$.
Here, we have used the fact that $\overline{\mathrm{x}}/2 \ge (\overline{\mathrm{x}}+\underbar{x})/2 \ge \gamma w$ and
$\overline{\mathrm{x}} - \gamma w \ge \overline{\mathrm{x}}/2 \ge (\alpha\cos\theta)/2$.  Hence,
setting
\[
   \frac{\min\{C_1,C_2\}\min\{1,C_1\}}{2\sqrt{2}} 
  \ge \frac{\min\left\{C_1,\frac{1}{4}\right\}
   \min\{1,C_1\}}{2\sqrt{2}} =: C_3,
\]
we obtain
\begin{align*}
  R_P \ge R_\theta \ge \frac{\sqrt{2}C_3}{t_2}\max\{\alpha,\gamma\}
   = \frac{\sqrt{2}C_3 \max\{\alpha,\gamma\}}{\sqrt{1-\mathbf{s}_2^2}}
   \ge C_3 \frac{\max\{\alpha,\gamma\}}{\sqrt{1-\mathbf{s}_2}},
\end{align*}
and the key inequality \eqref{RP-est} is shown.

Fix $\varphi$ such that
\[
   0 < \varphi < \frac{\pi}{6}, \qquad 
   \sin2\varphi\tan2\varphi \le \frac{1}{6}.
\]
In the following, we will show that, according to 
$\mathbf{x}_4 = (\gamma s_{21}, \gamma s_{22}, \gamma t_2)^\top$, we can
take an appropriate
$\theta \in \left[-\frac{\pi}{3},\frac{\pi}{3}\right]$ such that
conditions \eqref{key-inequality1} and \eqref{key-inequality2}
hold with $C_1 = \sin\varphi$.  

\vspace{0.3cm}
\noindent
\textbf{Case 1.} Suppose that $|s_{22}|\tan\varphi \le |s_{21}|$. \\
In this case, we set $\theta = 0$ and have $\underbar{x} = 0$, $\overline{\mathrm{x}} = \alpha$,
and $\gamma w = \gamma s_{21} \le \alpha/2 = (\underbar{x} + \overline{\mathrm{x}})/2$
because of \eqref{eq:tetra-param}.  Hence, \eqref{key-inequality1} holds.
For \eqref{key-inequality2}, we note that
\begin{align*}
   |\underbar{x} - \gamma w| = \gamma|s_{21}| & = \gamma
   \left(s_{21}^2\sin^2\varphi + s_{21}^2\cos^2\varphi\right)^{1/2} \\
  & \ge \gamma
   \left(s_{21}^2\sin^2\varphi + s_{22}^2\sin^2\varphi\right)^{1/2}
   = \gamma \mathbf{s}_2 \sin\varphi,
\end{align*}
and so \eqref{key-inequality2} holds with $C_1 := \sin\varphi$.

\vspace{0.3cm}
\noindent
\textbf{Case 2:} Suppose that $|s_{22}|\tan\varphi > |s_{21}|$
and $3\gamma s_{22} \tan2\varphi \le \alpha$. \\
In this case, set $\theta = -2\varphi$. We then have $\underbar{x} = 0$
and $\overline{\mathrm{x}} \ge \alpha\cos\theta = \alpha \cos2\varphi$.  If
$s_{22} > 0$, then
\begin{align*}
 \gamma w & = \gamma (s_{21}\cos2\varphi + s_{22}\sin2\varphi)
  \le \gamma (s_{22}\cos2\varphi\tan\varphi + s_{22}\sin2\varphi) \\
  & = \frac{\gamma s_{22}}{2}
   (3 \sin2\varphi - 2\sin^2\varphi\tan\varphi) \\
  & \le \frac{2\gamma s_{22}}{2}\sin2\varphi
   \le \frac{\alpha \cos2\varphi}{2} \le \frac{\underbar{x} + \overline{\mathrm{x}}}{2}.
\end{align*}
   If $s_{22} \le 0$, then
\begin{align*}
 \gamma w & = \gamma (s_{21}\cos2\varphi + s_{22}\sin2\varphi)
  \le \gamma (- s_{22}\cos2\varphi\tan\varphi + s_{22}\sin2\varphi) \\
  & = \gamma s_{22}\tan\varphi \le 0 \le \frac{\underbar{x} + \overline{\mathrm{x}}}{2}.
\end{align*}
Thus, in either case, \eqref{key-inequality1} holds.
For \eqref{key-inequality2}, we note that
\begin{align*}
   |\underbar{x} - \gamma w| & = \gamma
  |s_{21}\cos2\varphi + s_{22}\sin2\varphi| \ge \gamma \left(
  |s_{22}|\sin2\varphi - |s_{21}|\cos2\varphi\right) \\
  & \ge \gamma \left(
  |s_{22}|\sin2\varphi - |s_{22}|\cos2\varphi\tan\varphi\right)
  = \gamma|s_{22}|\tan\varphi \\
  & = \gamma \left(
    s_{22}^2 \sin^2\varphi\tan^2\varphi + s_{22}^2\sin^2\varphi
   \right)^{1/2} \\
  & = \gamma \left(s_{21}^2 \sin^2\varphi + s_{22}^2\sin^2\varphi
   \right)^{1/2} = \gamma\mathbf{s}_2\sin\varphi,
  \end{align*}
and so \eqref{key-inequality2} holds with $C_1 := \sin\varphi$.

\vspace{0.3cm}
\noindent
\textbf{Case 3:} Suppose that $|s_{22}|\tan\varphi > |s_{21}|$
and $3\gamma s_{22} \tan2\varphi > \alpha$. \\
In this case, set $\theta = 2\varphi$.  We then have
\[
   \underbar{x} = \min\{\eta\cos2\varphi - \xi\sin2\varphi,0\}, \qquad
    \overline{\mathrm{x}} = \alpha \cos2\varphi.
\]
If $\eta\cos2\varphi - \xi\sin2\varphi \le 0$, then
\begin{align*}
   \underbar{x} & = \eta\cos2\varphi - \xi\sin2\varphi
   = \alpha \cos2\varphi + (\eta - \alpha) \cos2\varphi
      - \xi\sin2\varphi \\
  & \ge \alpha \cos2\varphi -
   \left((\eta - \alpha)^2 + \xi^2\right)^{1/2} 
   \ge \alpha \cos2\varphi - \alpha = -2\alpha \sin^2\varphi,
\end{align*}
because the lengths of all edges of the base are less than $\alpha$.
Even if $\eta\cos2\varphi - \xi\sin2\varphi \ge 0 = \underbar{x}$, the above
inequality obviously holds.  Because 
\begin{align*}
  \gamma w & = \gamma (s_{21}\cos2\varphi - s_{22}\sin2\varphi) \\
  & \le \gamma (s_{22}\cos2\varphi\tan\varphi - s_{22}\sin2\varphi)
  = -\gamma s_{22}\tan\varphi,
\end{align*}
we have
\begin{align*}
  \underbar{x} - \gamma w & \ge -2\alpha\sin^2\varphi + \gamma s_{22}\tan\varphi \\
  & \ge -2\alpha\sin^2\varphi + \frac{\alpha\tan\varphi}{6\tan2\varphi}
 + \frac{\gamma}{2} s_{22}\tan\varphi \\
  & = \alpha\left(\frac{1}{6}- \sin2\varphi\tan2\varphi\right)
   \frac{\tan\varphi}{\tan2\varphi}
   + \frac{\gamma}{2} s_{22}\tan\varphi
    \ge \frac{\gamma}{2} s_{22}\tan\varphi > 0.
\end{align*}
Therefore,
\[
   \gamma w < \underbar{x} < \frac{\underbar{x} + \overline{\mathrm{x}}}{2},
\]
and \eqref{key-inequality1} holds.
For \eqref{key-inequality2}, we note that
\begin{align*}
   |\underbar{x} - \gamma w| & \ge \frac{\gamma}{2} s_{22}\tan\varphi \\
  & = \gamma \left(
    s_{22}^2 \sin^2\varphi\tan^2\varphi + s_{22}^2\sin^2\varphi
   \right)^{1/2} \\
  & \ge \gamma \left(s_{21}^2 \sin^2\varphi + s_{22}^2\sin^2\varphi
   \right)^{1/2} = \gamma\mathbf{s}_2\sin\varphi,
  \end{align*}
and so \eqref{key-inequality2} holds with $C_1 := \sin\varphi$.

\vspace{0.3cm}
Using inequalities \eqref{RB-est} and \eqref{RP-est}, we have shown the
following lemma.

\begin{lemma}\label{lemma:geometric}
Let $K$ be a tetrahedron with vertices given by \eqref{eq:position}
and the parameters in \eqref{eq:tetra-param}.
We then have
\[
  \prod_{i=1}^2(1-\mathbf{s}_i)^{-1/2} \le C
   \frac{R_BR_P}{h_B\max\{\alpha,\gamma\}},
\]
where $C$ is a constant independent of $K$.
\end{lemma}

Take a facet $B$ so that $R_BR_P/h_B = R_K$.
Combining Theorem~\ref{thmB} and Lemma~\ref{lemma:geometric} with the 
projected circumradius $R_K$, we have
\begin{align*}
 B_p^{m,k}(K) := \sup_{v \in \mathcal{T}_p^k}\frac{|v|_{m,p,K}}{|v|_{k+1,p,K}}
  & \le C  \frac{\left(\max\{\alpha,\gamma\}\right)^{k+1-m}}
  {\prod_{i=1}^2(1-\mathbf{s}_i)^{m/2}} \\
  & \le C \left(\frac{R_BR_P}{h_B}\right)^m
    \left(\max\{\alpha,\gamma\}\right)^{k+1-2m} \\
  & \le C R_K^m h_K^{k+1-2m},
\end{align*}
where the constant $C = C(k,m,p)$ is independent of $K$
with vertices \eqref{eq:position}.  Note that Sobolev (semi-)norms 
are affected by rotation up to a constant.  Therefore, we have
derived the following result, which is the main theorem of this
paper.

\begin{theorem}\label{thm:main}
Let $K$ be an arbitrary tetrahedron.  Let $h_K:=\mathrm{diam}K$ and
$R_K$ be the projected circumradius of $K$ defined by
\eqref{def:RK}.  Assume that $k$, $m$ are integers with $k \ge 1$,
$0 \le m \le k$, and $p$ $(1 \le p \le \infty)$ is taken as in \eqref{p-cond}.
Then, for arbitrary $v \in W^{k+1,p}(K)$,
there exists a constant $C = C(k,m,p)$ independent of $K$ such that
\begin{align}
  |v - \mathcal{I}_K^k v|_{m,p,K} \le C R_K^m h_K^{k+1-2m}|v|_{k+1,p,K}
   = C \left(\frac{R_K}{h_K}\right)^m h_K^{k+1-m}|v|_{k+1,p,K}.
  \label{main-est}
\end{align}
\end{theorem}
Note that the error estimation \eqref{main-est} is exactly same as
\eqref{eq:circumradius-est}.

\section{Concluding remarks}
In this final section, we present some remarks on the newly obtained
error estimation.

(1) For finite element error analysis, the most important
case is $p=2$ and $m=1$.  In this case, $k$ should be greater than or
equal to $2$
because of the restriction in \eqref{p-cond}.  In this case, the main
estimation is
\[
   |v - \mathcal{I}_K^k v|_{1,p,K} \le C R_K h_K^{k-1}|v|_{k+1,p,K}
   = C \left(\frac{R_K}{h_K}\right) h_K^{k}|v|_{k+1,p,K}.
\]

Let $\Omega \subset \mathbb R^3$ be a bounded polygonal domain.  
Suppose that we compute a numerical solution of the Poisson equation
\[
    -\Delta u = f \text{ in } \Omega, \quad
    u = 0 \text{ on } \partial\Omega
\]
by the finite element method with conforming simplicial elements.
One difficulty we may encounter in the numerical computation is the bad
triangulation of $\Omega$ with many slivers \cite{Sliver}.  Consider
the typical sliver $K$ mentioned in Section~\ref{intro}, whose vertices are
$(\pm h,0,0)^\top$ and $(0,\pm h,h^\alpha)^\top$.  We can see that
$R_K = \mathcal{O}(h_K^{2-\alpha})$ and
$R_Kh_K^{k-1} = \mathcal{O}(h_K^{k+1-\alpha})$.  Thus, if we use
a $k$th-order conforming Lagrange finite element method with
triangulation that contains many slivers like $K$, the theoretical
convergence rate may be $\mathcal{O}(h^{k+1-\alpha})$.  This is
worse than the expected rate $\mathcal{O}(h^{k})$, but we can still
expect convergence if $k+1-\alpha > 0$.   Therefore,
bad triangulation with many slivers can be remedied by using
higher-order Lagrange elements.

(2) Let $K$ be a tetrahedron with the inscribed sphere $S_K$.
If the regularity assumption is imposed on $K$, there exists
a constant $\sigma$ such that $h_K/\rho_K \le \sigma$, where
$\rho_K := \mathrm{diam}S_K$.
Suppose that $K$ is at the standard position.  Recall
$\delta_\theta$, the projection introduced in
Section~\ref{section:Circumradius}.  It is clear that
$\delta_\theta(S_K) \subset \delta_\theta(K)$ and 
$\mathrm{diam}\delta_\theta(S_K) = \rho_K$.  Therefore, recalling that $R_P$
is the maximum value of the circumradius of $\delta_\theta(K)$, there
exists a constant $C$ such that $R_P \le C h_K$.
The regularity assumption implies the maximum angle condition, and
as was pointed out in \cite[Section~4.1]{KobayashiTsuchiya3}, the
maximum angle condition of $B$ implies the boundedness of $R_B/h_B$.
Thus, Theorem~\ref{thm:main} is an extension of the standard estimation
\eqref{eq:standard-est}.

(3) For tetrahedrons, a (generalized) maximum angle condition was given
by K\v{r}\'{i}\v{z}ek \cite{Krizek2}, \cite{Duran}.
A slightly more general condition has been stated by Jamet \cite{Jamet}.
The authors conjecture that Theorem~\ref{thm:main} is an extension
of these results.

(4) In Theorem~\ref{thmB}, we obtained an error estimation for
Lagrange interpolation on tetrahedrons in terms of the singular values
${\prod_{i=1}^2(1-\mathbf{s}_i)^{-1/2}}$ of the linear transformation.
In Theorem~\ref{thm:main}, we showed that the projected circumradius
is a geometric interpretation of the singular values.  The authors,
however, are not completely sure whether the projected circumradius is
the \textit{best} interpretation. 

Further research on the geometry of tetrahedrons is required 
to ascertain the geometric properties of the singular values 
${\prod_{i=1}^2(1-\mathbf{s}_i)^{-1/2}}$ and the projected circumradius,
and the relationship between these singular values and prior results.



2016-06-14T02:18:15.000
1606
1606.03918
en
https://arxiv.org/abs/1606.03918
This paper describes the analysis of Lagrange interpolation errors on tetrahedrons. In many textbooks, the error analysis of Lagrange interpolation is conducted under geometric assumptions such as shape regularity or the (generalized) maximum angle condition. In this paper, we present a new estimation in which the error is bounded in terms of the diameter and projected circumradius of the tetrahedron. Because we do not impose any geometric restrictions on the tetrahedron itself, our error estimation may be applied to any tetrahedralizations of domains including very thin tetrahedrons.
Numerical Analysis (math.NA)
Error Analysis of Lagrange Interpolation on Tetrahedrons
Qwen/Qwen-72B
1. YES
2. YES