Add the simple reaction energy for heterogenous material

Author: ShiWeiHuH

src/structure/physics/ablation/ablation.jl

@@ -103,7 +103,13 @@ function update_new_bcs!(chunk::Peridynamics.AbstractBodyChunk, new_add_bcs_idx:
 
     for id in intersect(new_add_bcs_idx, chunk.system.chunk_handler.loc_points) 
         i_ch = chunk.system.chunk_handler.localizer[id]
-        chunk.storage.hsource[1, i_ch] = chunk.paramsetup.energy
+        paramsetup = chunk.paramsetup
+        if paramsetup isa Peridynamics.AbstractPointParameters
+            chunk.storage.hsource[1, i_ch] = chunk.paramsetup.energy
+        else
+            params_i = Peridynamics.get_params(paramsetup, i_ch)
+            chunk.storage.hsource[1, i_ch] = params_i.energy
+        end
     end
 end