def testProblemLoad(self): width = 10 height = 10 [nodes, boundary_nodes, tris] = generateRectangularMesh((width, height), (0,0), (1,1)) p = Problem(nodes, boundary_nodes, tris) f = p.getLoad(lambda x,y: x+y, lambda x,y: x*y, map(lambda x: 0, p.boundary_nodes))