For the proposed model, we implement a block gradient descent algorithm to optimize the absolute error loss. Except for the last four parameters in ζ, all the other parameters are updated within a subgradient descent algorithm. The last four parameters are pre-specified and the details are in Section S2 in the Supplementary Information.