Efficient WENO schemes for nonuniform grids

A set of arbitrarily high-order WENO schemes for reconstructions on nonuniform grids is presented. These non-linear interpolation methods use simple smoothness indicators with a linear cost with respect to the order, making them easy to implement and computationally efficient. The theoretical analysis to verify the accuracy and the essentially non-oscillatory properties are presented together with some numerical experiments involving algebraic problems in order to validate them. Also, these general schemes are applied for the solution of conservation laws and hyperbolic systems in the context of finite volume methods.

The idea behind the design of the WENO algorithm is quite powerful.It is based on the construction of a piecewise interpolator trying to avoid the points where the function is discontinuous.For this, a convex combination of interpolations of lower degree is proposed, using non-linear weights for the contribution of each interpolator.The final result is an interpolator of high-order accuracy if the data used are free of discontinuities and a lower order of accuracy interpolator if the utilized stencil crosses a singularity.The main ingredients of WENO algorithms are the optimal weights, which are the values needed to obtain the highest order of accuracy, and the smoothness indicators which indicate if a subset of points contains a discontinuity.There are certain situations where the optimal weights are negative and non-linear WENO weights do not satisfy the required conditions (see [27]).Some solutions for this issue are summarized in [27].Levy et al. proposed in [21] a central scheme based on a centered version of WENO, named CWENO, which overcomes the mentioned matter.
On the one hand, in [7], the authors introduced a central WENO scheme based on a global average weight independent of the optimal weights.This method has two key ingredients: first, the global average weight, which is defined using the smoothness indicators through elementary operations, and second, the evaluation of the interpolator polynomial using the whole stencil.This method was developed for structured grids.However, in some applications (see [12,25], mentioned in [27]), WENO reconstructions on non-uniform meshes are needed.
On the other hand, in [11], the authors presented an alternative family of smoothness indicators that are simpler and computationally cheaper than those originally proposed by Jiang and Shu in [19].They proved that the computational time needed to obtain the proposed smoothness indicators increased linearly with respect to the accuracy order of the WENO scheme considered while for Jiang and Shu's smoothness indicators the growth is quadratic.
In this paper, we combine and extend the ideas proposed in [7,11] for nonuniform grids and prove theoretically that the new proposed scheme satisfies the desired accuracy.In addition, we show that this scheme efficiently solves, in some cases, the loss of order of accuracy near smooth extrema associated with the original WENO and CWENO schemes.We will present some mimetic examples to check these theoretical results.Another key feature of this method is that the cost of the computation of the smoothness indicators for a nonuniform stencil does not have a significant impact on the whole computational cost of the scheme, compared with the cost that the computation of the smoothness indicators has on the original WENO and CWENO schemes.We will apply the new WENO reconstruction method in finite volume schemes to obtain numerical solutions in shock problems from hyperbolic conservation laws to illustrate its performance.
1.2.Outline of this paper.The remainder of the paper is divided as follows: in Section 2 we briefly review the classical Lagrange interpolation and the notation that will be used throughout this paper.We present linear reconstructions from point values and cell averages on nonuniform grids.The new WENO method will be introduced in Section 3 and some theoretical results will be provided.In particular, we will prove that the new WENO scheme satisfies the necessary properties to achieve maximum order of accuracy when the data do not cross any singularity and the resulting interpolator is essentially non-oscillatory.Section 4 is devoted to showing the computational aspects of the algorithm.In Section 5, we present some numerical experiments to validate the theoretical results and to test the performance of the new scheme in shock problems from hyperbolic conservation laws.Finally, some conclusions and remarks will be drawn in Section 6.

Preliminaries
We start introducing basic definitions which will be used throughout this work.We also review the definition of linear pointwise reconstructions on nonuniform grids from given point values and cell averages.Definition 2.1.In the remainder of the paper, we denote by Π k , k ∈ N 0 , the space of polynomials of maximal degree k, i.e.
and by Πk the space of polynomials of exact degree k, k ∈ N 0 , i.e.
2.1.Reconstructions on nonuniform grids.Since we want to design WENO schemes working on nonuniform grids, the first step is to describe the procedure to obtain a reconstruction from given data at the desired order.We will focus both on the pointwise reconstructions from point values (the traditional Lagrange interpolation/extrapolation) and the reconstructions from cell averages.
Reconstructions from point values.Let S = {x 0 , . . ., x R−1 } be a stencil, x i < x j , i < j, and f i := f (x i ), 0 ≤ i ≤ R − 1 for some function f .It is well-known that the polynomial p ∈ Π R−1 satisfying p(x i ) = f i is given by the Lagrange form: Reconstructions from cell averages.Now, let us consider a stencil consisting of cell interfaces given by S = {x 0 , . . ., x R }, x i < x j , i < j, being the associated cells

