Refactor RotM transformation out of WhitneyAV

[ettaka]

Original idea by @raback #429 (comment)

So the idea would be to refactor RotM functionality and its application to tensor transformation.
Use a test case, for example this and refactor the AV solver side such that the
conductivity and homogenization parameters are transformed correctly. The tensor would be treated like this:

TYPE(Rotate_t) :: Rotate
IF(Init) CALL ListInitRotation(Model, Rotate)

MagnetizationAtIp = ListGetElementReal3D(Mag_h,Basis,Element,...)
CALL ListGetElementRotate(Rotate,MagnetizationAtIp)

Instead of the old way:
Call to get the transformation matrix RotM

elmerfem/fem/src/modules/MagnetoDynamics/WhitneyAVHarmonicSolver.F90

Line 492 in 1f5be3a

CALL GetElementRotM(Element, RotM, n)

Then that is used to transform the conductivity tensor

elmerfem/fem/src/modules/MagnetoDynamics/WhitneyAVHarmonicSolver.F90

Line 1439 in 1f5be3a

C = MATMUL(MATMUL(RotMLoc, C),TRANSPOSE(RotMLoc))