function order = gl_level_to_order ( dim_num, level )

%*****************************************************************************80
%
%% gl_level_to_order() converts a GL level to a GL order.
%
%  Discussion:
%
%    Gauss-Legendre rules cannot be nested, but in order to facilitate
%    comparison with Clenshaw Curtis grids, we need to be able to
%    convert a "level" into an order.  Except for the
%    first case of LEVEL = 0, the relationship is
%
%      ORDER = 2^LEVEL + 1
%
%    Level      Order
%
%    0          1
%    1          3 =  2 + 1
%    2          5 =  4 + 1
%    3          9 =  8 + 1
%    4         17 = 16 + 1
%    5         33 = 32 + 1
%
%    In this routine, we assume that a vector of levels is given,!
%    and the corresponding orders are desired.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    06 November 2006
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer DIM_NUM, the spatial dimension.
%
%    integer LEVEL(DIM_NUM), the level.
%
%  Output:
%
%    integer ORDER(DIM_NUM), the order (number of points) of the
%    corresponding rule.
%
  for dim = 1 : dim_num

    if ( level(dim) < 0 )
      order(dim) = -1;
    elseif ( level(dim) == 0 )
      order(dim) = 1;
    else
      order(dim) = ( 2^level(dim) ) + 1;
    end

  end

  return
end
