Non-stationary two-stage relaxation based on the principle of aggregation multi-grid

We propose an automatization procedure to find an optimal value of the parameter used in some relaxation schemes. The basis of the method is formed by an aggregation-type procedure into a single aggregate and is meant to be applicable as a patch for existing codes.