function is_gaussian_prime_real_test ( )

%*****************************************************************************80
%
%% is_gaussian_prime_real_test() tests is_gaussian_prime_real().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    26 January 2023
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'is_gaussian_prime_real_test():\n' );
  fprintf ( 1, '  is_gaussian_prime_real() determines if a real integer\n' );
  fprintf ( 1, '  is a Gaussian prime.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '           p          is_gaussian_prime(p)\n' );
  fprintf ( 1, '\n' );

  n_data = 0;

  while ( true )

    [ n_data, p ] = gaussian_prime_real_values ( n_data );

    if ( n_data == 0 )
      break;
    end

    tf = is_gaussian_prime ( p );

    if ( tf )
      ch = 'true';
    else
      ch = 'false';
    end

    fprintf ( 1, '%12d  %s\n', p, ch );

  end

  return
end
