previous overview next
Via the Cubic
$
\def \hieruit {\quad \Longrightarrow \quad}
\def \slechts {\quad \Longleftrightarrow \quad}
\def \EN {\quad \mbox{and} \quad}
\def \OF {\quad \mbox{or} \quad}
$
The cubic equation for grid refinement with TripleGrid Calculus is:
$$
y_{k+1} (y_{k+1} - 3)^2 = y_k
$$ $$
y_{k+1}^3 - 6 y_{k+1}^2 + 9 y_{k+1} - y_k = 0
$$
Which is of the form:
$$
a_3 x^3 + a_2 x^2 + a_1 x + a_0 = 0
$$
The discriminant of such a cubic equation is far less well known than the one
of a quadratic equation. It may nevertheless be found on the Internet at
Wikipedia: Discriminant.
Where we can copy and paste the following formula:
$$
D = a_1^2 a_2^2 - 4a_0 a_2^3 - 4a_1^3 a_3 + 18a_0 a_1 a_2 a_3 - 27a_0^2 a_3^2
$$
The general theory about the discriminant of a polynomial is quite interesting
in itself. The abovementioned web page offers pointers to topics like i.e. the
resultant and via the resultant link the Sylvester matrix of two polynomials.
Whatever. Let's apply the general discriminant of the cubic to the specific
one we have at hand with TripleGrid:
$$
y_{k+1}^3 - 6 y_{k+1}^2 + 9 y_{k+1} - y_k = 0
$$
With $\left[\; a_3 = 1 \quad a_2 = - 6 \quad a_1 = 9 \quad a_0 = - y_k \;
\right]$ its discriminant $D$ becomes:
$$
D = 9^2 (-6)^2 - 4(-y_k)(-6)^3 - 4(9)^3 1 + 18(-y_k)9(-6)1 - 27(-y_k)^21^2
$$ $$
= 3^6 2^2 - 3^6 2^2 - 2^5 3^3 y_k + 3^5 2^2 y_k - 3^3 y_k^2
= 2^2 3^3 (-8+9) y_k - 3^3 y_k^2
$$ $$
\hieruit D = 27.y_k(4-y_k)
$$
It is known that the above cubic equation has three real solutions if and only
if this discriminant $D$ is positive. Meaning that the following is a necessary
and sufficient condition for it:
$$
0 \le y_k \le 4
$$
Thus, with TripleGrid coarsening, the Dangerous Interval becomes apparent
almost immediately, from the sign of the discriminant alone.
Key reference for the rest of this paragraph is:
The Geometry of the Cubic Formula.
Which is based on "A new approach to solving the cubic" by R.W.D. Nickalls
(not available anymore).
Recall the cubic equation which is associated with triple grid refinement:
$$
y_{k+1}^3 - 6 y_{k+1}^2 + 9 y_{k+1} - y_k = 0
$$
It is of the form:
$$
a x^3 + b x^2 + c x + d = 0
$$
With $\left[\quad a = 1 \quad b = -6 \quad c = 9 \quad d = -y_k \quad\right]$ .
According to the above reference, we can define a quantity $x_N$ by:
$$
x_N = \frac{-b}{3a} = \frac{6}{3} = 2
$$
This means that a substitution $z = y_{k+1}-2$ or $y_{k+1} = z+2$ will bring
our cubic into its reduced form:
$$
y_{k+1}^3 - 6 y_{k+1}^2 + 9 y_{k+1} - y_k =
(z+2)^3 - 6(z+2)^2 + 9(z+2) - y_k =
$$ $$
z^3 + 6 z^2 + 12 z + 8 - 6 z^2 - 24 z - 24 + 9 z + 18 - y_k =
z^3 - 3 z + 2 - y_k = 0
$$
Furthermore, a quantity $\delta$ is defined by:
$$
\delta^2 = \frac{b^2 - 3 a c}{9 a^2} = \frac{36 - 27}{9} = 1
$$
Let's consider the case where $0 \le y_k \le 4$. Then there are three real
solutions and we have the Casus Irreducibilis for the cubic equation
at hand. According to the above reference, we must use a trigonometric
substitution in order to be able to find any solutions:
$$
z = 2\delta \cos(\theta) = 2 \cos(\theta)
$$
Upon substitution it reads:
$$
8 \cos^3(\theta) - 6 \cos(\theta) + 2 - y_k = 0
$$
In the paragraph TripleGrid Product we found the following (well, eh)
formula:
$$
2 + 2\cos(3\theta) = 8\cos^3(\theta) - 6\cos(\theta) + 2
$$
Therefore we have to solve:
$$
2 + 2\cos(3\theta) = y_k \slechts y_{k+1} (y_{k+1} - 3)^2 = y_k
$$
Now define:
$$
y_k = 2 + 2\cos(\phi_k) \EN y_{k+1} = 2 + 2\cos(\phi_{k+1})
$$
Then we have that $\theta = \phi_{k+1}$ and:
$$
y_{k+1} (y_{k+1} - 3)^2 = y_k \slechts
2 + 2\cos(3\phi_{k+1}) = 2 + 2\cos(\phi_k) \slechts
$$ $$
\cos(\phi_k) = \cos(3\phi_{k+1})
$$
Compare this with:
$$
\cos(\phi_k) = \cos(2\phi_{k+1})
$$
And it seems that we have landed on solid ground: the method of angles.