I wanted to know guass point coordinates used by form.

mesh=UnitCubeMesh(1,1,1)
x = SpatialCoordinate(mesh)
x.get_local()
---> 3 x.get_local()
AttributeError: 'SpatialCoordinate' object has no attribute 'get_local'

Can you tell me how I can get the guass points location of mesh.

Thanks @dokken for your response. The symfem really helps.

How can we get the determinant of jacobian matrix which is multiplied during assemble operation?
If not values, at least, I could get the formulation used by FEniCS for tetrahedron elements.

from dolfin import *
mesh=UnitCubeMesh(1,1,1)
x = SpatialCoordinate(mesh)
assemble((x[1]*x[2])*dx)

Can we make the mesh to brick element from default tetrahedron element ?