stella PETSc framework for solving: \div dcoef \grad phi + bcoef phi = rhs [n . dcoef \grad phi] = jump phi = g_{dirichlet} on \Gamma_D n \grad phi = 0 on \Gamma_N Tools Single Block Testing