Let $f:\mathbb{R}^n\to \mathbb{R}$ that $f=\frac{1}{2}\|(Ax-b)_+\|^2$, where $A\in\mathbb{R}^{m\times n}$, $b \in \mathbb{R}^m$ and if $x\in \mathbb{R}^n$: $((x)_+)_i=max\{0,x_i\}$.
I think the gradient of this function is as bellow: \begin{equation} \nabla f(x)=A^T(Ax-b)_+ \end{equation} is It true?
Note that $f(x) = g(h(x))$, where $h(x) = Ax - b$ and $$ g(u) = \frac12 \| u_+ \|^2 = \frac12\max(u_1,0)^2 + \cdots + \frac12\max(u_m,0)^2. $$ The function $g$ is certainly differentiable, and $$ g'(u) = u_+^T. $$ Also, the derivative of $h$ is $h’(x) = A$. By the chain rule, $f$ is differentiable and $$ f'(x) = g'(h(x)) h'(x) = (Ax - b)_+^T A. $$ If we use the convention that $\nabla f(x)$ is a column vector, then $$ \nabla f(x) = f'(x)^T = A^T(Ax - b)_+. $$ This confirms that $f$ is a differentiable function and the formula you gave in your question is correct.