program main

!*****************************************************************************80
!
!! MAIN is the main program for CC_TO_ST_TEST.
!
!  Discussion:
!
!    CC_TO_ST_TEST tests the CC_TO_ST library.
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    23 July 2014
!
!  Author:
!
!    John Burkardt
!
  implicit none

  call timestamp ( )
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'CC_TO_ST_TEST'
  write ( *, '(a)' ) '  FORTRAN90 version'
  write ( *, '(a)' ) '  Test the CC_TO_ST library.'

  call test01 ( )
  call test02 ( )
!
!  Terminate.
!
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'CC_TO_ST_TEST'
  write ( *, '(a)' ) '  Normal end of execution.'
  write ( *, '(a)' ) ' '
  call timestamp ( )

  stop 0
end
subroutine test01 ( )

!*****************************************************************************80
!
!! TEST01 tests CC_TO_ST using a 1-based matrix.
!
!  Discussion:
!
!    This test uses a trivial matrix whose full representation is:
!
!          2  3  0  0  0
!          3  0  4  0  6
!      A = 0 -1 -3  2  0
!          0  0  1  0  0
!          0  4  2  0  1
!
!    The 1-based CC representation is
!
!      #  ICC  CCC  ACC
!     --  ---  ---  ---
!      1    1    1    2
!      2    2         3
!
!      3    1    3    3
!      4    3        -1
!      5    5         4
!
!      6    2    6    4
!      7    3        -3
!      8    4         1
!      9    5         2
!
!     10    3   10    2
!
!     11    2   11    6
!     12    5         1
!
!     13    *   13
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    18 July 2014
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer ( kind = 4 ), parameter :: m = 5
  integer ( kind = 4 ), parameter :: n = 5
  integer ( kind = 4 ), parameter :: ncc = 12

  real ( kind = 8 ), dimension ( ncc ) :: acc = (/ &
    2.0,  3.0, &
    3.0, -1.0,  4.0, &
    4.0, -3.0,  1.0, 2.0, &
    2.0, &
    6.0, 1.0 /)
  real ( kind = 8 ) ast(ncc)
  integer ( kind = 4 ), dimension ( n + 1 ) :: ccc = (/ &
    1, 3, 6, 10, 11, 13 /)
  integer ( kind = 4 ) i
  integer ( kind = 4 ), dimension ( ncc ) :: icc = (/ &
    1, 2, &
    1, 3, 5, &
    2, 3, 4, 5, &
    3, &
    2, 5 /)
  integer ( kind = 4 ) ist(ncc)
  integer ( kind = 4 ) jst(ncc)
  integer ( kind = 4 ) nst

  write ( *, '(a)' ) ''
  write ( *, '(a)' ) 'TEST01'
  write ( *, '(a)' ) '  Convert a 1-based CC matrix to ST format.'
!
!  Print the CC matrix.
!
  call cc_print ( m, n, ncc, icc, ccc, acc, '  The CC matrix:' )
!
!  Convert it.
!
  call cc_to_st ( m, n, ncc, icc, ccc, acc, nst, ist, jst, ast )
!
!  Print the ST matrix.
!
  call st_print ( m, n, nst, ist, jst, ast, '  The ST matrix:' )

  return
end
subroutine test02 ( )

!*****************************************************************************80
!
!! TEST02  tests CC_TO_ST using a 0-based matrix.
!
!  Discussion:
!
!    This test uses a trivial matrix whose full representation is:
!
!          2  3  0  0  0
!          3  0  4  0  6
!      A = 0 -1 -3  2  0
!          0  0  1  0  0
!          0  4  2  0  1
!
!    The 0-based CC representation is
!
!      #  ICC  CCC  ACC
!     --  ---  ---  ---
!      0    0    0    2
!      1    1         3
!
!      2    0    2    3
!      3    2        -1
!      4    4         4
!
!      5    1    5    4
!      6    2        -3
!      7    3         1
!      8    4         2
!
!      9    2    9    2
!
!     10    1   10    6
!     11    4         1
!
!     12    *   12
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    23 July 2014
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer ( kind = 4 ), parameter :: m = 5
  integer ( kind = 4 ), parameter :: n = 5
  integer ( kind = 4 ), parameter :: ncc = 12

  real ( kind = 8 ), dimension ( ncc ) :: acc = (/ &
    2.0,  3.0, &
    3.0, -1.0,  4.0, &
    4.0, -3.0,  1.0, 2.0, &
    2.0, &
    6.0, 1.0 /)
  real ( kind = 8 ) ast(ncc)
  integer ( kind = 4 ), dimension ( n + 1 ) :: ccc = (/ &
    0, 2, 5, 9, 10, 12 /)
  integer ( kind = 4 ) i
  integer ( kind = 4 ), dimension ( ncc ) :: icc = (/ &
    0, 1, &
    0, 2, 4, &
    1, 2, 3, 4, &
    2, &
    1, 4 /)
  integer ( kind = 4 ) ist(ncc)
  integer ( kind = 4 ) jst(ncc)
  integer ( kind = 4 ) nst

  write ( *, '(a)' ) ''
  write ( *, '(a)' ) 'TEST02'
  write ( *, '(a)' ) '  Convert a 0-based CC matrix to ST format.'
!
!  Print the CC matrix.
!
  call cc_print ( m, n, ncc, icc, ccc, acc, '  The CC matrix:' )
!
!  Convert it.
!
  call cc_to_st ( m, n, ncc, icc, ccc, acc, nst, ist, jst, ast )
!
!  Print the ST matrix.
!
  call st_print ( m, n, nst, ist, jst, ast, '  The ST matrix:' )

  return
end