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

Showing results for

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

Marcias_Martinez1

Beginner

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

01-18-2011
04:39 AM

99 Views

Polar Decomposition

Can anyone help in determining the function call in IMKL 10.1 (using C and Visual Studio 2008) to do a polar decomposition of a deformation gradient matrix F that is 3 by 3.

I need to determine the rotation matrix such that F = RU

If anyone knows please advice,

Thanks

Marcias

Link Copied

1 Reply

mecej4

Black Belt

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

01-18-2011
06:35 AM

99 Views

A = W S V*

you can obtain the left polar decomposition

A = U P

using

P = V S V*

U = W V*

Similarly, you can obtain the right polar decomposition (which is what you sought)

A = P U

using

P = W S W*

and

U = W V*

as before. Please note that the factors P in the left and right polar decompositions are different.

MKL/Lapack provides routines to perform the SVD. You may consider the intended application of the polar decomposition to see if you can use P and U in factorized form instead of forming them explicitly. In particular, if you are going to multiply some vector by P or U, you can accomplish that with two or three matrix-vector multiplications, avoiding the more costly matrix-matrix multiplications that are involved in explicitly forming P and U.

- 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.