<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Hello, <br><br>I want to obtain the inverse of a matrix, even though it's not advised by the FAQ ;-) and I have a problem that earlier responses didn't solve for me. Simply speaking : the inverse matrix that is a result of&nbsp; DI=Dinv as told by the FAQ, is empty for me. Here is my code(the significant part) :<br><br>&nbsp; ierr = MatCreate(PETSC_COMM_WORLD,&amp;D);CHKERRQ(ierr);<br>&nbsp; ierr = MatSetSizes(D,PETSC_DECIDE,PETSC_DECIDE,n,n);CHKERRQ(ierr);<br>&nbsp; ierr = MatSetFromOptions(D);CHKERRQ(ierr);<br><br>&nbsp; ierr = MatCreate(PETSC_COMM_WORLD,&amp;Dinv);CHKERRQ(ierr);<br>&nbsp; ierr = MatSetSizes(Dinv,PETSC_DECIDE,PETSC_DECIDE,n,n);CHKERRQ(ierr);<br>&nbsp; ierr = MatSetFromOptions(Dinv);CHKERRQ(ierr);<br><br>&nbsp; ierr = MatCreate(PETSC_COMM_WORLD,&amp;I);CHKERRQ(ierr);<br>&nbsp; ierr =
 MatSetSizes(I,PETSC_DECIDE,PETSC_DECIDE,n,n);CHKERRQ(ierr);<br>&nbsp; ierr = MatSetFromOptions(I);CHKERRQ(ierr);<br><br>(...)<br><br>&nbsp; ierr = MatAssemblyEnd(D,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr);<br>&nbsp; ierr = MatAssemblyEnd(I,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr);<br>&nbsp; ierr = MatAssemblyEnd(Dinv,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr);<br><br>&nbsp; MatView(D,PETSC_VIEWER_STDOUT_SELF);<br>&nbsp;&nbsp;&nbsp; <br>&nbsp; ierr = MatGetOrdering(D,&nbsp; MATORDERING_NATURAL,&nbsp; &amp;perm,&nbsp; &amp;iperm); CHKERRQ(ierr); &nbsp;&nbsp;&nbsp; <br>&nbsp; ierr = MatFactorInfoInitialize(&amp;info); CHKERRQ(ierr);<br>&nbsp; ierr = MatLUFactor(D, perm, iperm, &amp;info); CHKERRQ(ierr);<br><br>&nbsp; MatView(D,PETSC_VIEWER_STDOUT_SELF);<br>&nbsp; MatView(I,PETSC_VIEWER_STDOUT_SELF);<br>&nbsp;&nbsp;&nbsp; <br>&nbsp; ierr = MatMatSolve(D,I,Dinv); CHKERRQ(ierr);<br><br>&nbsp; MatView(Dinv,PETSC_VIEWER_STDOUT_SELF);<br><br>// Results in such output:<br><br>//====
 MATRIX D<br>row 0: (0, 2) <br>row 1: (1, 1) <br>row 2: (2, 1) <br>row 3: (3, 1) <br><br>//==== MATRIX D after LU decomposition<br>row 0: (0, 0.5) <br>row 1: (1, 1) <br>row 2: (2, 1) <br>row 3: (3, 1) <br><br>//==== Identity matrix<br>row 0: (0, 1) <br>row 1: (1, 1) <br>row 2: (2, 1) <br>row 3: (3, 1) <br><br>//==== inverse of D&nbsp; (result of MatMatSolve(D,I,Dinv))<br>row 0:<br>row 1:<br>row 2:<br>row 3:<br><br><br>The approach with :<br><br>&nbsp;&nbsp;&nbsp; ierr = MatGetOrdering(D,&nbsp; MATORDERING_NATURAL,&nbsp; &amp;perm,&nbsp; &amp;iperm); CHKERRQ(ierr); &nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; ierr = MatGetFactor(D, MAT_SOLVER_PETSC, MAT_FACTOR_LU,&amp;F);<br>&nbsp;&nbsp;&nbsp; ierr = MatFactorInfoInitialize(&amp;info); CHKERRQ(ierr);<br>&nbsp;&nbsp;&nbsp; ierr = MatLUFactorSymbolic(F,D,perm,iperm,&amp;info);CHKERRQ(ierr);<br>&nbsp;&nbsp;&nbsp; ierr = MatLUFactorNumeric(F,D,&amp;info);CHKERRQ(ierr);<br>&nbsp;&nbsp;&nbsp; ierr = MatLUFactor(D,
 perm, iperm, &amp;info); CHKERRQ(ierr);<br><br>Provides same output<br><br>Thank You for Your help.<br><br></td></tr></table>