Table of Contents Author Guidelines Submit a Manuscript
Abstract and Applied Analysis
VolumeΒ 2011, Article IDΒ 203453, 22 pages
http://dx.doi.org/10.1155/2011/203453
Research Article

Subdivision Depth Computation for Tensor Product 𝑛 -Ary Volumetric Models

Department of Mathematics, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan

Received 22 October 2010; Revised 27 January 2011; Accepted 28 February 2011

Academic Editor: YoshikazuΒ Giga

Copyright Β© 2011 Ghulam Mustafa and Muhammad Sadiq Hashmi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We offer computational formula of subdivision depth for tensor product 𝑛 -ary ( 𝑛 β©Ύ 2 ) volumetric models based on error bound evaluation technique. This formula provides and error control tool in subdivision schemes over regular hexahedron lattice in higher-dimensional spaces. Moreover, the error bounds of Mustafa et al. (2006) are special cases of our bounds.

1. Introduction

Subdivision is a simple and elegant method to describe smooth curves and surfaces. The approach of subdivision schemes is simple and efficient due to its mathematical formulation. Its application ranges from industrial design and animation to scientific visualization and simulation. Due to this, subdivision method is becoming a standard technique and now well-understood by both academic and industrial communities. It is an algorithm to generate smooth curves and surfaces as a sequence of successively refined control polygons. At each subdivision level, the subdivision scheme describe the source grid maps to the subdivided grid, which results in increase in the number of points. The number of points inserted at level π‘˜ + 1 between two consecutive points from level π‘˜ is called arity of the scheme. In the case when number of points inserted are 2 , 3 , … , 𝑛 the subdivision schemes are called binary, ternary, …, 𝑛 -ary, respectively. For more details on 𝑛 -ary subdivision schemes, we may refer to [14] thesis of Aspert [5] and Ko [6].

However tensor product trivariate schemes obtained from above schemes have been proven themselves an excellent tool for the modeling of largely regular volumetric/solid models over hexahedron lattice, for example, manufacturing of industrial regular block, font animation and garment pressures for the biomechanical design of functional apparel products, and so forth.

Although subdivision has many attractive advantages for modeling purposes, it has received far less attention in volumetric modeling. One of the reason is the topological complexity of domain meshes in higher-dimensional spaces. Another reason is lack of proper mathematical tools for extraordinary analysis. Even so, it is clear that subdivision is slowly gaining interest in solid modeling community. MacCracken and Joy [7] proposed a tensor product extension of the Catmull-Clark scheme in the volumetric setting over hexahedron lattice. Later on, Bajaj et al. [8] further extended the scheme with the analysis based on numerical experiments. In 2004, McDonnell et al. [9] present a volumetric subdivision scheme for interpolation of arbitrary hexahedral meshes. Mustafa and Liu [10] present subdivision scheme which exhibits control over shrink-age/size of volumetric models in 2005. The method here presented is much simpler and easier as compared to MacCracken and Joy [7]. Cheng and Yong [11] proposed subdivision scheme based on nonhexahedron lattice. Wang et al. [12] gave the application of the volumetric subdivision scheme in the simulation of elastic human body deformation and garment pressure.

Subdivision schemes have become important in recent years because they provide a precise and efficient way to describe smooth curves/surfaces/volumetric models, however the little have been done in the area of error control for tensor product 𝑛 -ary volumetric models. The investigation of error control for volumetric models arises two questions in mind.(i)How well the regular hexahedron lattice approximate to the limit volumetric model?(ii)How many subdivision steps are needed to satisfy a user-specified error tolerance?

For given error tolerance, the subdivision levels performed on the initial control polygon, so that the error/distance between the resulting control polygon and the limit volumetric models would be less than the error tolerance is called subdivision depth.

A subdivision depth and error bound based on forward differences of control points have been presented by [11, 1318], while the methods [1922] are based on eigenanalysis. But nothing in this area has been done for more general tensor product 𝑛 -ary volumetric models yet. In this paper, we will answer-above-said questions and present a subdivision depth computation technique based on error bounds for tensor product 𝑛 -ary volumetric models.

It is notified that the increase in arity offers greater freedom than offered by low arity subdivision volumetric scheme in term of coefficients. Higher arity volumetric schemes allow a range of different behaviors than the lower arity volumetric schemes. Ko [6] notified that subdivision curves/surfaces with higher arity results in higher smoothness and approximation order but smaller in support, which make it more practical in use. It is also noticed that higher arity volumetric models have slightly lower computational cost than lower arity volumetric models. This discussion motivate us to calculate error bound and depth for higher arity subdivision volumetric models, that is, in general for tensor product 𝑛 -ary subdivision volumetric models. Our method is generalization of Mustafa et al. [13, 1517]. The paper is arranged as follows.

Section 2 is devoted for basic definitions and notations. In Section 3, we have computed subdivision depth for tensor product 𝑛 -ary volumetric models. Section 4 presents applications of our results for tensor product 𝑛 -ary volumetric models. Conclusion and future research directions are given in Section 5. The typical mathematical proofs are placed in Appendices A and B for improved presentation of the paper.

2. Preliminaries

In this section, first we list all the basic facts about subdivision curve, surface and volumetric models needed in this paper. Then we settle some notations for fair reading and better understanding of Section 3.

2.1. Concepts
2.1.1. 𝑛 -Ary Subdivision Curve

Given a sequence of control points 𝑝 π‘˜ 𝑖 ∈ ℝ 𝑁 , 𝑖 ∈ β„€ , 𝑁 β©Ύ 2 , where the upper index π‘˜ β©Ύ 0 indicates the subdivision level, an 𝑛 -ary subdivision curve [5] is defined by 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 = π‘š  𝑗 = 0 π‘Ž 𝛼 , 𝑗 𝑝 π‘˜ 𝑖 + 𝑗 , 𝛼 = 0 , 1 , … , 𝑛 βˆ’ 1 , ( 2 . 1 ) where π‘š > 0 and π‘š  𝑗 = 0 π‘Ž 𝛼 , 𝑗 = 1 , 𝛼 = 0 , 1 , … , 𝑛 βˆ’ 1 . ( 2 . 2 ) The set of coefficients { π‘Ž 𝛼 , 𝑗 , 𝛼 = 0 , 1 , … , 𝑛 βˆ’ 1 } π‘š 𝑗 = 0 is called subdivision mask. Given initial values 𝑝 0 𝑖 ∈ ℝ 𝑁 , 𝑖 ∈ β„€ , then in the limit π‘˜ β†’ ∞ , the process (2.1) defines an infinite set of points in ℝ 𝑁 . The sequence of control points { 𝑝 π‘˜ 𝑖 } is related, in a natural way, with the dyadic mesh points 𝑑 π‘˜ 𝑖 = 𝑖 / 𝑛 π‘˜ , 𝑖 ∈ β„€ . The process then defines a scheme whereby 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 replaces the value 𝑝 π‘˜ 𝑖 + 𝛼 / 𝑛 for 𝛼 ∈ { 0 , 𝑛 } . Here 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 is inserted at the mesh point 𝑑 π‘˜ + 1 𝑛 𝑖 + 𝛼 = ( 1 / 𝑛 ) ( ( 𝑛 βˆ’ 𝛼 ) 𝑑 π‘˜ 𝑖 + 𝛼 𝑑 π‘˜ 𝑖 + 1 ) for 𝛼 = 0 , 1 , … , 𝑛 . Labelling of old and new points is shown in Figure 1 which illustrates subdivision scheme (2.1).

