# include <cmath>
# include <cstdlib>
# include <cstring>
# include <ctime>
# include <iomanip>
# include <iostream>
using namespace std;
# include "lobatto_polynomial.hpp"
int main ( );
void lobatto_polynomial_value_test ( );
void lobatto_polynomial_derivative_test ( );
void lobatto_polynomial_plot_test ( );
void timestamp ( );
//****************************************************************************80
int main ( )
//****************************************************************************80
//
// Purpose:
//
// lobatto_polynomial_test() tests lobatto_polynomial().
//
// Licensing:
//
// This code is distributed under the MIT license.
//
// Modified:
//
// 20 November 2014
//
// Author:
//
// John Burkardt
//
{
 timestamp ( );
 cout << "\n";
 cout << "LOBATTO_POLYNOMIAL_TEST():\n";
 cout << " C++ version.\n";
 cout << " Test LOBATTO_POLYNOMIAL().\n";
 lobatto_polynomial_value_test ( );
 lobatto_polynomial_derivative_test ( );
 lobatto_polynomial_plot_test ( );
//
// Terminate.
//
 cout << "\n";
 cout << "LOBATTO_POLYNOMIAL_TEST():\n";
 cout << " Normal end of execution.\n";
 cout << "\n";
 timestamp ( );
 return 0;
}
//****************************************************************************80
void lobatto_polynomial_value_test ( )
//****************************************************************************80
//
// Purpose:
//
// LOBATTO_POLYNOMIAL_VALUE_TEST tests LOBATTO_POLYNOMIAL_VALUE.
//
// Licensing:
//
// This code is distributed under the MIT license.
//
// Modified:
//
// 20 November 2014
//
// Author:
//
// John Burkardt
//
{
 double e;
 double fx1;
 double fx2;
 double *l;
 int m;
 int n;
 int n_data;
 double x;
 double xvec[1];
 m = 1;
 cout << "\n";
 cout << "LOBATTO_POLYNOMIAL_VALUE_TEST:\n";
 cout << " LOBATTO_POLYNOMIAL_VALUES stores values of\n";
 cout << " the completed Lobatto polynomial L(n,x).\n";
 cout << " LOBATTO_POLYNOMIAL_VALUE evaluates the completed Lobatto polynomial.\n";
 cout << "\n";
 cout << " Tabulated Computed\n";
 cout << " N X L(N,X) L(N,X) Error\n";
 cout << "\n";
 n_data = 0;
 for ( ; ; )
 {
 lobatto_polynomial_values ( n_data, n, x, fx1 );
 if ( n_data == 0 )
 {
 break;
 }
 xvec[0] = x;
 l = lobatto_polynomial_value ( m, n, xvec );
 fx2 = l[0+(n-1)*m];
 e = fx1 - fx2;
 cout << " " << setw(4) << n
 << " " << setw(12) << x
 << " " << setw(12) << fx1
 << " " << setw(12) << fx2
 << " " << setw(8) << e << "\n";
 delete [] l;
 }
 return;
}
//****************************************************************************80
void lobatto_polynomial_derivative_test ( )
//****************************************************************************80
//
// Purpose:
//
// LOBATTO_POLYNOMIAL_DERIVATIVE_TEST tests LOBATTO_POLYNOMIAL_DERIVATIVE.
//
// Licensing:
//
// This code is distributed under the MIT license.
//
// Modified:
//
// 20 November 2014
//
// Author:
//
// John Burkardt
//
{
 double e;
 double fx1;
 double fx2;
 double *lp;
 int m;
 int n;
 int n_data;
 double x;
 double xvec[1];
 m = 1;
 cout << "\n";
 cout << "LOBATTO_POLYNOMIAL_DERIVATIVE_TEST:\n";
 cout << " LOBATTO_POLYNOMIAL_DERIVATIVES stores derivatives of\n";
 cout << " the completed Lobatto polynomial L(n,x).\n";
 cout << " LOBATTO_POLYNOMIAL_DERIVATIVE evaluates the completed Lobatto polynomial.\n";
 cout << "\n";
 cout << " Tabulated Computed\n";
 cout << " N X L''(N,X) L''(N,X) Error\n";
 cout << "\n";
 n_data = 0;
 for ( ; ; )
 {
 lobatto_polynomial_derivatives ( n_data, n, x, fx1 );
 if ( n_data == 0 )
 {
 break;
 }
 xvec[0] = x;
 lp = lobatto_polynomial_derivative ( m, n, xvec );
 fx2 = lp[0+(n-1)*m];
 e = fx1 - fx2;
 cout << " " << setw(4) << n
 << " " << setw(12) << x
 << " " << setw(12) << fx1
 << " " << setw(12) << fx2
 << " " << setw(8) << e << "\n";
 delete [] lp;
 }
 return;
}
//****************************************************************************80
void lobatto_polynomial_plot_test ( )
//****************************************************************************80
//
// Purpose:
//
// LOBATTO_POLYNOMIAL_PLOT_TEST tests LOBATTO_POLYNOMIAL_PLOT.
//
// Licensing:
//
// This code is distributed under the MIT license.
//
// Modified:
//
// 20 November 2014
//
// Author:
//
// John Burkardt
//
{
 int ndx[7] = { 1, 2, 3, 4, 5, 6, 7 };
 int ndx_num = 7;
 string prefix = "test";
 cout << "\n";
 cout << "LOBATTO_POLYNOMIAL_PLOT_TEST:\n";
 cout << " LOBATTO_POLYNOMIAL_PLOT plots Lobatto polynomials.\n";
 lobatto_polynomial_plot ( ndx_num, ndx, prefix );
 return;
}
//****************************************************************************80
void timestamp ( )
//****************************************************************************80
//
// Purpose:
//
// TIMESTAMP prints the current YMDHMS date as a time stamp.
//
// Example:
//
// 31 May 2001 09:45:54 AM
//
// Licensing:
//
// This code is distributed under the MIT license. 
//
// Modified:
//
// 08 July 2009
//
// Author:
//
// John Burkardt
//
{
# define TIME_SIZE 40
 static char time_buffer[TIME_SIZE];
 const struct std::tm *tm_ptr;
 std::time_t now;
 now = std::time ( NULL );
 tm_ptr = std::localtime ( &now );
 std::strftime ( time_buffer, TIME_SIZE, "%d %B %Y %I:%M:%S %p", tm_ptr );
 std::cout << time_buffer << "\n";
 return;
# undef TIME_SIZE
}
</div><div class="naked_ctrl">
<form action="/index.cgi/larger-text" method="get" name="gate">
<p><a href="http://altstyle.alfasado.net">AltStyle</a> によって変換されたページ <a href="https://people.sc.fsu.edu/~jburkardt/cpp_src/lobatto_polynomial_test/lobatto_polynomial_test.cpp">(-&gt;オリジナル)</a>
/ <label>アドレス: <input type="text" name="naked_post_url" value="https://people.sc.fsu.edu/~jburkardt/cpp_src/lobatto_polynomial_test/lobatto_polynomial_test.cpp" size="22" /></label> <label>モード: <select name="naked_post_mode">
<option value="default">デフォルト</option>
<option value="speech">音声ブラウザ</option>
<option value="ruby">ルビ付き</option>
<option value="contrast">配色反転</option>
<option value="larger-text" selected="selected">文字拡大</option>
<option value="mobile">モバイル</option>
</select>
<input type="submit" value="表示" />
</p>
</form>
</div>