<div dir="ltr"><div><div><div>Dear Petsc user's<br><br>I am trying to use SuperLU_dist to compute LU decomposition of a non-square matrix. <br></div><br>My purpose is to build a basis of the nullspace of a matrix C . This matrix stores constraints (i.e. linear relationships between degrees of freedom): it is thus rectangular (size 9 x 33) and not square<br></div>I need to get the L U factors of C^T which is a tall skinny matrix. <br>It is possible when using directly SuperLu or SuperLU_dist<br></div>But when I use Petsc interface to SuperLU I get the following error : <br><div><div><div><br>[0]PETSC ERROR: Invalid argument<br>[0]PETSC ERROR: Must be square matrix, rows 9 columns 33<br>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Release Version 3.8.3, Dec, 09, 2017 <br>[0]PETSC ERROR: #1 MatGetOrdering() line 243 in /home/H03755/Librairies/petsc-3.8.3/src/mat/order/sorder.c<br>[0]PETSC ERROR: #2 PCSetUp_LU() line 84 in /home/H03755/Librairies/petsc-3.8.3/src/ksp/pc/impls/factor/lu/lu.c<br>[0]PETSC ERROR: #3 PCSetUp() line 924 in /home/H03755/Librairies/petsc-3.8.3/src/ksp/pc/interface/precon.c<br>[0]PETSC ERROR: #4 KSPSetUp() line 381 in /home/H03755/Librairies/petsc-3.8.3/src/ksp/ksp/interface/itfunc.c<br>[0]PETSC ERROR: #5 KSPSolve() line 612 in /home/H03755/Librairies/petsc-3.8.3/src/ksp/ksp/interface/itfunc.c<br>[0]PETSC ERROR: #6 User provided function() line 0 in User file<br><br></div><div>The use of non-square matrices is forbidden. <br></div><div><br>Is there a way to specify that I want to use a rectangular matrix ? How can I proceed to call SuperLU through PETSc with a rectangular matrix ? <br><br></div><div>Thanks a lot for your help <br></div><div>Best regards,<br><br></div><div>Natacha <br></div></div></div></div>