Is it possible to update Schur complement quickly with MKL/Pardiso?

Intel® oneAPI Math Kernel Library

Ask questions and share information with other developers who use Intel® Math Kernel Library.

Is it possible to update Schur complement quickly with MKL/Pardiso?

207 Views

I have a matrix

A = [A11 A12;A21 A22];

and I'd like to compute the Schur complement of A11 in A using pardiso: S = A22 - A21 * A11^-1 * A12. This seems well supported by the PARDISO-SCHUR routine IPARM(38).

However, I'm wondering if I can do some prefactorization on A11 so that if I change only A12,A21,A22 then I could compute S quickly (quicker than just recomputing everything).

Is there an interface to alter part of A?

Link Copied

0 Replies

Community support is provided during standard business hours (Monday to Friday 7AM - 5PM PST). Other contact methods are available here.

Intel does not verify all solutions, including but not limited to any file transfers that may appear in this community. Accordingly, Intel disclaims all express and implied warranties, including without limitation, the implied warranties of merchantability, fitness for a particular purpose, and non-infringement, as well as any warranty arising from course of performance, course of dealing, or usage in trade.

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