and fi the known values given by fi
Then, let us assume that we want to find p ∈ Π R−1 such that p(x) interpolates f (x).To do so, we impose (2.1) To find p, we use the strategy proposed in [3] for structured grids, which we next generalize to the nonuniform case.Let us consider with p i ∈ Π R−1 to be determined, 0 ≤ i ≤ R − 1.Then, imposing (2.1) on (2.2) and using linearity, we have for a fixed j, with δ i,j the Kronecker delta.Let P i ∈ Π R be a primitive of p i , namely, such that Hence, it suffices to impose to write P i in Lagrange form as Thus, taking derivatives at both sides of (2.5) we get and (2.2) then reads (2.7)

Generalized WENO schemes
The classical Lagrange methods introduced in Section 2.1 are useful to interpolate smooth data but some artifacts appear when they are used in applications with non-continuous inputs.To solve this problem, we present in this section a generalization of the WENO method for nonuniform grids.
3.1.Construction of the method and theoretical results.Let h > 0, R ≥ 3 and let S R,h = {x 0,h , . . ., x R−1,h } be a stencil, with x i,h = z + c i h, for some fixed z ∈ R and We denote by x * h ∈ R the point in which we desire to reconstruct the data from the stencil S R,h .This point is located between x r−1,h and x r,h if R is even, with r = R/2, or rather between x r−1,h and x r,h (left-biased) or between x r,h and x r+1,h (right-biased) if R is odd, with r = (R − 1)/2.Note that this is actually a generalization of the traditional centered or left/right-biased WENO/CWENO schemes on cartesian grids to nonuniform grids, respectively, (see e.g.[21,22]).
Let us first consider p R,i,h the reconstruction polynomials satisfying p R,i,h (x j,h ) = L[f ](x j,h ), where the operator L is defined depending on the framework chosen (point values or cell average reconstructions), i.e., ), for point value reconstructions, fj,h , for cell average reconstructions, with i ≤ j ≤ i + r, for 0 ≤ i ≤ r ′ , with r = ⌊(R − 1)/2⌋ and r ′ = ⌈(R − 1)/2⌉, and evaluate them at x * h , where ⌊•⌋ denotes the floor function and ⌈•⌉ the ceil function.Following the same idea based on the Yamaleev-Carpenter weight construction proposed in [31], we define d R,h as where f [x 0,h , . . ., x R−1,h ] is the undivided difference defined as The simplified smoothness indicators, akin to the idea proposed in [11], which in turn is the key to the simplicity of this algorithm, are defined by Since in this case, in the context of an arbitrary grid, there is no point in computing ideal weights in order to build the reconstruction with optimal accuracy from the lower order reconstruction, we simply consider a uniform convex combination based on the r ′ + 1 values 1 r ′ +1 and define with s an exponent to be determined (will be done so when the accuracy analysis is performed) and ε > 0 a small positive quantity in order to avoid divisions by zero.With this weight design, we can define the following weighted essentially non-oscillatory reconstruction: Since this reconstruction does not attain the desired optimal order, R, we use the reconstruction associated with the whole stencil S R,h , p R,h satisfying p R,h (x i,h ) = L[f ](x i,h ), which we evaluate at x * h , and the following global average weight, given by Then, the final reconstruction with essentially non-oscillatory properties is given by q R,h Now that we have introduced and defined all the components of the WENO method, we are in a position to prove the principal properties of the scheme in terms of the order of the weights and accuracy of the reconstructions obtained.Proposition 3.1.For the standard case, and keeping the same notation, if f has a critical point of order k at z, 0 where ).On the other hand, by all the considerations derived from [8, Section 2] and [11, Section 2] we conclude that, since r Therefore, we have Now, let us assume that a discontinuity crosses S R,h .Then, we define as A the set of indexes j such that the discontinuity does not cross S R,j,h and as B the set of indexes j such that the discontinuity crosses S R,j,h .Note that, by construction of r ′ , the stencils do not overlap, and thus A ̸ = ∅.Then, on the one hand, if j ∈ A: On the other hand, if j ∈ B: Consequently, if i ∈ A, we have: Now we have to analyze ω R,h .To do so, we first study the accuracy of J R,h,s : where the last equality holds since m then the optimal accuracy is attained both in the case of smoothness and the case of discontinuities, namely: Proof.Let us first assume that f ∈ C R .Then, by Proposition 3.1: with m = max 0≤i<R−1 m i , where Now, taking into account our hypothesis involving the critical points, we have that , and thus we can assume k < R − 2 for any R ≥ 3 (the case k > R − 2 being trivial, as we will point out at the end of the proof).Therefore, taking into account that and using that , where in the last equality we have used that R ≥ r + 1.Therefore, if one wants the optimal R-th order accuracy, one must impose On the other hand, since Let us consider now the case where a discontinuity crosses S R,h , then, again by Proposition 3.1, we have Combining in this case (3.4) and (3.5), we can conclude that Then, in order to attain the optimal accuracy, one must impose and therefore, it suffices with the choice to attain optimal accuracy both in the smooth and the discontinuous case.□ 3.2.Remarks regarding the computational cost.As stated in the introduction, one of the strongest points of the proposed WENO method is the ease of computation of the smoothness indicators, even taking into consideration the fact that it is defined in the context of a nonuniform grid.Unlike the case of the traditional Jiang-Shu smoothness indicators, which in turn are used in the method proposed in [7], the indicators used in this work, based on those proposed in [11], are considerably cheaper to compute, especially when the grid is non-uniform.
To be more specific, given R ∈ N, a stencil of the form ), for f a function to be approximated, and the reconstruction point x * h = z + c * h, the simplified expression for the classical Jiang-Shu smoothness indicators has the structure with α i,j,k to be determined for each 0 From (3.6), it can be clearly seen that the number of operations that are required to compute the Jiang-Shu smoothness indicators in terms of the accuracy order is O(R 3 ).Therefore, any WENO method using these smoothness indicators has, at least, a cubic cost with respect to the aforementioned order of the method.This includes the techniques proposed by Levy, Puppo and Russo in [21] or those proposed in [7].
In contrast, from (3.1) the smoothness indicators proposed in this paper have linear cost with respect to the accuracy order, namely, O(R).It is also important to mention that, depending on the coefficients c i and the location of the reconstruction point x * h , the ideal weights in the method proposed by Levy, Puppo and Russo can be negative.This issue was fixed in [7] and the method proposed here neither suffers from this phenomenon.
Since computing an interpolating polynomial under a non-uniform grid requires evaluating an expression of the form with β i,k to be determined, 0 ≤ i, k ≤ R − 1, the cost from that part of the algorithm requires a quadratic number of operations with respect to the accuracy order, namely, O(R 2 ).Thus, taking into account that ⌊ R 2 ⌋ polynomials and smoothness indicators have to be computed in the process, the computational cost of the WENO methods using the Jiang-Shu smoothness indicators have O(R 4 ) cost when they are applied on nonuniform grids, while the WENO method proposed in this paper has O(R 3 ) cost with respect to the accuracy order in a context of nonuniform grids.
In the particular case in which a static nonuniform grid is used during the whole time evolution, the terms α i,j,k and β i,k can be precomputed for each local stencil and stored to be used on each time iteration.This implies that the computational cost with respect to the order of the expression in (3.6) is reduced to quadratic, while the cost in (3.7) is reduced to linear.Therefore, the global computational cost of the CWENO methods [7,21] can be reduced to O(R 3 ), while the cost associated with the method presented in this paper is reduced to O(R 2 ).

