removed a bug affecting some restarts on the first FSI version of waveBEM

git-svn-id: svn://svn.sissa.it/openship@561 958231e9-8e66-0410-a4a0-f664109d2741

Author: andreamola

source/free_surface.cc

@@ -2068,12 +2068,16 @@ bool FreeSurface<dim>::solution_check(Vector<double> & solution,
           comp_dom.evaluate_ref_surf_distances(nodes_ref_surf_dist,true);
        differential_components();
        for (unsigned int i=0; i<comp_dom.vector_dh.n_dofs(); ++i)
-           if ( (comp_dom.vector_flags[i] & boat) &&
-                !(comp_dom.vector_flags[i] & near_water) )
+           if ( ((comp_dom.flags[i] & boat) &&
+                !(comp_dom.flags[i] & near_water) ) ||
+                (comp_dom.flags[i] & transom_on_water)  )
               {
               comp_dom.map_points(i) -= nodes_ref_surf_dist(i);
               comp_dom.old_map_points(i) = comp_dom.map_points(i);
-              solution(i) = comp_dom.map_points(i)-comp_dom.rigid_motion_map_points(i);
+              if (constraints.is_constrained(i))
+                 solution(i) = comp_dom.map_points(i);
+              else
+                 solution(i) = comp_dom.map_points(i)-comp_dom.rigid_motion_map_points(i);
               //cout<<i<<" "<<nodes_ref_surf_dist(i)<<endl;
               }
 

waveBem.prm

@@ -6,10 +6,10 @@ set Iges file name                  =  utilities/goteborgExperiments.iges
 set Boat displacement (Kg)          = 0.0 #674.077471967
 
 # se si conosce il valore di affondamento idrodinamico da aggiungere a quello idrostatico, inserirlo qui
-set Assigned sink (m)               = 0.02 #0.02
+set Assigned sink (m)               = 0.01 #0.02
 
 # se si conosce il valore di angolo di beccheggio idrodinamico da aggiungere a quello idrostatico, inserirlo qui
-set Assigned trim (rad)             = 0.02 #0.02
+set Assigned trim (rad)             = 0.0 #0.02
 
 # Per ora è l'utente che deve specificare, dopo avere guardato il file iges, da quanti edges
 # è composta una eventiale transom stern. Se la barca non ha transom stern, questo parametro è