Tightened epsilon, separated out pdesolve to make debugging easier