In Ciarlet the last factor is left unspecified. I could not follow the
choice used in this implementation. I guess it was intended to minimize
the energy under a hydrostatic stress F = t*I at the unit matrix, i.e.,
for t=1. A simple calculation shows that (independent on the dimension)

d = 2a + 4b + 2c

has to hold in this case. Values in the test were adjusted.