Summary of the algorithms
In this Section, we summarize the algorithm of the proposed WENO method, including the considerations involving the items to attain unconditionally the optimal accuracy.We divide this section into two parts: First, we introduce the description for the reconstruction from point values and after we present the necessary changes to adapt it for reconstructions from cell average values.

Algorithm for reconstructions from point values.
Input data: Then, the next steps are to be followed: (2) Compute the reconstruction polynomials at c * h of the corresponding substencils S R,i,h , 0 ≤ i ≤ r ′ , p R,i,h , and the one associated to the global stencil, S R,h , p R,h , which are: (3) Compute the smoothness indicators: Optionally, the denominators can be omitted in both expressions without any impact on the accuracy of the scheme.
(4) Obtain d R,h : ( Generate the weights: Obtain the half-optimal-order essentially non-oscillatory reconstruction: (8) Compute the optimal order essentially non-oscillatory reconstruction:

Algorithm for reconstructions from cell average values.
Input data: In step (2) reconstruction polynomials should be defined by: In step (3) the smoothness indicators are: and in step (4): Remark 4.1.Note that, if one works with fixed grids, then all the terms involving computations related to the grid spacings, namely, c i , as well as the relative position of the reconstructions points c * h , can be computed and stored only once, and then used to perform the computations, thus reducing to linear cost with respect to the order the computation of terms such as p R,i,h , p R,h and d R,h , yielding a comparable performance with respect to the classical WENO schemes with structured grids.

Numerical experiments
In this section, we perform some numerical examples to check the accuracy properties of our schemes.We divide the experiments into two sections: Firstly, we apply our algorithm in mimetic algebraic numerical tests to assess the order of accuracy.Secondly, we numerically solve several PDEs.In both cases, we implement the method to non-continuous data.
In all the experiments the C++ wrapper [18] of the GNU MPFR library [23] is used with a precision of 332 bits (≈ 100 digits) and taking ε = 10 −10 5 .5.1.Algebraic numerical experiments.The goal of this subsection is to check the behaviour of the proposed generalized WENO scheme in different scenarios.We present two different tests: one with a smooth function, and one for non-smooth functions computed both for point values and for cell average values for which the estimated order of accuracy is calculated.
In all the numerical experiments, we consider a mesh of the form x i,h = z + c i h n , 0 ≤ i < N , with h n+1 = hn 2 and provide the table of the resulting experiments based on taking h 0 = 0.2 and 0 ≤ n ≤ 20.The order values are thus computed as with x * hn the reconstruction node and u * hn the corresponding resulting WENO reconstruction.
5.1.1.Test 1: Reconstruction of a smooth function.Let us consider the function f (x) := xe x and the grid given in Figure 1 (a), defined by the values c i , i = 0, . . ., 11 showed in Table 1 (column 1), with N = 12 and z = 0. We consider the reconstruction point x * h = z + c * h, with c * = 0, hence in this case x * h is always located between the nodes x 5,h and x 6,h .
using the same function f and the grid given in Figure 1 (b) defined by the values c i , i = 0, . . ., 11 in Table 1 (column 2).We consider x * h = z + ch, with c * = 0.As can be seen in Table 2, the behaviour of the numerical solutions is the same as before: the expected order of accuracy is reached.

Test 2: Reconstruction of a discontinuous function. Let us now consider now the function
and the grid given by the data shown in column 3 of Table 1 with c * = 2.3251.Note that the reconstruction point is now located between x 4,h and x 5,h .In this case, the discontinuity divides the stencil into two halves, where the first three nodes contain the information corresponding to the left side of the discontinuity, namely, xe x , and the last eight nodes contain the information related to the right side of the discontinuity, that is, 2xe x + 1.The expected behaviour is therefore an accuracy drop to sixth order (since ⌊ 11+1 2 ⌋ = 6), which is confirmed in Table 3.
We repeat the experiment for cell averages with nodes given by the values c i , 0 ≤ i ≤ 11 in the fourth column of Table 1 and c * = 0.5041.In this case, the reconstruction point is located in the cell with bounds x 5,h and x 6,h as it can be seen in Figure 1 (d).As before, the discontinuity separates the cells into two parts.The first one is given by the first five cells, and the second one contains six cells with the information related to the right side of the discontinuity.As we expected, in Table 3 we observe that the accuracy order reaches sixth order, as we have studied in Theorem 3.1.

Numerical experiments involving conservation laws.
To illustrate the applicability of our proposed WENO interpolator in the context of numerical solutions of PDEs we next consider some experiments using one-dimensional hyperbolic conservation laws u t + f (u) x = 0. First, we will test our method using some academic examples such as the linear advection equation, with f (u) = u, and Burgers equation, which corresponds with f (u) = u 2 /2.Later, we will perform two more sophisticated experiments to compare the results obtained with the non-uniform approach and the uniform one in more challenging scenarios.
For this purpose, a finite volume scheme will be used (see e.g.[26]), which is suitable to attain a method with high order accuracy.Given a domain Ω = [a, b] and a partition of Ω, , the semidiscrete form of the conservation law reads where To approximate (5.2) for approximations u i (t) ≈ u t (x i , t), we use the conservative difference formula: where the numerical fluxes, f i+ 1 2 ≈ f (u(x i+ 1 2 , t)), are approximated using: For the autonomous linear advection equation: Left-biased upwind reconstructions are applied: 2 ), where q − 5,h is the WENO reconstruction using five cells left-biased, i.e., u i+l , −3 ≤ l ≤ 2. For the rest of the numerical experiments in this work: Local Lax-Friedrichs flux splitting is utilized: where q ± 5,h are the WENO reconstructions using five cells left and rightbiased and α = max u |f ′ (u)|.
The time-stepping solver that will be used along all the experiments is the secondor third-order Runge-Kutta TVD scheme [28].The spatial discretization will always have order 5.

