On modelling of anisotropic viscoelasticity for soft tissue simulation: numerical solution and GPU execution