Difference between revisions of "Resequence Command"

From OSUPDOCS
Jump to navigation Jump to search
 
Line 30: Line 30:
== Notes ==
== Notes ==


The resequencing is done using the "GPS Algorithm, named after Gibbs, Poole, and Stockmeyer (1976).<ref name='RS'>N. E. Gibbs, W. G. Poole, and P. K. Stockmeyer, "An Algorithm for Reducing the Bandwidth and Profile of a Sparse Matrix," <i>SIAM Journal of Numerical Analysis</i>, <b>13</b>, 236-250 (1976).</ref> Another bandwidth minimization method is the RCM method or the Reverse Cuthill-McKee method.<ref name="RCM">E. Cuthill and J. McKee. "Reducing the bandwidth of sparse symmetric matrices," In Proc. 24th Nat. Conf. ACM, 157–172 (1969).</ref> Now one method gets the absolute minimum. In testing, GPS and RCM get similar results but GPS is faster.
The resequencing is done using the "GPS Algorithm, named after Gibbs, Poole, and Stockmeyer (1976).<ref name='RS'>N. E. Gibbs, W. G. Poole, and P. K. Stockmeyer, "An Algorithm for Reducing the Bandwidth and Profile of a Sparse Matrix," <i>SIAM Journal of Numerical Analysis</i>, <b>13</b>, 236-250 (1976).</ref> Another bandwidth minimization method is the RCM method or the Reverse Cuthill-McKee method.<ref name="RCM">E. Cuthill and J. McKee. "Reducing the bandwidth of sparse symmetric matrices," In Proc. 24th Nat. Conf. ACM, 157–172 (1969).</ref> No method gets the absolute minimum. In testing, GPS and RCM get similar results but GPS is faster.


# It is best to start the resequencing at a node on the boundary of the object and probably on a corner. The final bandwidth may depend on the node selected for resequencing. The bandwidth is reported in FEA output results. You can vary the resequencing node to find the minimum value.
# It is best to start the resequencing at a node on the boundary of the object and probably on a corner. The final bandwidth may depend on the node selected for resequencing. The bandwidth is reported in FEA output results. You can vary the resequencing node to find the minimum value.

Latest revision as of 12:55, 4 February 2022

The Resequence command is used to renumber the numbers in an attempt to minimize the bandwidth of the problem. The smaller the bandwidth, the faster will be the FEA calculations.

Scripted Input Files

The two options in scripted input files are

Resequence (x),(y)

or

Resequence (id)

where

  • (x, y) define the coordinates (in length units) for a point (or (R,Z) coordinates if axisymmetric). The resequencing will start at the one node nearest to that point.
  • (id) is a previously defined keypoint. The resequencing will start at the node at that keypoint.

XML Input Files

In XML files, the two options are:

<Resequence x='(x)' y='(y)'/>

or

<Resequence keypt='(id)'/>

where (x), (y), and (id) are the same as defined above. Whichever method is used, it must be the last command in the single <GridBCs> block in the file.

Notes

The resequencing is done using the "GPS Algorithm, named after Gibbs, Poole, and Stockmeyer (1976).[1] Another bandwidth minimization method is the RCM method or the Reverse Cuthill-McKee method.[2] No method gets the absolute minimum. In testing, GPS and RCM get similar results but GPS is faster.

  1. It is best to start the resequencing at a node on the boundary of the object and probably on a corner. The final bandwidth may depend on the node selected for resequencing. The bandwidth is reported in FEA output results. You can vary the resequencing node to find the minimum value.
  2. Another use of this command is to verify mesh connectivity. Since disconnected sections of a static FEA mesh will cause a singular stiffness matrix, the calculations will fail. If you use a Resequence command on a disconnected mesh, it will detect the problem and abort the calculations.

References

  1. N. E. Gibbs, W. G. Poole, and P. K. Stockmeyer, "An Algorithm for Reducing the Bandwidth and Profile of a Sparse Matrix," SIAM Journal of Numerical Analysis, 13, 236-250 (1976).
  2. E. Cuthill and J. McKee. "Reducing the bandwidth of sparse symmetric matrices," In Proc. 24th Nat. Conf. ACM, 157–172 (1969).