5.2.1.
Construction of the nonuniform grid.For the sake of reproducibility, we provide all the details about the non-uniform grid construction in each experiment.Given n ∈ N, we define The parameter ξ controls the fluctuation of the cell interfaces with respect to a uniform grid, and the function r : N 3 → [0, 1] is the Wichmann-Hill pseudo-random number generator given by r(s with seed update s 1 (j) = 171s 1 (j 1) mod 30269, s 2 (j) = 172s 2 (j − 1) mod 30307, s 3 (j) = 170s 3 (j − 1) mod 30323.
In all the experiments, the selected starting seeds will be chosen as In the experiments analyzing the order under smooth solutions, namely, those involving several simulations in which the number of cells is duplicated, the seeds are not restarted upon the next doubled resolution; instead, the last updated seed from the previous experiment with halved resolution is used as starting seed for the next one.
5.2.2.Test 3: Linear advection equation.First of all, we compute the approximations of the solution of the linear advection equation: with initial and boundary conditions: which exact solution is given by: u(x, t) = 0.25 + 0.5 sin (π(x − t)) .
As the method used to solve the ODE is third-order accurate, for test purposes we choose ∆t = ( ∆x) 5/3 with ∆x = min i {∆x i }.With this selection, we guarantee that ∆t 3 = O(∆x 5 ), with ∆x = max i {∆x i } and ∆t/∆x ≤ 1.In this case, the fluctuation parameter is ξ = 0.1.We calculate the error for T = 1 and the numerical order in a similar way as we computed in Equation (5.1) for values of n = 20 • 2 j with j = 0, . . ., 4. The results are displayed in Table 4.It is clear that the expected order of accuracy is reached.We consider now the discontinuous initial data given by: u(x, 0) = −0.25,x ≤ 0, 1, x > 0, ∆t = 0.9 ∆x and the fluctuation parameter ξ = 0.25.The simulation is now run until T = 1.5 and with n = 40 and n = 100.We can see in Figure 2 that the behaviour of the numerical solution is oscillation-free and does not present Gibbs phenomenon close to the discontinuities present in the solution.5 and in Figure 3. Again, the expected approximation order, 5, is achieved for both norm 1 and norm infinity.The simulations in Figure 3 show that for small values of n some numerical oscillations appear near the boundaries.However, this problem is overcome when we refine the mesh.using the nonuniform approach significantly outperform those of the uniform one.We consider the equation (see [24]) with initial condition u(x, 0) = 1, x ∈ [0, 2π], and periodic boundary conditions.By the method of characteristics, the exact solution of (5.3) is given by For a fixed value of x ∈ [0, 2π], there holds and for fixed t ≥ 0 2π 0 u(x, t)dx = 2π, therefore, the solution of this equation behaves progressively like a Dirac delta function as time increases.This means that the solution will have a strong gradient near x = π, while it would remain close to zero in the rest of the domain.The structure of the solution suggests that the use of a nonuniform grid with decreasing distance between nodes as we approach the point x = π could be convenient.We set T = 8 and consider uniform grids with n = 100 • 2 j − 1, j = 1, . . ., 10 cells.We define non-uniform grids of 2m cells with cell boundaries defined by where u N j is the result of N time steps of length ∆t = T /N , subject to the following CFL restriction ∆t min The ODE solver used in this test is TVD-RK2.
In Figure 4 we show the numerical solutions obtained using both uniform and non-uniform meshes.Figures 4 (c) and (d) show that the accuracy of the approximated solutions obtained using non-uniform grids is higher than the ones obtained using uniform grids, even when a significantly lower resolution is considered.
Table 6 displays the errors and orders obtained when using uniform grids.For comparison, the error obtained using a non-uniform grid with 2m = 198 cells and q = 1.1 is 2.87e − 2, which is comparable to the error obtained with a uniform grid of 25599 cells.If we consider 2m = 398 cells and q = 1.04 the computed error is 2.67e − 3, which in this case is comparable to the error obtained with a uniform grid with 51199 cells.
It is worth mentioning that although the number of cells on the non-uniform grids considered is much smaller than the number of cells on the uniform ones, the smallest cell size for a non-uniform grid with 2m = 198 cells and q = 1.1 is 2.51e − 5 while for the grid with 2m = 398 and q = 1.04 is 5.13e − 5.Those sizes are smaller than the cell size for the uniform grids with 25599 and 51199 cells, given by 2.45e − 4 and 1.23e − 4, respectively.5.2.5.Test 6: Shu-Osher problem.In this test, we compare the performance of our method against the CWENO-type methods.Given that the method proposed by Levy, Puppo and Russo [21] can have negative weights depending on the distribution of the grid and the relative position of the reconstruction point, we will use the alternative proposed in [7], which was designed to avoid this issue as one of its purposes.The main goal of this test is to compare the impact of the computation of the smoothness indicators in the global cost of the method, which has been discussed theoretically in Subsection 3.2.
We simulate until T = 1.8 with a resolution of n = 256 cells, CFL = 0.5 and the same parameters as those stated in Section 5.2.1 to generate the non-uniform grid.Then, we compare the results against a reference solution computed with a resolution of N = 65536 cells.The results of the experiment can be seen in Figure 5.
In this simulation, our method is 1.13 times faster than CWENO for the fifthorder scheme and 8.24 times faster for the seventh-order scheme, mostly due to the considerably more complex computations required to obtain the classical Jiang-Shu smoothness indicators on a nonuniform grid, required by the chosen CWENO method.Furthermore, our method produces a better approximation than the CWENO method, as Figure 5 reflects.

Conclusions
In this paper, we have presented a novel WENO approach capable of performing essentially non-oscillatory reconstructions in a general context of nonuniform grids.Both the theoretical results and the numerical experiments show that the scheme has the desired accuracy.
The efficiency of the algorithm relies strongly on the simplification of the smoothness indicators that can be performed in weight designs akin to the one proposed by Yamaleev-Carpenter [31], as shown in [11].Another important factor is that the coefficients involving the grid spacing can be computed and stored previously if the numerical simulation uses a fixed grid, yielding in that case a similar efficiency to the uniform grid case.As for the cases in which the grid is not fixed, then Newton polynomials, together with the associated computation of the divided differences properly scaled, are recommended instead.
Our next goal is to tackle the challenge of redesigning the weights of the method so that it can be unconditionally optimal, regardless of the order of the critical point to which the stencil converges, for stencils with four or more nodes, similarly as done in [8] and [9] for the case of uniform stencils.
Another of our future research projects involves a generalization of the current approach into an interpolator capable of preserving the essentially non-oscillatory properties regardless of the relative position of the point to be interpolated or extrapolated with respect to the stencil.This approach can be especially interesting to handle boundary conditions in finite-difference schemes, among other possible applications.

4 x 0 x 1 x 2 − 2 x 3 x 4 x 5 0 4 x 0 x 1 x 2 x 3 − 2 x 4 x 5 0 = x * h x 6 x 7 x 8 2 x 9 x 10 x 11 4 (x 0 x 1 0x 2 x 3 x 4 2 x * h x 5 x 6 x 7 x 8 4 x 9 x 10 (d) − 4 x 0 x 1 x 2 x 3 − 2 x 4 x 5 0 x * h x 6 x 7 2 x 8 x 9 x 10 x 11 4 1Figure 1 .
Figure 1.Representation of the different nonuniform meshes considered in Tests 1 and 2 for h = 1, z = 0 and the values of c * indicated in each case.Note that to simplify the notation we have used x i to denote x i,h .(a) Mesh for Test 1 and smooth initial data, (b) mesh for Test 1 and discontinuous initial data, (c) mesh for Test 2 and smooth initial data and (d) mesh for Test 2 and discontinuous initial data.

5. 2 . 3 .
Test 4: Burgers equation.If we repeat the same experiments as in Test 3 but for the Burgers equation with T = 0.3 for smooth initial data and T = 1 for non-smooth initial data, we obtain the results shown in Table

Figure 5 .
Figure 5. Test 6: Numerical solution of the Shu-Osher problem with CWENO-type schemes with uniform (blue) and non-uniform (red) grids with n = 256 cells.

Table 2 .
Test 1: Numerical errors and orders of the reconstruction of the smooth function from point values with grid given in column 1 of Table1and from cell average with grid given in column 2 of Table1.