def helmert_inverse(n): #*****************************************************************************80 # ## HELMERT_INVERSE returns the inverse of the HELMERT matrix. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 26 March 2015 # # Author: # # John Burkardt # # Parameters: # # Input, integer N, the order of A. # # Output, real A(N,N), the inverse matrix. # import numpy as np a = helmert(n) a = np.transpose(a) return a
def helmert_inverse ( n ): #*****************************************************************************80 # ## HELMERT_INVERSE returns the inverse of the HELMERT matrix. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 26 March 2015 # # Author: # # John Burkardt # # Parameters: # # Input, integer N, the order of A. # # Output, real A(N,N), the inverse matrix. # import numpy as np a = helmert ( n ) a = np.transpose ( a ) return a
def helmert_determinant_test(): #*****************************************************************************80 # ## HELMERT_DETERMINANT_TEST tests HELMERT_DETERMINANT. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 11 February 2015 # # Author: # # John Burkardt # from helmert import helmert from r8mat_print import r8mat_print print '' print 'HELMERT_DETERMINANT_TEST' print ' HELMERT_DETERMINANT computes the HELMERT determinant.' m = 5 n = m a = helmert(n) r8mat_print(m, n, a, ' HELMERT matrix:') value = helmert_determinant(n) print '' print ' Value = %g' % (value) print '' print 'HELMERT_DETERMINANT_TEST' print ' Normal end of execution.' return
def helmert_determinant_test ( ): #*****************************************************************************80 # ## HELMERT_DETERMINANT_TEST tests HELMERT_DETERMINANT. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 11 February 2015 # # Author: # # John Burkardt # from helmert import helmert from r8mat_print import r8mat_print print '' print 'HELMERT_DETERMINANT_TEST' print ' HELMERT_DETERMINANT computes the HELMERT determinant.' m = 5 n = m a = helmert ( n ) r8mat_print ( m, n, a, ' HELMERT matrix:' ) value = helmert_determinant ( n ) print '' print ' Value = %g' % ( value ) print '' print 'HELMERT_DETERMINANT_TEST' print ' Normal end of execution.' return