Mercurial > hg > kwantix
changeset 49:7e0d6febb26d thesis-solved
Chose a better epsilon for SW problem. Thesis solved. w00t
author | Jordi Gutiérrez Hermoso <jordigh@gmail.com> |
---|---|
date | Sun, 25 Apr 2010 12:49:56 -0500 |
parents | 6f9071675ee7 |
children | 26111d9bfcb9 |
files | src/main-sw-rk4.cpp |
diffstat | 1 files changed, 3 insertions(+), 6 deletions(-) [+] |
line wrap: on
line diff
--- a/src/main-sw-rk4.cpp +++ b/src/main-sw-rk4.cpp @@ -125,7 +125,7 @@ //init the interps. conical::set_n(7); thin_plate_spline::set_n(6); - c_infty_rbf::set_epsilon(0.01); + c_infty_rbf::set_epsilon(20); cout << "Initialising... please wait." << endl; @@ -144,9 +144,6 @@ // of the bvp. k1(3,u0), k2(3,u0), k3(3,u0), k4(3,u0); - //debug - cout << "Size of interpolator is " << sizeof(k1[0]) << endl; - Fu<RBF_TYPE> f_u; Fv<RBF_TYPE> f_v; Fh<RBF_TYPE> f_h; @@ -158,12 +155,12 @@ f_v.set_dt(dt); f_h.set_dt(dt); - vtkplot plotter(h0,0,0.2); + vtkplot plotter(h0,0.08,0.12, true); plotter.set_view_direction(0,1,1); plotter.begin_interaction(); //main loop, timestepping with RK4. - size_t maxiter = 200; + size_t maxiter = 1000; for(size_t i = 1; i <= maxiter; i++){ cout << "Now on iteration #" << i << endl; plotter.update_values(h0.at());