Yes.  We're replicating vector calculus identities that leave us with Functions that are precisely 0, so it's nicer to obtain that with a direct solver rather than having a tolerance >> machine precision.