code working for full viscoelastic solution, full complement of test
cases