Function almost_identity 

pub fn almost_identity<T: Number + Float>(x: T, m: T, n: T) -> T

returns soft clipping (in a cubic fashion) of x; let m be the threshold (anything above m stays unchanged), and n the value things will take when the signal is zero