function triangulation_test127 ( )

%*****************************************************************************80
%
%% triangulation_test127 tests triangulation_order3_adjacency.
%
%  Discussion:
%
%    41--42--43--44  45--46--47--48
%     | \ | \ | \ |   | \ | \ | \ |
%    33--34--35--36  37--38--39--40
%     | \ |                   | \ |
%    29--30                  31--32
%     | \ |                   | \ |
%    25--26                  27--28
%     | \ |                   | \ |
%    21--22                  23--24
%     | \ |                   | \ |
%    17--18                  19--20
%     | \ |                   | \ |
%     9--10--11--12--13--14--15--16
%     | \ | \ | \ | \ | \ | \ | \ |
%     1---2---3---4---5---6---7---8
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    02 March 2014
%
%  Author:
%
%    John Burkardt
%
  element_num = 46;
  node_num = 48;
  
  element_node = [ ...
     1,  2,  9; ...
     2, 10,  9; ...
     2,  3, 10; ...
     3, 11, 10; ...
     3,  4, 11; ...
     4, 12, 11; ...
     4,  5, 12; ...
     5, 13, 12; ...
     5,  6, 13; ...
     6, 14, 13; ...
     6,  7, 14; ...
     7, 15, 14; ...
     7,  8, 15; ...
     8, 16, 15; ...
     9, 10, 17; ...
    10, 18, 17; ...
    15, 16, 19; ...
    16, 20, 19; ...
    17, 18, 21; ...
    18, 22, 21; ...
    19, 20, 23; ...
    20, 24, 23; ...
    21, 22, 25; ...
    22, 26, 25; ...
    23, 24, 27; ...
    24, 28, 27; ...
    25, 26, 29; ...
    26, 30, 29; ...
    27, 28, 31; ...
    28, 32, 31; ...
    29, 30, 33; ...
    30, 34, 33; ...
    31, 32, 39; ...
    32, 40, 39; ...
    33, 34, 41; ...
    34, 42, 41; ...
    34, 35, 42; ...
    35, 43, 42; ...
    35, 36, 43; ...
    36, 44, 43; ...
    37, 38, 45; ...
    38, 46, 45; ...
    38, 39, 46; ...
    39, 47, 46; ...
    39, 40, 47; ...
    40, 48, 47 ]';

  fprintf ( 1, '\n' );
  fprintf ( 1, 'triangulation_test127\n' );
  fprintf ( 1, '  For an order3 triangulation:\n' );
  fprintf ( 1, '  TRIANGULATION_ORDER3_ADJACENCY sets the full\n' );
  fprintf ( 1, '  adjacency matrix.\n' );

  adj = triangulation_order3_adjacency ( node_num, element_num, element_node );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Adjacency matrix:\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '               1         2         3         4        \n' );
  fprintf ( 1, '      123456789012345678901234567890123456789012345678\n' );
  fprintf ( 1, '\n' );
  for i = 1 : node_num
    fprintf ( 1, '  %2d  ', i );
    for j = 1 : node_num
      fprintf ( 1, '%d', adj(i,j) );
    end
    fprintf ( 1, '\n' );
  end

  return
end
