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());