Fixed bug in computeSources(). torus_2d still does not work. Ughh