Intel® Fortran Compiler
Build applications that can scale for the future with optimized code designed for Intel® Xeon® and compatible processors.

iteration observation - math

JohnNichols
Valued Contributor III
483 Views

I am using MKL to solve a nonlinear problem first with an inversion and then eigenvector solution, interestingly if the problem is well formed the solution converges quickly, 3 to 10 steps. But if the problem becomes not well formed the iterations can take 9000 or so steps for even a minor change in the loading vector.  I developed a problem that takes about 30 steps and the error reduces at an exponential rate - shown on the graph. 

I was thinking the non-linear solver is estimating a set of points, if I measure the change in the set of points in the first few iterations - I wonder if i can use newton ralphson to get a better estimate of the final solution points - the problem is there are 360 degrees of freedom, but likely only a few are non-linear -

Any ideas would be welcome. 

0 Kudos
1 Reply
avinashs
New Contributor I
483 Views

It is difficult to know the exact cause of the slowdown without seeing the full problem. However, the symptoms you describe are characteristic of successive substitution or fixed point iteration algorithms. For chemical engineering models that I deal with, a method known as Dominant Eigenvalue Acceleration has been successful in speeding up the convergence in difficult cases. It can be applied every few iterations and has many parameters that can be adjusted. If your system is dense and has 360 degrees of freedom, Newton's method may have other problems in its application.

0 Kudos
Reply