Beispiel #1
0
def milnes_determinant_test():

    #*****************************************************************************80
    #
    ## MILNES_DETERMINANT_TEST tests MILNES_DETERMINANT.
    #
    #  Licensing:
    #
    #    This code is distributed under the GNU LGPL license.
    #
    #  Modified:
    #
    #    19 February 2015
    #
    #  Author:
    #
    #    John Burkardt
    #
    from milnes import milnes
    from r8vec_uniform_ab import r8vec_uniform_ab
    from r8mat_print import r8mat_print

    print ''
    print 'MILNES_DETERMINANT_TEST'
    print '  MILNES_DETERMINANT computes the MILNES determinant.'

    m = 5
    n = m
    r8_lo = -5.0
    r8_hi = +5.0
    seed = 123456789
    x, seed = r8vec_uniform_ab(n, r8_lo, r8_hi, seed)

    a = milnes(m, n, x)

    r8mat_print(m, n, a, '  MILNES matrix:')

    value = milnes_determinant(n, x)

    print ''
    print '  Value =  %g' % (value)

    print ''
    print 'MILNES_DETERMINANT_TEST'
    print '  Normal end of execution.'

    return
Beispiel #2
0
def milnes_determinant_test ( ):

#*****************************************************************************80
#
## MILNES_DETERMINANT_TEST tests MILNES_DETERMINANT.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    19 February 2015
#
#  Author:
#
#    John Burkardt
#
  from milnes import milnes
  from r8vec_uniform_ab import r8vec_uniform_ab
  from r8mat_print import r8mat_print

  print ''
  print 'MILNES_DETERMINANT_TEST'
  print '  MILNES_DETERMINANT computes the MILNES determinant.'

  m = 5
  n = m
  r8_lo = -5.0
  r8_hi = +5.0
  seed = 123456789
  x, seed = r8vec_uniform_ab ( n, r8_lo, r8_hi, seed )

  a = milnes ( m, n, x )

  r8mat_print ( m, n, a, '  MILNES matrix:' )

  value = milnes_determinant ( n, x )

  print ''
  print '  Value =  %g' % ( value )

  print ''
  print 'MILNES_DETERMINANT_TEST'
  print '  Normal end of execution.'

  return
Beispiel #3
0
def milnes_test():

    #*****************************************************************************80
    #
    ## MILNES_TEST tests MILNES.
    #
    #  Licensing:
    #
    #    This code is distributed under the GNU LGPL license.
    #
    #  Modified:
    #
    #    18 February 2015
    #
    #  Author:
    #
    #    John Burkardt
    #
    from r8vec_uniform_ab import r8vec_uniform_ab
    from r8mat_print import r8mat_print

    print ''
    print 'MILNES_TEST'
    print '  MILNES computes the MILNES matrix.'

    m = 5
    n = m
    r8_lo = -5.0
    r8_hi = +5.0
    seed = 123456789
    x, seed = r8vec_uniform_ab(n, r8_lo, r8_hi, seed)

    a = milnes(m, n, x)

    r8mat_print(m, n, a, '  MILNES matrix:')

    print ''
    print 'MILNES_TEST'
    print '  Normal end of execution.'

    return
Beispiel #4
0
def milnes_test ( ):

#*****************************************************************************80
#
## MILNES_TEST tests MILNES.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    18 February 2015
#
#  Author:
#
#    John Burkardt
#
  from r8vec_uniform_ab import r8vec_uniform_ab
  from r8mat_print import r8mat_print

  print ''
  print 'MILNES_TEST'
  print '  MILNES computes the MILNES matrix.'

  m = 5
  n = m
  r8_lo = -5.0
  r8_hi = +5.0
  seed = 123456789
  x, seed = r8vec_uniform_ab ( n, r8_lo, r8_hi, seed )

  a = milnes ( m, n, x )
 
  r8mat_print ( m, n, a, '  MILNES matrix:' )

  print ''
  print 'MILNES_TEST'
  print '  Normal end of execution.'

  return