Error: Singular Matrix

Solution Number: 953
Title: Error: Singular Matrix
Platform: All Platforms
Applies to: All Products
Versions: All versions
Categories: Error Messages
Keywords:

Problem Description

What does the error message "Singular Matrix" mean?

Solution

The problem is that the stiffness matrix of the linear system is singular and the linear solver cannot invert it. Examples of practical modeling situations where this can occur are

  • Your problem (equation system combined with the boundary conditions) is over- or underspecified. Check the equations and boundary conditions. You will often see this error if you have an ill-posed initial condition. It can be a good start to look at the initial condition and set it to some educated guess. Another example can be if you apply periodic boundary conditions when the PDE problem you solve does not permit such a solution. Yet another example is a Navier-Stokes problem where the pressure is unconstrained on the boundaries. You will need to make a point constraint somewhere.
  • You have specified a dependency of an expression defined in terms of coordinate derivatives of the solution variable. As the default initial guess into nonlinear systems is a constant (making the initial guess for the solution-derivative dependent expression zero), this can cause the equation to become singular.  The cure is to specify an initial value with a non-zero derivative, such as 1e-6*sqrt(x^2+y^2+z^2). This problem can sometimes occur when working with nonlinear materials, for example nonlinear magnetic materials (magnetic saturation effects).
  • You are solving a nonlinear eigenvalue problem with a zero linearization point. Examples of models that often result in a nonlinear eigenvalue problem are eigenfrequency analyses with PMLs or lossy materials, in RF or acoustics. To avoid the error message, you need to specify a non-zero, possibly complex-valued eigenvalue linearization point in the Study->Solver Configurations->Solver->Eigenvalue Solver settings window. Once you have a first solution, use the Global Evaluation under Results->Derived Values to evaluate the expression lambda for the solution that you are interested in. Enter its value as the new eigenvalue linearization point, solve again, and repeat the process until the eigenvalue no longer changes. This should usually not take more than a few iterations.

Disclaimer

COMSOL makes every reasonable effort to verify the information you view on this page. Resources and documents are provided for your information only, and COMSOL makes no explicit or implied claims to their validity. COMSOL does not assume any legal liability for the accuracy of the data disclosed. Any trademarks referenced in this document are the property of their respective owners. Consult your product manuals for complete trademark details.