Radiation rework and subsystem

code/__HELPERS/maths.dm

@@ -43,6 +43,8 @@ GLOBAL_LIST_INIT(sqrtTable, list(1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4,
 /proc/Inverse(x)
 	return 1 / x
 
+#define InverseSquareLaw(initial_strength,cur_distance,initial_distance) (initial_strength*(initial_distance**2/cur_distance**2))
+
 /proc/IsAboutEqual(a, b, deviation = 0.1)
 	return abs(a - b) <= deviation
 

code/__HELPERS/radiation.dm

@@ -0,0 +1,35 @@
+/proc/get_rad_contents(atom/location, list/output=list()) // A special GetAllContents that doesn't search past things with rad insulation
+	. = output
+
+	if(!location)
+		return
+
+	output += location
+
+	var/datum/component/rad_insulation/insulation = location.GetComponent(/datum/component/rad_insulation)
+	if(insulation && insulation.protects)
+		return
+	
+	for(var/i in 1 to location.contents.len)
+		var/static/list/ignored_things = typecacheof(list(/mob/dead, /obj/effect, /obj/docking_port, /turf, /atom/movable/lighting_object))
+		var/atom/thing = location.contents[i]
+		if(ignored_things[thing.type])
+			continue
+		get_rad_contents(thing, output)
+
+/proc/radiation_pulse(turf/epicenter, intensity, range_modifier, log=FALSE, can_contaminate=TRUE)
+	if(!SSradiation.can_fire)
+		return
+	for(var/dir in GLOB.cardinals)
+		new /datum/radiation_wave(epicenter, dir, intensity, range_modifier, can_contaminate)
+
+	var/list/things = get_rad_contents(epicenter) //copypasta because I don't want to put special code in waves to handle their origin
+	for(var/k in 1 to things.len)
+		var/atom/thing = things[k]
+		if(!thing)
+			continue
+		thing.rad_act(intensity)
+
+	if(log)
+		log_game("Radiation pulse with intensity:[intensity] and range modifier:[range_modifier] in area [epicenter.loc.name] ")
+	return TRUE