Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- Schur with Sparse PARDISO Complex matrix

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

marcsolal

Beginner

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

01-27-2017
11:19 AM

31 Views

Schur with Sparse PARDISO Complex matrix

I am using PARDISO to compute a Schur complement matrix. I am testing with real matrix meaning that I am using complex numbers with 0 imaginary parts. When the matrix size is small, the imaginary part of the result is zero but when the size increases, I get a small imaginary part. For example, for a 110000X11000 matrix and a Schur matrix size of about 100x100, the imaginary part is about 1e-6 to 1e-5 smaller than the real part. I would assume the imaginary part should be always 0. I know that numeric error can exist but I would think it is impossible to have anything than 0 by combining zero numbers.Also, a 1e-6 accuracy in double precision does not look that good. Am I wrong? If it is needed I can provide an example. I will need to make a simplified code for that.

Thanks,

Marc

Link Copied

1 Reply

marcsolal

Beginner

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

01-27-2017
04:39 PM

31 Views

I am sorry. This was a mistake. I do not know how to remove this post.

Marc

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