On 14/07/14 10:28, Antonio Garcia-Uceda wrote: > Dear all, > > My name is Antonio Garcia-Uceda, and I'm working on implementing a FR > algorithm for hybrid meshes. > > Could you tell me where exactly in the pyFR sources I have to include > commands to print out the matrixes M1,....,M8 in which the FR equations > are casted? Many thanks.
The matrices are defined in shapes.py. You can place the relevant print/save statements inside the constructor of BaseShape. For example: print self.m0 will output some details about the m0 matrix. However, I would caution against the use of these matrices. The ordering of entries depends on the choice of face and flux point counting patterns along with how one decides to arrange values out in memory. Moreover, they also depend on the choice of reference element. Regards, Freddie.
signature.asc
Description: OpenPGP digital signature
