Fixed bug in InterpolateGrid that only causes problem when you compil…

…e with debug flags
plumed · Feb 3, 2024 · ac2b111 · ac2b111 · PlumedBot · Feb 3, 2024
1 parent 7a34fac
commit ac2b111
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/gridtools/InterpolateGrid.cpp b/src/gridtools/InterpolateGrid.cpp
@@ -166,7 +166,7 @@ void InterpolateGrid::gatherStoredValue( const unsigned& valindex, const unsigne
 
 void InterpolateGrid::gatherForcesOnStoredValue( const Value* myval, const unsigned& itask, const MultiValue& myvals, std::vector<double>& forces ) const {
   std::vector<double> pos(output_grid.getDimension()); double ff = myval->getForce(itask);
-  output_grid.getGridPointCoordinates( itask, pos ); input_grid.applyForce( 0, pos, ff, forces );
+  output_grid.getGridPointCoordinates( itask, pos ); input_grid.applyForce( this, pos, ff, forces );
 }