Efficient design matrix free procedure for solving a soft maximin problem for large scale array-tensor structured models. Currently Lasso and SCAD penalized estimation is implemented.