The ReLU was first used by Alston Householder in 1941 as a mathematical abstraction of biological neural networks.[9]
Kunihiko Fukushima in 1969 used ReLU in the context of visual feature extraction in hierarchical neural networks.[10][11] 30 years later, Hahnloser et al. argued that ReLU approximates the biological relationship between neural firing rates and input current, in addition to enabling recurrent neural network dynamics to stabilise under weaker criteria.[12][13]
Jarrett et al. (2009) noted that rectification by either absolute or ReLU (which they called "positive part") was critical for object recognition in convolutional neural networks (CNNs), specifically because it allows average pooling without neighboring filter outputs cancelling each other out. They hypothesized that the use of sigmoid or tanh was responsible for poor performance in previous CNNs.[15]
Nair and Hinton (2010) made a theoretical argument that the softplus activation function should be used, in that the softplus function numerically approximates the sum of an exponential number of linear models that share parameters. They then proposed ReLU as a good approximation to it. They noted that ReLU allows "intensity equivariance" in the sense that multiplying input image by a constant multiplies the output also, which is not true for sigmoid or tanh. They found that ReLU activation allowed good empirical performance in restricted Boltzmann machines.[16]
Glorot et al (2011) argued that ReLU has the following advantages over sigmoid or tanh. ReLU is more similar to biological neurons' responses in their main operating regime. ReLU avoids vanishing gradients. ReLU is cheaper to compute. ReLU creates sparse representation naturally, because many hidden units output exactly zero for a given input. They also found empirically that deep networks trained with ReLU can achieve strong performance without unsupervised pre-training, especially on large, purely supervised tasks.[4]
Better gradient propagation: fewer vanishing gradient problems compared to sigmoidal activation functions that saturate in both directions.[4]
Efficiency: only requires comparison and addition.
Scale-invariant (homogeneous, or "intensity equivariance"[16]):
.
Potential problems
Possible downsides can include:
Non-differentiability at zero (however, it is differentiable anywhere else, and the value of the derivative at zero can be chosen to be 0 or 1 arbitrarily).
Not zero-centered: ReLU outputs are always non-negative. This can make it harder for the network to learn during backpropagation, because gradient updates tend to push weights in one direction (positive or negative). Batch normalization can help address this.[citation needed]
ReLU is unbounded.
Redundancy of the parametrization: Because ReLU is scale-invariant, the network computes the exact same function by scaling the weights and biases in front of a ReLU activation by , and the weights after by .[4]
Dying ReLU: ReLU neurons can sometimes be pushed into states in which they become inactive for essentially all inputs. In this state, no gradients flow backward through the neuron, and so the neuron becomes stuck in a perpetually inactive state (it "dies"). This is a form of the vanishing gradient problem. In some cases, large numbers of neurons in a network can become stuck in dead states, effectively decreasing the model capacity and potentially even halting the learning process. This problem typically arises when the learning rate is set too high. It may be mitigated by using "leaky" ReLU instead, where a small positive slope is assigned for . However, depending on the task, performance may be reduced.
Variants
Piecewise-linear variants
Leaky ReLU allows a small, positive gradient when the unit is inactive,[6] helping to mitigate the vanishing gradient problem. This gradient is defined by a parameter , typically set to 0.01–0.3.[17][18]
The same function can also be expressed without the piecewise notation as:
Parametric ReLU (PReLU) takes this idea further by making a learnable parameter along with the other network parameters.[19]
Concatenated ReLU (CReLU) preserves positive and negative phase information:[20]
Other non-linear variants
DELU
ExtendeD Exponential Linear Unit (DELU) is an activation function which is smoother within the neighborhood of zero and sharper for bigger values, allowing better allocation of neurons in the learning process for higher performance. Thanks to its unique design, it has been shown that DELU may obtain higher classification accuracy than ReLU and ELU.[21]
In these formulas, , and are hyperparameter values which could be set as default constraints , and , as done in the original work.
This activation function is illustrated in the figure at the start of this article. It has a "bump" to the left of x < 0 and serves as the default activation for models such as BERT.[22]
The SiLU (sigmoid linear unit) or swish function[23] is another smooth approximation which uses the sigmoid function, first introduced in the GELU paper:[22]
which is called the softplus[24][4] or SmoothReLU function.[25] For large negative it is roughly , so just above 0, while for large positive it is roughly , so just above .
This function can be approximated as:
By making the change of variables , this is equivalent to
The multivariable generalization of single-variable softplus is the LogSumExp with the first argument set to zero:
The LogSumExp function is
and its gradient is the softmax; the softmax with the first argument set to zero is the multivariable generalization of the logistic function. Both LogSumExp and softmax are used in machine learning.
ELU
Exponential linear units try to make the mean activations closer to zero, which speeds up learning. It has been shown that ELUs can obtain higher classification accuracy than ReLUs.[26]
In these formulas, is a hyperparameter to be tuned with the constraint .
Given the same interpretation of , ELU can be viewed as a smoothed version of a shifted ReLU (SReLU), which has the form .
Mish
The mish function can also be used as a smooth approximation of the rectifier.[23] It is defined as
where is a hyperparameter that determines the "size" of the curved region near . (For example, letting yields ReLU, and letting yields the metallic mean function.)
Squareplus shares many properties with softplus: It is monotonic, strictly positive, approaches 0 as , approaches the identity as , and is smooth. However, squareplus can be computed using only algebraic functions, making it well-suited for settings where computational resources or instruction sets are limited. Additionally, squareplus requires no special consideration to ensure numerical stability when is large.
^Fukushima, K. (1969). "Visual feature extraction by a multilayered network of analog threshold elements". IEEE Transactions on Systems Science and Cybernetics. 5 (4): 322–333. doi:10.1109/TSSC.1969.300225.
^Fukushima, K.; Miyake, S. (1982). "Neocognitron: A Self-Organizing Neural Network Model for a Mechanism of Visual Pattern Recognition". Competition and Cooperation in Neural Nets. Lecture Notes in Biomathematics. Vol. 45. Springer. pp. 267–285. doi:10.1007/978-3-642-46466-9_18. ISBN978-3-540-11574-8. {{cite book}}: |journal= ignored (help)
^Dugas, Charles; Bengio, Yoshua; Bélisle, François; Nadeau, Claude; Garcia, René (2000-01-01). "Incorporating second-order functional knowledge for better option pricing"(PDF). Proceedings of the 13th International Conference on Neural Information Processing Systems (NIPS'00). MIT Press: 451–457. Since the sigmoid h has a positive first derivative, its primitive, which we call softplus, is convex.
^Clevert, Djork-Arné; Unterthiner, Thomas; Hochreiter, Sepp (2015). "Fast and Accurate Deep Network Learning by Exponential Linear Units (ELUs)". arXiv:1511.07289 [cs.LG].