<div dir="ltr"><div><div><div><div><div><div><div><div>Hello everyone,<br><br></div>I want to use PCApplyBAorAB <span id="result_box" class="" lang="en"><span class="">symmetrically with</span></span> a preconditioner from an ICC factorization.<br>
</div>However, I don't understand why the result is not the same if I use a natural or a RCM ordering during the factorization.<br></div>This problem does not occur if I use this routine with PC_LEFT or PC_RIGHT<br><br>
</div>I have put in attachment a code to illustrate the problem.<br></div>I compute a symmetric matix A in R^{n*n} and a full rank matrix B in R^{m*n}<br></div>I compute an ICC factorization of  the SPD matrix A+B^TB to precondition A.<br>
</div><div><br></div>Thank you!<br><br></div>Sylvain Mercier<br><div><div><br><div><div><div><div><br><br></div></div></div></div></div></div></div>