Thank you for your helpful response.
In my example, gfm is an array with length of 15 which each component is defined in the dual function space fesm. Let’s say gfm = (sigma, u) ; one can say the answer is gfm = [gfm1, gfm2, . . . , gfm15 ] where gfm1 = (sigma1, u1) ;
when the mulitidim =1, the following command works.
gfsigma = gfm.components[0]
gfu = gfm.components[1]
but when gfm is an array, how to call the sigma1? Thanks.