Old multi harmonic code working again + some additions.
Some performance issues are still there:
size of stack in Cal_WholeQuantity is to be modified if the number of harmonics is higher than 40 (NBR_MAX_HARMONIC);
assemble is very slow if -prealloc is not properly set;
for being able of increasing the size of the stack a bit more, I've added MAX_STACK_SIZE0 and fix its value to 2 => no tensors by default in the MH case;
real version of PetSC must still be used => feature not available in the pre-packaged version!
Working cases:
(nonlinear) ferromagnetic core with three windings (compumag2001, ieeetrmag 2002);