Surface of mesh element

Hello community,
Lately i want to compute a surface of each element of the mesh, for example :

pwc1 = Compress(L2(mesh, order=0, definedon=mesh.Materials(“Rotor001”)))
gf1 = GridFunction(pwc1)

I want to compute surfaces of all meshs in gf1. My actual method is to integrate each element ( it takes too much time)

So is there a way to know vertices of each element of gf1 like this i will do the calculus manually ?

Best,
LB

How do you integrate each element now? One easy way would be

areas = Integrate(1, mesh, element_wise=True)