function triangle_contains_ijk_test ( )

%*****************************************************************************80
%
%% triangle_contains_ijk_test tests triangle_contains_ijk().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    04 November 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'triangle_contains_ijk_test\n' );
  fprintf ( 1, '  triangle_contains_ijk() reports whether a\n' );
  fprintf ( 1, '  triangle is in a region specified by a triangle list.\n' );

  for k = 1 : 1
    t = triangle_hexiamond ( k );
    m = size ( t, 1 );
    name = name_hexiamond ( k );
    name = [ '  ', name ];
    triangle_print ( t, name );
    for i = 1 : m
      fprintf ( 1, '\n' );
      neighbors = ijk_neighbors ( t(i,1:3) );
      fprintf ( 1, '  triangle %2d           (%2d,%2d,%2d)\n', i, t(i,1:3) );
      for j = 1 : 3
        if ( triangle_contains_ijk ( t, neighbors(j,1:3) ) )
          fprintf ( 1, '  internal neighbor %2d  (%2d,%2d,%2d)\n', ...
            j, neighbors(j,1:3) );
        end
      end
    end
  end
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'triangle_contains_ijk_test\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end