fig1
Figure 1: Solid lines show coarse polygons whereas doted lines are refined polygons. (a)–(c) represent binary, ternary, and quaternary refinement of coarse polygon of scheme (2.1) for 𝑛 = 2 , 3 , 4 , respectively.
2.1.2. Tensor Product 𝑛 -Ary Subdivision Surface

Given a sequence of control points 𝑝 π‘˜ 𝑖 , 𝑗 ∈ ℝ 𝑁 , 𝑖 , 𝑗 ∈ β„€ , 𝑁 β©Ύ 2 , where the upper index π‘˜ β©Ύ 0 indicates the subdivision level, a tensor product 𝑛 -ary surface is a tensor product of (2.1) defined by 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 = π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛽 , 𝑠 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝛼 , 𝛽 = 0 , 1 , … , 𝑛 βˆ’ 1 , ( 2 . 3 ) where π‘Ž 𝛼 , π‘Ÿ satisfies (2.2). Given initial values 𝑝 0 𝑖 , 𝑗 ∈ ℝ 𝑁 , 𝑖 , 𝑗 ∈ β„€ , then in the limit π‘˜ β†’ ∞ , the process (2.3) defines an infinite set of points in ℝ β„• . The sequence of values { 𝑝 π‘˜ 𝑖 , 𝑗 } is related, in a natural way, with the dyadic mesh points ( 𝑖 / 𝑛 π‘˜ , 𝑗 / 𝑛 π‘˜ ) , 𝑖 , 𝑗 ∈ β„€ . The process then defines a scheme whereby 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 replaces the values 𝑝 π‘˜ 𝑖 + 𝛼 / 𝑛 , 𝑗 + 𝛽 / 𝑛 for 𝛼 , 𝛽 ∈ { 0 , 𝑛 } . Here the values 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 are inserted at the mesh points ( ( 𝑛 𝑖 + 𝛼 ) / 𝑛 π‘˜ + 1 , ( 𝑛 𝑗 + 𝛽 ) / 𝑛 π‘˜ + 1 ) for 𝛼 , 𝛽 = 0 , 1 , … , 𝑛 . Labelling of old and new points is shown in Figure 2 which illustrates subdivision scheme (2.3).

fig2
Figure 2: Solid lines show one face of coarse polygons whereas doted lines are refined polygons. (a)–(c) can be obtain by subdividing one face into four, nine and sixteen new faces by using (2.3) for 𝑛 = 2 , 3 , 4 (i.e., binary, ternary and quaternary), respectively.
2.1.3. Tensor Product 𝑛 -Ary Volumetric Model

Given a sequence of control points 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 ∈ ℝ 𝑁 , 𝑖 , 𝑗 , 𝑙 ∈ β„€ , 𝑁 β©Ύ 2 , where the upper index π‘˜ β©Ύ 0 indicates the subdivision level, a tensor product 𝑛 -ary volumetric model is tensor product of (2.3) defined by 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  π‘š 𝑠 = 0  𝑑 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛽 , 𝑠 π‘Ž 𝛾 , 𝑑 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 , 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 , ( 2 . 4 ) where π‘Ž 𝛼 , π‘Ÿ satisfies (2.2). Given initial values 𝑝 0 𝑖 , 𝑗 , 𝑙 ∈ ℝ 𝑁 , 𝑖 , 𝑗 , 𝑙 ∈ β„€ , then in the limit π‘˜ β†’ ∞ , the process (2.4) defines an infinite set of points in ℝ β„• . The sequence of values { 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 } is related, in a natural way, with the dyadic mesh points ( 𝑖 / 𝑛 π‘˜ , 𝑗 / 𝑛 π‘˜ , 𝑙 / 𝑛 π‘˜ ) , 𝑖 , 𝑗 , 𝑙 ∈ β„€ . The process then defines a scheme whereby 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 replaces the values 𝑝 π‘˜ 𝑖 + 𝛼 / 𝑛 , 𝑗 + 𝛽 / 𝑛 , 𝑙 + 𝛾 / 𝑛 for 𝛼 , 𝛽 , 𝛾 ∈ { 0 , 𝑛 } . Here the values 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 are inserted at the mesh points ( ( 𝑛 𝑖 + 𝛼 ) / 𝑛 π‘˜ + 1 , ( 𝑛 𝑗 + 𝛽 ) / 𝑛 π‘˜ + 1 , ( 𝑛 𝑙 + 𝛾 ) / 𝑛 π‘˜ + 1 ) for 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 . Labelling of old and new points is shown in Figure 3 which illustrates subdivision scheme (2.4).

fig3
Figure 3: Solid lines show one cube of coarse polygons whereas doted lines are refined polygons. (a)–(c) can be obtain by subdividing one cube into 23, 33, and 43 new cubes by using (2.4) for 𝑛 = 2 , 3 , 4 , respectively.
2.1.4. Subdivision Depth

Given control polygon of tensor product 𝑛 -ary volumetric model and an error tolerance πœ– , if we subdivide control polygon π‘˜ times so that the error between resulting polygon and volumetric model is smaller than πœ– , then π‘˜ is called subdivision depth of tensor product 𝑛 -ary volumetric model with respect to πœ– .

2.2. Notations

