New method HasStateBlocks::perturb()

include/core/state_block/has_state_blocks.h

@@ -68,6 +68,10 @@ class HasStateBlocks
         unsigned int    getSize() const;
         unsigned int    getLocalSize() const;
 
+        // Perturb state
+        void            perturb(double amplitude = 0.01);
+
+
     private:
         std::string structure_;
         std::unordered_map<std::string, StateBlockPtr> state_block_map_;
@@ -251,6 +255,5 @@ inline unsigned int HasStateBlocks::getLocalSize() const
     return size;
 }
 
-
 } // namespace wolf
 #endif /* STATE_BLOCK_HAS_STATE_BLOCKS_H_ */

src/state_block/has_state_blocks.cpp

@@ -44,5 +44,14 @@ void HasStateBlocks::removeStateBlocks(ProblemPtr _problem)
     }
 }
 
+void HasStateBlocks::perturb(double amplitude)
+{
+    for (const auto& pair_key_sb : state_block_map_)
+    {
+        auto& sb = pair_key_sb.second;
+        if (!sb->isFixed())
+            sb->perturb(amplitude);
+    }
+}
 
 }