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.
Best regards,
Antonio
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.