implemention of additional relaxation mechanisms
At the current state of the software only L=1 is implemented. Implementation of L>1 can be imitated from sofi2D. Maybe an f4tensor is needed to cover all 3 dimensions and L mechanisms. Performance has to be checked after the implementation (compared to old L=1).