core: virtual_sites: Removed default argument from virtual function and non-functional virtual COM feature

src/core/CMakeLists.txt

@@ -125,11 +125,10 @@ set(EspressoCore_SRC
   observables/LBVelocityProfile.cpp
   virtual_sites/lb_inertialess_tracers.cpp
   virtual_sites/lb_inertialess_tracers_cuda_interface.cpp
-  virtual_sites/virtual_sites_com.cpp
   virtual_sites/VirtualSitesInertialessTracers.cpp
   virtual_sites/VirtualSitesRelative.cpp
   accumulators/TimeSeries.cpp
-        observables/PidObservable.cpp)
+  observables/PidObservable.cpp)
 
 if(CUDA)
   set(EspressoCuda_SRC

src/core/event.cpp

@@ -466,7 +466,7 @@ void on_ghost_flags_change() {
 
 void update_dependent_particles() {
 #ifdef VIRTUAL_SITES
-  virtual_sites()->update();
+  virtual_sites()->update(true);
 #endif
   cells_update_ghosts();
 

src/core/integrate.cpp

@@ -189,7 +189,7 @@ void integrate_vv(int n_steps, int reuse_forces) {
     lb_lbcoupling_deactivate();
 
 #ifdef VIRTUAL_SITES
-    virtual_sites()->update();
+    virtual_sites()->update(true);
 #endif
 
     // Communication step: distribute ghost positions
@@ -248,7 +248,7 @@ void integrate_vv(int n_steps, int reuse_forces) {
 
 // VIRTUAL_SITES pos (and vel for DPD) update for security reason !!!
 #ifdef VIRTUAL_SITES
-    virtual_sites()->update();
+    virtual_sites()->update(true);
 #endif
 
     // Communication step: distribute ghost positions

src/core/virtual_sites/VirtualSites.hpp

@@ -43,7 +43,7 @@ class VirtualSites {
    *  Velocities are only updated if get_have_velocity() returns true.
    *  @param recalc_positions  Skip the recalculation of positions if false.
    */
-  virtual void update(bool recalc_positions = true) const = 0;
+  virtual void update(bool recalc_positions) const = 0;
   /** Back-transfer forces (and torques) to non-virtual particles. */
   virtual void back_transfer_forces_and_torques() const = 0;
   /** @brief Called after force calculation (and before rattle/shake) */