3.2.4.1. Norm Unbalance¶
This command is used to construct a convergence test which uses the norm of the right hand side of the matrix equation, i.e. \(b\) vector in \(Ax=b\), to determine if convergence has been reached. What the righthandside of the matrix equation is depends on integrator and constraint handler chosen. Usually, though not always, it is equal to the unbalanced forces in the system. The command to create a NormUnbalance test is the following:
 test NormUnbalance $tol $iter <$pFlag> <$nType>
Argument 
Type 
Description 

$tol 
float 
the tolerance criteria used to check for convergence 
$iter 
integer 
the max number of iterations to check before returning failure condition 
$pFlag 
integer 

$nType 
integer 
type of norm (optional: default is 2 (0 = maxnorm 1 = 1norm 2 = 2norm …)) 
Note
When using a penalty constraint handler, large forces (those necessary to enforce the constraint) are included in the \(b\) vector. Even for very small changes in the displacement, if user has selected overly large penalty factor, large forces can appear in the \(b\) vector.
Example:
The following examples demonstrate the command to create a NormUnbalance test which allows 10 iterations till failure with a 2norm in the \(b\) vector, i.e. \(\sqrt(b^T b)\) of 1.0e2.
Tcl Code
test NormUnbalance 1.0e2 10 2
Python Code
test('NormUnbalance', 1.0e2, 10, 2)
Code Developed by: fmk