private static void r8vec_sqct_test() //****************************************************************************80 // // Purpose: // // R8VEC_SQCT_TEST tests R8VEC_SQCTB and R8VEC_SQCTF. // // Licensing: // // This code is distributed under the GNU LGPL license. // // Modified: // // 27 February 2010 // // Author: // // John Burkardt // { const double ahi = 5.0; const double alo = 0.0; const int n = 256; Console.WriteLine(""); Console.WriteLine("R8VEC_SQCT_TEST"); Console.WriteLine(" R8VEC_SQCTF does a forward slow quarter wave cosine transform;"); Console.WriteLine(" R8VEC_SQCTB does a backward slow quarter wave cosine transform."); Console.WriteLine(""); Console.WriteLine(" The number of data items is N = " + n + ""); // // Set the data values. // int seed = 123456789; double[] x = UniformRNG.r8vec_uniform_ab_new(n, alo, ahi, ref seed); typeMethods.r8vec_print_part(n, x, 10, " The original data:"); // // Compute the coefficients. // double[] y = Slow.r8vec_sqctf(n, x); typeMethods.r8vec_print_part(n, y, 10, " The cosine coefficients:"); // // Now compute inverse transform of coefficients. Should get back the // original data. // double[] z = Slow.r8vec_sqctb(n, y); typeMethods.r8vec_print_part(n, z, 10, " The retrieved data:"); }