Inserted an extremely interesting development in NS R1rho 2site.
If one do a transpose of M0, one can calculate all the matrix evolutions in the start via numpy einsum.
Since M0 is in higher a dimensions, one should not do a numpy transpose, but swap/roll the outer M0 6x1 axis.
Task #7807 (https://gna.org/task/index.php?7807): Speed-up of dispersion models for Clustered analysis.