Hi Kiny,
As Freddie says, PyFR only deals with the divergence of these correction functions, and not the vector correction functions themselves.
The DG variants are the default coded into PyFR.
To try others you could adopt an approach similar to that which Brian adopted here for tensor product elements:
in which I believe he modified the code for calculating the m3 matrix to read from a file instead of calling gbasis_at (there should be a patch in the Supplementary Material).
One warning about this approach, is that the m3 coefficients you read from a file will have to have been generated in the context of a specific set of solution points, so if you run with anything but those solutions points it will lead to problems!
I guess another approach would be to generalise gbasis_at.
Thanks
Peter
Dr Peter Vincent MSci ARCS DIC PhD FRAeS
Reader in Aeronautics and EPSRC Fellow
Department of Aeronautics
Imperial College London
South Kensington
London
SW7 2AZ
UK