Here, we assume 𝛿 = m a x 𝛼 , 𝛽 , 𝛾 ⎧ βŽͺ ⎨ βŽͺ ⎩ | | | | | π‘š  π‘š π‘Ÿ = 0  π‘š 𝑠 = 0  𝑑 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛽 , 𝑠 𝑏 𝛾 , 𝑑 | | | | | ⎫ βŽͺ ⎬ βŽͺ ⎭ , , 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 ( 2 . 5 ) where 𝑏 𝛾 , 𝑗 = 𝑗  𝑑 = 0 ξ€· π‘Ž 𝛾 , 𝑑 βˆ’ π‘Ž 𝛾 + 1 , 𝑑 ξ€Έ 𝑏 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 2 , 𝑛 βˆ’ 1 , 𝑗 = π‘Ž 0 , 𝑗 βˆ’ 𝑛 βˆ’ 2  𝛾 = 0 𝑏 𝛾 , 𝑗 . ( 2 . 6 ) Suppose further for 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 , 𝑀 π‘˜ 𝛼 , 𝛽 , 𝛾 = β€– β€– 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 βˆ’ 1 𝑛 3  ( 𝑛 βˆ’ 𝛼 ) ( 𝑛 βˆ’ 𝛽 ) ( 𝑛 βˆ’ 𝛾 ) 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 𝛼 ( 𝑛 βˆ’ 𝛽 ) ( 𝑛 βˆ’ 𝛾 ) 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 + 𝛽 ( 𝑛 βˆ’ 𝛼 ) ( 𝑛 βˆ’ 𝛾 ) 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 𝛾 ( 𝑛 βˆ’ 𝛼 ) ( 𝑛 βˆ’ 𝛽 ) 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 + 𝛼 𝛽 ( 𝑛 βˆ’ 𝛾 ) 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 + 𝛼 𝛾 ( 𝑛 βˆ’ 𝛽 ) 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 + 1 + 𝛽 𝛾 ( 𝑛 βˆ’ 𝛼 ) 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1 + 𝛼 𝛽 𝛾 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 + 1  β€– β€– . ( 2 . 7 ) Assume Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 1 = 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 5 = 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 , Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 2 = 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 6 = 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 , Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 3 = 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 7 = 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1 , Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 4 = 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 , ( 2 . 8 ) πœ’ = m a x 𝑑 ξ‚» m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 𝑑 β€– β€– ξ‚Ό , , 𝑑 = 0 , 1 , … , 7 ( 2 . 9 ) πœ‚ 1 𝛼 , 𝛽 , 𝛾 = | | | | | π‘Ž 𝛽 , 0 π‘Ž π‘š 𝛾 , 0  𝑑 = 1 π‘Ž 𝛼 , 𝑑 βˆ’ 𝛼 ( 𝑛 βˆ’ 𝛽 ) ( 𝑛 βˆ’ 𝛾 ) 𝑛 3 | | | | | + | | | | | π‘Ž 𝛽 , 0 π‘Ž 𝛾 , 0 π‘š βˆ’ 1  𝑠 = 1 Μƒ π‘Ž 𝛼 , 𝑠 | | | | | , ( 2 . 1 0 ) πœ‚ 2 𝛼 , 𝛽 , 𝛾 = | | | | | π‘Ž π‘š 𝛾 , 0  𝑑 = 1 π‘Ž 𝛽 , 𝑑 βˆ’ 𝛽 ( 𝑛 βˆ’ 𝛾 ) 𝑛 2 | | | | | + | | | | | π‘Ž π‘š 𝛾 , 0  π‘Ÿ = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 𝛼 , π‘Ÿ Μƒ π‘Ž 𝛽 , 𝑠 | | | | | , ( 2 . 1 1 ) πœ‚ 3 𝛼 , 𝛽 , 𝛾 = | | | | | π‘š  𝑑 = 1 π‘Ž 𝛾 , 𝑑 βˆ’ 𝛾 𝑛 | | | | | + | | | | | π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛽 , 𝑑 Μƒ π‘Ž 𝛾 , 𝑠 | | | | | , ( 2 . 1 2 ) πœ‚ 4 𝛼 , 𝛽 , 𝛾 = | | | | | π‘Ž π‘š 𝛾 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 𝛼 , 𝑠 π‘Ž 𝛽 , 𝑑 βˆ’ 𝛼 𝛽 ( 𝑛 βˆ’ 𝛾 ) 𝑛 3 | | | | | + | | | | | π‘Ž π‘š 𝛾 , 0  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 𝛽 , 𝑑 Μƒ π‘Ž 𝛼 , 𝑠 | | | | | , ( 2 . 1 3 ) πœ‚ 5 𝛼 , 𝛽 , 𝛾 = | | | | | π‘Ž π‘š 𝛽 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 𝛼 , 𝑠 π‘Ž 𝛾 , 𝑑 βˆ’ 𝛼 𝛾 ( 𝑛 βˆ’ 𝛽 ) 𝑛 3 | | | | | + | | | | | π‘Ž π‘š 𝛽 , 0  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 𝛾 , 𝑑 Μƒ π‘Ž 𝛼 , 𝑠 | | | | | , ( 2 . 1 4 ) πœ‚ 6 𝛼 , 𝛽 , 𝛾 = | | | | | π‘š  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 𝛽 , 𝑠 π‘Ž 𝛾 , 𝑑 βˆ’ 𝛽 𝛾 𝑛 2 | | | | | + | | | | | π‘š  π‘š 𝑑 = 1  π‘Ÿ = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 𝛾 , 𝑑 π‘Ž 𝛼 , π‘Ÿ Μƒ π‘Ž 𝛽 , 𝑠 | | | | | , ( 2 . 1 5 ) πœ‚ 7 𝛼 , 𝛽 , 𝛾 = | | | | | π‘š  π‘š π‘Ÿ = 1  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛽 , 𝑠 π‘Ž 𝛾 , 𝑑 βˆ’ 𝛼 𝛽 𝛾 𝑛 3 | | | | | + | | | | | π‘š  π‘š π‘Ÿ = 1  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 𝛽 , π‘Ÿ π‘Ž 𝛾 , 𝑑 Μƒ π‘Ž 𝛼 , 𝑠 | | | | | , ( 2 . 1 6 ) where Μƒ π‘Ž 𝛼 , 0 = π‘š  𝑑 = 1 π‘Ž 𝛼 , 𝑑 βˆ’ 𝛼 𝑛 , 𝛼 = 0 , 1 , … , 𝑛 βˆ’ 1 , Μƒ π‘Ž 𝛼 , 𝑗 = π‘š  𝑑 = 𝑗 + 1 π‘Ž 𝛼 , 𝑑 , 𝑗 β©Ύ 1 , 𝛼 = 0 , 1 , … , 𝑛 βˆ’ 1 . ( 2 . 1 7 ) Furthermore, suppose πœ— = m a x 𝛼 , 𝛽 , 𝛾 ξƒ― ( πœ’ ) 7  𝑑 = 1 ξ‚€ πœ‚ 𝑑 𝛼 , 𝛽 , 𝛾  ξƒ° . , 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 ( 2 . 1 8 )

