The SEL code is described in [A.H.Glasser and X. Tang, Comput. Phys. Comm. 164 2004 p237-243].

The number of degrees of freedom (DOF) for a single scalar field is NX x NY x NP**2, where NP is the degree of the polynomial. Static condensation reduces this to NX x NY x (2 x NP - 1), although at the expense of increasing the fill of the matrix. SEL solves for 8 variables simultaneously, thus the rank of the matrix they are factoring is nrank = 8 x NX x NY x (2 x NP - 1).

The parameters used for the 3 test problems are as follows:

• Resistive MHD, hi-viscosity:
• Resistive MHD, low-viscosity:
• 2-fluid MHD: (for 1/4 domain) NX = NY = 24, NP = 8

return to GEM Code Comparisons home page