From RobustLight2.dm:
		set_color(red, green, blue, queued_run = 0)

			if (src.r == red && src.g == green && src.b == blue && !queued_run)
				return

			/*
			src.r_des = red
			src.g_des = green
			src.b_des = blue
			*/

			//hello yes now it's ZeWaka exporting my hellcode implementations across the code
			//scientific reasoning provided by Mokrzycki, Wojciech & Tatol, Maciej. (2011).
			/*
			var/R_sr = ((red + src.r*255) /2) //average value of R components in the two compared colors

			var/deltaR2 = abs(red   - (src.r*255))**2
			var/deltaG2 = abs(blue  - (src.b*255))**2
			var/deltaB2 = abs(green - (src.g*255))**2
			*/
			//this is our weighted euclidean distance function, weights based on red component
			//var/color_delta =( (2+(R_sr/256))*deltaR2 + (4*deltaG2) + (2+((255-R_sr)/256))*deltaB2 )

			//DEBUG_MESSAGE("[x],[y]:[temperature], d:[color_delta], [red]|[green]|[blue] vs [src.*255]|[src.*255]|[src.*255]")

			/*
			// This breaks everything if a light's value is 0, so, begone
			if (color_delta < 144) //determined via E'' sampling in science paper above, 144=12^2
				Z_LOG_DEBUG("Lighting", "Color update would be ignored due to color_delta ([color_delta]) under 144. ([R_sr], [deltaR2] [deltaG2] [deltaB2])")
				return
			*/

			if (src.enabled)
				if (SHOULD_QUEUE)
					light_update_queue.queue(src)
					dirty_flags |= D_COLOR
					return

				var/strip_gen = ++RL_Generation
				var/list/affected = src.strip(strip_gen)

				src.r = red
				src.g = green
				src.b = blue
				src.precalc()

				for (var/turf/T in src.apply())
					T.RL_UpdateLight()
				for (var/turf/T in affected)
					if (T.RL_UpdateGeneration <= strip_gen)
						T.RL_UpdateLight()
			else
				src.r = red
				src.g = green
				src.b = blue
				src.precalc()

From FEA_Fire.dm:
		// now this is ss13 level code
		proc/set_real_color()
			var/input = temperature / 100

			var/red
			if (input <= 66)
				red = 255
			else
				red = input - 60
				red = 329.698727446 * (red ** -0.1332047592)
				red = max(0, min(red, 255))

			var/green
			if (input <= 66)
				green = max(0.001, input)
				green = 99.4708025861 * log(green) - 161.1195681661
			else
				green = input - 60
				green = 288.1221695283 * (green ** -0.0755148492)
			green = max(0, min(green, 255))

			var/blue
			if (input >= 66)
				blue = 255
			else
				if (input <= 19)
					blue = 0
				else
					blue = input - 10
					blue = 138.5177312231 * log(blue) - 305.0447927307
					blue = max(0, min(blue, 255))

			color = rgb(red, green, blue) //changing obj.color is not expensive, and smooths imperfect light transitions


			// dear lord i apologuize for this conditional which i am about to write and commit. please be with the starving pygmies down in new guinea amen

			//hello yes now it's ZeWaka with an even more hellcode implementation that makes no sense
			//scientific reasoning provided by Mokrzycki, Wojciech & Tatol, Maciej. (2011).

			var/R_sr = ((red + light.r*255) /2) //average value of R components in the two compared colors

			var/deltaR2 = abs(red   - (light.r*255))**2
			var/deltaG2 = abs(blue  - (light.b*255))**2
			var/deltaB2 = abs(green - (light.g*255))**2

			//this is our weighted euclidean distance function, weights based on red component
			var/color_delta =( (2+(R_sr/256))*deltaR2 + (4*deltaG2) + (2+((255-R_sr)/256))*deltaB2 )

			//DEBUG_MESSAGE("[x],[y]:[temperature], d:[color_delta], [red]|[green]|[blue] vs [light.r*255]|[light.g*255]|[light.b*255]")

			if (color_delta > 144) //determined via E'' sampling in science paper above, 144=12^2

				red /= 255
				green /= 255
				blue /= 255

				light.set_color(red, green, blue, queued_run = 1)
				light.enable(queued_run = 1)