3. Depth for Tensor Product 𝑛 -Ary Volumetric Models

In this paragraph, we compute subdivision depth for tensor product 𝑛 -ary volumetric model. Moreover, we show that error bound for binary subdivision volumetric models [17] is special case of our bound. Here we need following lemmas for Theorem 3.5. The proof of first two lemmas are shown in Appendices A and B, respectively.

Lemma 3.1. Given initial control polygon 𝑝 0 𝑖 , 𝑗 , 𝑙 = 𝑝 𝑖 , 𝑗 , 𝑙 , 𝑖 , 𝑗 , 𝑙 ∈ β„€ , let the values 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , π‘˜ β©Ύ 1 be defined recursively by (2.4) together with (2.2), then m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 𝑑 β€– β€– β©½ ( 𝛿 ) π‘˜ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 𝑑 β€– β€– , ( 3 . 1 ) where 𝛿 and Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 𝑑 , 𝑑 = 1 , 2 , … , 7 , are defined by (2.5) and (2.8), respectively.

Lemma 3.2. Given initial control polygon 𝑝 0 𝑖 , 𝑗 , 𝑙 = 𝑝 𝑖 , 𝑗 , 𝑙 , 𝑖 , 𝑗 , 𝑙 ∈ β„€ , let the values 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , π‘˜ β©Ύ 1 be defined recursively by (2.4) together with (2.2), then 𝑀 π‘˜ 𝛼 , 𝛽 , 𝛾 β©½ ( πœ’ ) ( 𝛿 ) π‘˜ 7  𝑑 = 1 ξ‚€ πœ‚ 𝑑 𝛼 , 𝛽 , 𝛾  , ( 3 . 2 ) where 𝛿 , 𝑀 π‘˜ 𝛼 , 𝛽 , 𝛾 , πœ’ , πœ‚ 𝑑 𝛼 , 𝛽 , 𝛾 for 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 are defined by (2.5), (2.7), (2.9)–(2.16), respectively.

Lemma 3.3. Given initial control polygon 𝑝 0 𝑖 , 𝑗 , 𝑙 = 𝑝 𝑖 , 𝑗 , 𝑙 , 𝑖 , 𝑗 , 𝑙 ∈ β„€ , let the values 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , π‘˜ β©Ύ 1 be defined recursively by (2.4) together with (2.2). Suppose 𝑃 π‘˜ is the piecewise linear interpolant to the values 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 and 𝑃 ∞ is the limit volumetric model of (2.4). If 𝛿 < 1 , then error bound between tensor product 𝑛 -ary volumetric model and its control polygon after π‘˜ -fold subdivision is β€– β€– 𝑃 π‘˜ βˆ’ 𝑃 ∞ β€– β€– ∞ ξ‚΅ β©½ πœ— ( 𝛿 ) π‘˜ ξ‚Ά , 1 βˆ’ 𝛿 ( 3 . 3 ) where 𝛿 and πœ— are defined by (2.5) and (2.18), respectively.

Proof. Let β€– β‹… β€– ∞ denote the uniform norm. Since the maximum difference between 𝑃 π‘˜ + 1 and 𝑃 π‘˜ is attained at a point on the ( π‘˜ + 1 ) th mesh, we have β€– β€– 𝑃 π‘˜ + 1 βˆ’ 𝑃 π‘˜ β€– β€– ∞ β©½ m a x 𝛼 , 𝛽 , 𝛾  𝑀 π‘˜ 𝛼 , 𝛽 , 𝛾  , , 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 ( 3 . 4 ) where 𝑀 π‘˜ 𝛼 , 𝛽 , 𝛾 is defined by (2.7). By (3.2) and (3.4), we get β€– β€– 𝑃 π‘˜ + 1 βˆ’ 𝑃 π‘˜ β€– β€– ∞ β©½ πœ— ( 𝛿 ) π‘˜ , ( 3 . 5 ) where 𝛿 and πœ— are defined by (2.5) and (2.18), respectively.
By triangle inequality we get (3.3). This completes the proof.

Remark 3.4. Theorem 3.10 in [17] is designed to estimate error bound for binary subdivision volumetric model (i.e., each cube is divided in 8 subcubes). But for the higher arity subdivision schemes such as for 𝑛 = 3 , 4 , 5 , … (when each cube is divided in 3 3 , 4 3 , 5 3 , … subcubes) error estimates are not feasible by existing result. So estimation of error bounds for tensor product 𝑛 -ary volumetric model is quite necessary. Our Lemma 3.3 provides freedom to evaluate error bound for all arities. Here we also mention that Lemma 3.3 for 𝑛 = 2 reduces to [17, Theorem 3.10].

Now we offer the computational formula of subdivision depth for tensor product 𝑛 -ary volumetric model.

Theorem 3.5. Let π‘˜ be subdivision depth, and let 𝑑 π‘˜ be the error bound between tensor product 𝑛 -ary volumetric model 𝑃 ∞ and its π‘˜ -level control polygon 𝑃 π‘˜ . For arbitrary πœ– > 0 , if π‘˜ β©Ύ l o g 𝛿 βˆ’ 1 ξ‚΅ πœ— ξ‚Ά , πœ– ( 1 βˆ’ 𝛿 ) ( 3 . 6 ) then 𝑑 π‘˜ β©½ πœ– . ( 3 . 7 )

Proof. From (3.3), we have 𝑑 π‘˜ = β€– β€– 𝑃 π‘˜ βˆ’ 𝑃 ∞ β€– β€– ∞ ξ‚΅ β©½ πœ— ( 𝛿 ) π‘˜ ξ‚Ά . 1 βˆ’ 𝛿 ( 3 . 8 ) This implies, for arbitrary given πœ– > 0 , when subdivision depth π‘˜ satisfies the following inequality: π‘˜ β©Ύ l o g 𝛿 βˆ’ 1 ξ‚΅ πœ— ξ‚Ά , πœ– ( 1 βˆ’ 𝛿 ) ( 3 . 9 ) then 𝑑 π‘˜ β©½ πœ– . ( 3 . 1 0 ) This completes the proof.

4. Applications

4.1. Error Bound and Subdivision Depth of Tensor Product 𝑛 -Ary Interpolating Volumetric Models

