|[eigen] Adding zeroes to sparse matrix alters minimum degree ordering|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: [eigen] Adding zeroes to sparse matrix alters minimum degree ordering
- From: Avi Robinson-Mosher <avi.mosher@xxxxxxxxx>
- Date: Wed, 2 Mar 2016 14:57:51 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:from:date:message-id:subject:to; bh=plBadPCh7iOEEhAHvJNe82DDttZ44SF2GN3E5iADJ4U=; b=oOPD1Y7hFqY5vlXZhwNtnxv31ErnCeIHmUwLI936MQyofBgIip9hGJIkRiHK+Ki6Lz F06ocUuQApN/v17xQc96NLJLC9fnW0UVcWWJmNVIZbhPZ2xxLbm4NwuCUaxQF//9nbi7 qnAEakN3HgPLd2bWxuXEjJpfLpneVbeOq5MABt7+Ef05nHdx/EFhMT4EcD8X+xfV1dRE r46l0zHncsjSUh+mfjQJqwjyAV3yldl5S137EYKepE4ZfjqEc/gaVNv0ayfabILIMiFU vxW/7B4cNqFqNzrS1/ihHIJdig6TJgplDuyg5BPg2s2EU2LTFE23Acxs8NE9pDlxtzDi +AIA==
I'm encountering an issue with Incomplete Cholesky. I have a case where the factorization is successful for a particular matrix, but when I add zeroes to the matrix (that is, set entries in the matrix that were implicit zeroes to explicit zeroes) the factorization fails. I don't completely understand why this happens, but I have observed that the permutation produced by AMDOrdering is different for the two cases (and that's a level of magic too dark for me to dig into directly).
I haven't produced a minimum working example (I'm only observing this on medium-sized systems at the moment, 83x83), but I will attempt to do so if necessary. I'm hoping that someone may have an insight without it though.