Merge meca into threading

Keep main same as threading

Add class Material and CurrentCalculator

One change is that: calculate VCM as [theta * theta * J + phi * phi * J] (= r x J x r), or can't get result of same meaning as MECA