In this section, we estimate the error bound and subdivision depth of ( 2 𝑏 + 2 ) -point tensor product 𝑛 -ary interpolating volumetric models. By taking the tensor product of ( 2 𝑏 + 2 ) -point 𝑛 -ary scheme of [4], we get the following: 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 , 𝑛 𝑙 = 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝑠 1 , 𝑛 𝑗 , 𝑛 𝑙 = 𝑏 + 1  𝑑 1 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 1 , 𝑑 1 ξ€Έ 𝑝 π‘˜ 𝑑 1 + 𝑖 , 𝑗 , 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 + 𝑠 2 , 𝑛 𝑙 = 𝑏 + 1  𝑑 2 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 2 , 𝑑 2 ξ€Έ 𝑝 π‘˜ 𝑖 , 𝑑 2 + 𝑗 , 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 , 𝑛 𝑙 + 𝑠 3 = 𝑏 + 1  𝑑 3 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 3 , 𝑑 3 ξ€Έ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑑 3 + 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝑠 1 , 𝑛 𝑗 + 𝑠 2 , 𝑛 𝑙 = 𝑏 + 1  𝑑 1 = βˆ’ 𝑏 𝑏 + 1  𝑑 2 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 1 , 𝑑 1 𝔸 ξ€Έ ξ€· 𝑠 2 , 𝑑 2 ξ€Έ 𝑝 π‘˜ 𝑑 1 + 𝑖 , 𝑑 2 + 𝑗 , 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝑠 1 , 𝑛 𝑗 , 𝑛 𝑙 + 𝑠 3 = 𝑏 + 1  𝑑 1 = βˆ’ 𝑏 𝑏 + 1  𝑑 3 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 1 , 𝑑 1 𝔸 ξ€Έ ξ€· 𝑠 3 , 𝑑 3 ξ€Έ 𝑝 π‘˜ 𝑑 1 + 𝑖 , 𝑗 , 𝑑 3 + 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 + 𝑠 2 , 𝑛 𝑙 + 𝑠 3 = 𝑏 + 1  𝑑 2 = βˆ’ 𝑏 𝑏 + 1  𝑑 3 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 2 , 𝑑 2 𝔸 ξ€Έ ξ€· 𝑠 3 , 𝑑 3 ξ€Έ 𝑝 π‘˜ 𝑖 , 𝑑 2 + 𝑗 , 𝑑 3 + 𝑙 , 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝑠 1 , 𝑛 𝑗 + 𝑠 2 , 𝑛 𝑙 + 𝑠 3 = 𝑏 + 1  𝑑 1 = βˆ’ 𝑏 𝑏 + 1  𝑑 2 = βˆ’ 𝑏 𝑏 + 1  𝑑 3 = βˆ’ 𝑏 ξ€· 𝔸 𝑠 1 , 𝑑 1 𝔸 ξ€Έ ξ€· 𝑠 2 , 𝑑 2 𝔸 ξ€Έ ξ€· 𝑠 3 , 𝑑 3 ξ€Έ 𝑝 π‘˜ 𝑑 1 + 𝑖 , 𝑑 2 + 𝑗 , 𝑑 3 + 𝑙 , ( 4 . 1 ) where 𝔸 π‘₯ , 𝑦 = ∏ 𝑏 + 1 π‘š = βˆ’ 𝑏 ( π‘₯ βˆ’ 𝑛 π‘š ) ( π‘₯ βˆ’ 𝑛 𝑦 ) ( βˆ’ 1 ) ( 𝑏 βˆ’ 1 βˆ’ 𝑦 ) ( 𝑛 ) ( 2 𝑏 + 1 ) , ( 𝑏 + 𝑦 ) ! ( 𝑏 βˆ’ 𝑦 + 1 ) ! ( 4 . 2 ) 𝑠 1 , 𝑠 2 , 𝑠 3 = 1 , 2 , … , 𝑛 βˆ’ 1 , 𝑏 = 1 , 2 , 3 , … and 𝑛 stands for 𝑛 -ary interpolating subdivision volumetric model, that is, 𝑛 = 2 , 3 , 4 , … stands for binary, ternary, quaternary and so on, respectively.

The error bounds and subdivision depth of (4.1) are shown in Tables 1 and 2, respectively. In these tables, we have shown the error bounds and depth of different arity interpolating subdivision volumetric models by using Lemma 3.3 and Theorem 3.5 with πœ’ = 0 . 1 .

tab1
Table 1: Error bound of tensor product 𝑛 -ary interpolating volumetric models: here 𝑛 presents the arity of subdivision volumetric model and π‘˜ presents the subdivision level.
tab2
Table 2: Subdivision depth of tensor product 𝑛 -ary interpolating volumetric models: here 𝑛 presents the arity of subdivision volumetric model and πœ– presents error tolerance.
4.2. Error Bound and Subdivision Depth of Tensor Product 𝑛 -Ary Approximating Volumetric Models

In this section, we estimate the error bound and subdivision depth of tensor product ( 2 𝑏 + 2 ) -point 𝑛 -ary approximating volumetric models. By taking the tensor product of ( 2 𝑏 + 2 ) -point 𝑛 -ary scheme of [4], we get the following: 𝑝 π‘˜ + 1 𝑛 𝑖 + 𝑠 1 , 𝑛 𝑗 + 𝑠 2 , 𝑛 𝑙 + 𝑠 3 = 𝑏 + 1  𝑑 1 = βˆ’ 𝑏 𝑏 + 1  𝑑 2 = βˆ’ 𝑏 𝑏 + 1  𝑑 3 = βˆ’ 𝑏 ξ€· 𝔹 𝑠 1 , 𝑑 1 𝔹 ξ€Έ ξ€· 𝑠 2 , 𝑑 2 𝔹 ξ€Έ ξ€· 𝑠 3 , 𝑑 3 ξ€Έ 𝑝 π‘˜ 𝑑 1 + 𝑖 , 𝑑 2 + 𝑗 , 𝑑 3 + 𝑙 , ( 4 . 3 ) where 𝔹 π‘₯ , 𝑦 =  ∏ 𝑏 + 1 π‘š = βˆ’ 𝑏 ( 2 π‘₯ + 1 βˆ’ 2 𝑛 π‘š ) ( 2 π‘₯ + 1 βˆ’ 2 𝑛 𝑦 ) ( βˆ’ 1 ) ( 𝑏 βˆ’ 1 βˆ’ 𝑦 ) ( 2 𝑛 ) ( 2 𝑏 + 1 ) ξƒͺ , ( 𝑏 + 𝑦 ) ! ( 𝑏 βˆ’ 𝑦 + 1 ) ! ( 4 . 4 ) 𝑠 1 , 𝑠 2 , 𝑠 3 = 1 , 2 , … , 𝑛 βˆ’ 1 , 𝑏 = 1 , 2 , 3 , … , and 𝑛 stands for 𝑛 -ary approximating subdivision volumetric model, that is, 𝑛 = 2 , 3 , 4 , … stands for binary, ternary, quaternary, and so on, respectively.

