Peryzna viscodamage needs to override the hardening function if it wants to have unique hardening because there is currently an out-of-bounds memory access that results in Peryzna viscoplasticity's setStateVariableEvolutionRateDerivatives function. The hardening function is assumed to only be a single element long BUT peryzna viscodamage accesses more. I will try and prevent this in a future pull request.