%% TEST13 tests CC_ABSCISSA_LEVEL_ND.
%
%  Modified:
%
%    25 March 2007
%
%  Author:
%
%    John Burkardt
%
  clear

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST13\n' );
  fprintf ( 1, '  CC_ABSCISSA_LEVEL_ND reports the level on which\n' );
  fprintf ( 1, '  a Clenshaw Curtis abscissa of given index will first\n' );
  fprintf ( 1, '  be generated, assuming a series of grids that grow\n' );
  fprintf ( 1, '  in order as 2^LEVEL+1.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  This routine is applied for multidimensional cases.\n' );

  base = 3;
  order = 2^base + 1;
  dim_num = 2;
  test_num = order * order;

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Base B = %d\n', base );
  fprintf ( 1, '  ORDER 2^B+1 = %d\n', order );
  fprintf ( 1, '  DIM_NUM = %d\n', dim_num );

  test_val = zeros ( 2, test_num );

  k = 0;
  for i = 0 :order - 1
    for j = 0 : order-1
      k = k + 1;
      test_val(1:2,k) = [ i, j ]';
    end
  end

  test_level = cc_abscissa_level_nd ( base, dim_num, test_num, test_val );

  fprintf ( 1, '\n' );
  fprintf ( 1, '         I         J  Level(I,J)\n' );
  fprintf ( 1, '\n' );

  for i = 1 : test_num
    fprintf ( 1, '  %8d  %8d  %8d\n', test_val(1:2,i), test_level(i) );
  end

  test_level = [];
  test_val = [];