%% TEST04 tests CCOPY.
%
%  Modified:
%
%    15 April 2006
%
%  Author:
%
%    John Burkardt
%
  clear

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST04\n' );
  fprintf ( 1, '  CCOPY copies one complex vector into another.\n' );

  for j = 1 : 10
    x(j) = 10 * j + j * i;
  end

  for j = 1 : 10
    y(j) = 20 * j + 2 * j * i;
  end

  for j = 1 : 5
    for k = 1 : 5
      a(j,k) = 10 * j + k * i;
    end
  end

  fprintf ( 1, '\n' );
  fprintf ( 1, '  X = \n' );
  fprintf ( 1, '\n' );
  for j = 1 : 10
    fprintf ( 1, '  %6d  %12f  %12f\n', j, real ( x(j) ), imag ( x(j) ) );
  end
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Y =\n' );
  fprintf ( 1, '\n' );
  for j = 1 : 10
    fprintf ( 1, '  %6d  %12f  %12f\n', j, real ( y(j) ), imag ( y(j) ) );
  end
  fprintf ( 1, '\n' );
  fprintf ( 1, '  A = \n' );
  fprintf ( 1, '\n' );
  for j = 1 : 5
    for k = 1 : 5
      fprintf ( 1, '  %10f%10f', real ( a(j,k) ), imag ( a(j,k) ) );
    end
    fprintf ( 1, '\n' );
  end

  y = ccopy ( 5, x, 1, y, 1 );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Y = CCOPY ( 5, X, 1, Y, 1 )\n' );
  fprintf ( 1, '\n' );
  for j = 1 : 10
    fprintf ( 1, '  %6d  %12f  %12f\n', j, real ( y(j) ), imag ( y(j) ) );
  end

  for j = 1 : 10
    y(j) = 20 * j + 2 * j * i;
  end

  y = ccopy ( 3, x, 2, y, 3 );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Y = CCOPY ( 3, X, 2, Y, 3 )\n' );
  fprintf ( 1, '\n' );
  for j = 1 : 10
    fprintf ( 1, '  %6d  %12f  %12f\n', j, real ( y(j) ), imag ( y(j) ) );
  end

  a(1:5,1) = ccopy ( 5, x, 1, a(1:5,1), 1 );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  A(1:5,1) = CCOPY ( 5, X, 1, A(1:5,1), 1 )\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  A = \n' );
  fprintf ( 1, '\n' );
  for j = 1 : 5
    for k = 1 : 5
      fprintf ( 1, '  %10f%10f', real ( a(j,k) ), imag ( a(j,k) ) );
    end
    fprintf ( 1, '\n' );
  end

  for j = 1 : 5
    for k = 1 : 5
      a(j,k) = 10 * j + k * i;
    end
  end

  a(1,1:5) = ccopy ( 5, x, 2, a(1,1:5), 1 );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  A(1,1:5) = CCOPY ( 5, X, 2, A(1,1:5), 1 )\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  A =\n' );
  fprintf ( 1, '\n' );
  for j = 1 : 5
    for k = 1 : 5
      fprintf ( 1, '  %10f%10f', real ( a(j,k) ), imag ( a(j,k) ) );
    end
    fprintf ( 1, '\n' );
  end