Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- Pardiso w/ Quasi Newton method?

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

J_1

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-03-2017
11:41 AM

64 Views

Pardiso w/ Quasi Newton method?

I have read through the page describing using PARDISO to solve a nonlinear set of equations after they have been linearized, but I noticed the algorithm is just Newton's method applied with the PARDISO solving routing.

For faster and more reliable convergence, has anyone created a Quasi Newton solving algorithm using the PARDISO routine during the external iterations?

J

1 Solution

mecej4

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-09-2017
04:43 PM

64 Views

If the nonlinear equations being solved arise from discretization of nonlinear partial differential equations, application of a quasi-Newton method such as Broyden's scheme would result in a situation where a sparse solver such as Pardiso becomes worthy of consideration. See, for example, http://www.ams.org/journals/mcom/1970-24-109/S0025-5718-1970-0258276-9/ .

Many of the sparse matrices in the NIST Matrix Market collection (http://math.nist.gov/MatrixMarket/index.html) arose from FEA models of stress analysis, flow networks, etc.

Link Copied

3 Replies

Alexander_K_Intel2

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-09-2017
07:13 AM

64 Views

Hi,

I can be wrong, feel free to correct me, but Quasi Newton algorithm construct dense matrix, so it is not a point for Pardiso solver, is not it?

Thanks,

Alex

mecej4

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-09-2017
04:43 PM

65 Views

If the nonlinear equations being solved arise from discretization of nonlinear partial differential equations, application of a quasi-Newton method such as Broyden's scheme would result in a situation where a sparse solver such as Pardiso becomes worthy of consideration. See, for example, http://www.ams.org/journals/mcom/1970-24-109/S0025-5718-1970-0258276-9/ .

Many of the sparse matrices in the NIST Matrix Market collection (http://math.nist.gov/MatrixMarket/index.html) arose from FEA models of stress analysis, flow networks, etc.

J_1

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-10-2017
06:54 AM

64 Views

Alexander and Mecej4

You both are crazy smart, and both spot on. I only know enough to stir up trouble which is why I posted this. PARDISO is, without a doubt, needed to solve A*x = b where A is sparse, hands down no question. I am wanting to eventually update my solver to use a "quasi" Newton's method and was curious if anyone had an implementation posted on these forums or somewhere out there on the webs that utilized PARDISO in this type of iterative algorithm.

Mecej4 thanks for the links those are great resources!

J

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.