The error bounds and subdivision depth of (4.3) are shown in Tables 3 and 4, respectively. In these tables, we have shown the error bounds and depth of different arity approximating subdivision volumetric by using Lemma 3.3 and Theorem 3.5 with πœ’ = 0 . 1 .

tab3
Table 3: Error bound of tensor product 𝑛 -ary approximating volumetric models: here 𝑛 presents the arity of subdivision volumetric model and π‘˜ presents the subdivision level.
tab4
Table 4: Subdivision depth of tensor product 𝑛 -ary approximating volumetric models: here 𝑛 presents the arity of subdivision volumetric model and πœ– presents error tolerance.

5. Conclusion and Future Work

We have computed subdivision depth based on error bounds for tensor product 𝑛 -ary volumetric models. Furthermore, we have shown that error bounds for binary subdivision volumetric model [17] is special case of our bounds. It is noticed that the increase in arity results gradually decrease in error, which is shown in Tables 1, 3 and graphically in Figure 4. It is noticed from Tables 2 and 4 that higher arity subdivision volumetric models need less number of subdivision steps than lower arity to satisfy user-specified error tolerance.

fig4
Figure 4: (a) Presents comparison among the error bounds of different arity subdivision volumetric models (4.1), that is, for 𝑛 = 2 , 3 , 4 , 5 , 6 ; (b) presents comparison among the error bounds of different arity subdivision volumetric models (4.3), that is, for 𝑛 = 2 , 3 , 4 , 5 , 6 . Here π‘˜ presents subdivision level, 𝑛 presents the arity and πœ– presents the user-specified error tolerance.

The authors are looking, as a future work, to extend the computational techniques of subdivision depth for 𝑛 -ary arbitrary subdivision volumetric models over rectangular/triangular hexahedron lattice. we will discuss them elsewhere.

Appendices

A. Proof of Lemma 3.1

