Getting residual "nan" from the start

variational formulation
2

Perhaps Default absolute tolerance and relative tolerance - #4 by nate can help.