now two implementations:
1) kappa^alpha evolves separately for each slip system alpha
2) kappa evolves the same for each slip system alpha based on gamma^p
new classes to form shape function matricies for different interpolations
for different fields u and gamma^p