Proof. From (2.2), (2.4) for 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 , we obtain 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 + 1 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š 𝑠 = 0  𝑑 = 0 π‘Ž 𝛽 , 𝑠 π‘Ž 𝛾 , 𝑑  π‘š  π‘Ÿ = 0 𝑏 𝛼 , π‘Ÿ ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 , 𝑙 + 𝛾 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝛾  ξƒͺ , ( A . 1 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 + 1 , 𝑛 𝑙 + 𝛾 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛾 , 𝑑  π‘š  𝑠 = 0 𝑏 𝛽 , 𝑠 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑  ξƒͺ , ( A . 2 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝛽 + 1 , 𝑛 𝑙 + 𝛾 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘Ž 0 , π‘Ÿ π‘Ž 𝛾 , 𝑑  π‘š  𝑠 = 0 𝑏 𝛽 , 𝑠 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 , 𝑙 + 𝑑  ξƒͺ , ( A . 3 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝛽 + 1 , 𝑛 𝑙 + 𝑛 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝑛 = π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑑  π‘š  𝑠 = 0 𝑏 𝛽 , 𝑠 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 , 𝑙 + 𝑑 + 1  ξƒͺ , ( A . 4 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 + 1 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 𝛽 , 𝑠  π‘š  𝑑 = 0 𝑏 𝛾 , 𝑑 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑  ξƒͺ , ( A . 5 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 + 1 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘Ž 0 , π‘Ÿ π‘Ž 𝛽 , 𝑠  π‘š  𝑑 = 0 𝑏 𝛾 , 𝑑 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 , 𝑙 + 𝑑  ξƒͺ , ( A . 6 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝛾 + 1 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 0 , 𝑠  π‘š  𝑑 = 0 𝑏 𝛾 , 𝑑 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑  ξƒͺ , ( A . 7 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝛾 + 1 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝑛 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑠  π‘š  𝑑 = 0 𝑏 𝛾 , 𝑑 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑  ξƒͺ , ( A . 8 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 + 1 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝛾 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝛾 = π‘š  π‘š 𝑠 = 0  𝑑 = 0 π‘Ž 0 , 𝑠 π‘Ž 𝛾 , 𝑑  π‘š  π‘Ÿ = 0 𝑏 𝛼 , π‘Ÿ ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑  ξƒͺ , ( A . 9 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 + 1 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝑛 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝑛 = π‘š  π‘š 𝑠 = 0  𝑑 = 0 π‘Ž 𝛽 , 𝑠 π‘Ž 0 , 𝑑  π‘š  π‘Ÿ = 0 𝑏 𝛼 , π‘Ÿ ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 + 1  ξƒͺ , ( A . 1 0 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 + 1 , 𝑛 𝑙 + 𝑛 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝛽 , 𝑛 𝑙 + 𝑛 = π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘Ž 𝛼 , π‘Ÿ π‘Ž 0 , 𝑑  π‘š  𝑠 = 0 𝑏 𝛽 , 𝑠 ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 + 1  ξƒͺ , ( A . 1 1 ) 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 + 1 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝑛 βˆ’ 𝑝 π‘˜ 𝑛 𝑖 + 𝛼 , 𝑛 𝑗 + 𝑛 , 𝑛 𝑙 + 𝑛 = π‘š  π‘š 𝑠 = 0  𝑑 = 0 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑  π‘š  π‘Ÿ = 0 𝑏 𝛼 , π‘Ÿ ξ‚€ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ + 1 , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 + 1 βˆ’ 𝑝 π‘˜ βˆ’ 1 𝑖 + π‘Ÿ , 𝑗 + 𝑠 + 1 , 𝑙 + 𝑑 + 1  ξƒͺ , ( A . 1 2 ) where 𝑏 𝛾 , π‘Ÿ is defined by (2.6) and 𝛼 , 𝛽 , 𝛾 = 0 , 1 , … , 𝑛 βˆ’ 1 .
Now using (A.1) recursively together with notations defined by (2.8), we get m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 1 β€– β€– β©½ βŽ› ⎜ ⎜ ⎝ m a x 𝛼 , 𝛽 , 𝛾 | | | | | π‘š  π‘š 𝑠 = 0  π‘š 𝑑 = 0  π‘Ÿ = 0 π‘Ž 𝛽 , 𝑑 π‘Ž 𝛼 , 𝑠 𝑏 𝛾 , π‘Ÿ | | | | | ⎞ ⎟ ⎟ ⎠ π‘˜ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 1 β€– β€– . ( A . 1 3 ) From (2.5) and using the above inequality, we get m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 1 β€– β€– β©½ ( 𝛿 ) π‘˜ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 1 β€– β€– . ( A . 1 4 ) Again using (A.2)–(A.4) recursively and by utilizing (2.5) and (2.8), we have m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 2 β€– β€– β©½ ( 𝛿 ) π‘˜ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 2 β€– β€– . ( A . 1 5 ) Further using (A.5)–(A.8) recursively and by utilizing (2.5) and (2.8), we have m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 3 β€– β€– β©½ ( 𝛿 ) π‘˜ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 3 β€– β€– . ( A . 1 6 ) Similarly, using (A.9)–(A.12) recursively together with (2.5) and (2.8) separately for each 𝑑 = 4 , 5 , 6 , 7 , respectively, we have m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 𝑑 β€– β€– β©½ ( 𝛿 ) π‘˜ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” 0 𝑖 , 𝑗 , 𝑙 , 𝑑 β€– β€– . ( A . 1 7 ) This completes the proof.

B. Proof of Lemma 3.2

Proof. From (2.2) and (2.4), 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 , 𝑛 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 = π‘š βˆ‘ π‘Ÿ = 0 π‘Ž 0 , π‘Ÿ ξ‚΅ π‘š βˆ‘ 𝑠 = 0 π‘Ž 0 , 𝑠 ξ‚΅ π‘š βˆ‘ 𝑑 = 0 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  . ξ‚Ά ξ‚Ά ( B . 1 ) By expanding innermost summation, we get π‘š  𝑑 = 0 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž 0 , 1 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž 0 , 2 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 2 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 + 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + β‹― + π‘Ž 0 , π‘š ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘š βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘š βˆ’ 1 + β‹― + 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  , π‘š  𝑑 = 0 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  +  π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξƒͺ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž  , ( B . 2 ) where Μƒ π‘Ž 0 , π‘ž is defined by (2.16). Now π‘š  𝑠 = 0 π‘Ž 0 , 𝑠  π‘š  𝑑 = 0 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  ξƒͺ = π‘Ž π‘š 0 , 0  𝑠 = 0 π‘Ž 0 , 𝑠 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  π‘š 𝑝 = 1  𝑠 = 0 π‘Ž 0 , 𝑝 π‘Ž 0 , 𝑠 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑠 = 0 π‘Ž 0 , 𝑠 π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž  . ( B . 3 ) Since π‘š  𝑠 = 0 π‘Ž 0 , 𝑠 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙  , π‘š  𝑠 = 0 π‘Ž 0 , 𝑠 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙 + 1  , ( B . 4 ) then (B.3) implies π‘š  𝑠 = 0 π‘Ž 0 , 𝑠  π‘š  𝑑 = 0 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  ξƒͺ = π‘Ž 2 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž 0 , 0 π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  π‘š 𝑝 = 1  π‘ž = 1 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑝 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙 + 1  + π‘š  𝑠 = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑠 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž  . ( B . 5 ) Substituting it into (B.1), we get 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 , 𝑛 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 = ξ€· π‘Ž 0 , 0 ξ€Έ 2 π‘š  π‘Ÿ = 0 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  π‘š 𝑝 = 1  π‘Ÿ = 0 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  π‘š 𝑝 = 1  π‘š π‘ž = 1  π‘Ÿ = 0 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘ž π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  π‘š 𝑝 = 1  π‘Ÿ = 0 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  π‘Ÿ = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙  + π‘š  π‘š 𝑝 = 1  π‘Ÿ = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙 + 1  + π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑠 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž  . ( B . 6 ) Since π‘š  π‘Ÿ = 0 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 , 𝑙  , π‘š  π‘Ÿ = 0 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 + 1 , 𝑙  , π‘š  π‘Ÿ = 0 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 , 𝑙 + 1  , π‘š  π‘Ÿ = 0 π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  = π‘Ž 0 , 0 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 + 1 , 𝑙 + 1  , ( B . 7 ) then substituting these summations into (B.6) and rearranging, we get 𝑝 π‘˜ + 1 𝑛 𝑖 , 𝑛 𝑗 , 𝑛 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 = π‘Ž 2 π‘š 0 , 0  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  π‘š 𝑝 = 1  π‘ž = 1 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙  + π‘Ž π‘š 0 , 0  π‘š 𝑝 = 1  π‘ž = 1 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1  + π‘š  π‘š 𝑝 = 1  π‘ž = 1 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1  + π‘š  π‘š 𝑝 = 1  π‘š π‘ž = 1  π‘Ÿ = 1 π‘Ž 0 , 𝑝 π‘Ž 0 , π‘ž π‘Ž 0 , π‘Ÿ ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1  + 𝑁 2 π‘˜ , ( B . 8 ) where 𝑁 2 π‘˜ = π‘Ž 2 0 , 0 π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑝 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 + 1 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑝 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 , 𝑙 + 1  + π‘š  π‘š π‘Ÿ = 1  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑝 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 + 1 , 𝑙 + 1  + π‘Ž π‘š 0 , 0  π‘Ÿ = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙  + π‘š  π‘š π‘Ÿ = 0  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑝 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙 + 1  + π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑠 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž  . ( B . 9 ) This implies 𝑀 π‘˜ 0 , 0 , 0 β©½ βŽ› ⎜ ⎜ ⎝ | | | | | ξ€· π‘Ž 0 , 0 ξ€Έ 2 π‘š  𝑑 = 1 π‘Ž 0 , 𝑑 | | | | | + | | | | | ξ€· π‘Ž 0 , 0 ξ€Έ 2 π‘š βˆ’ 1  𝑠 = 1 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 1 β€– β€– + βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž π‘š 0 , 0  𝑑 = 1 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘Ž π‘š 0 , 0  π‘Ÿ = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 2 β€– β€– + βŽ› ⎜ ⎜ ⎝ | | | | | π‘š  𝑑 = 1 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 3 β€– β€– + βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž π‘š 0 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘Ž π‘š 0 , 0  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 4 β€– β€– + βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž π‘š 0 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘Ž π‘š 0 , 0  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 5 β€– β€– + βŽ› ⎜ ⎜ ⎝ | | | | | π‘š  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘š  π‘š 𝑑 = 1  π‘Ÿ = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , 𝑑 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 6 β€– β€– + βŽ› ⎜ ⎜ ⎝ | | | | | π‘š  π‘š π‘Ÿ = 1  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘š  π‘š π‘Ÿ = 1  𝑑 = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ m a x 𝑖 , 𝑗 , 𝑙 β€– β€– Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 7 β€– β€– , ( B . 1 0 ) where 𝑀 π‘˜ 0 , 0 , 0 and Ξ” π‘˜ 𝑖 , 𝑗 , 𝑙 , 𝑑 , 𝑑 = 0 , 1 , … , 7 are defined by (2.7) and (2.8). Now using notations (2.5) and (2.9), we have 𝑀 π‘˜ 0 , 0 , 0 β©½ ( 𝛿 ) π‘˜ πœ’ ⎧ βŽͺ ⎨ βŽͺ ⎩ βŽ› ⎜ ⎜ ⎝ | | | | | ξ€· π‘Ž 0 , 0 ξ€Έ 2 π‘š  𝑑 = 1 π‘Ž 0 , 𝑑 | | | | | + | | | | | ξ€· π‘Ž 0 , 0 ξ€Έ 2 π‘š βˆ’ 1  𝑠 = 1 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž π‘š 0 , 0  𝑑 = 1 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘Ž π‘š 0 , 0  π‘Ÿ = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝ | | | | | π‘š  𝑑 = 1 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘š  π‘š π‘Ÿ = 0  𝑑 = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž π‘š 0 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘Ž π‘š 0 , 0  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž π‘š 0 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘Ž π‘š 0 , 0  𝑑 = 1 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝ | | | | | π‘š  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘š  π‘š 𝑑 = 1  π‘Ÿ = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , 𝑑 π‘Ž 0 , π‘Ÿ Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝ | | | | | π‘š  π‘š π‘Ÿ = 1  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 | | | | | + | | | | | π‘š  π‘š π‘Ÿ = 1  𝑑 = 0 π‘š βˆ’ 1  𝑠 = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑑 Μƒ π‘Ž 0 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ ⎫ βŽͺ ⎬ βŽͺ ⎭ . ( B . 1 1 ) Using (2.10)–(2.16) for 𝛼 , 𝛽 , 𝛾 = 0 , we have 𝑀 π‘˜ 0 , 0 , 0 β©½ ( 𝛿 ) π‘˜ πœ’ ξ‚€ πœ‚ 1 0 , 0 , 0 + πœ‚ 2 0 , 0 , 0 + β‹― + πœ‚ 7 0 , 0 , 0  = ( πœ’ ) ( 𝛿 ) π‘˜ 7  𝑑 = 1 πœ‚ 𝑑 0 , 0 , 0 . ( B . 1 2 ) Similarly from (2.2) and (2.4) for 𝛼 = 1 , 𝛽 = 𝛾 = 0 , 𝑝 π‘˜ + 1 𝑛 𝑖 + 1 , 𝑛 𝑗 , 𝑛 𝑙 βˆ’ 1 𝑛  ( 𝑛 βˆ’ 1 ) 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙  = π‘š βˆ‘ π‘Ÿ = 0 π‘Ž 1 , π‘Ÿ ξ‚΅ π‘š βˆ‘ 𝑠 = 0 π‘Ž 0 , 𝑠 ξ‚΅ π‘š βˆ‘ 𝑑 = 0 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + 𝑑 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  . ξ‚Ά ξ‚Ά ( B . 1 3 ) Now after expanding and rearranging the above summation, we have 𝑝 π‘˜ + 1 𝑛 𝑖 + 1 , 𝑛 𝑗 , 𝑛 𝑙 βˆ’ 1 𝑛  ( 𝑛 βˆ’ 1 ) 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙  =  π‘Ž 2 π‘š 0 , 0  𝑝 = 1 π‘Ž 1 , 𝑝 βˆ’ 1 𝑛 ξƒͺ ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘Ž 1 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 1 , 𝑠 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙  + π‘š  𝑝 = 1 π‘Ž 0 , 𝑝 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙  + π‘š  π‘š π‘Ÿ = 1  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 1 , π‘Ÿ π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1  + π‘Ž π‘š 0 , 0  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 1 , 𝑠 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1  + π‘š  π‘š 𝑠 = 1  𝑑 = 1 π‘Ž 0 , 𝑠 π‘Ž 0 , 𝑑 ξ‚€ 𝑝 π‘˜ 𝑖 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 , 𝑗 , 𝑙 + 1  + 𝑁 3 π‘˜ , ( B . 1 4 ) where 𝑁 3 π‘˜ = π‘Ž 2 0 , 0 π‘š βˆ’ 1  π‘ž = 1 Μƒ π‘Ž 1 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑝 Μƒ π‘Ž 1 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 + 1 , 𝑙  + π‘Ž π‘š 0 , 0  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , 𝑝 Μƒ π‘Ž 1 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 , 𝑙 + 1  + π‘š  π‘š π‘Ÿ = 1  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ π‘Ž 0 , 𝑝 Μƒ π‘Ž 1 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘ž + 1 , 𝑗 + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘ž , 𝑗 + 1 , 𝑙 + 1  + π‘Ž π‘š 0 , 0  π‘Ÿ = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 1 , π‘Ÿ Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙  + π‘š  π‘š π‘Ÿ = 0  𝑝 = 1 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 0 , π‘Ÿ π‘Ž 1 , 𝑝 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž + 1 , 𝑙 + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + π‘ž , 𝑙 + 1  + π‘š  π‘š π‘Ÿ = 0  𝑠 = 0 π‘š βˆ’ 1  π‘ž = 1 π‘Ž 1 , π‘Ÿ π‘Ž 0 , 𝑠 Μƒ π‘Ž 0 , π‘ž ξ‚€ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž + 1 βˆ’ 𝑝 π‘˜ 𝑖 + π‘Ÿ , 𝑗 + 𝑠 , 𝑙 + π‘ž  . ( B . 1 5 ) Now using notations (2.5), (2.7)–(2.9) and using Lemma 3.1, we have 𝑀 π‘˜ 1 , 0 , 0 β©½ ( 𝛿 ) π‘˜ πœ’ ⎧ βŽͺ ⎨ βŽͺ ⎩ βŽ› ⎜ ⎜ ⎝ | | | | | π‘Ž 2 π‘š 0 , 0  𝑑 = 1 π‘Ž 1 , 𝑑 βˆ’ 1 𝑛 | | | | | + | | | | | π‘Ž 2 0 , 0 π‘š βˆ’ 1  𝑠 = 1 Μƒ π‘Ž 1 , 𝑠 | | | | | ⎞ ⎟ ⎟ ⎠ + βŽ› ⎜ ⎜ ⎝